[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 (875)

Search Parameters:
Keywords = MODIS NDVI

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
29 pages, 16636 KiB  
Article
An Explanation of the Differences in Grassland NDVI Change in the Eastern Route of the China–Mongolia–Russia Economic Corridor
by Zhengfei Wang, Jiayue Wang, Wenlong Wang, Chao Zhang, Urtnasan Mandakh, Danzanchadav Ganbat and Nyamkhuu Myanganbuu
Remote Sens. 2025, 17(5), 867; https://doi.org/10.3390/rs17050867 - 28 Feb 2025
Viewed by 156
Abstract
This study analyzed the spatiotemporal changes in grassland NDVI from 2000 to 2020 in the eastern route of the China–Mongolia–Russia Economic Corridor, a region with frequent ecological–economic interactions, and explained the main driving factors, influencing patterns, and degrees of grassland NDVI changes in [...] Read more.
This study analyzed the spatiotemporal changes in grassland NDVI from 2000 to 2020 in the eastern route of the China–Mongolia–Russia Economic Corridor, a region with frequent ecological–economic interactions, and explained the main driving factors, influencing patterns, and degrees of grassland NDVI changes in different regions. Based on MODIS NDVI data, the study employs emerging spatiotemporal hotspot analysis, Maximum Relevance Minimum Redundancy (mRMR) feature selection, and Gaussian Process Regression (GPR) to reveal the spatiotemporal variation characteristics of grassland NDVI, while identifying long-term stable trends, and to select the most relevant and non-redundant factors to analyze the main driving factors of grassland NDVI change. Partial dependence plots were used to visualize the response and sensitivity of grassland NDVI to various factors. The results show the following: (1) From 2000 to 2020, the NDVI of grassland in the study area showed an overall upward trend, from 0.61 to 0.65, with significant improvement observed in northeastern China and northeastern Russia. (2) Spatiotemporal hotspot analysis indicates that 51% of the area is classified as persistent hotspots for grassland NDVI, mainly distributed in Russia, whereas 12% of the area is identified as persistent cold spots, predominantly located in Mongolia. (3) The analysis of key drivers reveals that precipitation and land surface temperature are the dominant climatic factors shaping grassland NDVI trends, while the effects of soil conditions and human activity vary regionally. In China, NDVI is primarily driven by land surface temperature (LST), GDP, and population density; in Mongolia, precipitation, LST, and GDP exert the strongest influence; whereas in Russia, livestock density and soil organic carbon play the most significant roles. (4) For the whole study area, in persistent cold spot areas of grassland NDVI, the negative effects of rising land surface temperature were most pronounced, reducing NDVI by 36% in the 25–40 °C range. The positive effects of precipitation on NDVI were most evident under low to moderate precipitation conditions, with the effects diminishing as precipitation increased. Soil moisture and soil pH have stronger effects in persistent hotspot areas. Regarding human activity factors, the livestock factor in Mongolia shows an inverted U-shaped relationship with NDVI, and increasing population density contributed to grassland degradation in persistent cold spots. Proper grazing intensity regulation strategy is crucial in these areas with inappropriate grazing intensity, while social and economic activities promoted vegetation cover improvement in persistent hotspots in China and Russia. These findings provide practical insights to guide grassland ecosystem restoration and ensure sustainable development along the eastern route of the China–Mongolia–Russia Economic Corridor. China should prioritize ecological compensation policies. Mongolia needs to integrate traditional nomadic grazing with modern practices. Russia should focus on strengthening regulatory frameworks to prevent the over-exploitation of grasslands. Especially for persistent cold spot areas of grassland NDVI in Mongolia and Russia that are prone to grassland degradation, attention should be paid to the significant negative impact of livestock on grassland. Full article
(This article belongs to the Section Environmental Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>An overview of the areas along the main railway of the eastern line of the China–Mongolia–Russia Economic Corridor. (<b>a</b>) The elevation distribution within the study area and the geographical locations of China, Mongolia, and Russia; (<b>b</b>) the spatial distribution of average NDVI during growing seasons from 2000 to 2020 for grasslands in the study area. (<a href="https://maps.elie.ucl.ac.be/CCI" target="_blank">https://maps.elie.ucl.ac.be/CCI</a>, accessed on 20 September 2024) (<a href="https://www.webmap.cn/" target="_blank">https://www.webmap.cn/</a>, accessed on 21 August 2024).</p>
Full article ">Figure 2
<p>Average spatial distribution characteristics of major environmental and socioeconomic elements in study area: (<b>a</b>) downward shortwave radiation (DSR); (<b>b</b>) land surface temperature; (<b>c</b>) soil moisture; (<b>d</b>) livestock density (in sheep units); (<b>e</b>) population; (<b>f</b>) annual precipitation.</p>
Full article ">Figure 3
<p>(<b>a</b>) Spatial pattern of Theil–Sen slope analysis showing NDVI trends; (<b>b</b>) frequency distribution of Theil–Sen slopes.</p>
Full article ">Figure 4
<p>Spatiotemporal patterns of grassland NDVI in eastern route of China–Mongolia–Russia Economic Corridor from 2000 to 2020. (<b>a</b>) Spatial distribution of emerging hotspot analysis and its statistical composition in three countries; (<b>b</b>) percentage of different spot types in each country; (<b>c</b>) annual variations in NDVI in different spot types for China, Mongolia, and Russia.</p>
Full article ">Figure 5
<p>Importance of driving factors for persistent cold spots and hotspots in whole study area—China, Mongolia, and Russia—analyzed using MRMR method. Numbers in brackets represent ranking of importance.</p>
Full article ">Figure 6
<p>Comparison of model performance (R<sup>2</sup>) across different machine learning algorithms using test set (<b>a</b>) and validation set (<b>b</b>) in different grassland clusters. Black line represents mean R<sup>2</sup> across all clusters. Legend explains regions and clusters as follows: CPCS (China Persistent Cold Spot), MPCS (Mongolia Persistent Cold Spot), RPCS (Russia Persistent Cold Spot), CPHS (China Persistent Hotspot), MPHS (Mongolia Persistent Hotspot), RPHS (Russia Persistent Hotspot), PCS (persistent cold spot for entire study area), and PHS (persistent hotspot for entire study area).</p>
Full article ">Figure 7
<p>Partial dependence plots of driving factors affecting grassland NDVI in CPHS. (<b>a</b>) Land surface temperature (°C); (<b>b</b>) soil pH; (<b>c</b>) GDP (millions of 2017 US dollar/km<sup>2</sup>); (<b>d</b>) downward surface shortwave radiation (W/m<sup>2</sup>); (<b>e</b>) soil moisture (mm); (<b>f</b>) soil organic carbon (g/kg); (<b>g</b>) air temperature (°C); (<b>h</b>) annual precipitation (mm/year).</p>
Full article ">Figure 8
<p>Partial dependence plots of driving factors affecting grassland NDVI in CPCS. (<b>a</b>) Land surface temperature (°C); (<b>b</b>) population density (person/km<sup>2</sup>); (<b>c</b>) annual precipitation (mm/year); (<b>d</b>) soil organic carbon (g/kg); (<b>e</b>) downward surface shortwave radiation (W/m<sup>2</sup>); (<b>f</b>) GDP (millions of 2017 US dollar/km<sup>2</sup>); (<b>g</b>) soil moisture (mm); (<b>h</b>) air temperature (°C).</p>
Full article ">Figure 9
<p>Partial dependence plots of driving factors affecting grassland NDVI in MPHS. (<b>a</b>) Soil organic carbon (g/kg); (<b>b</b>) downward surface shortwave radiation (W/m<sup>2</sup>); (<b>c</b>) GDP (millions of 2017 US dollar/km<sup>2</sup>); (<b>d</b>) land surface temperature (°C); (<b>e</b>) livestock density (sheep units); (<b>f</b>) annual precipitation (mm/year); (<b>g</b>) soil moisture (mm); (<b>h</b>) soil pH.</p>
Full article ">Figure 10
<p>Partial dependence plots of driving factors affecting grassland NDVI in MPCS. (<b>a</b>) Annual precipitation (mm/year); (<b>b</b>) soil organic carbon (g/kg); (<b>c</b>) GDP (millions of 2017 US dollar/km<sup>2</sup>); (<b>d</b>) land surface temperature (°C); (<b>e</b>) soil moisture (mm); (<b>f</b>) downward surface shortwave radiation (W/m<sup>2</sup>); (<b>g</b>) livestock density (sheep units); (<b>h</b>) population density (person/km<sup>2</sup>).</p>
Full article ">Figure 11
<p>Partial dependence plots of driving factors affecting grassland NDVI in RPHS. (<b>a</b>) Land surface temperature (°C); (<b>b</b>) GDP (millions of 2017 US dollar/km<sup>2</sup>); (<b>c</b>) soil organic carbon (g/kg); (<b>d</b>) annual precipitation (mm/year); (<b>e</b>) air temperature (°C); (<b>f</b>) soil bulk density (10 kg/m<sup>3</sup>; (<b>g</b>) population density (person/km<sup>2</sup>); (<b>h</b>) livestock density (sheep units).</p>
Full article ">Figure 12
<p>Partial dependence plots of driving factors affecting grassland NDVI in RPCS. (<b>a</b>) Livestock density (sheep units); (<b>b</b>) GDP (millions of 2017 US dollar/km<sup>2</sup>); (<b>c</b>) soil organic carbon (g/kg); (<b>d</b>) population density (person/km<sup>2</sup>); (<b>e</b>) soil pH; (<b>f</b>) land surface temperature (°C); (<b>g</b>) soil bulk density (10 kg/m<sup>3</sup>); (<b>h</b>) soil moisture (mm).</p>
Full article ">Figure 13
<p>Response intensity of NDVI to key factors.</p>
Full article ">
14 pages, 1607 KiB  
Article
Global NDVI-LST Correlation: Temporal and Spatial Patterns from 2000 to 2024
by Ehsan Rahimi, Pinliang Dong and Chuleui Jung
Environments 2025, 12(2), 67; https://doi.org/10.3390/environments12020067 - 17 Feb 2025
Viewed by 265
Abstract
While numerous studies have investigated the NDVI-LST relationship at local or regional scales, existing global analyses are outdated and fail to incorporate recent environmental changes driven by climate change and human activity. This study aims to address this gap by conducting an extensive [...] Read more.
While numerous studies have investigated the NDVI-LST relationship at local or regional scales, existing global analyses are outdated and fail to incorporate recent environmental changes driven by climate change and human activity. This study aims to address this gap by conducting an extensive global analysis of NDVI-LST correlations from 2000 to 2024, utilizing multi-source satellite data to assess latitudinal and ecosystem-specific variability. The MODIS dataset, which provides global daily LST data at a 1 km resolution from 2000 to 2024, was used alongside MODIS-derived NDVI data, which offers global vegetation indices at a 1 km resolution and 16-day temporal intervals. A correlation analysis was performed by extracting NDVI and LST values for each raster cell. The analysis revealed significant negative correlations in regions such as the western United States, Brazil, southern Africa, and northern Australia, where increased temperatures suppress vegetation activity. A total of 38,281,647 pixels, or 20% of the global map, exhibited statistically significant correlations, with 80.4% showing negative correlations, indicating a reduction in vegetation activity as temperatures rise. The latitudinal distribution of significant correlations revealed two prominent peaks: one in the tropical and subtropical regions of the Southern Hemisphere and another in the temperate zones of the Northern Hemisphere. This study uncovers notable spatial and latitudinal patterns in the LST-NDVI relationship, with most regions exhibiting negative correlations, underscoring the cooling effects of vegetation. These findings emphasize the crucial role of vegetation in regulating surface temperatures, providing valuable insights into ecosystem health, and informing conservation strategies in response to climate change. Full article
Show Figures

Figure 1

Figure 1
<p>Correlation map of LST-NDVI in six classes (<b>a</b>), and significant and non-significant pixels (<b>b</b>).</p>
Full article ">Figure 2
<p>Latitudinal distribution of significant correlations (<b>a</b>), and proportions of positive and negative significant correlations (<b>b</b>).</p>
Full article ">
13 pages, 1915 KiB  
Article
Temporal Relationships of Breeding Landbirds and Productivity on a Working Landscape
by Janel L. Ortiz, April A. T. Conkey, Maia L. Lipschutz, Leonard A. Brennan, David B. Wester, Tyler A. Campbell and Humberto L. Perotto-Baldivieso
Wild 2025, 2(1), 4; https://doi.org/10.3390/wild2010004 - 17 Feb 2025
Viewed by 293
Abstract
The Normalized Difference Vegetation Index (NDVI) is a measurement of landscape “greenness” and is used as a proxy for productivity to assess species distributions and habitats. Seasonal levels of productivity have been strongly related to avian population dynamics, suggesting dependence upon biomass production [...] Read more.
The Normalized Difference Vegetation Index (NDVI) is a measurement of landscape “greenness” and is used as a proxy for productivity to assess species distributions and habitats. Seasonal levels of productivity have been strongly related to avian population dynamics, suggesting dependence upon biomass production for completing annual life cycle events. The breeding season is a critical component of the avian life cycle that involves higher nutritional requirements to feed young, avoiding predators, and attracting mates. Our objective was to determine how the NDVI affects avian abundance and richness across breeding seasons with varied rainfall in South Texas, USA. Breeding bird point-count surveys were conducted, and MODIS Terra NDVI data were collected. We observed both positive and negative effects between May and June avian abundance, richness, and the NDVI depending upon the year (i.e., wet or average rainfall) and NDVI values in the months prior to (i.e., April) and during the peak of breeding season (May), with no significant effect of the NDVI in June, suggesting the months prior to peak breeding season may be most influential. This information can aid land management recommendations and better predict how environmental changes like rainfall may affect avian dynamics on a landscape for both wildlife and domestic animals. Full article
Show Figures

Figure 1

Figure 1
<p>East Foundation’s El Sauz Ranch is located within Willacy and Kenedy counties, Texas, USA, along the south Texas coast. Breeding bird survey (BBS) point-count stations (37) are illustrated with black dots and were used for the length of this study from 2014 to 2016.</p>
Full article ">Figure 2
<p>The total annual precipitation received by Port Mansfield, Texas, USA, from 2006 to 2016. Precipitation is shown for a ten-year span to display fluctuations in annual precipitation and drought beginning in 2010. This study took place from 2014 to 2016 and is indicated by the red arrows, with 2014 considered a recovery, post-drought year that received average precipitation. The year 2015 was above average, and 2016 was an average year. Average annual precipitation is depicted by the dashed line.</p>
Full article ">
18 pages, 2266 KiB  
Article
Soybean Yield Modeling and Analysis with Weather Dynamics in the Greater Mississippi River Basin
by Weiwei Xie, Yanbo Huang and Qingmin Meng
Climate 2025, 13(2), 33; https://doi.org/10.3390/cli13020033 - 6 Feb 2025
Viewed by 479
Abstract
Accurate crop yield prediction and modeling are essential for ensuring food security, optimizing resource allocation, and guiding policy decisions in agriculture, ultimately benefiting society at large. With the increasing threat of weather change, it is important to understand the impacts of weather dynamics [...] Read more.
Accurate crop yield prediction and modeling are essential for ensuring food security, optimizing resource allocation, and guiding policy decisions in agriculture, ultimately benefiting society at large. With the increasing threat of weather change, it is important to understand the impacts of weather dynamics on agricultural productivity, particularly for crucial crops like soybeans. This study considers the study area of the Greater Mississippi River Basin, where most soybeans are typically planted, with a large variety of weather across from the North to the South in the US. Leveraging the greenness and density measured by the normalized difference vegetation index (NDVI) from NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) satellite images, along with weather variables including mean precipitation, minimum temperature, and maximum temperature, we aim to uncover the relationships between these variables and soybean yield for different geographical and weather regions. Our analysis focuses on the four weather regions within the US: Very Cold, Cold, Mixed Humid, and Hot Humid, where most soybeans are planted in the Mississippi River Basin. The findings reveal that soybean yield in the Cold and Very Cold regions is positively correlated with minimum temperatures, whereas in the Mixed Humid and Hot Humid regions, negative correlations between maximum temperatures and yields are found. We identify a significant positive correlation between precipitation and soybean yield across all regions. In addition, the NDVI shows significant positive correlations with the soybean yield. Both linear and nonlinear regression models, including support vector machine and random forest models, are trained with remotely sensed data and weather data, showing a reliable and improved crop yield prediction. The findings of this study contribute to a better understanding of how soybean yield responds to climatic variations and will help the national agricultural management system in better monitoring and predicting crop yield when facing the increasing challenge of weather dynamics. Full article
Show Figures

Figure 1

Figure 1
<p>The study area of the Greater Mississippi Rivier Basin for soybean yield modeling.</p>
Full article ">Figure 2
<p>The schematic diagram for geographic data processing and yield modeling.</p>
Full article ">Figure 3
<p>Average soybean yield of four environmental regions over the years.</p>
Full article ">Figure 4
<p>The temporal pattern of VGM<sub>max</sub> (<b>a</b>) and three corresponding climatic factors: mean precipitation (<b>b</b>), minimum temperature (<b>c</b>), and maximum temperature (<b>d</b>).</p>
Full article ">Figure 5
<p>The relationship between the yield and four explanatory variables: VGM<sub>max</sub>, mean precipitation, maximum temperature, and minimum temperature at the county level in (<b>a</b>) Cold, (<b>b</b>) Hot Humid, (<b>c</b>) Mixed Humid, and (<b>d</b>) Very Cold regions.</p>
Full article ">Figure 6
<p>The temporal relationship between VGM<sub>max</sub> and yield, both normalized, over time in Cold (<b>a</b>), Hot Humid (<b>b</b>), Mixed Humid (<b>c</b>), and Very Cold (<b>d</b>) regions.</p>
Full article ">
20 pages, 4669 KiB  
Article
Monitoring Mangrove Phenology Based on Gap Filling and Spatiotemporal Fusion: An Optimized Mangrove Phenology Extraction Approach (OMPEA)
by Yu Hong, Runfa Zhou, Jinfu Liu, Xiang Que, Bo Chen, Ke Chen, Zhongsheng He and Guanmin Huang
Remote Sens. 2025, 17(3), 549; https://doi.org/10.3390/rs17030549 - 6 Feb 2025
Viewed by 517
Abstract
Monitoring mangrove phenology requires frequent, high-resolution remote sensing data, yet satellite imagery often suffers from coarse resolution and cloud interference. Traditional methods, such as denoising and spatiotemporal fusion, faced limitations: denoising algorithms usually enhance temporal resolution without improving spatial quality, while spatiotemporal fusion [...] Read more.
Monitoring mangrove phenology requires frequent, high-resolution remote sensing data, yet satellite imagery often suffers from coarse resolution and cloud interference. Traditional methods, such as denoising and spatiotemporal fusion, faced limitations: denoising algorithms usually enhance temporal resolution without improving spatial quality, while spatiotemporal fusion models struggle with prolonged data gaps and heavy noise. This study proposes an optimized mangrove phenology extraction approach (OMPEA), which integrates Landsat and MODIS data with a denoising algorithm (e.g., Gap Filling and Savitzky–Golay filtering, GF–SG) and a spatiotemporal fusion model (e.g., Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model, ESTARFM). The key of OMPEA is that GF–SG algorithm filled data gaps from cloud cover and satellite transit gaps, providing high-quality input to ESTARFM and improving its accuracy of NDVI imagery reconstruction in mangrove phenology extraction. By conducting experiments on the GEE platform, OMPEA generates 1-day, 30 m NDVI imagery, from which phenological parameters (i.e., the start (SoS), end (EoS), length (LoS), and peak (PoS) of the growing season) are derived using the maximum separation (MS) method. Validation in four mangrove areas along the coastal China shows that OMPEA significantly improves the potential to capture mangrove phenology in the presence of incomplete data. The OMPEA significantly increased usable data, adding 7–33 Landsat images and 318–415 MODIS images per region. The generated NDVI series exhibits strong spatiotemporal consistency with original data (R2: 0.788–0.998, RMSE: 0.007–0.253) and revealed earlier SoS and longer LoS at lower latitudes. Cross-correlation analysis showed a 2–3 month lagged effects of temperature on mangroves’ growth, with precipitation having minimal impact. The proposed OMPEA improves the possibility of capturing mangrove phenology under non-continuous and low-resolution data, providing valuable insights for large-scale and long-term mangrove conservation and management. Full article
Show Figures

Figure 1

Figure 1
<p>Map of study area: (<b>a</b>) the overall distribution of study area; (<b>b1</b>–<b>b4</b>) Zhangjiangkou National Mangrove Nature Reserve (ZNR) in Fujian Province, Qi’ao Island Provincial Nature Reserve (QPR) in Guangdong Province, Beilun Estuary National Nature Reserve (BNR) in Guangxi Province, and Dongzhaigang National Mangrove Nature Reserve (DNR) in Hainan Province.</p>
Full article ">Figure 2
<p>Workflow of mangrove phenology extraction based on OMPEA.</p>
Full article ">Figure 3
<p>Landsat 8 NDVI (16-day 30 m) and denoised Landsat NDVI (16-day 30 m) generated by OMPEA. Gray pixel indicates pixel with no data.</p>
Full article ">Figure 4
<p>MODIS NDVI (1-day 500 m) and denoised MODIS NDVI (1-day 30 m) generated by OMPEA. Gray pixel indicates pixel with no data.</p>
Full article ">Figure 5
<p>The OMPEA-generated fused NDVI imagery. Gray pixel indicates pixel with no data.</p>
Full article ">Figure 6
<p>Scatter density plots and marginal histograms of fused NDVI and denoised Landsat NDVI.</p>
Full article ">Figure 7
<p>Composite scatter plots and line plots of various NDVI time series.</p>
Full article ">Figure 8
<p>Fused NDVI time-series curve and phenological parameters.</p>
Full article ">Figure 9
<p>Boxplots of mangrove phenological parameters.</p>
Full article ">Figure 10
<p>The time-series curves for fused NDVI, precipitation, temperature, and their lagged time-series curves with corresponding lag days.</p>
Full article ">Figure 11
<p>The OMPEA-generated fused NDVI in QPR from 17 January 2020 to 24 March 2021. (<b>a</b>) Description of denoised Landsat 8 NDVI in a full-time range. (<b>b</b>) Description of denoised Landsat 8 NDVI across three different time ranges, (<b>c</b>,<b>d</b>) is fused NDVI that using (<b>a</b>,<b>b</b>) as inputs, respectively. Gray pixel indicates pixel with no data.</p>
Full article ">
24 pages, 8969 KiB  
Article
Integrating Climate Data and Remote Sensing for Maize and Wheat Yield Modelling in Ethiopia’s Key Agricultural Region
by Asfaw Kebede Kassa, Hongwei Zeng, Bingfang Wu, Miao Zhang, Kibebew Kibret Tsehai, Xingli Qin and Tesfay G. Gebremicael
Remote Sens. 2025, 17(3), 491; https://doi.org/10.3390/rs17030491 - 30 Jan 2025
Viewed by 694
Abstract
Traditional methods for crop data collection are labor-intensive, inefficient and, more costly compared to remote sensing (RS) techniques. This study aims to identify key climatic variables influencing maize and wheat yields and develop predictive models while also evaluating the performance of the CropWatch [...] Read more.
Traditional methods for crop data collection are labor-intensive, inefficient and, more costly compared to remote sensing (RS) techniques. This study aims to identify key climatic variables influencing maize and wheat yields and develop predictive models while also evaluating the performance of the CropWatch cloud yield prediction model (CW_YPM) in major agricultural regions of Ethiopia. Climate data from 54 meteorological stations spanning 2000–2021 were analyzed. RS data, including NDVI from MODIS at 250 m resolution, agroecological zones, and observed crop yield data, were utilized for model prediction and validation. Correlation analysis and a stepwise modeling approach with multiple regression models were applied. The results revealed regional variations in the effects of climatic parameters on yields, with vapor pressure deficits showing negative correlations and rainfall exhibiting positive correlations. Non-linear models generally outperformed linear models in yield prediction—using both climate-only (CO) and combined climate-NDVI data. The best CO model for maize in the Horo Guduru area achieved an RMSE of 0.392 tons/ha, an R2 of 0.94, and an index of agreement (d) of 0.984. Incorporating NDVI improved accuracy, with the best maize model in the Illu Ababor area achieving an RMSE of 0.477 tons/ha, an R2 of 0.91, and d of 0.976. CW_YPM also performed effectively across the study area. This research highlights the value of integrating critical climatic variables with the NDVI to enhance crop yield forecasting in Ethiopia, thereby-supporting agricultural planning and food security initiatives. Full article
Show Figures

Figure 1

Figure 1
<p>Location map of the study area, highlighting 13 selected administrative zones known for wheat and maize cultivation, along with the distribution of meteorological stations used in the study.</p>
Full article ">Figure 2
<p>Historical (2000–2021) climate variables (seasonal areal rainfall, average temperature (Tmean), and vapor pressure deficit) of selected administrative zones: (<b>a</b>) Arsi from wheat growing area (June to October) and (<b>b</b>) Illu Ababora from maize growing area (May to September).</p>
Full article ">Figure 3
<p>Historical grain yield: (<b>a</b>) total grain production in regional states, Ethiopia (2019/2020 and 2020/2021; (<b>b</b>) maize and wheat yield data at selected administrative zones in Oromia region (2000 to 2021), Ethiopia [<a href="#B15-remotesensing-17-00491" class="html-bibr">15</a>].</p>
Full article ">Figure 3 Cont.
<p>Historical grain yield: (<b>a</b>) total grain production in regional states, Ethiopia (2019/2020 and 2020/2021; (<b>b</b>) maize and wheat yield data at selected administrative zones in Oromia region (2000 to 2021), Ethiopia [<a href="#B15-remotesensing-17-00491" class="html-bibr">15</a>].</p>
Full article ">Figure 4
<p>General methodology flow chart. (RF = rainfall).</p>
Full article ">Figure 5
<p>Climate variables and Normalized Difference Vegetation Index (NDVI) correlation analysis with (<b>a</b>) maize and (<b>b</b>) wheat yield at selected administrative zones in the study area.</p>
Full article ">Figure 6
<p>Residual plot for observed and model-fitted crop yield (maize and wheat) for all the study areas in the zonal administrations.</p>
Full article ">Figure 6 Cont.
<p>Residual plot for observed and model-fitted crop yield (maize and wheat) for all the study areas in the zonal administrations.</p>
Full article ">Figure 7
<p>Comparison between observed maize and wheat yield and their corresponding predicted yields generated by the top-performing “Climate only” and “Climate and NDVI” models across the study region.</p>
Full article ">Figure 8
<p>Scatter plots for predicted versus observed maize yield, “CO = Climate only and CaNDVI = Climate-NDVI variables” across the study area.</p>
Full article ">Figure 8 Cont.
<p>Scatter plots for predicted versus observed maize yield, “CO = Climate only and CaNDVI = Climate-NDVI variables” across the study area.</p>
Full article ">Figure 9
<p>Scatter plots for predicted versus observed wheat yield, “CO = Climate only and CaNDVI = Climate-NDVI variables” across the study area.</p>
Full article ">Figure 9 Cont.
<p>Scatter plots for predicted versus observed wheat yield, “CO = Climate only and CaNDVI = Climate-NDVI variables” across the study area.</p>
Full article ">Figure 10
<p>Spatial distribution of maize and wheat crop yield in 2021 in two zones predicted using CropWatch yield prediction model.</p>
Full article ">Figure 11
<p>Comparison of observed and predicted (CropWatch crop yield prediction model) yield for the period 2013 to 2021: (<b>a</b>) maize, Illu Ababora zone; (<b>b</b>) wheat, Bale zone.</p>
Full article ">
1 pages, 118 KiB  
Retraction
RETRACTED: Ji et al. Temporal-Difference Graph-Based Optimization for High-Quality Reconstruction of MODIS NDVI Data. Remote Sens. 2024, 16, 2713
by Shengtai Ji, Shuxin Han, Jiaxin Hu, Yuguang Li and Jing-Cheng Han
Remote Sens. 2025, 17(3), 349; https://doi.org/10.3390/rs17030349 - 21 Jan 2025
Viewed by 417
Abstract
The Remote Sensing Editorial Office retracts the article titled “Temporal-Difference Graph-Based Optimization for High-Quality Reconstruction of MODIS NDVI Data” by Ji et al [...] Full article
(This article belongs to the Section Environmental Remote Sensing)
23 pages, 6752 KiB  
Article
Development of Fractional Vegetation Cover Change and Driving Forces in the Min River Basin on the Eastern Margin of the Tibetan Plateau
by Shuyuan Liu, Li Zhou, Huan Wang, Jin Lin, Yuduo Huang, Peng Zhuo and Tianqi Ao
Forests 2025, 16(1), 142; https://doi.org/10.3390/f16010142 - 14 Jan 2025
Viewed by 621
Abstract
Fractional vegetation cover (FVC) is an important indicator of regional ecological environment change, and quantitative research on the spatial and temporal distribution of FVC and the trend of change is of great significance to the monitoring, evaluation, protection, and restoration of regional ecology. [...] Read more.
Fractional vegetation cover (FVC) is an important indicator of regional ecological environment change, and quantitative research on the spatial and temporal distribution of FVC and the trend of change is of great significance to the monitoring, evaluation, protection, and restoration of regional ecology. This study estimates the FVC of the eastern Tibetan Plateau margin from 2000 to 2020 using the image element dichotomous model based on the Google Earth Engine platform using MODIS-NDVI images. It also investigates the temporal and spatial changes of the FVC in this region and its drivers using the Theil–Sen and Mann–Kendall trend tests, spatial autocorrelation analysis, geodetector, and machine learning approaches impact. The results of this study indicated a generally erratic rising tendency, with the Min River Basin (MRB) near the eastern tip of the Tibetan Plateau having an annual average FVC of 0.67 and an annual growth rate of 0.16%. The percentage of places with better vegetation reached 60.37%. The regional FVC showed significant positive spatial autocorrelation and was clustered. Driver analyses showed that soil type, DEM, temperature, potential evapotranspiration, and land use type were the main drivers influencing FVC on the eastern margin of the Tibetan Plateau. In addition, the random forest (RF) model outperformed the support vector machine (SVM), backpropagation neural network (BP), and long short-term memory network (LSTM) in FVC regression fitting. In summary, this study shows that the overall FVC in the eastern margin of the Tibetan Plateau is on an upward trend, and the regional ecological environment has improved significantly over the past two decades. Full article
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Location of the MRB study area, (<b>a</b>) Specific location on the Tibetan Plateau, (<b>b</b>) DEM, (<b>c</b>) Land use.</p>
Full article ">Figure 2
<p>Spatial distribution of the drivers in 2015.</p>
Full article ">Figure 3
<p>(<b>a</b>) Proportion of each FVC type from 2000 to 2020; (<b>b</b>) temporal trend of FVC variation.</p>
Full article ">Figure 4
<p>Spatial pattern of different classes of FVC on the eastern margin of the Tibetan Plateau, 2000–2020.</p>
Full article ">Figure 5
<p>FVC spatial transfer area distribution (<b>a</b>) from 2000 to 2010; (<b>b</b>) from 2010 to 2020.</p>
Full article ">Figure 6
<p>Trends in FVC and their significance.</p>
Full article ">Figure 7
<p>FVC global spatial autocorrelation.</p>
Full article ">Figure 8
<p>FVC localized spatial autocorrelation LISA aggregation distribution.</p>
Full article ">Figure 9
<p>FVC factor detection results.</p>
Full article ">Figure 10
<p>Interaction test results of vegetation cover drivers in different years (NE indicates nonlinear enhancement, BE indicates two-factor enhancement).</p>
Full article ">Figure 11
<p>Significance statistics for differences in the impact of each driver.</p>
Full article ">Figure 12
<p>Statistical findings for various FVC types or ranges for every factor.</p>
Full article ">Figure 13
<p>Comparison of true and regression values of FVC: (<b>a</b>) SVM, (<b>b</b>) BP, (<b>c</b>) LSTM, (<b>d</b>) RF.</p>
Full article ">
27 pages, 5909 KiB  
Article
A Phenologically Simplified Two-Stage Clumping Index Product Derived from the 8-Day Global MODIS-CI Product Suite
by Ge Gao, Ziti Jiao, Zhilong Li, Chenxia Wang, Jing Guo, Xiaoning Zhang, Anxin Ding, Zheyou Tan, Sizhe Chen, Fangwen Yang and Xin Dong
Remote Sens. 2025, 17(2), 233; https://doi.org/10.3390/rs17020233 - 10 Jan 2025
Viewed by 545
Abstract
The clumping index (CI) is a key structural parameter that quantifies the nonrandomness of the spatial distribution of vegetation canopy leaves. Investigating seasonal variations in the CI is crucial, especially for estimating the leaf area index (LAI) and studying global carbon and water [...] Read more.
The clumping index (CI) is a key structural parameter that quantifies the nonrandomness of the spatial distribution of vegetation canopy leaves. Investigating seasonal variations in the CI is crucial, especially for estimating the leaf area index (LAI) and studying global carbon and water cycles. However, accurate estimations of the seasonal CI have substantial challenges, e.g., from the need for accurate hot spot measurements, i.e., the typical feature of the bidirectional reflectance distribution function (BRDF) shape in the current CI algorithm framework. Therefore, deriving a phenologically simplified stable CI product from a high-frequency CI product (e.g., 8 days) to reduce the uncertainty of CI seasonality and simplify CI applications remains important. In this study, we applied the discrete Fourier transform and an improved dynamic threshold method to estimate the start of season (SOS) and end of season (EOS) from the CI time series and indicated that the CI exhibits significant seasonal variation characteristics that are generally consistent with the MODIS land surface phenology (LSP) product (MCD12Q2), although seasonal differences between them probably exist. Second, we divided the vegetation cycle into two phenological stages based on the MODIS LSP product, ignoring the differences mentioned above, i.e., the leaf-on season (LOS, from greenup to dormancy) and the leaf-off season (LFS, after dormancy and before greenup of the next vegetation cycle), and developed the phenologically simplified two-stage CI product for the years 2001–2020 using the MODIS 8-day CI product suite. Finally, we assessed the accuracy of this CI product (RMSE = 0.06, bias = 0.01) via 95 datasets from 14 field-measured sites globally. This study revealed that the CI exhibited an approximately inverse trend in terms of phenological variation compared with the NDVI. Globally, based on the phenologically simplified two-stage CI product, the CILOS is smaller than the CILFS across all land cover types. Compared with the LFS stage, the quality for this CI product is better in the LOS stage, where the QA is basically identified as 0 and 1, accounting for more than ~90% of the total quality flag, which is significantly higher than that in the LFS stage (~60%). This study provides relatively reliable CI datasets that capture the general trend of seasonal CI variations and simplify potential applications in modeling ecological, meteorological, and other surface processes at both global and regional scales. Therefore, this study provides both new perspectives and datasets for future research in relation to CI and other biophysical parameters, e.g., the LAI. Full article
(This article belongs to the Section Atmospheric Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Distribution of collected field CI measurements (red dots) and typical pixels (white triangles) for all IGBP classes at the global scale. These datasets are mainly distributed along the mid-latitudes where vegetation seasonality tends to be easily identified.</p>
Full article ">Figure 2
<p>Seasonal variation analysis of the CI and flow chart of the MODIS time-share two-stage CI product.</p>
Full article ">Figure 3
<p>Diagram of the phenometrics retrieved for a single hypothetical vegetation cycle of the MODIS LSP product (MCD12Q2, V061).</p>
Full article ">Figure 4
<p>Accuracy evaluation results of the estimation of typical pixel vegetation phenology parameters from the CI time series. The accuracy evaluation metrics are the root mean square error (RMSE) and bias, with the minimum errors for each IGBP class highlighted by a bold square, and the smaller the error is, the smaller the square. Red indicates overestimated SOS and EOS values, whereas blue represents underestimated values. ENF: evergreen needleleaf forests; EBF: evergreen broadleaf forests; DNF: deciduous needleleaf forests; DBF: deciduous broadleaf forests; MF: mixed forests; Csh: closed shrublands; Osh: open shrublands; Wsa: woody savannas; Sav: savannas; GL: grasslands; PWe: permanent wetlands; CL<sup>1</sup>: annual croplands; CL<sup>2-1</sup>: the first vegetation cycles of biannual cropland; CL<sup>2-2</sup>: the second vegetation cycles of biannual cropland; CVM: cropland/natural vegetation mosaics. (<b>a</b>) RMSE of the estimated SOS from the CI. (<b>b</b>) RMSE of the estimated EOS from the CI. (<b>c</b>) Bias of the estimated SOS from the CI. (<b>d</b>) Bias of the estimated EOS from the CI.</p>
Full article ">Figure 5
<p>Accuracy evaluation results of the estimation of typical pixel vegetation phenology parameters from the NDVI time series. (<b>a</b>) RMSE of the estimated SOS from the NDVI. (<b>b</b>) RMSE of the estimated EOS from the NDVI. (<b>c</b>) Bias of the estimated SOS from the NDVI. (<b>d</b>) Bias of the estimated EOS from the NDVI.</p>
Full article ">Figure 6
<p>The accuracy evaluation results compared the time-share two-stage CIs with field-measured CIs, marking data overestimated or underestimated by more than 0.1 with gray dots (kinds of outliers). The values in parentheses represent the accuracy evaluation results after removing the gray-dotted data.</p>
Full article ">Figure 7
<p>Temporal variation in the time-share two-stage CIs at the (<b>a</b>) mixed forest, (<b>b</b>) woody savanna, (<b>c</b>) deciduous broadleaf forest, and (<b>d</b>) savanna field sites. Details of the field-measured CI data are shown in <a href="#app1-remotesensing-17-00233" class="html-app">Appendix A</a>. Red pentagons indicate the field-measured CI data. The black and blue dots indicate the CI<sub>LOS</sub> and CI<sub>LFS</sub>, respectively.</p>
Full article ">Figure 8
<p>Global distribution of multiyear average time-share two-stage CIs for (<b>a</b>) LOS and (<b>b</b>) LFS in the first vegetation cycle from 2001 to 2020.</p>
Full article ">Figure 9
<p>Distribution of the multiyear average CI<sub>LOS</sub> and CI<sub>LFS</sub> across different land cover types from 2001 to 2020.</p>
Full article ">Figure 10
<p>Global distribution of the mode of QA for (<b>a</b>) LOS and (<b>b</b>) LFS, and histogram distribution of QA for the MODIS time-share two-stage CI product for (<b>c</b>) LOS and (<b>d</b>) LFS from 2001 to 2020.</p>
Full article ">
20 pages, 14318 KiB  
Article
Multi-Feature Driver Variable Fusion Downscaling TROPOMI Solar-Induced Chlorophyll Fluorescence Approach
by Jinrui Fan, Xiaoping Lu, Guosheng Cai, Zhengfang Lou and Jing Wen
Agronomy 2025, 15(1), 133; https://doi.org/10.3390/agronomy15010133 - 8 Jan 2025
Viewed by 573
Abstract
Solar-induced chlorophyll fluorescence (SIF), as a direct indicator of vegetation photosynthesis, offers a more accurate measure of plant photosynthetic dynamics than traditional vegetation indices. However, the current SIF satellite products have low spatial resolution, limiting their application in fine-scale agricultural research. To address [...] Read more.
Solar-induced chlorophyll fluorescence (SIF), as a direct indicator of vegetation photosynthesis, offers a more accurate measure of plant photosynthetic dynamics than traditional vegetation indices. However, the current SIF satellite products have low spatial resolution, limiting their application in fine-scale agricultural research. To address this, we leveraged MODIS data at a 1 km resolution, including bands b1, b2, b3, and b4, alongside indices such as the NDVI, EVI, NIRv, OSAVI, SAVI, LAI, FPAR, and LST, covering October 2018 to May 2020 for Shandong Province, China. Using the Random Forest (RF) model, we downscaled SIF data from 0.05° to 1 km based on invariant spatial scaling theory, focusing on the winter wheat growth cycle. Various machine learning models, including CNN, Stacking, Extreme Random Trees, AdaBoost, and GBDT, were compared, with Random Forest yielding the best performance, achieving R2 = 0.931, RMSE = 0.052 mW/m2/nm/sr, and MAE = 0.031 mW/m2/nm/sr for 2018–2019 and R2 = 0.926, RMSE = 0.058 mW/m2/nm/sr, and MAE = 0.034 mW/m2/nm/sr for 2019–2020. The downscaled SIF products showed a strong correlation with TanSIF and GOSIF products (R2 > 0.8), and consistent trends with GPP further confirmed the reliability of the 1 km SIF product. Additionally, a time series analysis of Shandong Province’s wheat-growing areas revealed a strong correlation (R2 > 0.8) between SIF and multiple vegetation indices, underscoring its utility for regional crop monitoring. Full article
Show Figures

Figure 1

Figure 1
<p>Overview map of the Shandong study area.</p>
Full article ">Figure 2
<p>Downscaling method flow.</p>
Full article ">Figure 3
<p>Correlation analysis heatmap.</p>
Full article ">Figure 4
<p>Density fit plots for model projections from October 2018 to May 2019 (seven variables).</p>
Full article ">Figure 5
<p>Density fit plots for model projections from October 2018 to May 2019(welve variables).</p>
Full article ">Figure 6
<p>Density fitting map of RF model predictions from October 2019 to May 2020.</p>
Full article ">Figure 7
<p>Monthly eSIF from October 2018 to May 2019.</p>
Full article ">Figure 8
<p>Monthly 1 km SIF from October 2018 to May 2019.</p>
Full article ">Figure 9
<p>Monthly eSIF from October 2019 to May 2020.</p>
Full article ">Figure 10
<p>Monthly 1 km SIF from October 2019 to May 2020.</p>
Full article ">Figure 11
<p>Validation at the 0.05° spatial scale from October 2018 to May 2019.</p>
Full article ">Figure 12
<p>Validation at the 0.05° spatial scale from October 2019 to May 2020.</p>
Full article ">Figure 13
<p>Density fitting map of 1 km SIF and other SIF satellite products verified. (<b>a</b>) 2018–2019 TanSIF. (<b>b</b>) 2018–2019 GOSIF. (<b>c</b>) 2019–2020 GOSIF.</p>
Full article ">Figure 14
<p>Correlation analysis of SIF with a GPP density plot. (<b>a</b>) 2018–2019. (<b>b</b>) 2019–2020.</p>
Full article ">Figure 15
<p>Line chart of the SIF and vegetation indices of wheat.</p>
Full article ">Figure 16
<p>Correlation analysis between SIF and vegetation index.</p>
Full article ">
25 pages, 19410 KiB  
Article
Calibration and Validation of MODIS-Derived Ground-Level Air Temperature Models by Means of Ground Measurements
by Marica Teresa Rocca, Marica Franzini and Vittorio Marco Casella
Appl. Sci. 2025, 15(1), 184; https://doi.org/10.3390/app15010184 - 28 Dec 2024
Viewed by 669
Abstract
The research initiatives envisaged by the PNRR (Italian National Recovery and Resilience Plan) include the creation of innovation ecosystems to promote collaboration between universities, research centers, and local institutions with a focus on territorial integration and sustainability. The NODES Project (Nord-Ovest Digitale E [...] Read more.
The research initiatives envisaged by the PNRR (Italian National Recovery and Resilience Plan) include the creation of innovation ecosystems to promote collaboration between universities, research centers, and local institutions with a focus on territorial integration and sustainability. The NODES Project (Nord-Ovest Digitale E Sostenibile) is part of this research. In this context, the Laboratory of Geomatics of the University of Pavia, in collaboration with other partners, deals with the study of the suitability maps for the renowned Pinot Noir wine. To achieve this, we considered different thematic input layers: elevation, slope, aspect, soil depth and type, Land Use Land Cover maps, NDVI, and current and forecast climatic aspects. An important thematic layer is concerned with the air temperature, which requires high spatial and temporal resolution. In the selected study area, the Lombardy Region has some accurate and reliable weather stations with high temporal resolution but low spatial resolution (7 stations in 648.5 square kilometers, i.e., one every 92 square kilometers). In addition, we considered Land Surface Temperature (LST) MODIS maps: these maps have good spatial resolution but present some voids and low temporal resolution. From the first evaluations made, the temperatures reported by MODIS are not always in excellent agreement with the ones from monitoring stations. To evaluate MODIS as a data source, we decided to use Kriging spatio-temporal interpolation. Starting from multitemporal MODIS data matrices, we interpolate them to estimate the temperature of the weather stations, in order to compare the estimation with the real weather station data, thus allowing the validation of MODIS data. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the Study Area in the province of Pavia and the Lombardy Region (<b>a</b>), with a focus on the Digital Terrain Model (DTM) to highlight altitude variations across the Oltrepò Pavese (<b>b</b>).</p>
Full article ">Figure 2
<p>ARPA Weather Stations distribution in the Lombardy Region.</p>
Full article ">Figure 3
<p>ARPA weather stations within the study area and its vicinity, along with the corresponding sensor IDs.</p>
Full article ">Figure 4
<p>Series of the average hourly air temperature measured by the Fortunago Monitoring Station—ID 8007.</p>
Full article ">Figure 5
<p>MODIS data in the study area: raster with 1 km<sup>2</sup> tiles depicting the daily average daytime temperature, captured by the Aqua satellite on 30 March 2018.</p>
Full article ">Figure 6
<p>MODIS data in the study area: raster with 1 km<sup>2</sup> tiles depicting the 8-day average daytime temperature, captured by the Aqua satellite from 30 March to 6 April 2018.</p>
Full article ">Figure 7
<p>Workflow diagram illustrating the sequential process for the analysis of MODIS Land Surface Temperature (LST) and ARPA Monitoring Stations’ air temperature (Tair) data, to estimate Tair from satellite data.</p>
Full article ">Figure 8
<p>Excerpt of the stacked structure of the MODIS raster images, the <span class="html-italic">white</span> pixels represent an example of missing values (NaN).</p>
Full article ">Figure 9
<p>Example of horizontal fill for missing values. In order to retrieve the value of the selected no-value pixel, we performed an interpolation by computing the mean of the pixels within a 5 × 5 window centered on the red pixel.</p>
Full article ">Figure 10
<p>Example of vertical fill for missing values: after performing the horizontal fill, the selected pixel still has no value. To retrieve its value, we perform an interpolation considering the column of 3 layers of pixels aligned with the red one and computing the mean.</p>
Full article ">Figure 11
<p>Comparison between daytime and nighttime temperatures detected from MODIS Aqua and Terra satellites, interpolated at the Fortunago ARPA Monitoring Station coordinates (ID 8007).</p>
Full article ">Figure 12
<p>Time-series of the 8-day average air temperature of the ARPA Monitoring Station in Fortunago, and of the LST detected by the MODIS sensors, interpolated at the location coordinates.</p>
Full article ">Figure 13
<p>Comparison between MODIS-derived (dotted line) and ARPA-measured (continuous line) Tair at the coordinates of Fortunago Monitoring Station (ID 8007), while the gray line depicts their difference (ΔT).</p>
Full article ">Figure 14
<p>Data images related to the 8-day data 20/07/2022–28/07/2022. (<b>a</b>) From the combination of the four 8-day MODIS LST images (Aqua day, Aqua night, Terra day and Terra night) with their respective coefficients—deduced from the linear regression—we obtained (<b>b</b>) the estimated air temperature.</p>
Full article ">Figure A1
<p>Comparison between daytime and nighttime temperatures detected from MODIS Aqua and Terra satellites, interpolated to the Varzi Nivione ARPA Monitoring Station coordinates (ID 2082).</p>
Full article ">Figure A2
<p>Comparison between daytime and nighttime temperatures detected from MODIS Aqua and Terra satellites, interpolated to the Varzi v. Mazzini ARPA Monitoring Station coordinates (ID 8002).</p>
Full article ">Figure A3
<p>Comparison between daytime and nighttime temperatures detected from MODIS Aqua and Terra satellites, interpolated to the Voghera v. Cambiaso ARPA Monitoring Station coordinates (ID 8191).</p>
Full article ">Figure A4
<p>Comparison between daytime and nighttime temperatures detected from MODIS Aqua and Terra satellites, interpolated to the Santa Margherita di Staffora Casanova ARPA Monitoring Station coordinates (ID 8202).</p>
Full article ">Figure A5
<p>Comparison between daytime and nighttime temperatures detected from MODIS Aqua and Terra satellites, interpolated to the Canevino ARPA Monitoring Station coordinates (ID 9019).</p>
Full article ">Figure A6
<p>Comparison between daytime and nighttime temperatures detected from MODIS Aqua and Terra satellites, interpolated to the Broni ARPA Monitoring Station coordinates (ID 17432).</p>
Full article ">Figure A7
<p>Time-series of the 8-day average air temperature of the ARPA Monitoring Station in Varzi Nivione, and of the LST detected by the MODIS sensors, interpolated at the location coordinates.</p>
Full article ">Figure A8
<p>Time-series of the 8-day average air temperature of the ARPA Monitoring Station in Varzi v. Mazzini, and of the LST detected by the MODIS sensors, interpolated at the location coordinates.</p>
Full article ">Figure A9
<p>Time-series of the 8-day average air temperature of the ARPA Monitoring Station in Voghera v. Cambiaso, and of the LST detected by the MODIS sensors, interpolated at the location coordinates.</p>
Full article ">Figure A10
<p>Time-series of the 8-day average air temperature of the ARPA Monitoring Station in Santa Margherita di Staffora Casanova, and of the LST detected by the MODIS sensors, interpolated at the location coordinates.</p>
Full article ">Figure A11
<p>Time-series of the 8-day average air temperature of the ARPA Monitoring Station in Canevino, and of the LST detected by the MODIS sensors, interpolated at the location coordinates.</p>
Full article ">Figure A12
<p>Time-series of the 8-day average air temperature of the ARPA Monitoring Station in Broni, and of the LST detected by the MODIS sensors, interpolated at the location coordinates.</p>
Full article ">Figure A13
<p>Comparison between MODIS-derived (dotted line) and ARPA-measured (continuous line) Tair at the coordinates of Varzi Nivione Monitoring Station (ID 2082), while the gray line depicts their difference (ΔT).</p>
Full article ">Figure A14
<p>Comparison between MODIS-derived (dotted line) and ARPA-measured (continuous line) Tair at the coordinates of Varzi via Mazzini Monitoring Station (ID 8002), while the gray line depicts their difference (ΔT).</p>
Full article ">Figure A15
<p>Comparison between MODIS-derived (dotted line) and ARPA-measured (continuous line) Tair at the coordinates of Voghera via Cambiaso Monitoring Station (ID 8191), while the gray line depicts their difference (ΔT).</p>
Full article ">Figure A16
<p>Comparison between MODIS-derived (dotted line) and ARPA-measured (continuous line) Tair at the coordinates of Santa Margherita di Staffora Casanova Monitoring Station (ID 8202), while the gray line depicts their difference (ΔT).</p>
Full article ">Figure A17
<p>Comparison between MODIS-derived (dotted line) and ARPA-measured (continuous line) Tair at the coordinates of Canevino Monitoring Station (ID 9019), while the gray line depicts their difference (ΔT).</p>
Full article ">Figure A18
<p>Comparison between MODIS-derived (dotted line) and ARPA-measured (continuous line) Tair at the coordinates of Broni Monitoring Station (ID 17432), while the gray line depicts their difference (ΔT).</p>
Full article ">
16 pages, 8606 KiB  
Article
Annual Cropping Intensity Dynamics in China from 2001 to 2023
by Jie Ren, Yang Shao and Yufei Wang
Remote Sens. 2024, 16(24), 4801; https://doi.org/10.3390/rs16244801 - 23 Dec 2024
Viewed by 621
Abstract
Spatial and temporal information about cropping patterns of single and multiple crops is important for monitoring crop production and land-use intensity. We used time-series MODIS NDVI 8-day composite data to develop annual cropping pattern products at a 250 m spatial resolution for China, [...] Read more.
Spatial and temporal information about cropping patterns of single and multiple crops is important for monitoring crop production and land-use intensity. We used time-series MODIS NDVI 8-day composite data to develop annual cropping pattern products at a 250 m spatial resolution for China, covering the period from 2001 to 2023. To address the potential impacts of varying parameters in both data pre-processing and the peak detection algorithm on the accuracy of cropping pattern mapping, we employed a grid-search method to fine-tune these parameters. This process focused on optimizing the Savitzky–Golay smoothing window size and the peak width parameters using a calibration dataset. The results highlighted that an optimal combination of a five to seven MODIS composite window size in Savitzky–Golay smoothing and a peak width of four MODIS composites achieved good overall mapping accuracy. Pixel-wise accuracy assessments were conducted for the selected mapping years of 2001, 2011, and 2021. Overall accuracies were between 89.7% and 92.0%, with F1 scores ranging from 0.921 to 0.943. Nationally, this study observed a fluctuating trend in multiple cropping percentages, with a notable increase after 2013, suggesting shifts toward more intensive agricultural practices in recent years. At a finer spatial scale, the combination of Mann–Kendall and Sen’s slope analyses revealed that approximately 12.9% of 3 km analytical windows exhibited significant changes in cropping intensity. We observed spatial clusters of increasing and decreasing crop intensity trends across provinces such as Hebei, Shandong, Shaanxi, and Gansu. This study underscores the importance of data smoothing and peak detection methods in analyzing high temporal resolution remote sensing data. The generation of annual single/multiple cropping pattern maps at a 250 m spatial resolution enhances our comprehension of agricultural dynamics through time and across different regions. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>MODIS tiles (represented by dashed-line polygons) spanning horizontal zones 23 to 29 and vertical zones three to seven. Cropland distribution at a 250 m resolution. The 10 m land cover map products were employed to determine the percentage of croplands.</p>
Full article ">Figure 2
<p>Illustration of time-series NDVI data in 2023 for single crop (<b>a</b>) and double crop (<b>b</b>). Peak width at half-prominence is highlighted in red.</p>
Full article ">Figure 3
<p>Data and workflow for cropping intensity mapping.</p>
Full article ">Figure 4
<p>Comparison of F1 scores across different peak detection methods, examining the impact of various Savitzky–Golay (SG) smoothing combinations, SG window sizes, and peak width parameters.</p>
Full article ">Figure 5
<p>MODIS-derived map of single and multiple cropping patterns for the year 2023. The subplots show the details of the cropping intensity within major agricultural regions in China and their locations are indicated by red dots: (<b>a</b>) Northeast China Plain, (<b>b</b>) Qinghai-Tibet Plateau, (<b>c</b>) North China Plain, (<b>d</b>) Yangtze Plain, and (<b>e</b>) Southern China. Only the 2023 map is presented here for simplicity.</p>
Full article ">Figure 6
<p>The percentages of multiple crops within all cropland from the year 2001 to 2023.</p>
Full article ">Figure 7
<p>Slope coefficient (change rate) of cropping intensity trend model (2001–2023). Note only 3 km windows showing significant (<span class="html-italic">p</span> &lt; 0.05) upward/downward trends based on the Mann–Kendall test were used for trend model development.</p>
Full article ">Figure 8
<p>NDVI time series from 2001 to 2023 for two selected 3 km analytical windows: (<b>a</b>) areas transitioning from multiple to single cropping practices, and (<b>b</b>) areas transitioning from single to multiple cropping practices. NDVI values were averaged for all cropland pixels within each 3 km window for every MODIS composite.</p>
Full article ">
21 pages, 4999 KiB  
Article
Assessment of Spatio-Temporal Dynamic Vegetation Evolution and Its Driving Mechanism on the Kubuqi Desert Using Multi-Source Satellite Remote Sensing
by Linjiang Nan, Mingxiang Yang, Hejia Wang, Ping Miao, Hongli Ma, Hao Wang and Xinhua Zhang
Remote Sens. 2024, 16(24), 4769; https://doi.org/10.3390/rs16244769 - 21 Dec 2024
Viewed by 505
Abstract
Desert vegetation is undergoing complex and diverse changes due to global climate change and human activities. To effectively utilize water resources and promote ecological recovery in desert areas, it is necessary to clarify the main driving mechanisms of vegetation growth in these regions. [...] Read more.
Desert vegetation is undergoing complex and diverse changes due to global climate change and human activities. To effectively utilize water resources and promote ecological recovery in desert areas, it is necessary to clarify the main driving mechanisms of vegetation growth in these regions. In this study, based on MODIS and Landsat 8 remote sensing image data, the vegetation changes and driving mechanisms before and after water diversion in the Kubuqi Desert from 2001 to 2020 were quantitatively analyzed using multiple linear regression, random forest, support vector machine, and deep neural network. The results show that the average NDVI in the study area has increased from 0.08 to 0.13 over the past 20 years, and the year of NDVI mutation corresponded with the lowest precipitation, which occurred in 2010. After the water diversion, under the combined influence of human and natural factors, NDVI increased steadily without any abrupt changes, indicating that water is the main limiting factor for vegetation growth. The change of NDVI also showed obvious spatial heterogeneity, among which the improvement of the southwest irrigation area was the most significant, and the area with NDVI above 0.1 showed an expanding trend, and the maximum value exceeded 0.4. This demonstrates that moderate water diversion can reduce desert areas, expand lake areas, and promote vegetation growth, yielding positive ecological effects. The integration of multiple linear regression, support vector machines, random forests, and deep neural network methods effectively reveals the driving mechanisms of NDVI and indirectly informs future water diversion intervals. Overall, these research results can provide a reliable reference for the efficient development of water diversion projects and have high application value. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of the study area.</p>
Full article ">Figure 2
<p>Flowchart of the methods.</p>
Full article ">Figure 3
<p>Interannual variation trend of meteorological elements in the ecological management area from 2001 to 2020. ((<b>A</b>) indicates the variation trend of Pre and SSD; (<b>B</b>) indicates the variation trend of Tem and Wind; (<b>C</b>) indicates the variation trend of Tem-max and Tem-min; (<b>D</b>) indicates the variation trend of RHU).</p>
Full article ">Figure 4
<p>NDVI time series based on the Sen trend line.</p>
Full article ">Figure 5
<p>Results of the Pettitt mutation test.</p>
Full article ">Figure 6
<p>NDVI spatial distribution before and after water diversion in the ecological management area.</p>
Full article ">Figure 7
<p>Image display and spectral attribute ((<b>A</b>–<b>C</b>) represent true and false color composition and the spectral attribute of the lake in 2013, respectively; (<b>D</b>–<b>F</b>) are for Grassland in 2014, and (<b>G</b>–<b>I</b>) are for desert in 2019).</p>
Full article ">Figure 8
<p>Spatial change map of land use in the ecological management area.</p>
Full article ">
20 pages, 4146 KiB  
Article
Prospects for Drought Detection and Monitoring Using Long-Term Vegetation Indices Series from Satellite Data in Kazakhstan
by Irina Vitkovskaya, Madina Batyrbayeva, Nurmaganbet Berdigulov and Damira Mombekova
Land 2024, 13(12), 2225; https://doi.org/10.3390/land13122225 - 19 Dec 2024
Viewed by 498
Abstract
The rainfed cereal growing regions of Northern Kazakhstan experience significant yield fluctuations due to dependence on weather conditions. Early detection and monitoring of droughts is crucial for effective mitigation strategies in this region. This study emphasises the following objectives: (1) description of the [...] Read more.
The rainfed cereal growing regions of Northern Kazakhstan experience significant yield fluctuations due to dependence on weather conditions. Early detection and monitoring of droughts is crucial for effective mitigation strategies in this region. This study emphasises the following objectives: (1) description of the current vegetation condition with a possible separation of short-term weather effects and (2) analysing trends of changes with their directionality and quantification. Terra MODIS satellite images from 2000 to 2023 are used. Differential indices—Normalised Difference Vegetation Index (NDVI) and Vegetation Condition Index (VCI)—are used to determine the characteristics of each current season. A key component is the comparison of the current NDVI values with historical maximum, minimum, and average values to identify early indicators of drought. NDVI deviations from multiyear norms and VCI values below 0.3 visually reflect changing vegetation conditions influenced by seasonal weather patterns. The results show that the algorithm effectively detects early signs of drought through observed deviations in NDVI values, showing a trend towards increasing drought frequency and intensity in Northern Kazakhstan. The algorithm was particularly effective in detecting severe drought seasons in advance, as was the case in June 2010 and May 2012, thus supporting early recognition of drought onset. The Integrated Vegetation Index (IVI) and Integrated Vegetation Condition Index (IVCI) time series are used for integrated multiyear assessments, in analysing temporal changes in vegetation cover, determining trends in these changes, and ranking the weather conditions of each growing season in the multiyear series. Areas with high probability of drought based on low IVCI values are mapped. The present study emphasises the value of remote sensing as a tool for drought monitoring, offering timely and spatially detailed information on vulnerable areas. This approach provides critical information for agricultural planning, environmental management and policy making, especially in arid and semi-arid regions. The study emphasises the importance of multiyear data series for accurate drought forecasting and suggests that this methodology can be adapted to other drought-sensitive regions. Emphasising the socio-economic benefits, this study suggests that the early detection of drought using satellite data can reduce material losses and facilitate targeted responses. Full article
Show Figures

Figure 1

Figure 1
<p>Study area.</p>
Full article ">Figure 2
<p>Technological scheme for the formation of a series of vegetation indices.</p>
Full article ">Figure 3
<p>Changes in the HTI coefficient and the IVCI (2000–2023).</p>
Full article ">Figure 4
<p>Changes in the IVCI and average grain yield (2000–2023).</p>
Full article ">Figure 5
<p>NDVI distributions for different weather years, Akmola region.</p>
Full article ">Figure 6
<p>NDVI distributions in different weather conditions of vegetation seasons (Akkol district, Akmola region).</p>
Full article ">Figure 7
<p>Dynamics of the changes in the areas of zones of different productivity determined by IVInorm values for 2000–2023. (<b>A</b>) Location of zones of different productivity determined by IVInorm values. (<b>B</b>) Areas of zones with IVInorm values 0–0.1 and 0.1–0.2. (<b>C</b>) Areas of zones with IVInorm values 0.2–0.3. (<b>D</b>) Areas of zones with IVInorm values 0.3–0.4 and 0.4–1.</p>
Full article ">Figure 8
<p>Changes in the IVCI for the northern regions of Kazakhstan.</p>
Full article ">Figure 9
<p>Frequency of droughts on the territory of Kazakhstan, calculated from remote sensing data for April–September in 2000–2023.</p>
Full article ">
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 820
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 ">
Back to TopTop