[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (861)

Search Parameters:
Keywords = MODIS NDVI

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
17 pages, 3410 KiB  
Article
The Aerosol Optical Depth Retrieval from Wide-Swath Imaging of DaQi-1 over Beijing
by Zhongting Wang, Ruijie Zhang, Ruizhi Chen and Hui Chen
Atmosphere 2024, 15(12), 1476; https://doi.org/10.3390/atmos15121476 - 10 Dec 2024
Viewed by 447
Abstract
The Wide-Swath Imaging (WSI) sensor is a Chinese satellite launched in 2022, capable of providing data at resolutions ranging from 75 to 600 m for monitoring aerosols, fire points, and dust, among other uses. In this study, we developed a Dark Dense Vegetation [...] Read more.
The Wide-Swath Imaging (WSI) sensor is a Chinese satellite launched in 2022, capable of providing data at resolutions ranging from 75 to 600 m for monitoring aerosols, fire points, and dust, among other uses. In this study, we developed a Dark Dense Vegetation method to retrieve the Aerosol Optical Depth (AOD) quickly from WSI 600 m data. First, after splitting into three types according to the Normalized Difference Vegetation Index (NDVI), we calculated the empirical parameters of land reflectance between the red (0.65 μm) and blue (0.47 μm) channels using Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance products over the Beijing area. Second, the decrease in the NDVI was simulated and analyzed under different AODs and solar zenith angles, and we introduced an iterative inversion approach to account for it. The simulation retrievals demonstrated that the iterative inversion produced accurate results after less than four iterations. Thirdly, we utilized the atmospherically corrected NDVI for dark target identification and output the AOD result. Finally, retrieval experiments were conducted using WSI 600 m data collected over Beijing in 2023. The retrieved AOD images highlighted two air pollution events occurring during 3–8 March and 27–31 October 2023. The inversion results in 2023 showed a strong correlation with Aerosol Robotic Network station data (the correlation coefficient was greater than 0.9). Our method exhibited greater accuracy than the MODIS aerosol product, though it was less accurate than the Multi-Angle Implementation of Atmospheric Correction product. Full article
Show Figures

Figure 1

Figure 1
<p>Filter response functions of WSI and MODIS in the channels which range from 380 nm to 900 nm.</p>
Full article ">Figure 2
<p>Flow chart of AOD retrieval method for WSI.</p>
Full article ">Figure 3
<p>The percentage histogram of pixel count over the Beijing area: (<b>a</b>) is surface reflectance in blue, green, red and NIR channels, and (<b>b</b>) is NDVI.</p>
Full article ">Figure 4
<p>The comparison of surface reflectance between red and blue: (<b>a</b>) is low vegetation, (<b>b</b>) is medium vegetation, (<b>c</b>) is high vegetation, and (<b>d</b>) is all of the vegetation. Dashed line represents the linear fitting line. The color represents the percentage of pixel numbers.</p>
Full article ">Figure 5
<p>The decreased NDVI under different AODs. The left is SZA = 21 degrees (<b>a</b>,<b>c</b>,<b>e</b>), while the right is SZA = 63 degrees (<b>b</b>,<b>d</b>,<b>f</b>). The top is low vegetation (<b>a</b>,<b>b</b>), the middle is medium vegetation (<b>c</b>,<b>d</b>), and the bottom is high vegetation (<b>e</b>,<b>f</b>).</p>
Full article ">Figure 6
<p>The maximum number of iterations changing with AOD and errors from measurements: (<b>a</b>) is low vegetation, (<b>b</b>) is medium vegetation, and (<b>c</b>) is high vegetation.</p>
Full article ">Figure 7
<p>Daily PM<sub>2.5</sub> concentration on 3–8 March 2023.</p>
Full article ">Figure 8
<p>WSI AOD images during a pollution event on 3–8 March 2023.</p>
Full article ">Figure 9
<p>Daily PM<sub>2.5</sub> concentration on 27–31 October 2023.</p>
Full article ">Figure 10
<p>WSI AOD images during a pollution event on 27–31 October 2023. (<b>a</b>) is October 27, (<b>b</b>) is October 28, (<b>c</b>) is October 29, (<b>d</b>) is October 30, and (<b>e</b>) is October 31.</p>
Full article ">Figure 10 Cont.
<p>WSI AOD images during a pollution event on 27–31 October 2023. (<b>a</b>) is October 27, (<b>b</b>) is October 28, (<b>c</b>) is October 29, (<b>d</b>) is October 30, and (<b>e</b>) is October 31.</p>
Full article ">Figure 11
<p>AOD in 2023 over AERONET Beijing station. The yellow is AERONET, the green is WSI, the red is MYD04, and the black is MAIAC.</p>
Full article ">Figure 12
<p>AOD validation of WSI (<b>a</b>) MAIAC, (<b>b</b>) and MYD04, (<b>c</b>) with AERONET data. N represents the number of matched pixels. (<b>d</b>) All matched points, where the green is WSI, the red is MYD04, and the black is MAIAC.</p>
Full article ">Figure 12 Cont.
<p>AOD validation of WSI (<b>a</b>) MAIAC, (<b>b</b>) and MYD04, (<b>c</b>) with AERONET data. N represents the number of matched pixels. (<b>d</b>) All matched points, where the green is WSI, the red is MYD04, and the black is MAIAC.</p>
Full article ">
21 pages, 8184 KiB  
Article
Estimation of Vegetation Carbon Sinks and Their Response to Land Use Intensity in the Example of the Beijing–Tianjin–Hebei Region
by Qing Yao, Junping Zhang, Huayang Song, Rongxia Yu, Nina Xiong, Jia Wang and Liu Cui
Forests 2024, 15(12), 2158; https://doi.org/10.3390/f15122158 - 6 Dec 2024
Viewed by 335
Abstract
Accurate regional carbon sequestration estimates are essential for China’s emission reduction and carbon sink enhancement efforts to address climate change. Enhancing the spatial precision of vegetation carbon sink estimates is crucial for a deeper understanding of the underlying response mechanisms, yet this remains [...] Read more.
Accurate regional carbon sequestration estimates are essential for China’s emission reduction and carbon sink enhancement efforts to address climate change. Enhancing the spatial precision of vegetation carbon sink estimates is crucial for a deeper understanding of the underlying response mechanisms, yet this remains a significant challenge. In this study, the Beijing–Tianjin–Hebei (BTH) region was selected as the study area. We employed the GF-SG (Gap filling and Savitzky–Golay filtering) model to fuse Landsat and MODIS data, generating high-resolution imagery to enhance the accuracy of NPP (Net Primary Productivity) and NEP (Net Ecosystem Productivity) estimates for this region. Subsequently, the Sen+MK model was used to analyze the spatiotemporal variations in carbon sinks. Finally, the land use intensity index, which reflects human activity disturbances, was applied, and the bivariate Moran’s spatial autocorrelation method was used to analyze the response mechanisms of carbon sinks. The results indicate that the fused GF-SG NDVI (Normalized Difference Vegetation Index) data provided highly accurate 30 m resolution imagery for estimating NPP and NEP. The spatial distribution of carbon sinks in the study area showed higher values in the northeastern forest regions, relatively high values in the southeastern plains, and lower values in the northwestern plateau and central urban areas. Additionally, 58.71% of the area exhibited an increasing trend, with 11.73% showing significant or strongly significant growth. A generally negative spatial correlation was observed between land use intensity and carbon sinks, with the impact of land use intensity on carbon sinks exceeding 0.3 in 2010. This study provides methodological insights for obtaining vegetation monitoring data and estimating carbon sinks in large urban agglomerations and offers scientific support for developing ecological and carbon reduction strategies in the BTH region. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of research area.</p>
Full article ">Figure 2
<p>NDVI fusion results: a general overview.</p>
Full article ">Figure 3
<p>NDVI fusion results: local view. (<b>a</b>–<b>c</b>) are MODIS, Landsat, and GF-SG NDVI in cloud-free zones; (<b>d</b>–<b>f</b>) are MODIS, Landsat, and GF-SG NDVI in cloud-covered zones.</p>
Full article ">Figure 4
<p>Sample points. (<b>a</b>) Spatial distribution of sample points in the BTH region. (<b>b</b>–<b>d</b>) are scatter plots of actual vs. simulated NPP in 2005, 2010, and 2015.</p>
Full article ">Figure 5
<p>Time series of mean NEP in the BTH region from 2000 to 2020.</p>
Full article ">Figure 6
<p>Spatial distribution of NEP in the BTH region from 2000 to 2020.</p>
Full article ">Figure 7
<p>NEP for the BTH region from 2000 to 2020.</p>
Full article ">Figure 8
<p>LISA clustering of NEP and LUI in the BTH region from 2000 to 2020.</p>
Full article ">Figure 9
<p>Spatial correlation between LUI and NEP in each city in the BTH region (2000–2020).</p>
Full article ">
28 pages, 7193 KiB  
Article
Country-Scale Crop-Specific Phenology from Disaggregated PROBA-V
by Henry Rivas, Nicolas Delbart, Fabienne Maignan, Emmanuelle Vaudour and Catherine Ottlé
Remote Sens. 2024, 16(23), 4521; https://doi.org/10.3390/rs16234521 - 2 Dec 2024
Viewed by 355
Abstract
Large-scale crop phenology monitoring is essential for agro-ecosystem policy. Remote sensing helps track crop development but requires high-temporal and spatial resolutions. While datasets with both attributes are now available, their large-scale applications require significant resources. Medium-resolution data offer daily observations but lack detail [...] Read more.
Large-scale crop phenology monitoring is essential for agro-ecosystem policy. Remote sensing helps track crop development but requires high-temporal and spatial resolutions. While datasets with both attributes are now available, their large-scale applications require significant resources. Medium-resolution data offer daily observations but lack detail for smaller plots. This study generated crop-specific phenomaps for mainland France (2016–2020) using PROBA-V data. A spatial disaggregation method reconstructed NDVI time series for individual crops within mixed pixels. Then, phenometrics were extracted from disaggregated PROBA-V and Sentinel-2 separately and compared to observed phenological stages. Results showed that PROBA-V-based phenomaps closely matched observations at regional level, with moderate accuracy at municipal level. PROBA-V demonstrated a higher detection rate than Sentinel-2, especially in cloudy periods, and successfully generated phenomaps before Sentinel-2B’s launch. The study highlights PROBA-V’s potential for operational crop monitoring, i.e., wheat heading and oilseed rape flowering, with performance comparable to Sentinel-2. PROBA-V outputs complement Sentinel-2: phenometrics cannot be generated at plot level but are efficiently produced at regional or national scales to study phenological gradients more easily than with Sentinel-2 and with similar accuracy. This approach could be extended to MODIS or SPOT-VGT, to generate historical phenological data, providing that a crop map is available. Full article
(This article belongs to the Special Issue Remote Sensing of Land Surface Phenology II)
Show Figures

Figure 1

Figure 1
<p>General flowchart of this study. Crop-specific phenological mapping procedure using (<b>A</b>) PROBA-V (300 m) and (<b>B</b>) Sentinel-2 (10 m) data. PROBA-V-based phenometrics were extracted at pixel level, while those from Sentinel-2 were extracted at plot level. Both are inter-compared and compared with TEMPO data at the municipal level. PROBA-V-based phenometrics are then compared with Céré’Obs data at the regional level. For phenometrics extraction, thresholds were calibrated in Block A using disaggregated PROBA-V NDVI time series and TEMPO data, and were then applied identically to both disaggregated PROBA-V and Sentinel-2 NDVI time series.</p>
Full article ">Figure 2
<p>Study area location. The validation of the phenomaps was carried out at regional and municipal scales. (<b>A</b>) The regional comparison between the phenometrics (from PROBA-V) and the observed (from Céré’Obs) median phenological dates was made across administrative regions, highlighted in both yellow and orange (n = 14), while the comparison between phenometrics and observed phenological progression (i.e., the percentage of area reaching a given phenological stage, as a function of date) focused on regions highlighted in orange only (n = 3). Regional ground data were not available for the three regions highlighted in gray. (<b>B</b>,<b>C</b>) The comparison at municipal level between the phenometrics (from PROBA-V and Sentinel-2) and the observed (from TEMPO) median phenological stage dates was made across points highlighted in blue (winter wheat) and red (oilseed rape) within inter-comparison sites. Inter-comparison sites were delimited by seven Sentinel-2 tiles: 30UWU, 31UCQ, 31UDP, 31TCN, 31UGQ, 31TFN, and 31TCJ. Additionally, PROBA-V-based phenomaps were also compared to TEMPO data available outside these inter-comparison sites. Finally, the green mask shows the winter wheat and oilseed rape areas declared in 2019 according to LPIS.</p>
Full article ">Figure 3
<p>NDVI time profile for (<b>A</b>) winter wheat and (<b>B</b>) oilseed rape. Values represent the national average derived from disaggregated PROBA-V (300 m) data in 2019, which were fitted using the Whittaker smoother model. Vertical lines indicate observed median dates across France in 2019, of all phenological stages of interest. Vertical line labels represent winter wheat tillering (BBCH29), stem elongation (BBCH31), heading (BBCH51), development of fruits (BBCH75), and senescence (BBCH99); and oilseed rape stem elongation (BBCH31), flowering (BBCH65), and development of fruits (BBCH73). These dates were obtained from the TEMPO dataset, except for those of winter wheat stem elongation (BBCH31) and senescence (BBCH99), which were obtained from the Céré’Obs dataset. Panel (A) details the amplitude definition for each side of the curve and the calibrated threshold value for each phenometric associated with each phenological stage of interest. Curve sides are relative to the maximum of the growing season. In panel (B), calibrated threshold values are shown for phenometrics associated with stem elongation (BBCH31) and development of fruits (BBCH73), which were obtained using the same amplitude definition process detailed in panel (A). For flowering (BBCH65), the associated phenometric (NDVI<sub>local_min</sub>) is shown within the corresponding temporal window (gray band). A local minimum occurs when the first derivative (dashed line) is zero at time <span class="html-italic">t</span>, negative at <span class="html-italic">t</span> − 1, and positive at <span class="html-italic">t</span> + 1.</p>
Full article ">Figure 4
<p>PROBA-V-based phenomaps of winter wheat. Each column represents a phenometric associated with a phenological stage available in the Céré’Obs database, i.e., SOS<sub>54–60</sub> with stem elongation (BBCH31), SOS<sub>95–98</sub> with heading (BBCH51) and EOS<sub>10–15</sub> with senescence (BBCH99), respectively. Each row represents a year in our study period. The color palette represents the day of the year (DoY) on which the phenometric was detected.</p>
Full article ">Figure 5
<p>Regional comparison of phenometrics (from PROBA-V) and phenological stages (from Céré’Obs) for winter wheat. (<b>A</b>) SOS<sub>54–60</sub> vs. stem elongation, (<b>B</b>) SOS<sub>95–98</sub> vs. Heading, and (<b>C</b>) EOS<sub>10–15</sub> vs. senescence. Each point represents a region median date and its color the year of interest. Median dates are expressed in the day of the year (DoY).</p>
Full article ">Figure 6
<p>Intra-regional comparison of phenometrics and associated phenological stages in terms of phenological progression in 2019. Phenometrics progression (in % of area) was obtained from the PROBA-V-based phenomaps (dashed lines), while phenological stage progression was obtained from Céré’Obs (solid lines). Each column represents a specific winter wheat phenometric and its associated phenological stage: (<b>A</b>–<b>C</b>) SOS<sub>54–60</sub> vs. stem elongation and (<b>D</b>–<b>F</b>) SOS<sub>95–98</sub> vs. heading, while each row represents a region of interest. The curves represent the values fitted with a logistic function.</p>
Full article ">Figure 7
<p>Municipal comparison of phenometrics and associated phenological stages in terms of median dates. Phenometrics dates were obtained from both PROBA-V- and Sentinel-2-based phenomaps, while phenological stages were obtained from the TEMPO dataset. Each column represents a specific winter wheat phenometric compared to associated phenological stage: (<b>A</b>–<b>D</b>) SOS<sub>42–53</sub> vs. tillering, (<b>E</b>–<b>H</b>) SOS<sub>95–98</sub> vs. Heading, and (<b>I</b>–<b>L</b>) EOS<sub>60–89</sub> vs. development of fruits. First row shows the locations where ground observations of each stage were made. Second row shows scatter-plots from all these observed municipalities, while third and fourth rows show scatter-plots from only the observed municipalities within the inter-comparison sites (dashed boxes), where we were able to extract phenometrics from Sentinel-2: third row represents estimated dates from PROBA-V, while fourth row represents those from Sentinel-2. In these last two rows, assessed municipalities were identical for both sensors. Scatter-plots and maps share the same color code representing the year of interest. Median dates are expressed in the day of the year (DoY).</p>
Full article ">Figure 8
<p>Municipal comparison of phenometrics and associated phenological stages in terms of median dates. Phenometrics dates were obtained from both PROBA-V and Sentinel-2-based phenomaps, while phenological stages came from the TEMPO dataset. Each column represents a specific oilseed rape phenometric compared to associated phenological stage: (<b>A</b>–<b>D</b>) SOS<sub>30–45</sub> vs. stem elongation, (<b>E</b>–<b>H</b>) NDVI<sub>local_min</sub> vs. Flowering, and (<b>I</b>–<b>L</b>) EOS<sub>97–99</sub> vs. development of fruits. First row shows the locations where ground observations of each stage were made. Second row shows scatter-plots from all these observed municipalities, while third and fourth rows show scatter-plots from only the observed municipalities within the inter-comparison sites (dashed boxes), where we were able to extract phenometrics from Sentinel-2: third row represents estimated dates from PROBA-V, while fourth row represents those from Sentinel-2. In these last two rows, assessed municipalities were identical for both sensors. Scatter-plots and maps share the same color code representing the year of interest. Median dates are expressed in the day of the year (DoY).</p>
Full article ">Figure 9
<p>Phenological detection rate from both PROBA-V and Sentinel-2 data. Each row represents a specific crop, and each column a specific phenometric: (<b>A</b>–<b>C</b>) winter wheat SOS<sub>42–53</sub>, SOS<sub>95–98</sub> and EOS<sub>60–89</sub>, (<b>D</b>–<b>F</b>) oilseed rape SOS<sub>30–45</sub>, NDVI<sub>local_min</sub>, and EOS<sub>97–99</sub>. The assessment was conducted based on all municipalities in the intercomparison sites, not only on those in which ground observations were carried out. On average between 2016 and 2020, 1943 municipalities were assessed for winter wheat and 1293 for oilseed rape.</p>
Full article ">
19 pages, 2495 KiB  
Article
Spatial–Temporal Differentiation and Driving Factors of Vegetation Landscape Pattern in Beijing–Tianjin–Hebei Region Based on the ESTARFM Model
by Yilin Wang, Ao Zhang, Xintong Gao, Wei Zhang, Xiaohong Wang and Linlin Jiao
Sustainability 2024, 16(23), 10498; https://doi.org/10.3390/su162310498 - 29 Nov 2024
Viewed by 445
Abstract
Urbanization and industrialization have led to obvious changes in the ecological environment and landscape pattern in the Beijing–Tianjin–Hebei region. Therefore, it is crucial to clarify the spatial–temporal changes in vegetation cover and its landscape pattern and conduct its analysis with the driving factors [...] Read more.
Urbanization and industrialization have led to obvious changes in the ecological environment and landscape pattern in the Beijing–Tianjin–Hebei region. Therefore, it is crucial to clarify the spatial–temporal changes in vegetation cover and its landscape pattern and conduct its analysis with the driving factors for ecological preservation in the Beijing–Tianjin–Hebei region. This study combined AVHRR GIMMS NDVI and MODIS NDVI data based on the ESTARFM model to obtain a high spatial–temporal resolution for vegetation cover; it then analyzed the vegetation cover changes at the type and landscape scales using a landscape index and explored the driving factors of the landscape pattern through principal component analysis. The results show that (1) the vegetation is mainly of medium and higher coverage and is distributed in the northeast, the western part of the Taihang Mountains and the central plains in the study area. From 1985 to 2022, there was no statistically significant difference in the overall change in its coverage. (2) From 1985 to 2022, at the landscape level, the vegetation cover landscape exhibited the following characteristics: increased fragmentation, an increase in the complexity of the landscape shape, a decrease in connectivity, a discrete landscape and a decrease in species diversity. At the type level, the medium vegetation demonstrated the most significant degree of fragmentation. The high-vegetation-cover areas exhibited a more concentrated distribution. Additionally, the low, lower and higher vegetation types displayed an increase in complexity, shape, discreteness and heterogeneity within the landscape. (3) Meanwhile, the principal component analysis showed that the changes in the landscape pattern of vegetation cover were mainly the result of the combined effects of climatic and anthropogenic factors in the Beijing–Tianjin–Hebei region. The human factor played the dominant role; this was followed by larger contributions from climatic factors. In addition to offering pertinent scientific insights for the maximization of the ecological environment and the fostering of regional ecological and sustainable development in the Beijing–Tianjin–Hebei region, the aforementioned analysis and research could serve as the foundation for the sustainable management and planning of vegetation cover. Full article
Show Figures

Figure 1

Figure 1
<p>Overview of the study area.</p>
Full article ">Figure 2
<p>Scatter plot of GIMMS NDVI3g data after real MODIS NDVI fusion.</p>
Full article ">Figure 3
<p>Comparison of the images before and after fusion in 2000.</p>
Full article ">Figure 4
<p>Comparison of the images before and after fusion in 2005.</p>
Full article ">Figure 5
<p>FVC in 1985–2022.</p>
Full article ">Figure 6
<p>Sliding <span class="html-italic">t</span>−test for mutation point detection.</p>
Full article ">Figure 7
<p>Vegetation coverage distribution map.</p>
Full article ">Figure 8
<p>Vegetation cover change map.</p>
Full article ">Figure 9
<p>Trend of vegetation cover from 1985 to 2022.</p>
Full article ">Figure 10
<p>Landscape index changes during 1985–2022 at the Beijing–Tianjin–Hebei type level.</p>
Full article ">Figure 11
<p>Changes in socioeconomic and demographic factors in the study area from 1985 to 2022.</p>
Full article ">Figure 12
<p>Changes in climate factors.</p>
Full article ">
18 pages, 20080 KiB  
Article
Driving Factors for Vegetation NDVI Changes in a Cold Temperate Zone: Climate, Topography, and Land Use
by Dandan Zhao, Weijia Hu, Jianmiao Wang and Jiping Liu
Forests 2024, 15(12), 2098; https://doi.org/10.3390/f15122098 - 27 Nov 2024
Viewed by 663
Abstract
Exploring the spatio-temporal evolution and driving mechanism of the NDVI (Normalized Difference Vegetation Index) is important in order to understand the operating forces of the ecosystem and the response process of environmental change. We analyzed spatio-temporal vegetation changes by using the trend analysis [...] Read more.
Exploring the spatio-temporal evolution and driving mechanism of the NDVI (Normalized Difference Vegetation Index) is important in order to understand the operating forces of the ecosystem and the response process of environmental change. We analyzed spatio-temporal vegetation changes by using the trend analysis method during 2001–2020 based on the MODIS NDVI, the meteorological data, the DEM (Digital Elevation Model) and land use types data. We quantitatively revealed the influence degree and mechanism of each detection factor and their interaction on the spatial differentiation of vegetation by using the geographical detector model. Results showed that the vegetation NDVI showed an increasing trend with an increasing rate of 0.021/10 a during 2001–2020 and mainly distributed in the northwest and southwest of the Greater Khingan Mountains. The explanatory power values of each driving factor are as follows: land use (0.384) > elevation (0.193) > slope (0.159) > annual precipitation (0.104) > aspect (0.069) > average annual temperature (0.056). The explanatory power of interaction between driving factors were relatively high, as follows: Land use ∩ Aspect (0.490) > Land use ∩ Slope (0.471) > Land use ∩ Annual precipitation (0.460) > Land use ∩ elevation (0.443) > Land use ∩ Annual temperature (0.421) > Aspect ∩ elevation (0.408). Our research was of great significance for understanding the growth law of vegetation, protecting the ecological environment, and sustainable development in cold temperate zones. Full article
(This article belongs to the Section Forest Ecology and Management)
Show Figures

Figure 1

Figure 1
<p>The Greater Khingan Mountains area of China.</p>
Full article ">Figure 2
<p>(<b>a</b>) NDVI variation trends (<b>b</b>) Annual average NDVI (<b>c</b>) Slope variation trends of NDVI during 2001–2020 in the Greater Khingan Mountains.</p>
Full article ">Figure 3
<p>Annual NDVI changes of different land use in the Greater Khingan Mountains.</p>
Full article ">Figure 4
<p>Spatial distribution of land use transformation in the Greater Khingan Mountains during 2001–2020.</p>
Full article ">Figure 5
<p>Interannual changes of mean NDVI and average temperature (<b>a</b>) and interannual changes of mean NDVI and annual precipitation (<b>b</b>) during 2001–2020.</p>
Full article ">Figure 6
<p>Spatial distribution of correlation between annual NDVI and average temperature (<b>a</b>) Spatial distribution of correlation between annual NDVI and annual precipitation (<b>b</b>).</p>
Full article ">Figure 7
<p>Variation trend of NDVI with elevation (<b>a</b>), slope (<b>b</b>), and aspect (<b>c</b>) in the Greater Khingan Mountains during 2001–2020.</p>
Full article ">Figure 8
<p>Explanatory power of driving factors of the NDVI spatial differentiation in the Greater Khingan Mountains.</p>
Full article ">Figure 9
<p>Explanatory power of interactive detection of driving factors in the Greater Khingan Mountains.</p>
Full article ">
26 pages, 13283 KiB  
Article
Reconstruction of 30 m Land Cover in the Qilian Mountains from 1980 to 1990 Based on Super-Resolution Generative Adversarial Networks
by Xiaoya Wang, Bo Zhong, Kai Ao, Bailin Du, Longfei Hu, He Cai, Yang Qiao, Junjun Wu, Aixia Yang, Shanlong Wu and Qinhuo Liu
Remote Sens. 2024, 16(22), 4252; https://doi.org/10.3390/rs16224252 - 14 Nov 2024
Viewed by 678
Abstract
Long time series of annual land cover with fine spatio-temporal resolutions play a crucial role in studying environmental climate change, biophysical modeling, carbon cycling models, and land management. Despite a strong consistency exhibited by several publicly available medium to fine resolution global land [...] Read more.
Long time series of annual land cover with fine spatio-temporal resolutions play a crucial role in studying environmental climate change, biophysical modeling, carbon cycling models, and land management. Despite a strong consistency exhibited by several publicly available medium to fine resolution global land cover datasets, significant discrepancies exist at the regional scale; moreover, only every 5/10 year land cover were available. Consequently, high-quality annual land cover datasets before 2000 are unavailable in China. In this study, we proposed a deep learning-based method by integrating multiple remote sensing data from different platforms with historical high spatial resolution land cover datasets (CNLUCC) to derive the 30 m annual land cover maps from 1980 to 1990 for Qilian Mountain. First, the super-resolution generative adversarial network models for upscaling the 5.5 km AVHRR NDVI to 250 m were established by employing the AVHRR and MODIS NDVI data with the same year as input, and the early time series AVHRR NDVI data were subsequently upscaled to 250 m through the above models. Second, the breaks for the additive seasonal and trend (BFAST) change detection algorithm was applied to the upscaled time series NDVI data to detect the change time of different land cover types. Third, the CNLUCC data in 1980 and 1990 were updated to annual land cover datasets from 1980 to 1990 and the annual mapping results provided insights into the dynamic processes of urbanization, deforestation, water bodies, and farmland from 1980 to 1990. Finally, comprehensive analysis and validation were carried out for evaluation and an overall accuracy of 77.26% for the land cover product in 1986 was achieved. Full article
(This article belongs to the Special Issue Surface Radiative Transfer: Modeling, Inversion, and Applications)
Show Figures

Figure 1

Figure 1
<p>The location of major river basins in the Qilian Mountains (<b>left</b>) and the visualization of geographical characteristics including color composite from remote sensing image (<b>right</b>).</p>
Full article ">Figure 2
<p>Multiple satellite images’ time range and availability.</p>
Full article ">Figure 3
<p>Preprocessing rendering. (<b>a</b>) Quality check and cloud mask; (<b>b</b>) fill in missing values by temporal filter; (<b>c</b>) fill in missing values by spatial filter.</p>
Full article ">Figure 4
<p>Workflow of the annual land cover mapping process.</p>
Full article ">Figure 5
<p>The process of making SR model training dataset (LR-HR image pairs). (<b>a</b>) Creating a square grid; (<b>b</b>) creating centroids of the grids.</p>
Full article ">Figure 6
<p>The examples of training dataset (LR-HR image pairs). (<b>a</b>) Gobi; (<b>b</b>) Lake; (<b>c</b>) Forest; (<b>d</b>) River.</p>
Full article ">Figure 7
<p>The architecture of generator and discriminator network.</p>
Full article ">Figure 8
<p>Minibatch statistic layer.</p>
Full article ">Figure 9
<p>Example of change pixels for 3 × 3 grids (C: Cropland, F: Forest, G: Grassland, W: Water body, B: Built-up land, and U: Unused land).</p>
Full article ">Figure 10
<p>Breakpoint test (example for one pixel).</p>
Full article ">Figure 11
<p>The location of Jiayuguan City and Suzhou District (indicated by the red line) and the Heihe River Basin (indicated by the blue line).</p>
Full article ">Figure 12
<p>Annual land cover maps from 1980 to 1990.</p>
Full article ">Figure 12 Cont.
<p>Annual land cover maps from 1980 to 1990.</p>
Full article ">Figure 13
<p>The expansion of built-up in the study area.</p>
Full article ">Figure 14
<p>June 1986 NDVI super-resolution results (The left column (<b>a</b>,<b>c</b>,<b>e</b>) displays the results of SR NDVI data in the Qilian Mountains from original resolution to 1 km and further to 250 m; the right column (<b>b</b>,<b>d</b>,<b>f</b>) shows enlarged details of the black box area and (<b>g</b>,<b>h</b>) show 30 m SR data and resampled data, respectively).</p>
Full article ">Figure 14 Cont.
<p>June 1986 NDVI super-resolution results (The left column (<b>a</b>,<b>c</b>,<b>e</b>) displays the results of SR NDVI data in the Qilian Mountains from original resolution to 1 km and further to 250 m; the right column (<b>b</b>,<b>d</b>,<b>f</b>) shows enlarged details of the black box area and (<b>g</b>,<b>h</b>) show 30 m SR data and resampled data, respectively).</p>
Full article ">Figure 15
<p>The spatial distribution of validation samples in 1986.</p>
Full article ">Figure 16
<p>Confusion matrix for our land cover map (<b>left</b>) and CLUD-A (<b>right</b>) in 1986.</p>
Full article ">Figure 17
<p>Spatial distribution of 1000 random sample points.</p>
Full article ">Figure 18
<p>Area change in different land cover classes from 1995 to 2005. The fluctuations for each land cover type are enlarged using the different ranges of y-axis.</p>
Full article ">Figure 19
<p>Visual comparison of the detected change years (the first column) with the images from Google Earth and Landsat (the second and third column). The highlighted areas with blue shapes were the change regions. (<b>a</b>) Change from unused land to water body, (<b>b</b>) the conversion from unused land to urban area, (<b>c</b>) change from unused land to cropland.</p>
Full article ">Figure 20
<p>Deviation of detected change years from the change samples.</p>
Full article ">
18 pages, 11426 KiB  
Article
Spring Phenological Responses of Diverse Vegetation Types to Extreme Climatic Events in Mongolia
by Qier Mu, Sainbuyan Bayarsaikhan, Gang Bao, Battsengel Vandansambuu, Siqin Tong, Byambakhuu Gantumur, Byambabayar Ganbold and Yuhai Bao
Sustainability 2024, 16(22), 9931; https://doi.org/10.3390/su16229931 - 14 Nov 2024
Viewed by 468
Abstract
The increasing frequency of extreme climate events may significantly alter the species composition, structure, and functionality of ecosystems, thereby diminishing their stability and resilience. This study draws on temperature and precipitation data from 53 meteorological stations across Mongolia, covering the period from 1983 [...] Read more.
The increasing frequency of extreme climate events may significantly alter the species composition, structure, and functionality of ecosystems, thereby diminishing their stability and resilience. This study draws on temperature and precipitation data from 53 meteorological stations across Mongolia, covering the period from 1983 to 2016, along with MODIS normalized difference vegetation index (NDVI) data from 2001 to 2016. The climate anomaly method and the curvature method of cumulative NDVI logistic curves were employed to identify years of extreme climate events and to extract the start of the growing season (SOS) in Mongolia. Furthermore, the study assessed the impact of extreme climate events on the SOS across different vegetation types and evaluated the sensitivity of the SOS to extreme climate indices. The study results show that, compared to the multi-year average green-up period from 2001 to 2016, extreme climate events significantly impact the SOS. Extreme dryness advanced the SOS by 6.9 days, extreme wetness by 2.5 days, and extreme warmth by 13.2 days, while extreme cold delayed the SOS by 1.2 days. During extreme drought events, the sensitivity of SOS to TN90p (warm nights) was the highest; in extremely wet years, the sensitivity of SOS to TX10p (cool days) was the strongest; in extreme warm events, SOS was most sensitive to TX90p (warm days); and during extreme cold events, SOS was most sensitive to TNx (maximum night temperature). Overall, the SOS was most sensitive to extreme temperature indices during extreme climate events, with a predominantly negative sensitivity. The response and sensitivity of SOS to extreme climate events varied across different vegetation types. This is crucial for understanding the dynamic changes of ecosystems and assessing potential ecological risks. Full article
Show Figures

Figure 1

Figure 1
<p>Location of Mongolia and spatial distribution of meteorological stations (<b>a</b>) and vegetation types (<b>b</b>).</p>
Full article ">Figure 2
<p>Plotted mean climatological departures of Mongolia from 1983 to 2016 for (<b>a</b>) maximum temperature, (<b>b</b>) mean temperature, (<b>c</b>) minimum temperature, and (<b>d</b>) precipitation. The pink, red, blue, and green lines correspond to extremely dry, warm, cold, and wet years, respectively. The red and black circles represent extreme and normal values of climate observations, respectively.</p>
Full article ">Figure 3
<p>Spatial distribution of SOS anomalies on extreme climate events (<b>a</b>–<b>d</b>). Plot of the relative frequency of SOS anomalies for extreme climate events (<b>e</b>–<b>h</b>).</p>
Full article ">Figure 4
<p>(<b>a</b>) The changes in the average SOS anomalies (black point) across the entire study area in Mongolia. (<b>b</b>) SOS anomalies among the four vegetation types during extreme climate events.</p>
Full article ">Figure 5
<p>The importance of 10 extreme climate indices to the diverse vegetation types of spring phenology is shown in each year of extreme climate events: (<b>a</b>–<b>e</b>) extremely dry; (<b>f</b>–<b>j</b>) extremely warm; (<b>k</b>–<b>o</b>) extremely cold; and (<b>p</b>–<b>t</b>) extremely wet years.</p>
Full article ">Figure 5 Cont.
<p>The importance of 10 extreme climate indices to the diverse vegetation types of spring phenology is shown in each year of extreme climate events: (<b>a</b>–<b>e</b>) extremely dry; (<b>f</b>–<b>j</b>) extremely warm; (<b>k</b>–<b>o</b>) extremely cold; and (<b>p</b>–<b>t</b>) extremely wet years.</p>
Full article ">Figure 6
<p>The sensitivity of four extreme climate indices to the entire study area and diverse vegetation types of spring phenology is shown in each year of extreme climate events: (<b>a</b>) extreme dry, (<b>b</b>) extreme warm, (<b>c</b>) extreme cold, and (<b>d</b>) extreme wet years. ** represents 0.01 significance level.</p>
Full article ">Figure 7
<p>The correlation of the SOS with precipitation and temperature across the entire study area during extreme climate events (<b>a</b>–<b>d</b>). ** represents 0.01 significance level.</p>
Full article ">
22 pages, 4548 KiB  
Article
MODIS Evapotranspiration Forecasting Using ARIMA and ANN Approach at a Water-Stressed Irrigation Scheme in South Africa
by Mbulelo Phesa, Nkanyiso Mbatha and Akinola Ikudayisi
Hydrology 2024, 11(10), 176; https://doi.org/10.3390/hydrology11100176 - 21 Oct 2024
Viewed by 911
Abstract
The forecasting of evapotranspiration (ET) in some water-stressed regions remains a major challenge due to the lack of reliable and sufficient historical datasets. For efficient water balance, ET remains the major component and its proper forecasting and quantifying is of the utmost importance. [...] Read more.
The forecasting of evapotranspiration (ET) in some water-stressed regions remains a major challenge due to the lack of reliable and sufficient historical datasets. For efficient water balance, ET remains the major component and its proper forecasting and quantifying is of the utmost importance. This study utilises the 18-year (2001 to 2018) MODIS ET obtained from a drought-affected irrigation scheme in the Eastern Cape Province of South Africa. This study conducts a teleconnection evaluation between the satellite-derived evapotranspiration (ET) time series and other related remotely sensed parameters such as the Normalised Difference Vegetation Index (NDVI), Normalised Difference Water Index (NDWI), Normalised Difference Drought Index (NDDI), and precipitation (P). This comparative analysis was performed by adopting the Mann–Kendall (MK) test, Sequential Mann–Kendall (SQ-MK) test, and Multiple Linear Regression methods. Additionally, the ET detailed time-series analysis with the Keiskamma River streamflow (SF) and monthly volumes of the Sandile Dam, which are water supply sources close to the study area, was performed using the Wavelet Analysis, Breaks for Additive Seasonal and Trend (BFAST), Theil–Sen statistic, and Correlation statistics. The MODIS-obtained ET was then forecasted using the Autoregressive Integrated Moving Average (ARIMA) and Artificial Neural Networks (ANNs) for a period of 5 years and four modelling performance evaluations such as the Root Mean Square Error (RMSE), Mean Absolute Percentage Error (MAPE), Mean Absolute Error (MAE), and the Pearson Correlation Coefficient (R) were used to evaluate the model performances. The results of this study proved that ET could be forecasted using these two time-series modeling tools; however, the ARIMA modelling technique achieved lesser values according to the four statistical modelling techniques employed with the RMSE for the ARIMA = 37.58, over the ANN = 44.18; the MAE for the ARIMA = 32.37, over the ANN = 35.88; the MAPE for the ARIMA = 17.26, over the ANN = 24.26; and for the R ARIMA = 0.94 with the ANN = 0.86. These results are interesting as they give hope to water managers at the irrigation scheme and equally serve as a tool to effectively manage the irrigation scheme. Full article
Show Figures

Figure 1

Figure 1
<p>Keiskammahoek Irrigation Scheme map marked with gray colour.</p>
Full article ">Figure 2
<p>Shows the BAST graphs for ET (<b>a</b>), NDWI (<b>b</b>), P (<b>c</b>), NDVI (<b>d</b>), MV (<b>e</b>), and SF (<b>f</b>).</p>
Full article ">Figure 3
<p>ET (<b>a</b>) and P (<b>b</b>) normalised wavelet power spectra of yearly mean variability during 2001 to 2018 at KIS.</p>
Full article ">Figure 4
<p>Wavelet coherence between ET and precipitation for period of 2001 to 2018 with phase relationship depicted by the arrows.</p>
Full article ">Figure 5
<p>Pearson Correlation Coefficient for NDDI, ET, NDVI, NDWI, SF, and P.</p>
Full article ">Figure 6
<p>Long-term trend of monthly SF (<b>a</b>), MV (<b>b</b>), ET (<b>c</b>), and P (<b>d</b>).</p>
Full article ">Figure 7
<p>SQ-MK test for ET and P data for trend period in years 2001–2018. The sequent statistics values of progressive (Prog) u(t) (red solid line) and retrograde u’(t) (solid black line) attained by the Sequential Man-Kendal (SQ-) test for ET data (<b>a</b>) and Precipitation data (<b>b</b>) for 18 years period.</p>
Full article ">Figure 8
<p>Autocorrelation Function (ACF) (<b>a</b>) and the histogram of residuals (<b>b</b>) for KIS for the 18 years (2001 to 2018) showing a Keiskammahoek Irrigation Scheme best fitted model for ET data series.</p>
Full article ">Figure 9
<p>(<b>a</b>) ET dataset from 2001 to 2018 and 3-year forecast from 2015 to 2028, and (<b>b</b>) ARIMA scatter plot with correlation between forecasted ET and obtained ET.</p>
Full article ">Figure 10
<p>(<b>a</b>) ET dataset from 2001 to 2015 and the remaining forecasted from 2015 to 2018, and (<b>b</b>) ANN scatter plot with correlation between forecasted ET and obtained ET.</p>
Full article ">
24 pages, 5753 KiB  
Article
Phenological Monitoring of Irrigated Sugarcane Using Google Earth Engine, Time Series, and TIMESAT in the Brazilian Semi-Arid
by Diego Rosyur Castro Manrique, Pabrício Marcos Oliveira Lopes, Cristina Rodrigues Nascimento, Eberson Pessoa Ribeiro and Anderson Santos da Silva
AgriEngineering 2024, 6(4), 3799-3822; https://doi.org/10.3390/agriengineering6040217 - 18 Oct 2024
Viewed by 763
Abstract
Monitoring sugarcane phenology is essential since the globalized market requires reliable information on the quantity of raw materials for the industrial production of sugar and alcohol. In this context, the general objective of this study was to evaluate the phenological seasonality of the [...] Read more.
Monitoring sugarcane phenology is essential since the globalized market requires reliable information on the quantity of raw materials for the industrial production of sugar and alcohol. In this context, the general objective of this study was to evaluate the phenological seasonality of the sugarcane varieties SP 79-1011 and VAP 90-212 observed from the NDVI time series over 19 years (2001–2020) from global databases. In addition, this research had the following specific objectives: (i) to estimate phenological parameters (Start of Season (SOS), End of Season (EOS), Length of Season (LOS), and Peak of Season (POS)) using TIMESAT software in version 3.3 applied to the NDVI time series over 19 years; (ii) to characterize the land use and land cover obtained from the MapBiomas project; (iii) to analyze rainfall variability; and (iv) to validate the sugarcane harvest date (SP 79-1011). This study was carried out in sugarcane growing areas in Juazeiro, Bahia, Brazil. The results showed that the NDVI time series did not follow the rainfall in the region. The sugarcane areas advanced over the savanna formation (Caatinga), reducing them to remnants along the irrigation channels. The comparison of the observed harvest dates of the SP 79-1011 variety to the values estimated with the TIMESAT software showed an excellent fit of 0.99. The mean absolute error in estimating the sugarcane harvest date was approximately ten days, with a performance index of 0.99 and a correlation coefficient of 0.99, significant at a 5% confidence level. The TIMESAT software was able to estimate the phenological parameters of sugarcane using MODIS sensor images processed on the Google Earth Engine platform during the evaluated period (2001 to 2020). Full article
Show Figures

Figure 1

Figure 1
<p>Map of the sugarcane with the physical boundaries in RGB (red, green and blue) color composite Landsat-8 and the location under study.</p>
Full article ">Figure 2
<p>Graphical abstract steps for obtaining the phenological metrics. where: SOS = Start of Season, EOS = End Of Season, LOS = Length of the Season and POS = Peak of Season. Source: Adapted of Rodigheri et al. [<a href="#B40-agriengineering-06-00217" class="html-bibr">40</a>].</p>
Full article ">Figure 3
<p>TIMESAT software modules for processing NDVI time series in the TIMESAT software module.</p>
Full article ">Figure 4
<p>Application of the Savitsky–Golay filter (red curve) for a time series of NDVI scaled (black curve) as a function of time (days) to estimate phenological parameters: points (a) and (b) mark, respectively, start and end of the season, points (c) and (d) give the 80% levels, (e) displays the point with the maximum value, (f) displays the seasonal amplitude, (g) the seasonal length, and (h) and (i) are integrals showing the cumulative effect of vegetation during the season. Source: Jönsson and Eklundh [<a href="#B48-agriengineering-06-00217" class="html-bibr">48</a>].</p>
Full article ">Figure 5
<p>Meteorological data from the Meteorology Laboratory (LabMet) automatic weather station for the period 2008 to 2012, Juazeiro, Bahia, Brazil.</p>
Full article ">Figure 6
<p>Classification of the land use and land cover of the watershed using MapBiomas in its Collection 6 (2006–2012).</p>
Full article ">Figure 7
<p>Temporal distribution of the area cultivated with sugarcane in the watershed from 2001 to 2020 in Juazeiro, Bahia, Brazil.</p>
Full article ">Figure 8
<p>MODIS NDVI (2001–2020) for the total area and rainfall of Labmet Juazeiro (2008–2020) time series.</p>
Full article ">Figure 9
<p>NDVI time series and Savitsky–Golay filter for the sugarcane total area. The dots represent the start (in blue) and the end (in yellow) of the sugarcane phenological cycles.</p>
Full article ">Figure 10
<p>Sugarcane agricultural calendar in the test area from the SP 79-1011 and VAP 90-212. In blue, months referring to the phenological phases of sugarcane.</p>
Full article ">Figure 11
<p>Comparison of variety SP 79-1011 harvest dates observed with estimated values with TIMESAT software for the test area between 2006 to 2012.</p>
Full article ">
19 pages, 11653 KiB  
Article
Influence of Vegetation Phenology on the Temporal Effect of Crop Fractional Vegetation Cover Derived from Moderate-Resolution Imaging Spectroradiometer Nadir Bidirectional Reflectance Distribution Function–Adjusted Reflectance
by Yinghao Lin, Tingshun Fan, Dong Wang, Kun Cai, Yang Liu, Yuye Wang, Tao Yu and Nianxu Xu
Agriculture 2024, 14(10), 1759; https://doi.org/10.3390/agriculture14101759 - 5 Oct 2024
Viewed by 731
Abstract
Moderate-Resolution Imaging Spectroradiometer (MODIS) Nadir Bidirectional Reflectance Distribution Function (BRDF)-Adjusted Reflectance (NBAR) products are being increasingly used for the quantitative remote sensing of vegetation. However, the assumption underlying the MODIS NBAR product’s inversion model—that surface anisotropy remains unchanged over the 16-day retrieval period—may [...] Read more.
Moderate-Resolution Imaging Spectroradiometer (MODIS) Nadir Bidirectional Reflectance Distribution Function (BRDF)-Adjusted Reflectance (NBAR) products are being increasingly used for the quantitative remote sensing of vegetation. However, the assumption underlying the MODIS NBAR product’s inversion model—that surface anisotropy remains unchanged over the 16-day retrieval period—may be unreliable, especially since the canopy structure of vegetation undergoes stark changes at the start of season (SOS) and the end of season (EOS). Therefore, to investigate the MODIS NBAR product’s temporal effect on the quantitative remote sensing of crops at different stages of the growing seasons, this study selected typical phenological parameters, namely SOS, EOS, and the intervening stable growth of season (SGOS). The PROBA-V bioGEOphysical product Version 3 (GEOV3) Fractional Vegetation Cover (FVC) served as verification data, and the Pearson correlation coefficient (PCC) was used to compare and analyze the retrieval accuracy of FVC derived from the MODIS NBAR product and MODIS Surface Reflectance product. The Anisotropic Flat Index (AFX) was further employed to explore the influence of vegetation type and mixed pixel distribution characteristics on the BRDF shape under different stages of the growing seasons and different FVC; that was then combined with an NDVI spatial distribution map to assess the feasibility of using the reflectance of other characteristic directions besides NBAR for FVC correction. The results revealed the following: (1) Generally, at the SOSs and EOSs, the differences in PCCs before vs. after the NBAR correction mainly ranged from 0 to 0.1. This implies that the accuracy of FVC derived from MODIS NBAR is lower than that derived from MODIS Surface Reflectance. Conversely, during the SGOSs, the differences in PCCs before vs. after the NBAR correction ranged between –0.2 and 0, suggesting the accuracy of FVC derived from MODIS NBAR surpasses that derived from MODIS Surface Reflectance. (2) As vegetation phenology shifts, the ensuing differences in NDVI patterning and AFX can offer auxiliary information for enhanced vegetation classification and interpretation of mixed pixel distribution characteristics, which, when combined with NDVI at characteristic directional reflectance, could enable the accurate retrieval of FVC. Our results provide data support for the BRDF correction timescale effect of various stages of the growing seasons, highlighting the potential importance of considering how they differentially influence the temporal effect of NBAR corrections prior to monitoring vegetation when using the MODIS NBAR product. Full article
Show Figures

Figure 1

Figure 1
<p>Spatial extent of the Wancheng District study area (in Henan Province, China). (<b>a</b>) Map of land cover types showing the location of sampling points across the study area. This map came from MCD12Q1 (v061). (<b>b</b>–<b>d</b>) True-color images of the three mixed pixels, obtained from Sentinel-2. The distribution characteristics are as follows: crops above with buildings below (<b>b</b>); crops below with buildings above (<b>c</b>); and buildings in the upper-left corner, crops in the remainder (<b>d</b>).</p>
Full article ">Figure 2
<p>Monthly average temperature and monthly total precipitation in the study area, from 2017 to 2021.</p>
Full article ">Figure 3
<p>Data processing flow chart. The green rectangles from top to the bottom represent three steps: crop phenological parameters extraction with TIMESAT; Fractional Vegetation Cover (FVC) derived from MOD09GA and MCD43A4; and accuracy evaluation, respectively. Blue solid rectangles refer to a used product or derived results, while blue dashed rectangles refer to the software or model used in this study. NDVI<sub>MOD09GA</sub>: NDVI derived from MOD09GA, NDVI<sub>MCD43A4</sub>: NDVI derived from MCD43A4, FVC<sub>MOD09GA</sub>: FVC derived from MOD09GA, FVC<sub>MCD43A4</sub>: FVC derived from MCD43A4. PCC<sub>MOD09GA</sub>: Pearson correlation coefficient (PCC) calculated for FVC<sub>MOD09GA</sub> and GEOV3 FVC, PCC<sub>MCD43A4</sub>: PCC calculated for FVC<sub>MCD43A4</sub> and GEOV3 FVC.</p>
Full article ">Figure 4
<p>NDVI and EVI time series fitted curves and phenological parameters of crops. SOS: start of season; EOS: end of season; SGOS: stable growth of season.</p>
Full article ">Figure 5
<p>Spatial distribution of Fractional Vegetation Cover (FVC) derived from MOD09GA and MCD43A4, and the difference images of FVC. FVC<sub>MOD09GA</sub>: FVC derived from MOD09GA, FVC<sub>MCD43A4</sub>: FVC derived from MCD43A4. (<b>a</b>–<b>c</b>) FVC derived from MOD09GA, MCD43A4, and the difference between FVC<sub>MOD09GA</sub> and FVC<sub>MCD43A4</sub> on 15 November 2020, respectively; (<b>d</b>–<b>f</b>) FVC derived from MOD09GA, MCD43A4, and the difference between FVC<sub>MOD09GA</sub> and FVC<sub>MCD43A4</sub> on 10 February 2021, respectively; (<b>g</b>–<b>i</b>) FVC derived from MOD09GA, MCD43A4, and the difference between FVC<sub>MOD09GA</sub> and FVC<sub>MCD43A4</sub> on 30 September 2021, respectively.</p>
Full article ">Figure 6
<p>Pearson correlation coefficients (PCCs) of Fractional Vegetation Cover (FVC) derived before and after the NBAR correction with GEOV3 FVC at different stages of the growing seasons. FVC<sub>MOD09GA</sub>: FVC derived from MOD09GA. FVC<sub>MCD43A4</sub>: FVC derived from MCD43A4. PCC<sub>MOD09GA</sub>: PCC calculated for FVC<sub>MOD09GA</sub> and GEOV3 FVC, PCC<sub>MCD43A4</sub>: PCC calculated for FVC<sub>MCD43A4</sub> and GEOV3 FVC. (<b>a</b>) PCC<sub>MOD09GA</sub> and PCC<sub>MCD43A4</sub> in 2018–2021; (<b>b</b>) Scatterplot of numerical differences between PCC<sub>MOD09GA</sub> and PCC<sub>MCD43A4</sub>. SOS: start of season; EOS: end of season; SGOS: stable growth of season.</p>
Full article ">Figure 7
<p>NDVI spatial distribution maps of crop pixel, savanna pixel, and grassland pixel in different stages of the growing seasons. (<b>a</b>–<b>d</b>) Crop. (<b>e</b>–<b>h</b>) Savanna. (<b>i</b>–<b>l</b>) Grassland. SZA: Solar Zenith Angle, FVC: Fractional Vegetation Cover, AFX_RED: Anisotropic Flat Index (AFX) in the red band, AFX_NIR: AFX in the near-infrared band.</p>
Full article ">Figure 8
<p>NDVI spatial distribution maps of mixed pixels in different stages of the growing seasons. (<b>a</b>–<b>d</b>) Crops above and buildings below. (<b>e</b>–<b>h</b>) Crops below and buildings above. (<b>i</b>–<b>l</b>) Buildings in the upper-left corner and crops in the remainder. SZA: Solar Zenith Angle, FVC: Fractional Vegetation Cover, AFX_RED: Anisotropic Flat Index (AFX) in the red band, AFX_NIR: AFX in the near-infrared band.</p>
Full article ">
19 pages, 6249 KiB  
Article
Carbon and Energy Balance in a Primary Amazonian Forest and Its Relationship with Remote Sensing Estimates
by Mailson P. Alves, Rommel B. C. da Silva, Cláudio M. Santos e Silva, Bergson G. Bezerra, Keila Rêgo Mendes, Larice A. Marinho, Melahel L. Barbosa, Hildo Giuseppe Garcia Caldas Nunes, José Guilherme Martins Dos Santos, Theomar Trindade de Araújo Tiburtino Neves, Raoni A. Santana, Lucas Vaz Peres, Alex Santos da Silva, Petia Oliveira, Victor Hugo Pereira Moutinho, Wilderclay B. Machado, Iolanda M. S. Reis, Marcos Cesar da Rocha Seruffo, Avner Brasileiro dos Santos Gaspar, Waldeir Pereira and Gabriel Brito-Costaadd Show full author list remove Hide full author list
Remote Sens. 2024, 16(19), 3606; https://doi.org/10.3390/rs16193606 - 27 Sep 2024
Viewed by 1045
Abstract
With few measurement sites and a great need to validate satellite data to characterize the exchange of energy and carbon fluxes in tropical forest areas, quantified by the Net Ecosystem Exchange (NEE) and associated with phenological measurements, there is an increasing need for [...] Read more.
With few measurement sites and a great need to validate satellite data to characterize the exchange of energy and carbon fluxes in tropical forest areas, quantified by the Net Ecosystem Exchange (NEE) and associated with phenological measurements, there is an increasing need for studies aimed at characterizing the Amazonian environment in its biosphere–atmosphere interaction, considering the accelerated deforestation in recent years. Using data from a flux measurement tower in the Caxiuanã-PA forest (2005–2008), climatic data, CO2 exchange estimated by eddy covariance, as well as Gross Primary Productivity (GPP) data and satellite vegetation indices (from MODIS), this work aimed to describe the site’s energy, climatic and carbon cycle flux patterns, correlating its gross primary productivity with satellite vegetation indices. The results found were: (1) marked seasonality of climatic variables and energy flows, with evapotranspiration and air temperature on the site following the annual march of solar radiation and precipitation; (2) energy fluxes in phase and dependent on available energy; (3) the site as a carbon sink (−569.7 ± 444.9 gC m−2 year−1), with intensity varying according to the site’s annual water availability; (4) low correlation between productivity data and vegetation indices, corroborating data in the literature on these variables in this type of ecosystem. The results show the importance of preserving this type of environment for the mitigation of global warming and the need to improve satellite estimates for this region. NDVI and EVI patterns follow radiative availability, as does LAI, but without direct capture related to GPP data, which correlates better with satellite data only in the months with the highest LAI. The results show the significant difference at a point measurement to a satellite interpolation, presenting how important preserving any type of environment is, even related to its size, for the global climate balance, and also the need to improve satellite estimates for smaller areas. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the micrometeorological tower in the Caxiuanã-PA forest, with a view of the base of the installed tower and the temperature and humidity sensor at the top of the tower.</p>
Full article ">Figure 2
<p>Data density distribution of monthly air temperature (°C) by precipitation intensity (mm/day) in Caxiuanã-PA.</p>
Full article ">Figure 3
<p>Annual wind patterns at the Caxiuanã site.</p>
Full article ">Figure 4
<p>Monthly averages of (<b>a</b>) air temperature (°C), (<b>b</b>) Evapotranspiration (mm), (<b>c</b>) photosynthetic photon flux density (μmol m<sup>−2</sup> s<sup>−1</sup>) and (<b>d</b>) precipitation (mm month<sup>−1</sup>) at Caxiuanã-PA site. The hatched area indicates the 95% confidence interval.</p>
Full article ">Figure 5
<p>Monthly (<b>a</b>) and hourly (<b>b</b>) averages of the energy balance fluxes: radiation balance (Rn), sensible heat flux (H) and latent heat flux (LE), all in W m<sup>−2</sup>. The hatched area indicates the 95% confidence interval.</p>
Full article ">Figure 6
<p>NEE response to PPFD in different seasons ((<b>a</b>) = dry; (<b>b</b>) = wet) and day periods.</p>
Full article ">Figure 7
<p>Average diurnal cycle of NEE at Caxiuanã, representing the CO<sub>2</sub> flux above the canopy, for the dry and rainy seasons. The hatched area indicates the 95% confidence interval.</p>
Full article ">Figure 8
<p>Monthly variation of NEE, GPP and RECO fluxes (<b>a</b>) and vegetation indices (<b>b</b>) at Caxiuanã. The hatched area indicates the 95% confidence interval.</p>
Full article ">Figure 9
<p>Pearson’s linear correlation between the GPP data measured in the tower and the vegetation indices: NDVI (<b>a</b>) and EVI (<b>b</b>). The hatched area indicates the 95% confidence interval.</p>
Full article ">Figure 10
<p>Comparison between measured gross primary production (Tower) and MODIS gross primary production (MODIS) in (<b>a</b>,<b>b</b>) Monthly GPP boxplot (gC m<sup>−2</sup> day<sup>−1</sup>) for eddy covariance observed data (Tower) and MODIS—derived data (Satellite). The points refer to the outliers of the data, while the central lines refer to the median of each month.</p>
Full article ">Figure 11
<p>Monthly variation of LAI in the Caxiuanã site. The hatched area indicates the 95% confidence interval.</p>
Full article ">Figure 12
<p>Annual carbon balance (<b>a</b>) and accumulated annual precipitation (<b>b</b>) at the Caxiuanã site for the years 2005 (aquamarine line), 2006 (blue line), 2007 (green line) and 2008 (red line).</p>
Full article ">
20 pages, 16133 KiB  
Article
Changes in Vegetation Cover and the Relationship with Surface Temperature in the Cananéia–Iguape Coastal System, São Paulo, Brazil
by Jakeline Baratto, Paulo Miguel de Bodas Terassi and Emerson Galvani
Remote Sens. 2024, 16(18), 3460; https://doi.org/10.3390/rs16183460 - 18 Sep 2024
Viewed by 881
Abstract
The objective of this article is to investigate the possible correlations between vegetation indices and surface temperature in the Cananéia–Iguape Coastal System (CICS), in São Paulo (Brazil). Vegetation index data from MODIS orbital products were used to carry out this work. The Normalized [...] Read more.
The objective of this article is to investigate the possible correlations between vegetation indices and surface temperature in the Cananéia–Iguape Coastal System (CICS), in São Paulo (Brazil). Vegetation index data from MODIS orbital products were used to carry out this work. The Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) were acquired from the MODIS/Aqua sensor (MYD13Q1) and the leaf area index (LAI) from the MODIS/Terra (MOD15A2H). Surface temperature data were acquired from MODIS/Aqua (MYD11A2). The data were processed using Google Earth Engine and Google Colab. The data were collected, and spatial and temporal correlations were applied. Correlations were applied in the annual and seasonal period. The annual temporal correlation between vegetation indices and surface temperature was positive, but statistically significant for the LAI, with r = 0.43 (90% significance). In the seasonal period, positive correlations occurred in JFM for all indices (95% significance). Spatially, the results of this research indicate that the largest area showed a positive correlation between VI and LST. The hottest and rainiest periods (OND and JFM) had clearer and more significant correlations. In some regions, significant and clear correlations were observed, such as in some areas in the north, south and close to the city of Iguape. This highlights the complexity of the interactions between vegetation indices and climatic attributes, and highlights the importance of considering other environmental variables and processes when interpreting changes in vegetation. However, this research has significantly progressed the field, by establishing new correlations and demonstrating the importance of considering climate variability, for a more accurate understanding of the impacts on vegetation indices. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area (<b>A</b>,<b>B</b>) and land use mapping (<b>C</b>).</p>
Full article ">Figure 2
<p>Variation in surface temperature and monthly (<b>A</b>) and annual (<b>B</b>) rainfall for the Cananéia-Iguape Coastal System for the 20032022 period. Source: MODIS/Aqua and CHIRPS, 2024.</p>
Full article ">Figure 3
<p>Annual variation in vegetation indices for the 2003–2022 period in the Cananéia–Iguape Coastal System.</p>
Full article ">Figure 4
<p>Scatter plot of annual NDVI (<b>a</b>), EVI (<b>b</b>) and LAI (<b>c</b>) values and surface temperature from 2003 to 2022.</p>
Full article ">Figure 5
<p>Scatter plot of seasonal values of VI–NDVI (<b>a</b>,<b>d</b>), EVI (<b>b</b>,<b>e</b>) and LAI (<b>c</b>,<b>f</b>)—and surface temperature for the JFM (<b>a</b>–<b>c</b>) and AMJ (<b>d</b>–<b>f</b>) quarter.</p>
Full article ">Figure 6
<p>Scatter plot of seasonal values of VI–NDVI (<b>a</b>,<b>d</b>), EVI (<b>b</b>,<b>e</b>) and LAI (<b>c</b>,<b>f</b>)—and climate variables for the JAS (<b>a</b>–<b>c</b>) and OND (<b>d</b>–<b>f</b>) quarter.</p>
Full article ">Figure 7
<p>Annual linear correlation between surface temperature and NDVI (<b>A</b>), EVI (<b>B</b>) and LAI (<b>C</b>) between 2003 and 2022.</p>
Full article ">Figure 8
<p>Seasonal linear correlation between surface temperature and VI between 2004 and 2022 for the JFM (<b>A</b>–<b>C</b>) and AMJ (<b>D</b>–<b>F</b>) periods.</p>
Full article ">Figure 9
<p>Seasonal linear correlation between surface temperature and VI between 2004 and 2022 for the JAS (<b>A</b>–<b>C</b>) and OND (<b>D</b>–<b>F</b>) periods.</p>
Full article ">
19 pages, 7218 KiB  
Article
Relationship between Vegetation and Soil Moisture Anomalies Based on Remote Sensing Data: A Semiarid Rangeland Case
by Juan José Martín-Sotoca, Ernesto Sanz, Antonio Saa-Requejo, Rubén Moratiel, Andrés F. Almeida-Ñauñay and Ana M. Tarquis
Remote Sens. 2024, 16(18), 3369; https://doi.org/10.3390/rs16183369 - 11 Sep 2024
Viewed by 753
Abstract
The dynamic of rangelands results from complex interactions between vegetation, soil, climate, and human activity. This scenario makes rangeland’s condition challenging to monitor, and degradation assessment should be carefully considered when studying grazing pressures. In the present work, we study the interaction of [...] Read more.
The dynamic of rangelands results from complex interactions between vegetation, soil, climate, and human activity. This scenario makes rangeland’s condition challenging to monitor, and degradation assessment should be carefully considered when studying grazing pressures. In the present work, we study the interaction of vegetation and soil moisture in semiarid rangelands using vegetation and soil moisture indices. We aim to study the feasibility of using soil moisture negative anomalies as a warning index for vegetation or agricultural drought. Two semiarid agricultural regions were selected in Spain for this study: Los Vélez (Almería) and Bajo Aragón (Teruel). MODIS images, with 250 m and 500 m spatial resolution, from 2002 to 2019, were acquired to calculate the Vegetation Condition Index (VCI) and the Water Condition Index (WCI) based on the Normalised Difference Vegetation Index (NDVI) and soil moisture component (W), respectively. The Optical Trapezoid Model (OPTRAM) estimated this latter W index. From them, the anomaly (Z-score) for each index was calculated, being ZVCI and ZWCI, respectively. The probability of coincidence of their negative anomalies was calculated every 10 days (10-day periods). The results show that for specific months, the ZWCI had a strong probability of informing in advance, where the negative ZVCI will decrease. Soil moisture content and vegetation indices show more similar dynamics in the months with lower temperatures (from autumn to spring). In these months, given the low temperatures, precipitation leads to vegetation growth. In the following months, water availability depends on evapotranspiration and vegetation type as the temperature rises and the precipitation falls. The stronger relationship between vegetation and precipitation from autumn to the beginning of spring is reflected in the feasibility of ZWCI to aid the prediction of ZVCI. During these months, using ZWCI as a warning index is possible for both areas studied. Notably, November to the beginning of February showed an average increase of 20–30% in the predictability of vegetation anomalies, knowing moisture soil anomalies four lags in advance. We found other periods of relevant increment in the predictability, such as March and April for Los Vélez, and from July to September for Bajo Aragón. Full article
(This article belongs to the Special Issue Advances in Remote Sensing for Regional Soil Moisture Monitoring)
Show Figures

Figure 1

Figure 1
<p>Map representing the selected pixels. In purple, the selected pixels of Los Vélez (<b>A</b>), and in red, the pixels of Bajo Aragón (<b>B</b>).</p>
Full article ">Figure 2
<p>Sketch illustrating parameters of the OPTRAM model used in Equations (3)–(5) to estimate parameters. Adapted from [<a href="#B25-remotesensing-16-03369" class="html-bibr">25</a>].</p>
Full article ">Figure 3
<p>Indices series (NDVI and WCI) for the average of the selected pixels of Los Velez (<b>A</b>) and Bajo Aragón (<b>B</b>).</p>
Full article ">Figure 4
<p>Time lagged cross-correlation of Indices series (NDVI and WCI) for the average of the selected pixels of Los Velez (<b>A</b>) and Bajo Aragón (<b>B</b>). Confident bounds at 95% are represented in red lines.</p>
Full article ">Figure 5
<p>Boxplots for vegetation condition index (VCI) and water condition index (WCI) for Los Vélez (<b>A</b>) and Bajo Aragón (<b>B</b>) every 10 days of the year. The blue vertical lines represent the phases split based on the VCI dynamics.</p>
Full article ">Figure 6
<p>Boxplots for temperature (orange) and precipitation (blue) for Los Vélez (<b>A</b>) and Bajo Aragón (<b>B</b>) every 10 days of the year. The blue vertical lines represent the phases split based on the VCI dynamics.</p>
Full article ">Figure 7
<p>Z-score series (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>Z</mi> </mrow> <mrow> <mi>V</mi> <mi>C</mi> <mi>I</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>Z</mi> </mrow> <mrow> <mi>W</mi> <mi>C</mi> <mi>I</mi> </mrow> </msub> </mrow> </semantics></math>) for the average of the selected pixels of Los Velez (<b>A</b>) and Bajo Aragón (<b>B</b>).</p>
Full article ">Figure 8
<p>Base probabilities for the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>Z</mi> </mrow> <mrow> <mi>V</mi> <mi>C</mi> <mi>I</mi> </mrow> </msub> </mrow> </semantics></math> to pass below the three thresholds every 10 days for Los Velez (<b>A</b>) and Bajo-Aragón (<b>B</b>). Thresholds −0.5 in blue, −0.7 in orange, and −1.0 in grey.</p>
Full article ">Figure 9
<p>Base probabilities for the <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>Z</mi> </mrow> <mrow> <mi>W</mi> <mi>C</mi> <mi>I</mi> </mrow> </msub> </mrow> </semantics></math> to pass below the three thresholds for each 10-day period for Los Velez (<b>A</b>) and Bajo-Aragón (<b>B</b>). Thresholds −0.5 in blue, −0.7 in orange, and −1.0 in grey.</p>
Full article ">Figure 10
<p>Base probability for <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>Z</mi> </mrow> <mrow> <mi>V</mi> <mi>C</mi> <mi>I</mi> </mrow> </msub> </mrow> </semantics></math> to be below −0.5 (grey), conditional probability without lag (dotted blue), and condition probability for lag-4 (red) in Los Velez (<b>A</b>) and Bajo-Aragón (<b>B</b>).</p>
Full article ">Figure 11
<p>Base probability for <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>Z</mi> </mrow> <mrow> <mi>V</mi> <mi>C</mi> <mi>I</mi> </mrow> </msub> </mrow> </semantics></math> to be below −0.7 (grey), conditional probability without lag (dotted blue), and conditional probability for lag-4 (red) in Los Velez (<b>A</b>) and Bajo-Aragón (<b>B</b>).</p>
Full article ">Figure 12
<p>Base probability for <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>Z</mi> </mrow> <mrow> <mi>V</mi> <mi>C</mi> <mi>I</mi> </mrow> </msub> </mrow> </semantics></math> to be below −1.0 (grey), conditional probability without lag (dotted blue), and conditional probability for lag-4 (red) in Los Velez (<b>A</b>) and Bajo-Aragón (<b>B</b>).</p>
Full article ">Figure 13
<p>Comparison of positive differences (lag-4 conditional probability minus base probability) for both areas of study and the thresholds −0.5 (<b>A</b>), −0.7 (<b>B</b>), and −1.0 (<b>C</b>).</p>
Full article ">
16 pages, 10159 KiB  
Article
Contribution of Climatic Change and Human Activities to Vegetation Dynamics over Southwest China during 2000–2020
by Gang Qi, Nan Cong, Man Luo, Tangzhen Qiu, Lei Rong, Ping Ren and Jiangtao Xiao
Remote Sens. 2024, 16(18), 3361; https://doi.org/10.3390/rs16183361 - 10 Sep 2024
Cited by 1 | Viewed by 791
Abstract
Southwest China is an important carbon sink area in China. It is critical to track and assess how human activity (HA) and climate change (CC) affect plant alterations in order to create effective and sustainable vegetation restoration techniques. This study used MODIS NDVI [...] Read more.
Southwest China is an important carbon sink area in China. It is critical to track and assess how human activity (HA) and climate change (CC) affect plant alterations in order to create effective and sustainable vegetation restoration techniques. This study used MODIS NDVI data, vegetation type data, and meteorological data to examine the regional and temporal variations in the normalized difference vegetation index (NDVI) in Southwest China from 2000 to 2020. Using trend analysis, the study looks at the temporal and geographical variability in the NDVI. Partial correlation analysis was also used to assess the effects of precipitation, extreme climate indicators, and mean temperature on the dynamics of the vegetation. A new residual analysis technique was created to categorize the effects of CC and HA on NDVI changes while taking extreme climate into consideration. The findings showed that the NDVI in Southwest China grew at a rate of 0.02 per decade between 2000 and 2020. According to the annual NDVI, there was a regional rise in around 85.59% of the vegetative areas, with notable increases in 36.34% of these regions. Temperature had a major influence on the northern half of the research region, but precipitation and extreme climate had a notable effect on the southern half. The rates at which climatic variables and human activity contributed to changes in the NDVI were 0.0008/10a and 0.0034/10a, respectively. These rates accounted for 19.1% and 80.9% of the variances, respectively. The findings demonstrate that most areas displayed greater HA-induced NDVI increases, with the exception of the western Sichuan Plateau. This result suggests that when formulating vegetation restoration and conservation strategies, special attention should be paid to the impact of human activities on vegetation to ensure the sustainable development of ecosystems. Full article
Show Figures

Figure 1

Figure 1
<p>The elevation map (<b>a</b>) and vegetation type map (<b>b</b>) of Southwest China. The terms “Needleleaf Forest (NLF)”, “Broadleaf Forest (BLF)”, “Shrublands (SHR)”, “Grasslands (GRA)”, and “Croplands (CRO)” refer to different kinds of vegetation.</p>
Full article ">Figure 2
<p>The spatial distribution of the multi-year average NDVI from 2000 to 2010 (<b>a</b>), the multi-year average NDVI from 2010 to 2020 (<b>b</b>), the annual trends in the NDVI for the years 2000–2020 (<b>c</b>), and the associated significances (<b>d</b>).</p>
Full article ">Figure 3
<p>Partial correlation coefficients’ spatial distributions from 2000 to 2020 between the NDVI of Southwest China and (<b>a</b>) temperature, (<b>b</b>) precipitation, (<b>c</b>) extreme climate, and (<b>d</b>) spatial distribution of the largest correlation factor. Correlation coefficients with significant values of 5% and 1% are represented by values of 0.46 (–0.46) and 0.53 (–0.53). The unit is the correlation coefficient.</p>
Full article ">Figure 4
<p>HA and CC’s contributions to NDVI variations in the study area. The spatial distribution of two contributions: (<b>a</b>) those resulting from CC, and (<b>b</b>) those resulting from HA. Significant upward trend, non-significant upward trend, unrelated trend, non-significant downward trend, and significant downward trend are denoted by the letters SUT, NUT, U, UNT, and SDT, respectively.</p>
Full article ">Figure 5
<p>Contributions of CC and HA to NDVI fluctuations under various vegetation types in the studied area. (<b>a</b>) Data about the impact of CC; (<b>b</b>) data regarding the impact of HA. Significant upward trend, non-significant upward trend, unrelated trend, non-significant downward trend, and significant downward trend are denoted by the letters SUT, NUT, U, UNT, and SDT, respectively. “Needleleaf Forest”, “Broadleaf Forest”, “Shrublands”, “Grasslands”, and “Croplands” are represented by the acronyms NLF, BLF, SHR, GRA, and CRO, in that order.</p>
Full article ">Figure 6
<p>Proportional contributions of HA and CC to NDVI fluctuations in the studied area. (<b>a</b>) The relative impact of climate change; (<b>b</b>) the relative impact of human activity.</p>
Full article ">Figure 7
<p>Spatial distribution of driving variables for plant cover changes in Southwest China from 2000 to 2020 (CC and HA denote climate change and human activities, respectively).</p>
Full article ">Figure 8
<p>Percentage of pixels under various plant species in Southwest China, where variations in vegetation cover are mostly caused by human activity and climate change. The symbols CC, HA, NDF, BDF, SHR, GRA, and CRO stand for “climate change”, “human activity”, “needleleaf forest”, “broadleaf forest”, “shrublands”, “grasslands”, and “croplands”, respectively.</p>
Full article ">
16 pages, 6543 KiB  
Article
Climate Warming Has Contributed to the Rise of Timberlines on the Eastern Tibetan Plateau but Slowed in Recent Years
by Xuefeng Peng, Yu Feng, Han Zang, Dan Zhao, Shiqi Zhang, Ziang Cai, Juan Wang and Peihao Peng
Atmosphere 2024, 15(9), 1083; https://doi.org/10.3390/atmos15091083 - 6 Sep 2024
Viewed by 791
Abstract
The alpine timberline is a component of terrestrial ecosystems and is highly susceptible to climate change. Since 2000, the Tibetan Plateau’s high-altitude zone has been experiencing a persistent warming, clarifying that the response of the alpine timberline to climate warming is important for [...] Read more.
The alpine timberline is a component of terrestrial ecosystems and is highly susceptible to climate change. Since 2000, the Tibetan Plateau’s high-altitude zone has been experiencing a persistent warming, clarifying that the response of the alpine timberline to climate warming is important for mitigating the negative impacts of global warming. However, it is difficult for traditional field surveys to clarify changes in the alpine timberline over a wide range of historical periods. Therefore, alpine timberline sites were extracted from 2000–2021, based on remote sensing data sources (LANDSAT, MODIS), to quantify the timberline vegetation growth in the Gexigou National Nature Reserve and to explore the impacts of climate change on timberline vegetation growth. The results show that the mean temperature increased significantly from 2000 to 2021 (R2 = 0.35, p = 0.0036) at a rate of +0.03 °C/year. The alpine timberline continued to shift upwards, but at a slower rate, by +22.87 m, +23.23 m, and +2.73 m in 2000–2007, 2007–2014, and 2014–2021, respectively. The sample plots of the timberline showing an upward shift experienced a decreasing trend. The timberline NDVI increased significantly from 2000 to 2021 (R2 = 0.2678, p = 0.0136) with an improvement in its vegetation. The timberline NDVI is positively correlated with the annual mean temperature (p < 0.05), February mean temperature (p < 0.05), June minimum temperature (p < 0.05), February maximum temperature (p < 0.01), June maximum temperature (p < 0.01), and June mean temperature (p < 0.01). It was also found to be negatively correlated with annual precipitation (p < 0.01). The study showcases the practicality of using remote sensing techniques to investigate the alpine timberline shifts and timberline vegetation. The findings are valuable in developing approaches to the sustainable management of timberline ecosystems. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the Gexigou National Nature Reserve in eastern Tibet Plateau. The map is reproduced from Tianditu (<a href="https://www.tianditu.gov.cn/" target="_blank">https://www.tianditu.gov.cn/</a>, accessed on 05 July 2024).</p>
Full article ">Figure 2
<p>Flowchart for the study of the impacts of climate change on timberline change on the eastern Tibetan Plateau (GNNR).</p>
Full article ">Figure 3
<p>Climate change trends at the timberline: (<b>a</b>) precipitation, (<b>b</b>) minimum temperature, (<b>c</b>) mean temperature, and (<b>d</b>) maximum temperature. The red line represents the linear trend of climate variables.</p>
Full article ">Figure 4
<p>Results of timberline extraction based on Landsat imagery. The background image is a Landsat-8 satellite image (LC08_131039_20211001), date: 01 October 2021.</p>
Full article ">Figure 5
<p>Spatial variations in timberline altitude: (<b>a</b>) 2000, (<b>b</b>) 2007, (<b>c</b>) 2014, and (<b>d</b>) 2021.</p>
Full article ">Figure 6
<p>Timberline altitude shift in (<b>a</b>) 2000–2021, (<b>b</b>) 2000–2007, (<b>c</b>) 2007–2014, and (<b>d</b>) 2014–2021.</p>
Full article ">Figure 7
<p>Trends in annual NDVI changes in the timberline area. The red line represents the linear trend of the annual NDVI.</p>
Full article ">Figure 8
<p>Spatial variations in timberline NDVI: (<b>a</b>) NDVI changes (NDVI<sub>diff</sub>) and (<b>b</b>) NDVI change trends (NDVI<sub>slope</sub>).</p>
Full article ">Figure 9
<p>Pearson correlation between timberline annual NDVI and climatic variables, (<b>a</b>) precipitation, (<b>b</b>) minimum temperature, (<b>c</b>) mean temperature, and (<b>d</b>) maximum temperature. * Indicates a significant correlation (<span class="html-italic">p</span> &lt; 0.05), and ** indicates an extremely significant correlation (<span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">
Back to TopTop