Nothing Special   »   [go: up one dir, main page]

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
remove_circle_outline

Search Results (392)

Search Parameters:
Keywords = LandTrendR

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
24 pages, 7921 KiB  
Article
Comprehensive Comparison and Validation of Forest Disturbance Monitoring Algorithms Based on Landsat Time Series in China
by Yunjian Liang, Rong Shang, Jing M. Chen, Xudong Lin, Peng Li, Ziyi Yang, Lingyun Fan, Shengwei Xu, Yingzheng Lin and Yao Chen
Remote Sens. 2025, 17(4), 680; https://doi.org/10.3390/rs17040680 - 17 Feb 2025
Viewed by 130
Abstract
Accurate long-term and high-resolution forest disturbance monitoring are pivotal for forest carbon modeling and forest management. Many algorithms have been developed for this purpose based on the Landsat time series, but their nationwide performance across different regions and disturbance types remains unexplored. Here, [...] Read more.
Accurate long-term and high-resolution forest disturbance monitoring are pivotal for forest carbon modeling and forest management. Many algorithms have been developed for this purpose based on the Landsat time series, but their nationwide performance across different regions and disturbance types remains unexplored. Here, we conducted a comprehensive comparison and validation of six widely used forest disturbance- monitoring algorithms using 12,328 reference samples in China. The algorithms included three annual-scale (VCT, LandTrendr, mLandTrendr) and three daily-scale (BFAST, CCDC, COLD) algorithms. Results indicated that COLD achieved the highest accuracy, with F1 and F2 scores of 81.81% and 81.25%, respectively. Among annual-scale algorithms, mLandTrendr exhibited the best performance, with F1 and F2 scores of 73.04% and 72.71%, and even outperformed the daily-scale BFAST algorithm. Across China’s six regions, COLD consistently achieved the highest F1 and F2 scores, showcasing its robustness and adaptability. However, regional variations in accuracy were observed, with the northern region exhibiting the highest accuracy and the southwestern region the lowest. When considering different forest disturbance types, COLD achieved the highest accuracies for Fire, Harvest, and Other disturbances, while CCDC was most accurate for Forestation. These findings highlight the necessity of region-specific calibration and parameter optimization tailored to specific disturbance types to improve forest disturbance monitoring accuracy, and also provide a solid foundation for future studies on algorithm modifications and ensembles. Full article
Show Figures

Figure 1

Figure 1
<p>The study area of 12,328 reference forest disturbance samples and two regions. (<b>a</b>) The frequency of forest disturbance; (<b>b</b>) The type of forest disturbance. NE: Northeast China; N: North China; NW: Northwest China; E: East China; S: South China; SW: Southwest China.</p>
Full article ">Figure 2
<p>The statistics of reference samples with different forest disturbance types in the six regions of China. NE: Northeast China; N: North China; NW: Northwest China; E: East China; S: South China; SW: Southwest China.</p>
Full article ">Figure 3
<p>Sensitivity analysis for determining the optimal thresholds of key parameters in the VCT algorithm across the entire China: (<b>a</b>) Compositing periods; (<b>b</b>) forThrMax; (<b>c</b>) minNdvi. Five points in each line from left to right represent using the maximum Z-score of 2, 3, 4, 5, and 6, respectively. The dotted line represents the 1:1 line.</p>
Full article ">Figure 4
<p>Example of forest disturbance monitoring using the six algorithms. (<b>a</b>) VCT; (<b>b</b>) LandTrendr; (<b>c</b>) mLandTrendr; (<b>d</b>) BFAST; (<b>e</b>) CCDC; (<b>f</b>) COLD. IFZ: Integrated Forest Z-score; NBR: Normalized Burn Ratio. The vertical dotted line represents the disturbance date.</p>
Full article ">Figure 5
<p>Sensitivity analysis for determining the optimal thresholds of key parameters in the LandTrendr algorithm across the entire China: (<b>a</b>) bestModelProportion; (<b>b</b>) Compositing periods; (<b>c</b>) Indices; (<b>d</b>) maxSegments; (<b>e</b>) pvalThreshold; (<b>f</b>) recoveryThreshold. Five points in each line from left to right represent using the spikeThreshold of 0.6, 0.75, 0.85, 0.9, and 1, respectively. The dotted line represents the 1:1 line.</p>
Full article ">Figure 6
<p>Sensitivity analysis for determining the optimal thresholds of the key parameters in the mLandTrendr algorithm across the entire China: (<b>a</b>,<b>b</b>) Index Combinations, 1–5 represent NBR, NDMI, TCW, NDVI, and TCA; (<b>c</b>) Tn. Five points in each line from left to right represent using the Tc of 0.9, 0.95, 0.99, 0.999, and 0.9999, respectively. The dotted line represents the 1:1 line.</p>
Full article ">Figure 7
<p>Sensitivity analysis for determining the optimal thresholds of the key parameters in the BFAST algorithm across the entire China: (<b>a</b>) h; (<b>b</b>) harmonics; (<b>c</b>) Indices; (<b>d</b>) period. Five points in each line from left to right represent using the alpha of 0.05, 0.025, 0.01, 0.005, and 0.001, respectively. The dotted line represents the 1:1 line.</p>
Full article ">Figure 8
<p>Sensitivity analysis for determining the optimal thresholds of the key parameters in the CCDC algorithm across the entire China: (<b>a</b>) breakpointBands; (<b>b</b>) lambda; (<b>c</b>) maxIterations; (<b>d</b>) minNumOfYearsScaler; (<b>e</b>) minObservations; (<b>f</b>) tmaskBands. Five points in each line from left to right represent using the chi-square distribution thresholds of 0.9, 0.95, 0.99, 0.999, and 0.9999, respectively. The dotted line represents the 1:1 line.</p>
Full article ">Figure 9
<p>Sensitivity analysis for determining the optimal thresholds of the five key parameters in the COLD algorithm across the entire China: (<b>a</b>) conObservations; (<b>b</b>) detectBands; (<b>c</b>) minNumOfYearsScaler; (<b>d</b>) nsign; (<b>e</b>) tmaskBands. Five points in each line from left to right represent using the chi-square distribution thresholds of 0.9, 0.95, 0.99, 0.999, and 0.9999, respectively. The dotted line represents the 1:1 line.</p>
Full article ">Figure 10
<p>The parameter calibration results for the combination of all optimal thresholds determined through sensitivity analyses for the six algorithms across the entire China. For VCT, five points in each line from left to right represent using the maximum Z-score of 2, 3, 4, 5, and 6, respectively. For LandTrendr, five points in each line from left to right represent using the spikeThreshold of 0.6, 0.75, 0.85, 0.9, and 1, respectively. For mLandTrendr, five points in each line from left to right represent using the Tc of 0.9, 0.95, 0.99, 0.999, and 0.9999, respectively. For BFAST, five points in each line from left to right represent using the alpha of 0.05, 0.025, 0.01, 0.005, and 0.001, respectively. For CCDC and COLD, five points in each line from left to right represent using the chi-square distribution thresholds of 0.9, 0.95, 0.99, 0.999, and 0.9999, respectively. The dotted line represents the 1:1 line.</p>
Full article ">Figure 11
<p>Validations of six forest disturbance monitoring algorithms in China.</p>
Full article ">Figure 12
<p>Validations of six forest disturbance monitoring algorithms in the six regions of China.</p>
Full article ">Figure 13
<p>Validations of six algorithms in monitoring different types of forest disturbance in China.</p>
Full article ">Figure 14
<p>Regional example of forest disturbance monitoring using the six algorithms in northeastern China. The first column is Landsat 8 False Color Composited images (R: SWIR1, G: NIR, B: R) on 22 May 2015 and 8 May 2016. The second to last columns are maps of monitored forest disturbance (red colors) by the VCT, LandTrendr, mLandTrendr, BFAST, CCDC, and COLD algorithms. The second and fourth rows are the enlarged views of the blue rectangle in the first and third rows.</p>
Full article ">Figure 15
<p>Regional example of forest disturbance monitoring using the six algorithms in southern China. The first column is Landsat 8 False Color Composited images (R: SWIR1, G: NIR, B: R) on 10 November 2010 and 13 November 2011. The second to last columns are maps of monitored forest disturbance (red colors) by the VCT, LandTrendr, mLandTrendr, BFAST, CCDC, and COLD algorithms. The second and fourth rows are the enlarged views of the blue rectangle in the first and third rows.</p>
Full article ">
30 pages, 3465 KiB  
Article
Weather-Driven Cycling: Developing a Predictive Model for Urban Bicycle Usage Based on Five Key Weather Factors
by Nahid Falah, Nadia Falah and Jaime Solis-Guzman
Urban Sci. 2025, 9(2), 41; https://doi.org/10.3390/urbansci9020041 - 11 Feb 2025
Viewed by 416
Abstract
Weather conditions significantly influence urban cycling, shaping both its frequency and intensity. This study develops a predictive model to evaluate the impact of five key meteorological factors, namely temperature, humidity, precipitation, wind speed, and daylight duration, on urban cycling trends. Using non-linear regression [...] Read more.
Weather conditions significantly influence urban cycling, shaping both its frequency and intensity. This study develops a predictive model to evaluate the impact of five key meteorological factors, namely temperature, humidity, precipitation, wind speed, and daylight duration, on urban cycling trends. Using non-linear regression analysis, the research examines cycling data from 2017 to 2019 in Hamburg, Germany, comparing predicted values for 2019 with actual data to assess model accuracy. The statistical analyses reveal strong correlations between weather parameters and cycling activity, highlighting each factor’s unique influence. The model achieved high accuracy, with R2 values of 0.942 and 0.924 for 2017 and 2019, respectively. To further validate its robustness, the model is applied to data from 2021 and 2023—years not included in its initial development—yielding R2 values of 0.893 and 0.919. These results underscore the model’s reliability and adaptability across different timeframes. This study not only confirms the critical influence of weather on urban cycling patterns, but also provides a scalable framework for broader urban planning applications. Beyond the immediate findings, this research proposes expanding the model to incorporate urban factors, such as land use, population density, and socioeconomic conditions, offering a comprehensive tool for urban planners and policymakers to enhance sustainable transportation systems. Full article
Show Figures

Figure 1

Figure 1
<p>Topography of Hamburg [<a href="#B98-urbansci-09-00041" class="html-bibr">98</a>].</p>
Full article ">Figure 2
<p>The current state of Hamburg’s bicycle infrastructure [<a href="#B98-urbansci-09-00041" class="html-bibr">98</a>].</p>
Full article ">Figure 3
<p>Comparison of the weather parameters (mean air temperature °C, mean humidity %, precipitation mm, daylight minute, mean wind speed kph) in 2017, 2019, 2021, and 2023 and the bicycle trip volume during the same years.</p>
Full article ">Figure 4
<p>Methodological framework illustrating data collection and model validation steps.</p>
Full article ">Figure 5
<p>Average monthly variability in the bicycle volume, monthly mean daylight hours, monthly mean precipitation, monthly mean humidity, and monthly mean air temperature in 2017.</p>
Full article ">Figure 6
<p>The significance of the relationship between the number of real cycling trips and the number of predicted cycling trips in 2017.</p>
Full article ">Figure 7
<p>Comparison of predicted vs. actual bicycle trips in 2019.</p>
Full article ">
17 pages, 9559 KiB  
Article
Vegetation Carbon Source/Sink Dynamics and Extreme Climate Response in the Yangtze River Delta Coastal Zone
by Yuhang Han and Zhen Han
Sustainability 2025, 17(4), 1456; https://doi.org/10.3390/su17041456 - 11 Feb 2025
Viewed by 394
Abstract
Coastal zones, as transition areas for sea/land interaction, have substantial carbon sequestration potential while also being particularly vulnerable to extreme climate. Consequently, it has become essential to evaluate the vegetation carbon sinks in coastal zone areas under extreme climate conditions. In this study, [...] Read more.
Coastal zones, as transition areas for sea/land interaction, have substantial carbon sequestration potential while also being particularly vulnerable to extreme climate. Consequently, it has become essential to evaluate the vegetation carbon sinks in coastal zone areas under extreme climate conditions. In this study, we evaluated the vegetation net ecosystem productivity (NEP) in typical regions within the Yangtze River Delta coastal zone from 2000 to 2020. We studied the regional and chronological properties of NEP and its response to extreme climate. The results revealed the following: (1) Vegetation NEP demonstrated a fluctuating rising trend over the past 21 years, with an interannual change rate of 1.96 gC·m−2·a−1, and the 21-year average was 249.22 gC·m−2·a−1. Spatially, the southern part of the region had a higher NEP than the northern part, and the northern part had a higher NEP than the central part. (2) The overall area showed characteristics of a vegetation carbon sink, with carbon sink areas accounting for 82.41%. Among the ecosystems, forest ecosystems exhibited the strongest carbon sink capacity, followed by cropland ecosystems, while wetland ecosystems, urban ecosystems, and grassland ecosystems had relatively weaker carbon sink capacities. (3) The overall spatial change trend showed an upward trend, consistent with the temporal trend. There is also a high risk of vegetation NEP degradation in the future. (4) The NEP’s response to extreme temperature was more pronounced. The largest explanatory power was observed with SU25 and TMAX during single-factor analysis. The strongest explanatory power in the interaction analysis was found in the following three factor groups: R99p∩TMAX, SU25∩TNx, and TXx∩LST. The results highlight a complex synergistic interplay among these influences on NEP. The findings offer a scientific basis for ecological protection and the attainment of dual-carbon goals in the coastal zone of the Yangtze River Delta. Full article
(This article belongs to the Special Issue Ecology, Environment, and Watershed Management)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Overview of the study area; (<b>b</b>) Ecosystem distribution.</p>
Full article ">Figure 2
<p>Validation of NPP estimation results.</p>
Full article ">Figure 3
<p>Interannual variation in vegetation NEP in different ecosystems: (<b>a</b>) Study area; (<b>b</b>) Cropland; (<b>c</b>) Forest; (<b>d</b>) Grassland; (<b>e</b>) Wetland; (<b>f</b>) Urban.</p>
Full article ">Figure 4
<p>(<b>a</b>) Spatial distribution of vegetation NEP multi-year averages; (<b>b</b>) Carbon source/sink areas.</p>
Full article ">Figure 5
<p>(<b>a</b>) Trend of vegetation carbon source/sink; (<b>b</b>) Significance of trend.</p>
Full article ">Figure 6
<p>(<b>a</b>) Hurst index; (<b>b</b>) Future trend of vegetation carbon source/sink.</p>
Full article ">Figure 7
<p>Interaction detection results. Note: Lower triangles: q-value magnitude; Upper triangles: blue slashes indicate bivariable enhancement, and red circles indicate nonlinear enhancement.</p>
Full article ">Figure 8
<p>Comparison of NPP results in 2020 (<b>a</b>) Shanghai NPP and NPP(Sentinel-2); (<b>b</b>) Nantong NPP and NPP(Sentinel-2); (<b>c</b>) Shanghai NPP and NPP(Sentinel-2) modification; (<b>d</b>) Nantong NPP and NPP(Sentinel-2) modification.</p>
Full article ">
24 pages, 12137 KiB  
Article
Spatiotemporal Changes of Vegetation Growth and Its Influencing Factors in the Huojitu Mining Area from 1999 to 2023 Based on kNDVI
by Zhichao Chen, Yiqiang Cheng, Xufei Zhang, Zhenyao Zhu, Shidong Wang, Hebing Zhang, Youfeng Zou and Chengyuan Hao
Remote Sens. 2025, 17(3), 536; https://doi.org/10.3390/rs17030536 - 5 Feb 2025
Viewed by 394
Abstract
Vegetation indices are important representatives of plant growth. Climate change and human activities seriously affect vegetation. This study focuses on the Huojitu mining area in the Shendong region, utilizing the kNDVI index calculated via the Google Earth Engine (GEE) cloud platform. The Mann–Kendall [...] Read more.
Vegetation indices are important representatives of plant growth. Climate change and human activities seriously affect vegetation. This study focuses on the Huojitu mining area in the Shendong region, utilizing the kNDVI index calculated via the Google Earth Engine (GEE) cloud platform. The Mann–Kendall mutation test and linear regression analysis were employed to examine the spatiotemporal changes in vegetation growth over a 25-year period from 1999 to 2023. Through correlation analysis, geographic detector models, and land use map fusion, combined with climate, topography, soil, mining, and land use data, this study investigates the influencing factors of vegetation growth evolution. The key findings are as follows: (1) kNDVI is more suitable for analyzing vegetation growth in this study compared to NDVI. (2) Over the past 25 years, vegetation growth has exhibited an overall fluctuating upward trend, with an annual growth rate of 0.0041/a. The annual average kNDVI value in the mining area is 0.121. Specifically, kNDVI initially increased gradually, then rapidly increased, and subsequently declined rapidly. (3) Vegetation growth in the study area has significantly improved, with areas of improved vegetation accounting for 89.08% of the total mining area, while degraded areas account for 11.02%. (4) Precipitation and air temperature are the primary natural factors influencing vegetation growth fluctuations in the mining area, with precipitation being the dominant factor (r = 0.81, p < 0.01). The spatial heterogeneity of vegetation growth is influenced by land use, topography, soil nutrients, and mining activities, with land use having the greatest impact (q = 0.43). Major land use changes contribute 46.45% to vegetation improvement and 13.43% to vegetation degradation. The findings of this study provide a scientific basis for ecological planning and the development of the Huojitu mining area. Full article
Show Figures

Figure 1

Figure 1
<p>The map of the geographic location of study area.</p>
Full article ">Figure 2
<p>Layout of the sampling points in the study area.</p>
Full article ">Figure 3
<p>Spatial interpolation results of soil data.</p>
Full article ">Figure 4
<p>Spatial distribution of topographic data: (<b>a</b>) DEM, (<b>b</b>) slope, and (<b>c</b>) topographic position index (TPI).</p>
Full article ">Figure 5
<p>Spatial distribution of mining thickness.</p>
Full article ">Figure 6
<p>The spatial distribution of land use types in the Huojitu mining area in 1999 and 2023: (<b>a</b>) land use types in 1999, and (<b>b</b>) land use types in 2023.</p>
Full article ">Figure 7
<p>kNDVI and NDVI histograms.</p>
Full article ">Figure 8
<p>Spatial distribution of changes in NDVI and kNDVI from 1999 to 2023.</p>
Full article ">Figure 9
<p>Remote sensing image data acquired in the study area: (<b>a</b>) RGB Orthophoto Map, (<b>b</b>) kNDVI, and (<b>c</b>) NDVI.</p>
Full article ">Figure 10
<p>kNDVI seasonal changes in the Huojitu mining area from 1999 to 2023.</p>
Full article ">Figure 11
<p>Time variation of kNDVI in the Huojitu mining area from 1999 to 2023: (<b>a</b>) MK mutation test, and (<b>b</b>) interannual variation.</p>
Full article ">Figure 12
<p>kNDVI spatial distribution in different time periods in the Huojitu mining area.</p>
Full article ">Figure 13
<p>Interannual variation trend and fluctuation of vegetation growth in the Huojitu mining area from 1999 to 2023: (<b>a</b>) slope value of kNDVI change, (<b>b</b>) change trend of vegetation growth, and (<b>c</b>) Significant interannual change trend of vegetation growth.</p>
Full article ">Figure 14
<p>Changes in kNDVI, mean temperature, and precipitation in the Huojitu mining area.</p>
Full article ">Figure 15
<p>Explanatory power of each influence factor.</p>
Full article ">Figure 16
<p>Chord diagram of land use type Changes in Huojitu mining area from 1999 to 2023.</p>
Full article ">Figure 17
<p>Land use change and its impact on vegetation growth between 1999 and 2023: (<b>a</b>) map of main land use transformation types, and (<b>b</b>) effects of land use change on vegetation growth.</p>
Full article ">
20 pages, 12491 KiB  
Article
Forest Disturbance and Restoration in China's North-South Transition Zone: A Case from the Funiu Mountains
by Qifan Wu, Jiacheng Hou, Shiwen Wu, Fuyuan Su, Shilong Hao, Tailai Yin, Haoyuan Chen, Yunpeng Xu and Hailong He
Forests 2025, 16(2), 269; https://doi.org/10.3390/f16020269 - 4 Feb 2025
Viewed by 516
Abstract
Accurate monitoring and assessment of forest disturbance and recovery dynamics are essential for sustainable forest management, particularly in ecological transition zones. This study analyzed forest disturbance and recovery patterns in China’s Funiu Mountains from 1991 to 2020 by integrating the LandTrendr algorithm with [...] Read more.
Accurate monitoring and assessment of forest disturbance and recovery dynamics are essential for sustainable forest management, particularly in ecological transition zones. This study analyzed forest disturbance and recovery patterns in China’s Funiu Mountains from 1991 to 2020 by integrating the LandTrendr algorithm with space-time cube analysis. Using Landsat time series data and the Geodetector method, we examined both the spatiotemporal characteristics and driving factors of forest change across three periods. The results showed that (1) between 1991 and 2020, the study area experienced 131.19 km2 of forest disturbance and 495.88 km2 of recovery, with both processes most active during the 1990s; (2) spatiotemporal analysis revealed that both disturbance and recovery patterns were predominantly characterized by cold spots, suggesting relatively stable forest conditions despite localized changes; (3) human activities were the primary drivers of forest disturbance in the early period, while forest recovery was consistently influenced by the combined effects of topographic conditions and precipitation. Additionally, forest fires emerged as an important factor affecting both disturbance and recovery patterns after 2010. These findings enhance our understanding of forest dynamics in transition zones and provide empirical support for regional forest management strategies. The results also highlight the importance of considering both spatial and temporal dimensions when monitoring long-term forest changes. Full article
(This article belongs to the Special Issue Mapping and Modeling Forests Using Geospatial Technologies)
Show Figures

Figure 1

Figure 1
<p>Geographical location of the study area: (<b>a</b>) Location of the study area in China. (<b>b</b>) Location in Henan Province. (<b>c</b>) Topographic map of the Funiu Mountains.</p>
Full article ">Figure 2
<p>Methodological framework for analyzing the spatiotemporal patterns and driving factors of forest disturbance and recovery: (<b>a</b>) data preprocessing and forest area extraction, (<b>b</b>) forest change detection, (<b>c</b>) spatiotemporal pattern analysis, and (<b>d</b>) driving factor identification.</p>
Full article ">Figure 3
<p>Recognition of forest change by different spectral indices.</p>
Full article ">Figure 4
<p>Space-time cube of forest disturbance and recovery: (<b>a</b>) is the forest disturbance cube, and (<b>b</b>) is the forest recovery space-time cube.</p>
Full article ">Figure 5
<p>Area of forest disturbance and recovery per year and cumulative area.</p>
Full article ">Figure 6
<p>Spatial distribution of local outliers: (<b>a</b>) forest disturbance and (<b>b</b>) forest recovery.</p>
Full article ">Figure 7
<p>Statistical distribution of local outlier patterns for forest disturbance and recovery.</p>
Full article ">Figure 8
<p>Map of an emerging spatiotemporal hotspot analysis of forest disturbance and recovery, with (<b>a</b>) showing cold spots and hotspot patterns of forest disturbance and (<b>b</b>) showing cold spots and hotspot patterns of forest recovery.</p>
Full article ">Figure 9
<p>Factor detection for different periods: (<b>a</b>, <b>c</b>, <b>e</b>) are the factor detections for forest disturbance, while (<b>b</b>, <b>d</b>, <b>f</b>) are the factor detections for forest recovery.</p>
Full article ">Figure 10
<p>Heatmap of the factor interactions for forest disturbances (<b>a</b>–<b>c</b>) and recovery (<b>d</b>–<b>f</b>) across three periods (1991–2000, 2001–2010, and 2011–2020). Abbreviations: ELEV (elevation), SLP (slope), MAP (mean annual precipitation), MAT (mean annual temperature), ST (soil type), LUI (land use intensity), RD (road density), PD (population density), GDP (gross domestic product), and FO (forest fires).</p>
Full article ">Figure 11
<p>Distribution of forest disturbance at different elevations (<b>a</b>) and slopes (<b>b</b>) from 1991 to 2000.</p>
Full article ">Figure 12
<p>Distribution of forest recovery at different elevations (<b>a</b>) and slopes (<b>b</b>) from 1991 to 2020.</p>
Full article ">Figure 13
<p>Spatial distribution of forest disturbances: (<b>a</b>,<b>c</b>) represent disturbances after 2000, while (<b>b</b>,<b>d</b>) represent disturbances before 2000.</p>
Full article ">
22 pages, 11614 KiB  
Article
Analysis of the Spatial–Temporal Characteristics of Vegetation Cover Changes in the Loess Plateau from 1995 to 2020
by Zhihong Yao, Yichao Huang, Yiwen Zhang, Qinke Yang, Peng Jiao and Menghao Yang
Land 2025, 14(2), 303; https://doi.org/10.3390/land14020303 - 1 Feb 2025
Viewed by 542
Abstract
The Loess Plateau is one of the most severely affected regions by soil erosion in the world, with a fragile ecological environment. Vegetation plays a key role in the region’s ecological restoration and protection. This study employs the Geographical Detector (Geodetector) model to [...] Read more.
The Loess Plateau is one of the most severely affected regions by soil erosion in the world, with a fragile ecological environment. Vegetation plays a key role in the region’s ecological restoration and protection. This study employs the Geographical Detector (Geodetector) model to quantitatively assess the impact of natural and human factors, such as temperature, precipitation, soil type, and land use, on vegetation growth. It aims to reveal the characteristics and driving mechanisms of vegetation cover changes on the Loess Plateau over the past 26 years. The results indicate that from 1995 to 2020, the vegetation coverage on the Loess Plateau shows an increasing trend, with a fitted slope of 0.01021 and an R2 of 0.96466. The Geodetector indicates that the factors with the greatest impact on vegetation cover in the Loess Plateau are temperature, precipitation, soil type, and land use. The highest average vegetation coverage is achieved when the temperature is between −4.8 and 2 °C or 12 and 16 °C, precipitation is between 630.64 and 935.51 mm, the soil type is leaching soil, and the land use type is forest. And the interaction between all factors has a greater effect on the vegetation cover than any single factor alone. This study reveals the factors influencing vegetation growth on the Loess Plateau, as well as their types and ranges, providing a scientific basis and guidance for improving vegetation coverage in this region. Full article
Show Figures

Figure 1

Figure 1
<p>Map of the Loess Plateau geographic location.</p>
Full article ">Figure 2
<p>The long−term average precipitation and temperature values of the Loess Plateau: (<b>a</b>) temperature; (<b>b</b>) precipitation.</p>
Full article ">Figure 3
<p>Monthly average NDVI from 2001 to 2015.</p>
Full article ">Figure 4
<p>Spatial distributions of natural and human factors in 2020: (<b>a</b>) slope; (<b>b</b>) aspect; (<b>c</b>) temperature; (<b>d</b>) precipitation; (<b>e</b>) soil type; (<b>f</b>) land use type; (<b>g</b>) population density; and (<b>h</b>) GDP.</p>
Full article ">Figure 5
<p>The principle of geographical detector.</p>
Full article ">Figure 6
<p>Annual mean FVC changes in the Loess Plateau from 1995 to 2020.</p>
Full article ">Figure 7
<p>Trend of vegetation coverage change from 1995 to 2020, using the Mann–Kendall test.</p>
Full article ">Figure 8
<p>Average FVC value for each precipitation zone in 1995, 2000, 2010, and 2020.</p>
Full article ">Figure 9
<p>Average FVC value for each temperature zone in 1995, 2000, 2010, and 2020.</p>
Full article ">Figure 10
<p>Average FVC under different soil types in 1995, 2000, 2010, and 2020.</p>
Full article ">Figure 11
<p>Average FVC under different land use types in 1995, 2000, 2010, and 2020.</p>
Full article ">Figure 12
<p><span class="html-italic">q</span> value for detection of interaction effects of various factors in 1995, 2000, 2010, and 2020: (<b>a</b>) 1995; (<b>b</b>) 2000; (<b>c</b>) 2010; (<b>d</b>) 2020.</p>
Full article ">Figure 12 Cont.
<p><span class="html-italic">q</span> value for detection of interaction effects of various factors in 1995, 2000, 2010, and 2020: (<b>a</b>) 1995; (<b>b</b>) 2000; (<b>c</b>) 2010; (<b>d</b>) 2020.</p>
Full article ">
25 pages, 4735 KiB  
Article
Remote Sensing Detection of Forest Changes in the South Ridge Corridor and an Attribution Analysis
by Nan Wu, Linghui Huang, Meng Zhang, Yaqing Dou, Kehan Mo and Junang Liu
Forests 2025, 16(2), 205; https://doi.org/10.3390/f16020205 - 23 Jan 2025
Viewed by 673
Abstract
As the largest mountain range in Southern China, the natural vegetation of Nanling plays an irreplaceable role in maintaining the stability of the ecosystem and exerting its functions. The forested area of the Nanling Corridor encompasses 168,633 km2, with a forest [...] Read more.
As the largest mountain range in Southern China, the natural vegetation of Nanling plays an irreplaceable role in maintaining the stability of the ecosystem and exerting its functions. The forested area of the Nanling Corridor encompasses 168,633 km2, with a forest coverage rate exceeding 60% of all cities together. Long-term analysis of the temporal and spatial evolution of this forest and the disturbance factors in this region is of great importance for realizing the “dual carbon” goals, sustainable forest management, and protecting biodiversity. In this study, remote sensing images from a Landsat time series with a resolution of 30 m were obtained from the GEE (Google Earth Engine) cloud processing platform, and forest disturbance data were obtained using the LandTrendr algorithm. Using a machine learning random forest algorithm, the forest disturbance status and disturbance factors were explored from 2001 to 2020. The results show that the estimated disturbed forest area from 2001 to 2020 was 11,904.3 km2, accounting for 7.06% of the total area of the 11 cities in the Nanling Corridor, and the average annual disturbed area was 595.22 km2. From 2001 to 2016, the overall disturbed area increased, reaching a peak value of 1553.36 km2 in 2008, with a low value of 37.71 km2 in 2002. After 2016, the disturbed area showed a downward trend. In this study, an attribution analysis of forest disturbance factors was carried out. The results showed that the overall accuracy of forest disturbance factor attribution was as high as 82.48%, and the Kappa coefficient was 0.70. Among the disturbance factors, deforestation factors accounted for 58.45% of the total area of forest disturbance, followed by fire factors (28.69%) and building or road factors (12.85%). The regional distribution of each factor also had significant characteristics, and the Cutdown factors were mostly distributed in the lower elevations of the mountain margin, with most of them distributed in sheets. The fire factors were spatially distributed in the center of the mountains, and their distribution was loose. Building or road factors were mostly distributed in clusters or lines. These research results are expected to provide technical and data support for the study of the large-scale spatiotemporal evolution of forests and its driving mechanisms. Full article
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The location of the study area.</p>
Full article ">Figure 2
<p>Flow chart of forest disturbance monitoring based on LandTrendr.</p>
Full article ">Figure 3
<p>The corresponding schematic diagram of the NDVI values and the changes in the high-resolution Google images.</p>
Full article ">Figure 4
<p>Spatiotemporal distribution pattern of forest disturbance in Nanling Corridor from 2001 to 2020. A, B, and C are comparisons between Google imagery of 2020 and forest disturbance monitoring results from 2001 to 2020 in three detailed areas of the study area, respectively.</p>
Full article ">Figure 5
<p>Plot of disturbed forest area from 2001 to 2020.</p>
Full article ">Figure 6
<p>Columnar map of forest disturbance area from 2001 to 2020. (<b>a</b>) shows an image of Google Earth in 2020 for a typical area of the study area; (<b>b</b>,<b>c</b>) are the results of forest disturbances monitored in 2001 and 2020, respectively, for the corresponding locations.</p>
Full article ">Figure 7
<p>Percentages of forest disturbance factor areas from 2001 to 2020.</p>
Full article ">Figure 8
<p>Distribution of forest disturbance factors. A, B and C are comparisons between Google Images and the results of the 2020 forest disturbance attribution monitoring in three detailed regions of the study area, respectively.</p>
Full article ">Figure 9
<p>Disturbance due to forest Cutdown factors from 2001 to 2020. A, B, and C, are enlarged images of three typical areas, which are marked with corresponding letters in the panel below.</p>
Full article ">Figure 10
<p>Disturbance due to fire factors from 2001 to 2020. A, B, and C are enlarged images of three typical areas, which are marked with corresponding letters in the panel below.</p>
Full article ">Figure 11
<p>Disturbance due to building or road factors from 2001 to 2020. A, B, and C are enlarged images of three typical areas, which are marked with corresponding letters in the panel below.</p>
Full article ">
18 pages, 4891 KiB  
Article
Monitoring Forest Disturbances and Associated Driving Forces in Guangdong Province Using Long-Term Landsat Time Series Images
by Lin Qiu, Zhongbing Chang, Xiaomei Luo, Songjia Chen, Jun Jiang and Li Lei
Forests 2025, 16(1), 189; https://doi.org/10.3390/f16010189 - 20 Jan 2025
Viewed by 731
Abstract
Research on monitoring forest disturbances and analyzing its driving factors is crucial for the sustainable management of forest ecosystems. To quantitatively identify the spatial distribution and dynamic changes of forest disturbance and its driving factors in Guangdong Province from 1990 to 2019, the [...] Read more.
Research on monitoring forest disturbances and analyzing its driving factors is crucial for the sustainable management of forest ecosystems. To quantitatively identify the spatial distribution and dynamic changes of forest disturbance and its driving factors in Guangdong Province from 1990 to 2019, the long-term Landsat time series imagery and the LandTrendr change detection algorithm were utilized. The impact of forest disturbances on four types of landscape fragmentation (attrition, perforation, shrinkage, and subdivision) was analyzed using the Forman index. The Geodetector model was used to analyze the driving factors of forest disturbance from human activity and the natural environment. The results showed that the LandTrendr algorithm achieved a Kappa coefficient of 0.79, with an overall accuracy of approximately 82.59%. The findings indicate a consistent increase in shrinkage patches, both in quantity and area. Spatially, the centroids of forest fragmentation processes exhibited a clear inland migration trend, reflecting the growing ecological pressures faced by inland forest ecosystems. Furthermore, interactions among driving factors, particularly between population density and economic factors, significantly amplified their combined impacts. The correlation between forest disturbances and socio-economic factors revealed distinct regional variations, highlighting significant differences in forest disturbance dynamics across cities with varying levels of economic development. This study provides critical insights into the spatiotemporal dynamics of forest disturbances under rapid urbanization and economic development. It lays the groundwork for sustainable forest management strategies in Guangdong Province and may contribute to global discussions on managing forest ecosystems during periods of rapid socio-economic transformation. Full article
(This article belongs to the Section Urban Forestry)
Show Figures

Figure 1

Figure 1
<p>Location of the study area and spatial distribution of forests from GlobeLand30 (2020).</p>
Full article ">Figure 2
<p>Research technology roadmap.</p>
Full article ">Figure 3
<p>Construction process of forest subdivision process model (The red square is an example of an eight-neighborhood).</p>
Full article ">Figure 4
<p>Disturbance results for three typical areas ((<b>a</b>–<b>c</b>) were three representative areas).</p>
Full article ">Figure 5
<p>Result of centroid analysis in the spatial process of forest subdivision.</p>
Full article ">Figure 6
<p>Correlation coefficients between the area of forest disturbance and various factors.</p>
Full article ">
19 pages, 8430 KiB  
Article
Spatiotemporal Variation of Water Use Efficiency and Its Responses to Climate Change in the Yellow River Basin from 1982 to 2018
by Jie Li, Fen Qin, Yingping Wang, Xiuyan Zhao, Mengxiao Yu, Songjia Chen, Jun Jiang, Linhua Wang and Junhua Yan
Remote Sens. 2025, 17(2), 316; https://doi.org/10.3390/rs17020316 - 17 Jan 2025
Viewed by 562
Abstract
The ecosystem water use efficiency (WUE) plays a critical role in many aspects of the global carbon cycle, water management, and ecological services. However, the response mechanisms and driving processes of WUE need to be further studied. This research was conducted based on [...] Read more.
The ecosystem water use efficiency (WUE) plays a critical role in many aspects of the global carbon cycle, water management, and ecological services. However, the response mechanisms and driving processes of WUE need to be further studied. This research was conducted based on Gross Primary Productivity (GPP), Evapotranspiration (ET), meteorological station data, and land use/cover data, and the methods of Ensemble Empirical Mode Decomposition (EEMD), trend variation analysis, the Mann–Kendall Significant Test (M-K test), and Partial Correlation Analysis (PCA) methods. Our study revealed the spatio-temporal trend of WUE and its influencing mechanism in the Yellow River Basin (YRB) and compared the differences in WUE change before and after the implementation of the Returned Farmland to Forestry and Grassland Project in 2000. The results show that (1) the WUE of the YRB showed a significant increase trend at a rate of 0.56 × 10−2 gC·kg−1·H2O·a−1 (p < 0.05) from 1982 to 2018. The area showing a significant increase in WUE (47.07%, Slope > 0, p < 0.05) was higher than the area with a significant decrease (14.64%, Slope < 0, p < 0.05). The region of significant increase in WUE in 2000–2018 (45.35%, Slope > 0, p < 0.05) was higher than that of 1982–2000 (8.23%, Slope > 0, p < 0.05), which was 37.12% higher in comparison. (2) Forest WUE (1.267 gC·kg−1·H2O) > Cropland WUE (0.972 gC·kg−1·H2O) > Grassland WUE (0.805 gC·kg−1·H2O) under different land cover types. Forest ecosystem WUE has the highest rate of increase (0.79 × 10−2 gC·kg−1·H2O·a−1) from 2000 to 2018. Forest ecosystem WUE increased by 0.082 gC·kg−1·H2O after 2000. (3) precipitation (37.98%, R > 0, p < 0.05) and SM (10.30%, R > 0, p < 0.05) are the main climatic factors affecting WUE in the YRB. A total of 70.39% of the WUE exhibited an increasing trend, which is mainly attributed to the simultaneous increase in GPP and ET, and the rate of increasing GPP is higher than the rate of increasing ET. This study could provide a scientific reference for policy decision-making on the terrestrial carbon cycle and biodiversity conservation. Full article
(This article belongs to the Special Issue Remote Sensing for Terrestrial Hydrologic Variables)
Show Figures

Figure 1

Figure 1
<p>Study area, vegetation type, basin boundary, and elevation.</p>
Full article ">Figure 2
<p>Temporal trends of WUE in 1982–2018. (<b>a</b>) Annual; (<b>b</b>) Grow.</p>
Full article ">Figure 3
<p>Spatial variation characteristics of WUE in the YRB. (<b>a</b>) Annual WUE in 1982–2018; (<b>b</b>) Annual WUE in 1982–2000; (<b>c</b>) Annual WUE in 2000–2018; (<b>d</b>) Grow WUE in 1982–2018; (<b>e</b>) Annual WUE in 1982–2000; (<b>f</b>) Annual WUE in 2000–2018.</p>
Full article ">Figure 4
<p>Spatial characteristics of significant variation trend of the WUE in different time periods.</p>
Full article ">Figure 5
<p>Variation in WUE in different land cover types.</p>
Full article ">Figure 6
<p>The trends of WUE for different ecosystem types in the YRB. (<b>a</b>) Farmland; (<b>b</b>) Forest; (<b>c</b>) Grassland; (<b>d</b>) Other.</p>
Full article ">Figure 7
<p>Spatial distribution of partial correlation coefficient between WUE and climate factors in the YRB from 1982 to 2018. (<b>a</b>) WUE and temperature; (<b>b</b>) WUE and precipitation; (<b>c</b>) WUE and vapor pressure deficit; (<b>d</b>) WUE and soil moisture.</p>
Full article ">Figure 8
<p>Spatial distribution of partial correlation coefficient between WUE and climate factors in the YRB from 1982 to 2018 (significance test <span class="html-italic">p</span> &lt; 0.05). (<b>a</b>) WUE and temperature; (<b>b</b>) WUE and precipitation; (<b>c</b>) WUE and vapor pressure deficit; (<b>d</b>) WUE and soil moisture.</p>
Full article ">Figure 9
<p>WUE changes in response to GPP and ET across different time periods.</p>
Full article ">Figure 10
<p>WUE significant changes in response to GPP significant changes and ET significant changes across different time periods (significant test <span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">
16 pages, 5513 KiB  
Technical Note
Identifying Optimal Variables to Predict Soil Organic Carbon in Sandy, Saline, and Black Soil Regions: Remote Sensing, Terrain, or Climate Factors?
by Liping Wang, Huanjun Liu, Xiang Wang, Xiaofeng Xu, Liyuan He, Chong Luo, Yong Li, Xinle Zhang, Deqiang Zang, Shufeng Zheng and Xiaodan Mei
Remote Sens. 2025, 17(2), 237; https://doi.org/10.3390/rs17020237 - 10 Jan 2025
Viewed by 606
Abstract
Environmental variables have a substantial effect on the reliability of soil organic carbon (SOC) mapping. However, it is still challenging to identify which environmental variables are effective in cropland SOC prediction in sandy, saline, and black soil regions. To address this issue, we [...] Read more.
Environmental variables have a substantial effect on the reliability of soil organic carbon (SOC) mapping. However, it is still challenging to identify which environmental variables are effective in cropland SOC prediction in sandy, saline, and black soil regions. To address this issue, we used the principal component analysis (PCA) method for the feature selection of bands, spectral indexes, and terrain factors for each region. Based on the selection feature, we used global RF and local RF for SOC prediction for these three regions. Our results indicated that (1) climate factors, particularly mean annual precipitation and mean annual temperature, were the most effective predictors in SOC mapping across sandy, saline, and black soil regions, as indicated by their significant contribution to RF model performance (R2 > 0.63); (2) followed by climate factors, the Transformed Vegetation Index (TVI) was consistently identified as the most influential variable for SOC prediction among spectral indexes in all three regions; (3) a local regression method based on RF models showed good performance compared to a global model; (4) desertification and salinization were the main reasons for the spatial differences in AH and DM&LD, respectively. The SOC of HL in black soil regions was consistent with the climate change trend because of the latitude difference. This study provides valuable information for constructing a more precise soil prediction strategy for cultivated land in sandy, saline, and black soil regions. Full article
(This article belongs to the Section Remote Sensing for Geospatial Science)
Show Figures

Figure 1

Figure 1
<p>Locations of the study area: (<b>a</b>) location of the study area; (<b>b</b>) map of sandy region and sampling sites in AH; (<b>c</b>) image of sandy soil in AH; (<b>d</b>) map of saline soil and sampling sites in DM&amp;LD; (<b>e</b>) image of soil salinization in DM&amp;LD; (<b>f</b>) map of black soil region and sampling sites in HL; (<b>g</b>) image of black soil in HL.</p>
Full article ">Figure 2
<p>Original reflectance (OR, at the top of each figure) and continuum removal (CR, at the bottom of each figure) with different SOC contents. (<b>a</b>) AH, sandy region; (<b>b</b>) DM&amp;LD, saline region; (<b>c</b>) HL, black soil region.</p>
Full article ">Figure 3
<p>The weighting values of bands, spectral indexes, and terrain factors ((<b>a</b>): AH; (<b>b</b>): DM&amp;LD; (<b>c</b>): HL; (<b>d</b>): all regions). The yellow points mean two selected variables with the highest values in each part of the study areas (Bn (The nth bands of Sentinel-2), Green Normalized Difference Vegetation Index (GNDVI), Enhanced Vegetation Index (EVI), Soil Adjusted Total Vegetation Index (SATVI), Transformed Vegetation Index (TVI), Ratio Vegetation Index (RVI), Green Ratio Vegetation Index (GRVI), Land Surface Water Index (LSWI), Moisture Stress Index (MSI), Soil Adjusted Vegetation Index (SAVI), Normalized Differences Vegetation Index (NDVI), slope (S), aspect (A), plan curvatures (PlC), profile curvatures (PrC), topographic wetness index (TWI), roughness (Rn), relief (RL), slope length (SL), and hillshade (HS)).</p>
Full article ">Figure 4
<p>Spatial map of precipitation and temperature in the three different regions.</p>
Full article ">Figure 5
<p>Training results for each region and all regions based on RF model. ((<b>a</b>). Sandy soil area in AH. (<b>b</b>). Saline soil area in DM&amp;LD. (<b>c</b>). Black soil region in HL. (<b>d</b>). All regions.)</p>
Full article ">Figure 6
<p>Validation results of each region and all regions based on RF model. ((<b>a</b>). Sandy soil area in AH. (<b>b</b>). Saline soil area in DM&amp;LD. (<b>c</b>). Black soil region in HL. (<b>d</b>). All regions.).</p>
Full article ">Figure 7
<p>Training and validation results using the local regression method based on RF.</p>
Full article ">Figure 8
<p>Spatial distribution of predicted SOC content in cultivated land in RF model ((<b>a</b>): SOC of AH; (<b>b</b>): SOC of DM&amp;LD; (<b>c</b>): SOC of HL).</p>
Full article ">
19 pages, 7445 KiB  
Article
An Interpretable Model for Salinity Inversion Assessment of the South Bank of the Yellow River Based on Optuna Hyperparameter Optimization and XGBoost
by Xia Liu, Yu Hu, Xiang Li, Ruiqi Du, Youzhen Xiang and Fucang Zhang
Agronomy 2025, 15(1), 18; https://doi.org/10.3390/agronomy15010018 - 26 Dec 2024
Viewed by 457
Abstract
Soil salinization is a serious land degradation phenomenon, posing a severe threat to regional agricultural resource utilization and sustainable development. It has been a mainstream trend to use machine-learning methods to achieve monitoring of large-scale salinized soil quickly. However, machine learning model training [...] Read more.
Soil salinization is a serious land degradation phenomenon, posing a severe threat to regional agricultural resource utilization and sustainable development. It has been a mainstream trend to use machine-learning methods to achieve monitoring of large-scale salinized soil quickly. However, machine learning model training requires many samples and hyper-parameter optimization and lacks solvability. To compare the performance of different machine-learning models, this study conducted a soil sampling experiment on saline soils along the south bank of the Yellow River in Dalate Banner. The experiment lasted two years (2022 and 2023) during the spring bare soil period, collecting 304 soil samples. The soil salinity was estimated with the multi-source remote sensing satellite data by combining the extreme gradient boosting model (XGBoost), Optuna hyper-parameter optimization, and Shapley addition (SHAP) interpretable model. Correlation analysis and continuous variable projection were employed to identify key inversion factors. The regression effects of partial least squares regression (PLSR), geographically weighted regression (GWR), long short-term memory networks (LSTM), and extreme gradient boosting (XGBoost) were compared. The optimal model was selected to estimate soil salinity in the study area from 2019 to 2023. The results showed that the XGBoost model fitted optimally, the test set had high R2 (0.76) and the ratio of performance to deviation (2.05), and the estimation results were consistent with the measured salinity values. SHAP analysis revealed that the salinity index and topographic factors were the primary inversion factors. Notably, the same inversion factor influenced varying soil salinity estimates at different locations. The saline soils of the study area in 2019 and 2023 were 65% and 44%, respectively, and the overall trend of soil salinization decreased. From the viewpoint of spatial distribution, the degree of soil salinization showed a gradually increasing trend from south to north, and it was most serious on the side near the Yellow River. This study is of great significance for the quantitative estimation of salinized soil in the irrigated area on the south bank of the Yellow River, the prevention and control of soil salinization, and the sustainable development of agriculture. Full article
(This article belongs to the Section Precision and Digital Agriculture)
Show Figures

Figure 1

Figure 1
<p>Overview of the study area and distribution of sampling sites ((<b>a</b>) geographic location map; (<b>b</b>,<b>c</b>) elevation map, distribution of sampling sites).</p>
Full article ">Figure 2
<p>(<b>a</b>) Descriptive statistics box plot of measured SSC (CV: coefficient of variation); (<b>b</b>) SSC distribution map of different types of saline–alkali soils.</p>
Full article ">Figure 3
<p>Correlation analysis between soil salinity and inversion factors.</p>
Full article ">Figure 4
<p>Performance comparison of different model training sets.</p>
Full article ">Figure 5
<p>Performance comparison of different model test sets.</p>
Full article ">Figure 6
<p>Soil content grading chart 2019–2023.</p>
Full article ">Figure 7
<p>Area of different classes of saline land from 2019 to 2023 ((<b>a</b>) change in area; (<b>b</b>) rate of change).</p>
Full article ">Figure 8
<p>Transfer matrix between areas of different types of saline soils.</p>
Full article ">Figure 9
<p>(<b>a</b>) SHAP global interpretation map: feature summary map for SHAP; (<b>b</b>) heat map of SHAP-based features.</p>
Full article ">Figure 10
<p>SSC inversion data-processing procedure. Table (<b>a</b>) in the figure shows the 8th data point, and table (<b>b</b>) shows the 15th.</p>
Full article ">
20 pages, 8060 KiB  
Article
The Wheel of Vegetation: A Spatial and Temporal Story of Vegetation Evolution in the Shennongjia Forest District
by Xueli Wang, Xiaolong Du, Chunyan Zhao, An Luo, Hui Chen, Shaobin Li and Hewei Du
Forests 2024, 15(12), 2252; https://doi.org/10.3390/f15122252 - 22 Dec 2024
Viewed by 517
Abstract
As one of the most well-preserved areas in the vertical band spectrum of vegetation in central China and even in the northern hemisphere at the same latitude, the vegetation in Shennongjia Forest District is vital to global ecological balance. In order to fully [...] Read more.
As one of the most well-preserved areas in the vertical band spectrum of vegetation in central China and even in the northern hemisphere at the same latitude, the vegetation in Shennongjia Forest District is vital to global ecological balance. In order to fully understand the vegetation change in the study area, remotely sensed data since 1990, combined with the Sen-MK test, Geo detector, and LandTrendr algorithm, were used to analyze the vegetation distribution characteristics and change trends. The results showed that: (1) the overall NDVI in the study area displayed an upward trend. (2) Vegetation disturbance occurred frequently before 2000 and decreased significantly after 2000. The most severely disturbed year was 1991 when the disturbed area amounted to 4.0851 km2, accounting for 16.76% of the total disturbed area. The analysis of the topographic environment reveals that most of the vegetation disturbances occur in areas with slopes of 15–25° and elevations of 1500–2000 m, which indicates that these areas have frequent human activity. (3) The explanatory power of different influences on vegetation changes varied, with altitude having the most significant effect and the superposition of two influences increasing the effect on vegetation change. Over the past 30 years, vegetation in the Shennongjia Forest District has shown a general trend of recovery, with natural forest protection initiatives playing a critical role in mitigating disturbance. This comprehensive study of vegetation changes in Shennongjia offers a valuable research paradigm for forest conservation and sustainable development in temperate forests at similar latitudes, providing significant insights into the protection and management of similar ecosystems. Full article
(This article belongs to the Section Forest Ecology and Management)
Show Figures

Figure 1

Figure 1
<p>Location of the study area. (<b>a</b>) The location of Hubei Province in China, (<b>b</b>) the location of the Shennongjia Forest District in Hubei Province, and (<b>c</b>) the elevation of the Shennongjia Forest District.</p>
Full article ">Figure 2
<p>Methodology flowchart.</p>
Full article ">Figure 3
<p>Time series fitting curve. (<b>a</b>) NDVI; (<b>b</b>) NBR.</p>
Full article ">Figure 4
<p>The LandTrendr segmentation process, adapted with permission from ref. [<a href="#B23-forests-15-02252" class="html-bibr">23</a>].</p>
Full article ">Figure 5
<p>Annual variation of mean NDVI in the study area from 1990 to 2020.</p>
Full article ">Figure 6
<p>(<b>a</b>) NDVI change trends; (<b>b</b>) types of changes in NDVI.</p>
Full article ">Figure 7
<p>(<b>a</b>–<b>c</b>) are Landsat images; (<b>d</b>–<b>f</b>) are the disturbance years detected by the LandTrendr algorithm; (<b>g</b>–<b>i</b>) are for fitting NBR disturbance trajectories. The red markers on the left correspond to the green markers in the center, and each row is a set of images showing sampling points for low, medium, and high vegetation cover, respectively. Note: the sampling points marked in the figure are also the points used to fit the perturbation trajectories.</p>
Full article ">Figure 8
<p>(<b>a</b>) Year of disturbance; (<b>b</b>) inter-annual variation characteristics of forest disturbance in Shennongjia Forest District from 1990 to 2020.</p>
Full article ">Figure 9
<p>Magnitude of disturbance.</p>
Full article ">Figure 10
<p>Area of forest disturbance at different elevations and slopes in Shennongjia Forest District: (<b>a</b>) slope and (<b>b</b>) elevation.</p>
Full article ">Figure 11
<p>Average annual temperature and average annual precipitation in Shennongjia Forest District.</p>
Full article ">Figure 12
<p>Factor detector. Note: from left to right, 2010, 2011, and 2020. X1, X2, X3, X4, and X5 represent annual average temperature, average annual precipitation, slope, elevation, and population density, respectively.</p>
Full article ">Figure 13
<p>Interaction detector. Note: from left to right, 2010, 2011, and 2020. X1, X2, X3, X4, and X5 represent annual average temperature, average annual precipitation, slope, elevation, and population density, respectively.</p>
Full article ">
32 pages, 10269 KiB  
Article
Impact of Ridge Tillage and Mulching on Water Dynamics of Summer Maize Fields Under Climate Change in the Semi-Arid Region of Northwestern Liaoning, China
by Yao Li, Wanting Zhang, Mengxi Bai, Jiayu Wu, Chenmengyuan Zhu and Yujuan Fu
Agronomy 2024, 14(12), 3032; https://doi.org/10.3390/agronomy14123032 - 19 Dec 2024
Viewed by 578
Abstract
The ridge–furrow plastic mulching technique has been widely applied due to its benefits of increasing temperature, conserving moisture, reducing evaporation, and boosting yields. Hydrus-2D is a computer model designed to simulate the two-dimensional movement of water in soil characterized by a low cost [...] Read more.
The ridge–furrow plastic mulching technique has been widely applied due to its benefits of increasing temperature, conserving moisture, reducing evaporation, and boosting yields. Hydrus-2D is a computer model designed to simulate the two-dimensional movement of water in soil characterized by a low cost and high flexibility compared to field experiments. This study, based on field experiment data from Jianping County, Liaoning Province, China, during 2017–2018, developed Hydrus-2D models for two distinct field management practices: non-mulched flat cultivation (NM-FC) and mulched ridge tillage (M-RT). Furthermore, it simulated the dynamic changes in farmland water variations during the summer maize growth period (2021–2100) under climate change scenarios, specifically medium and high emission pathways (SSP2-4.5 and SSP5-8.5), based on the FGOALS-g3 model, which exhibits the highest similarity to the climate pattern of Jianping County in the Coupled Model Intercomparison Project Phase 6 (CMIP6) global climate models and the Shared Socioeconomic Pathways (SSPs). The results showed that in the future FGOALS-g3 model, net radiation exhibited a significant upward trend under the SSP2-4.5 scenario (Z = 2.38), while the average air temperature showed a highly significant increase under both SSP2-4.5 and SSP5-8.5 scenarios, with Z-values of 6.48 and 8.90, respectively. The Hydrus-2D model demonstrated high simulation accuracy in both NM-FC and M-RT treatments (R2 ranging from 0.86 to 0.96, with RMSE not exceeding 0.011), accurately simulating the dynamic changes in soil water content (SWC) under future climate change. Compared to NM-FC, M-RT reduced evaporation, increased transpiration, and effectively decreased the leakage caused by increased future precipitation, resulting in a 0.04 and 0.01 cm3/cm3 increase in surface and deep soil SWC, respectively, during the summer maize growing season, significantly improving water use efficiency. Moreover, M-RT treatment reduced the impact coefficients of climate change on various water balance parameters, stabilizing changes in these parameters and SWC under future climate conditions. This study demonstrates the significant advantages of M-RT in coping with climate change, providing key scientific evidence for future agricultural water resource management. These findings offer valuable insights for policymakers and farmers, particularly in developing adaptive land management and irrigation strategies, helping to improve water use efficiency and promote sustainable agricultural practices. Full article
(This article belongs to the Special Issue Advances in Tillage Methods to Improve the Yield and Quality of Crops)
Show Figures

Figure 1

Figure 1
<p>Overview map of the study area’s geographic location.</p>
Full article ">Figure 2
<p>Schematic diagram of the field experiment. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage.</p>
Full article ">Figure 3
<p>Schematic diagram of boundary conditions and finite element mesh division. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage.</p>
Full article ">Figure 4
<p>Changes in future meteorological data under SSP2-4.5 and SSP5-8.5 emission scenarios for the FGOALS-g3 model. Notes: <span class="html-italic">Tair</span>, <span class="html-italic">PRE</span>, <span class="html-italic">RH</span>, and <span class="html-italic">Rn</span> represent daily mean temperature, precipitation, daily mean relative humidity, and net radiation, respectively. SSP2-4.5 and SSP5-8.5 represent Shared Socioeconomic Pathways 2–4.5 (medium forcing scenario) and Shared Socioeconomic Pathways 5–8.5 (high forcing scenario), respectively.</p>
Full article ">Figure 5
<p>Measured and simulated soil water content values at different soil depths for each treatment in 2017 and 2018. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage.</p>
Full article ">Figure 5 Cont.
<p>Measured and simulated soil water content values at different soil depths for each treatment in 2017 and 2018. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage.</p>
Full article ">Figure 6
<p>Changes in water balance under future climate conditions for the NM-FC and M-RT treatments.</p>
Full article ">Figure 6 Cont.
<p>Changes in water balance under future climate conditions for the NM-FC and M-RT treatments.</p>
Full article ">Figure 6 Cont.
<p>Changes in water balance under future climate conditions for the NM-FC and M-RT treatments.</p>
Full article ">Figure 7
<p>Changes in SWC at various depths under future climate conditions for the NM-FC and M-RT treatments. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage. SSP2-4.5 and SSP5-8.5 represent Shared Socioeconomic Pathways 2–4.5 (medium forcing scenario) and Shared Socioeconomic Pathways 5–8.5 (high forcing scenario), respectively.</p>
Full article ">Figure 7 Cont.
<p>Changes in SWC at various depths under future climate conditions for the NM-FC and M-RT treatments. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage. SSP2-4.5 and SSP5-8.5 represent Shared Socioeconomic Pathways 2–4.5 (medium forcing scenario) and Shared Socioeconomic Pathways 5–8.5 (high forcing scenario), respectively.</p>
Full article ">Figure 8
<p>Path analysis between future meteorological variables and various factors of farmland water balance. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage. SSP2-4.5 and SSP5-8.5 represent Shared Socioeconomic Pathways 2–4.5 (medium forcing scenario) and Shared Socioeconomic Pathways 5–8.5 (high forcing scenario), respectively. Tair, <span class="html-italic">PRE</span>, <span class="html-italic">RH,</span> and <span class="html-italic">Rn</span> represent daily mean temperature, precipitation, daily mean relative humidity, and net radiation, respectively. The symbols *, **, and *** indicate the significance levels of one factor’s effect on another, where * represents <span class="html-italic">p</span> &lt; 0.05 (statistically significant), ** represents <span class="html-italic">p</span> &lt; 0.01 (highly significant), and *** represents <span class="html-italic">p</span> &lt; 0.001 (extremely significant).</p>
Full article ">Figure 8 Cont.
<p>Path analysis between future meteorological variables and various factors of farmland water balance. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage. SSP2-4.5 and SSP5-8.5 represent Shared Socioeconomic Pathways 2–4.5 (medium forcing scenario) and Shared Socioeconomic Pathways 5–8.5 (high forcing scenario), respectively. Tair, <span class="html-italic">PRE</span>, <span class="html-italic">RH,</span> and <span class="html-italic">Rn</span> represent daily mean temperature, precipitation, daily mean relative humidity, and net radiation, respectively. The symbols *, **, and *** indicate the significance levels of one factor’s effect on another, where * represents <span class="html-italic">p</span> &lt; 0.05 (statistically significant), ** represents <span class="html-italic">p</span> &lt; 0.01 (highly significant), and *** represents <span class="html-italic">p</span> &lt; 0.001 (extremely significant).</p>
Full article ">Figure 8 Cont.
<p>Path analysis between future meteorological variables and various factors of farmland water balance. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage. SSP2-4.5 and SSP5-8.5 represent Shared Socioeconomic Pathways 2–4.5 (medium forcing scenario) and Shared Socioeconomic Pathways 5–8.5 (high forcing scenario), respectively. Tair, <span class="html-italic">PRE</span>, <span class="html-italic">RH,</span> and <span class="html-italic">Rn</span> represent daily mean temperature, precipitation, daily mean relative humidity, and net radiation, respectively. The symbols *, **, and *** indicate the significance levels of one factor’s effect on another, where * represents <span class="html-italic">p</span> &lt; 0.05 (statistically significant), ** represents <span class="html-italic">p</span> &lt; 0.01 (highly significant), and *** represents <span class="html-italic">p</span> &lt; 0.001 (extremely significant).</p>
Full article ">Figure A1
<p>Measured and simulated soil water content values at different soil depths for each treatment in 2017 and 2018. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage.</p>
Full article ">Figure A1 Cont.
<p>Measured and simulated soil water content values at different soil depths for each treatment in 2017 and 2018. Note: NM-FC represents non-mulched flat cultivation, while M-RT represents mulched ridge tillage.</p>
Full article ">
17 pages, 3872 KiB  
Article
Impact of Land Use Types on Soil Physico-Chemical Properties, Microbial Communities, and Their Fungistatic Effects
by Giuseppina Iacomino, Mohamed Idbella, Salvatore Gaglione, Ahmed M. Abd-ElGawad and Giuliano Bonanomi
Soil Syst. 2024, 8(4), 131; https://doi.org/10.3390/soilsystems8040131 - 16 Dec 2024
Viewed by 2221
Abstract
Soilborne plant pathogens significantly impact agroecosystem productivity, emphasizing the need for effective control methods to ensure sustainable agriculture. Soil fungistasis, the soil’s ability to inhibit fungal spore germination under optimal conditions, is pivotal for biological control. This study explores soil fungistasis variability across [...] Read more.
Soilborne plant pathogens significantly impact agroecosystem productivity, emphasizing the need for effective control methods to ensure sustainable agriculture. Soil fungistasis, the soil’s ability to inhibit fungal spore germination under optimal conditions, is pivotal for biological control. This study explores soil fungistasis variability across land-use intensities, spanning deciduous and evergreen forests, grasslands, shrublands, and horticultural cultivations in both open fields and greenhouses. Soil characterization encompassed organic matter, pH, total nitrogen, C/N ratio, key cations (Ca2+, Mg2+, K+, Na+), enzymatic activities, microbial biomass, and soil microbiota analyzed through high-throughput sequencing of 16s rRNA genes. Fungistasis was evaluated against the pathogenic fungi Botrytis cinerea and the beneficial microbe Trichoderma harzianum. Fungistasis exhibited similar trends across the two fungi. Specifically, the application of glucose to soil temporarily annulled soil fungistasis for both B. cinerea and T. harzianum. In fact, a substantial fungal growth, i.e., fungistasis relief, was observed immediately (48 h) after the pulse application with glucose. In all cases, the fungistasis relief was proportional to the glucose application rate, i.e., fungal growth was higher when the concentration of glucose was higher. However, the intensity of fungistasis relief largely varied across soil types. Our principal component analysis (PCA) demonstrated that the growth of both Trichoderma and Botrytis fungi was positively and significantly correlated with organic carbon content, total nitrogen, iron, magnesium, calcium, and sodium while negatively correlated with fluorescein diacetate (FDA) hydrolysis. Additionally, bacterial diversity and composition across different ecosystems exhibited a positive correlation with FDA hydrolysis and a negative correlation with phosphoric anhydride and soil pH. Analysis of bacterial microbiomes revealed significant differences along the land use intensity gradient, with higher fungistasis in soils dominated by Pseudoarthrobacter. Soils under intensive horticultural cultivation exhibited a prevalence of Acidobacteria and Cyanobacteria, along with reduced fungistasis. This study sheds light on soil fungistasis variability in diverse ecosystems, underscoring the roles of soil texture rather than soil organic matter and microbial biomass to explain the variability of fungistasis across landscapes. Full article
Show Figures

Figure 1

Figure 1
<p>Images of the selected ecosystems across a climatic and land use intensity gradient in terms of organic amendment input, synthetic fertilizers, and pesticide application in the Campania Region (Southern Italy). All pictures by Giuliano Bonanomi.</p>
Full article ">Figure 2
<p>Box plots illustrating the variation in species richness (<b>A</b>) and the Shannon diversity index (<b>B</b>) for bacterial communities across the ecosystem soils. The boxes represent the interquartile range (IQR), with the lower and upper bounds indicating the 25th and 75th percentiles, respectively. The horizontal line within each box marks the median, while the whiskers extend to the range of data within 1.5 times the IQR. Different letters indicate significant differences (<span class="html-italic">p</span> &lt; 0.05). (<b>C</b>) Non-metric multidimensional scaling (NMDS) plots depict bacterial community composition in the different soils. The MDS axis1 and MDS axis2 correspond to the two axes of the two-dimensional ordination space, with each point representing a replicate sample. The stress level, shown on each plot, indicates how well the distances between objects are preserved (values closer to 0 indicate a better representation of the data in the ordination space). The <span class="html-italic">p</span>- and F-values represent the results of the PERMANOVA test conducted with 999 permutations on the bacterial data. (<b>D</b>) Bar charts display the relative abundance of bacterial phyla in the different soils.</p>
Full article ">Figure 3
<p>Heatmap showing the relative abundance of the 100 most frequent Amplicon Sequence Variants in the bacterial community in the soil of each ecosystem. The grouping of variables is based on Whittaker’s association index.</p>
Full article ">Figure 4
<p>Fungal growth of <span class="html-italic">B. cinerea</span> conidia (expressed as a percentage compared to the control (0%)) on soil watery extracts from the selected ecosystems during a 168 h incubation period that followed a single application of glucose at four application rates (0.10%, 0.30%, 1%, and 3%). Values are averages ± standard deviation.</p>
Full article ">Figure 5
<p>Fungal growth of <span class="html-italic">T. harzianum</span> conidia (expressed as a percentage compared to the control (0%)) on soil watery extracts from the selected ecosystems during a 168 h incubation period that followed a single application of glucose at four application rates (0.10%, 0.30%, 1%, and 3%). Values are averages ± standard deviation.</p>
Full article ">Figure 6
<p>Principal component analysis (PCA) based on soil physico-chemical characteristics (<b>A</b>) and SIMPER resulting taxa (<b>B</b>) as variables. Microbial biomass, fungal growth, and bacterial diversity and composition were fitted as factors with significance &lt;0.05 onto the ordination.</p>
Full article ">
27 pages, 14009 KiB  
Article
Model Development for Estimating Sub-Daily Urban Air Temperature Patterns in China Using Land Surface Temperature and Auxiliary Data from 2013 to 2023
by Yuchen Guo, János Unger and Tamás Gál
Remote Sens. 2024, 16(24), 4675; https://doi.org/10.3390/rs16244675 - 14 Dec 2024
Viewed by 823
Abstract
Near-surface air temperature (Tair) is critical for addressing urban challenges in China, particularly in the context of rapid urbanization and climate change. While many studies estimate Tair at a national scale, they typically provide only daily data (e.g., maximum and minimum Tair), with [...] Read more.
Near-surface air temperature (Tair) is critical for addressing urban challenges in China, particularly in the context of rapid urbanization and climate change. While many studies estimate Tair at a national scale, they typically provide only daily data (e.g., maximum and minimum Tair), with few focusing on sub-daily urban Tair at high spatial resolution. In this study, we integrated MODIS-based land surface temperature (LST) data with 18 auxiliary data from 2013 to 2023 to develop a Tair estimation model for major Chinese cities, using random forest algorithms across four diurnal and seasonal conditions: warm daytime, warm nighttime, cold daytime, and cold nighttime. Four model schemes were constructed and compared by combining different auxiliary data (time-related and space-related) with LST. Cross-validation results were found to show that space-related and time-related variables significantly affected the model performance. When all auxiliary data were used, the model performed best, with an average RMSE of 1.6 °C (R2 = 0.96). The best performance was observed on warm nights with an RMSE of 1.47 °C (R2 = 0.97). The importance assessment indicated that LST was the most important variable across all conditions, followed by specific humidity, and convective available potential energy. Space-related variables were more important under cold conditions (or nighttime) compared with warm conditions (or daytime), while time-related variables exhibited the opposite trend and were key to improving model accuracy in summer. Finally, two samples of Tair patterns in Beijing and the Pearl River Delta region were effectively estimated. Our study offered a novel method for estimating sub-daily Tair patterns using open-source data and revealed the impacts of predictive variables on Tair estimation, which has important implications for urban thermal environment research. Full article
Show Figures

Figure 1

Figure 1
<p>Research area and the location of meteorological stations. As the reference line for the study area, the Hu Line (Heihe–Tengchong Line) is represented by the dashed line. Beijing and the PRD were chosen as sample areas for estimated Tair illustration. Their satellite images and LCZ maps are shown in detail. The LCZ type codes refer to the specific LCZ types, as follows: 1 (compact high-rise), 2 (compact mid-rise), 3 (compact low-rise), 4 (open high-rise), 5 (open mid-rise), 6 (open low-rise), 7 (lightweight low-rise), 8 (large low-rise), 9 (sparsely built), 10 (heavy industry), A (dense trees), B (scattered trees), C (bush, scrub), D (low plants), E (bare rock or paved), F (bare soil or sand), and G (water).</p>
Full article ">Figure 2
<p>The overall framework of this study. The main steps are highlighted in blue.</p>
Full article ">Figure 3
<p>The RMSE (<b>a</b>) between the predicted and measured Tair based on the tenfold cross-validation of four model schemes under four diurnal and seasonal conditions and the RMSE gaps (ΔRMSE) between Model 1 and the other three model schemes (<b>b</b>).</p>
Full article ">Figure 4
<p>Scatter plots and fitting results between observed and estimated Tair of the final model scheme (Model 4) under four diurnal and seasonal conditions based on the tenfold cross-validation. The panels (<b>a</b>–<b>d</b>) represent warm daytime, warm nighttime, cold daytime, and cold nighttime, respectively. The color of the scatter plot represents the point density.</p>
Full article ">Figure 5
<p>The relative values of VIMs were calculated for all predictor variables under four diurnal and seasonal conditions, based on the impurity-corrected method. The subfigure (<b>a</b>–<b>d</b>) represent warm daytime, warm nighttime, cold daytime, and cold nighttime, respectively.</p>
Full article ">Figure 6
<p>The annual variation in daily RMSE based on cross-validation using the entire dataset, under warm (<b>a</b>) and cold (<b>b</b>) conditions. The gray shades mask the time period with less data. Under the warm condition (<b>a</b>) less than 70% stations have usable data in gray-shaded period and this proportion is 30% during cold condition (<b>b</b>).</p>
Full article ">Figure 7
<p>The spatial distribution of RMSE at each station based on cross-validation under four diurnal and seasonal conditions. The panels (<b>a</b>–<b>d</b>) represent warm daytime, warm nighttime, cold daytime, and cold nighttime, respectively.</p>
Full article ">Figure 8
<p>Spatiotemporal patterns of the estimated air temperature under the warm (25 July 2023) condition in Beijing. The black oval on the 14:00 map highlights the heat island at the airport.</p>
Full article ">Figure 9
<p>Spatiotemporal patterns of the estimated air temperature under the cold (9 January 2018) condition in Beijing. The black ovals on the 02:00 map highlight the heat spots in the northwest mountain regions surrounding Beijing.</p>
Full article ">Figure 10
<p>Spatiotemporal patterns of the estimated air temperature under warm conditions (30 September 2019) in the PRD. The black ovals with numbers on the map highlight the major cities within the region. The cities corresponding to the numbers are as follows: 1. Dongguan (coastal areas), 2. Shenzhen (coastal areas), 3. Hong Kong, 4. Macau, 5. Guangzhou and Foshan, and 6. Dongguan (inland areas).</p>
Full article ">Figure 11
<p>The annual variation in daily RMSE based on the cross-validation of Model 4, using the limited dataset under warm conditions, with 300 samples randomly selected each day.</p>
Full article ">Figure 12
<p>The annual variation in daily RMSE based on cross-validation of three model schemes, Model 1 (<b>a</b>), Model 2 (<b>b</b>), and Model 3 (<b>c</b>). All RMSEs are computed under warm conditions using the same limited dataset as <a href="#remotesensing-16-04675-f011" class="html-fig">Figure 11</a>.</p>
Full article ">Figure 13
<p>The spatial patterns of RMSE at each station based on the cross-validation of Model 2 (<b>a</b>) and Model 3 (<b>b</b>) under warm nighttime conditions.</p>
Full article ">
Back to TopTop