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 (373)

Search Parameters:
Keywords = LandTrendR

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
18 pages, 5546 KiB  
Article
Climatological Evaluation of Three Assimilation and Reanalysis Datasets on Soil Moisture over the Tibetan Plateau
by Yinghan Sang, Hong-Li Ren and Mei Li
Remote Sens. 2024, 16(22), 4198; https://doi.org/10.3390/rs16224198 - 11 Nov 2024
Viewed by 408
Abstract
Soil moisture is critical in the linkage between the land and atmosphere of energy and water exchange, especially over the Tibetan Plateau (TP). However, due to the lack of in situ plateau soil moisture measurements, the reanalyzed and assimilated data are the major [...] Read more.
Soil moisture is critical in the linkage between the land and atmosphere of energy and water exchange, especially over the Tibetan Plateau (TP). However, due to the lack of in situ plateau soil moisture measurements, the reanalyzed and assimilated data are the major supplements for TP climate research. Based on observations from 1992 to 2013, this study provides a comprehensive evaluation of three sets of assimilation and reanalysis products (GLDAS, ERA5-Land, and MERRA-2) on the climatic mean and variability of soil moisture over the Tibetan Plateau (TPSM). For the climatic mean, GLDAS captures the spatial distribution and annual cycle of TPSM better than other datasets in terms of lower spatial RMSE (0.07 m3×m-3) and bias (0.06 m3×m-3). In terms of the climatic variability of TPSM, the multi-data average (MDA) highlights its advantages in reducing the bias relative to any single data product. MDA describes the TPSM anomalies more stably and accurately in terms of temporal trend and variation (r = 0.94), as well as the dipole spatial pattern in EOF1. When considering both the climatic mean and spatial variability, the performance of MDA is more accurate and balanced than that of a single data product. This study overcomes the deficiency of limited time and space in previous evaluations of TPSM and indicates that multi-data averaging may be a more effective approach in the climate investigation of TPSM. Full article
Show Figures

Figure 1

Figure 1
<p>Locations and names of 17 soil moisture stations in the Tibetan Plateau selected by this study. The contours indicate the altitude (units: m).</p>
Full article ">Figure 2
<p>Annual and seasonal mean TPSM (units: m3×m-3): (<b>a</b>,<b>f</b>,<b>k</b>,<b>p</b>,<b>u</b>) observations as the reference; (<b>b</b>,<b>g</b>,<b>l</b>,<b>q</b>,<b>v</b>) GLDAS; (<b>c</b>,<b>h</b>,<b>m</b>,<b>r</b>,<b>w</b>) ERA5-Land; (<b>d</b>,<b>i</b>,<b>n</b>,<b>s</b>,<b>x</b>) MERRA-2; and (<b>e</b>,<b>j</b>,<b>o</b>,<b>t</b>,<b>y</b>) MDA. The pattern correlation coefficients are shown in the bottom-left corner of each panel.</p>
Full article ">Figure 3
<p>Differences in annual and seasonal mean TPSM (units: m3×m-3) between observations and data products: (<b>a</b>,<b>e</b>,<b>i</b>,<b>m</b>,<b>q</b>) differences between GLDAS and observations; (<b>b</b>,<b>f</b>,<b>j</b>,<b>n</b>,<b>r</b>) differences between ERA5-Land and observations; (<b>c</b>,<b>g</b>,<b>k</b>,<b>o</b>,<b>s</b>) differences between MERRA-2 and observations; and (<b>d</b>,<b>h</b>,<b>l</b>,<b>p</b>,<b>t</b>) differences between MDA and observations. The spatial RMSEs are shown in the bottom-left corner of each panel.</p>
Full article ">Figure 4
<p>Annual cycles of monthly mean TPSM (units: m3×m-3) on the regional average of TP (<b>a</b>) and separate 17 stations (<b>b</b>–<b>r</b>). Names and location information of each station are given on the top right corner of corresponding panels.</p>
Full article ">Figure 5
<p>Linear trends of annual and seasonal mean TPSM temporal variation (units: m3×m-3/yr): (<b>a</b>,<b>f</b>,<b>k</b>,<b>p</b>,<b>u</b>) observations as the reference; (<b>b</b>,<b>g</b>,<b>l</b>,<b>q</b>,<b>v</b>) GLDAS; (<b>c</b>,<b>h</b>,<b>m</b>,<b>r</b>,<b>w</b>) ERA5-Land; (<b>d</b>,<b>i</b>,<b>n</b>,<b>s</b>,<b>x</b>) MERRA-2; and (<b>e</b>,<b>j</b>,<b>o</b>,<b>t</b>,<b>y</b>) MDA. The pattern correlation coefficients are shown in the bottom-left corner of each panel.</p>
Full article ">Figure 6
<p>Temporal standard deviations of annual and seasonal mean TPSM (units: m3×m-3): (<b>a</b>,<b>f</b>,<b>k</b>,<b>p</b>,<b>u</b>) observations as the reference; (<b>b</b>,<b>g</b>,<b>l</b>,<b>q</b>,<b>v</b>) GLDAS; (<b>c</b>,<b>h</b>,<b>m</b>,<b>r</b>,<b>w</b>) ERA5-Land; (<b>d</b>,<b>i</b>,<b>n</b>,<b>s</b>,<b>x</b>) MERRA-2; and (<b>e</b>,<b>j</b>,<b>o</b>,<b>t</b>,<b>y</b>) MDA. The pattern correlation coefficients are shown in the bottom-left corner of each panel.</p>
Full article ">Figure 7
<p>Time series of soil moisture annual anomalies on the regional average of TP (<b>a</b>) and 17 separate stations (<b>b</b>–<b>r</b>). Names and location information of each station are given in the top right corner of the corresponding panels.</p>
Full article ">Figure 8
<p>TCC spatial distribution of (<b>a</b>) GLDAS; (<b>b</b>) ERA5-Land; (<b>c</b>) MERRA-2; and (<b>d</b>) MDA.</p>
Full article ">Figure 9
<p>EOF1s of annual mean TPSM anomalies and corresponding PC1 time series. (<b>a</b>) Observations; (<b>b</b>) MDA; (<b>c</b>) GLDAS; (<b>d</b>) ERA5-Land; and (<b>e</b>) MERRA-2.</p>
Full article ">Figure 10
<p>Taylor diagram for soil moisture data products based on observations over the Tibetan Plateau.</p>
Full article ">
21 pages, 14797 KiB  
Article
A Parameter Optimized Method for InVEST Model in Sub-Pixel Scale Integrating Machine Learning Algorithm and Vegetation–Impervious Surface–Soil Model
by Linlin Wu and Fenglei Fan
Land 2024, 13(11), 1876; https://doi.org/10.3390/land13111876 - 10 Nov 2024
Viewed by 343
Abstract
The InVEST model, with its ability to perform spatial visualization and quantification, is an important tool for mapping ecosystem services. However, the spatial accuracy and simulating performance of the model are deeply influenced by the land use parameter, which often relies on the [...] Read more.
The InVEST model, with its ability to perform spatial visualization and quantification, is an important tool for mapping ecosystem services. However, the spatial accuracy and simulating performance of the model are deeply influenced by the land use parameter, which often relies on the accuracy of land use/cover data. To address this issue, we propose a novel method for optimizing the land use parameter of the InVEST model based on the vegetation–impervious surface–soil (V–I–S) model and a machine learning algorithm. The optimized model is called Sub-InVEST, and it improves the performance of assessing ecosystem services on a sub-pixel scale. The conceptual steps are (i) extracting the V–I–S fraction of remote sensing images based on the spectral unmixing method; (ii) determining the mapping relationship of the V–I–S fraction between land use/cover type using a machine learning algorithm and field observation data; (iii) inputting the V–I–S fraction into the original model instead of the land use/cover parameter of the InVEST model. To evaluate the performance and spatial accuracy of the Sub-InVEST model, we employed the habitat quality module of InVEST and multi-source remote sensing data, which were applied to acquire Sub-InVEST and estimate the habitat quality of central Guangzhou city from 2000 to 2020 with the help of the LSMA and ISODATA methods. The experimental results showed that the Sub-InVEST model is robust in assessing ecosystem services in sets of complex ground scenes. The spatial distribution of the habitat quality of both models revealed a consistent increasing trend from the southwest to the northeast. Meanwhile, linear regression analyses observed a robust correlation and consistent linear trends, with R2 values of 0.41, 0.35, 0.42, 0.39, and 0.47 for the years 2000, 2005, 2010, 2015, and 2020, respectively. Compared with the original model, Sub-InVEST had a more favorable performance in estimating habitat quality in central Guangzhou. The spatial depictions and numerical distribution of the results of the Sub-InVSET model manifest greater detail and better concordance with remote sensing imagery and show a more seamless density curve and a substantially enhanced probability distribution across interval ranges. Full article
(This article belongs to the Section Land Environmental and Policy Impact Assessment)
Show Figures

Figure 1

Figure 1
<p>The V–I–S fraction combination for a mixed pixel.</p>
Full article ">Figure 2
<p>A flowchart of optimizing the land use parameter of the InVEST model based on the V–I–S model.</p>
Full article ">Figure 3
<p>The location of the study area in (<b>a</b>) China, (<b>b</b>) Guangdong Province, and Guangzhou City, and (<b>c</b>) remote sensing imagery of the study area in 2020.</p>
Full article ">Figure 4
<p>The spatial distribution of V–I–S fractions in central Guangzhou from 2000 to 2020.</p>
Full article ">Figure 5
<p>(<b>a</b>) The spatial distribution of habitat quality based on Sub-InVEST and InVEST. (<b>b</b>) The numerical distribution of habitat quality based on Sub-InVEST and InVEST.</p>
Full article ">Figure 6
<p>The location of sample points and sample regions for comparative assessment.</p>
Full article ">Figure 7
<p>The linear fitting of the InVEST and Sub-InVEST habitat quality results in (<b>a</b>) 2000, (<b>b</b>) 2005, (<b>c</b>) 2010, (<b>d</b>) 2015, (<b>e</b>) 2020.</p>
Full article ">Figure 8
<p>The habitat quality results based on Sub-InVEST, InVEST, and Landsat imagery for (<b>a</b>) 2000 and (<b>b</b>) 2020.</p>
Full article ">Figure 9
<p>The habitat quality results based on Sub-InVEST, InVEST, and remote sensing imagery in 2020.</p>
Full article ">
16 pages, 7229 KiB  
Article
How Has the Source Apportionment of Heavy Metals in Soil and Water Evolved over the Past 20 Years? A Bibliometric Perspective
by Huading Shi, Zexin He, Chenning Deng, Anfu Liu, Yao Feng, Li Li, Guohua Ji, Minghui Xie and Xu Liu
Water 2024, 16(22), 3171; https://doi.org/10.3390/w16223171 - 6 Nov 2024
Viewed by 436
Abstract
Exploring soil heavy metal sources is of great significance for ensuring the safety of ecological environments and agricultural product safety, as well as for guiding pollution control and management policies. This paper retrieved 452 research papers on soil heavy metal source analysis published [...] Read more.
Exploring soil heavy metal sources is of great significance for ensuring the safety of ecological environments and agricultural product safety, as well as for guiding pollution control and management policies. This paper retrieved 452 research papers on soil heavy metal source analysis published over the 2004–2024 period from the Web of Science database. The collected literature was subjected to multidimensional bibliometric analysis using the CiteSpace 6.3.R1. The results showed significantly increasing trends in the scientific outputs and the number of papers on heavy metal source analysis in soils and water over the study period. In addition, related research topics have expanded from single to multiple heavy metal elements in environmental media and have increasingly recognized the impact of water pollution on soil contamination. Research methods have also evolved from basic statistical analysis to complex spatial analysis techniques, covering agricultural and urban soils. Previous related studies have focused on heavy metal pollution in different areas, and related research on heavy metal source analysis has now extended from ecological environments to associated human health risks. The present study provides directions for future related research and guidance for ensuring effective source control of heavy metal pollution and safe utilization of land and water resources. Full article
Show Figures

Figure 1

Figure 1
<p>Statistical results of publication volume on soil heavy metal source analysis.</p>
Full article ">Figure 2
<p>Keyword co-occurrence map.</p>
Full article ">Figure 3
<p>Keyword clustering map.</p>
Full article ">Figure 4
<p>Keyword timeline map.</p>
Full article ">Figure 5
<p>Keyword emergence map.</p>
Full article ">Figure 6
<p>Map of the publishing institutions.</p>
Full article ">Figure 7
<p>Author collaboration network map.</p>
Full article ">
23 pages, 7056 KiB  
Article
Land Subsidence Predictions Based on a Multi-Component Temporal Convolutional Gated Recurrent Unit Model in Kunming City
by Tao Chen, Di Ning and Yuhang Liu
Appl. Sci. 2024, 14(21), 10021; https://doi.org/10.3390/app142110021 - 2 Nov 2024
Viewed by 529
Abstract
Land subsidence (LS) is a geological hazard driven by both natural conditions and human activities. Traditional LS time-series prediction models often struggle to accurately capture nonlinear data characteristics, leading to suboptimal predictions. To address this issue, this paper introduces a multi-component temporal convolutional [...] Read more.
Land subsidence (LS) is a geological hazard driven by both natural conditions and human activities. Traditional LS time-series prediction models often struggle to accurately capture nonlinear data characteristics, leading to suboptimal predictions. To address this issue, this paper introduces a multi-component temporal convolutional gate recurrent unit (MC-TCGRU) model, which integrates a fully adaptive noise-ensemble empirical-mode decomposition algorithm with a deep neural network to account for the complexity of time-series data. The model was validated using typical InSAR subsidence data from Kunming, analyzing the impact of each component on the prediction performance. A comparative analysis with the TCGRU model and models based on seasonal-trend decomposition using LOESS (STL) and empirical-mode decomposition (EMD) revealed that the MC-TCGRU model significantly enhanced the prediction accuracy by reducing the complexity of the original data. The model achieved R² values of 0.90, 0.93, 0.51, 0.93, and 0.96 across five points, outperforming the compared models. Full article
(This article belongs to the Special Issue Advanced Remote Sensing Technologies and Their Applications)
Show Figures

Figure 1

Figure 1
<p>Location of the study area, with an indication of the administrative boundaries of Kunming City.</p>
Full article ">Figure 2
<p>Structure of the GRU network unit.</p>
Full article ">Figure 3
<p>Structure of the TCN residual block.</p>
Full article ">Figure 4
<p>Structure of the TCGRU model.</p>
Full article ">Figure 5
<p>Structure of the proposed MC-TCGRU model.</p>
Full article ">Figure 6
<p>Deformation velocity map of the study area, with indication of five distinct LS regions and points.</p>
Full article ">Figure 7
<p>CEEMDAN decomposition results for 5 LS points. (<b>a</b>) P1; (<b>b</b>) P2; (<b>c</b>) P3; (<b>d</b>) P4; and (<b>e</b>) P5.</p>
Full article ">Figure 8
<p>Prediction results of each component based on CEEMDAN for point P1; (<b>a</b>) IMF<sub>1</sub> component; (<b>b</b>) IMF<sub>2</sub> component; (<b>c</b>) IMF<sub>3</sub> component; and (<b>d</b>) residual component.</p>
Full article ">Figure 8 Cont.
<p>Prediction results of each component based on CEEMDAN for point P1; (<b>a</b>) IMF<sub>1</sub> component; (<b>b</b>) IMF<sub>2</sub> component; (<b>c</b>) IMF<sub>3</sub> component; and (<b>d</b>) residual component.</p>
Full article ">Figure 9
<p>Point P1 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 10
<p>Point P2 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 11
<p>Point P3 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 12
<p>Point P4 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 13
<p>Point P5 LS time-series prediction results (<b>left</b>) and enlarged view of test set (<b>right</b>).</p>
Full article ">Figure 14
<p>Comparison of five AMEs on five LS points. (<b>a</b>) P1. (<b>b</b>) P2. (<b>c</b>) P3. (<b>d</b>) P4. (<b>e</b>) P5.</p>
Full article ">Figure 14 Cont.
<p>Comparison of five AMEs on five LS points. (<b>a</b>) P1. (<b>b</b>) P2. (<b>c</b>) P3. (<b>d</b>) P4. (<b>e</b>) P5.</p>
Full article ">
26 pages, 9980 KiB  
Article
Detecting Trends in Post-Fire Forest Recovery in Middle Volga from 2000 to 2023
by Eldar Kurbanov, Ludmila Tarasova, Aydin Yakhyayev, Oleg Vorobev, Siyavush Gozalov, Sergei Lezhnin, Jinliang Wang, Jinming Sha, Denis Dergunov and Anna Yastrebova
Forests 2024, 15(11), 1919; https://doi.org/10.3390/f15111919 - 31 Oct 2024
Viewed by 528
Abstract
Increased wildfire activity is the most significant natural disturbance affecting forest ecosystems as it has a strong impact on their natural recovery. This study aimed to investigate how burn severity (BS) levels and climate factors, including land surface temperature (LST) and precipitation variability [...] Read more.
Increased wildfire activity is the most significant natural disturbance affecting forest ecosystems as it has a strong impact on their natural recovery. This study aimed to investigate how burn severity (BS) levels and climate factors, including land surface temperature (LST) and precipitation variability (Pr), affect forest recovery in the Middle Volga region of the Russian Federation. It provides a comprehensive analysis of post-fire forest recovery using Landsat time-series data from 2000 to 2023. The analysis utilized the LandTrendr algorithm in the Google Earth Engine (GEE) cloud computing platform to examine Normalized Burn Ratio (NBR) spectral metrics and to quantify the forest recovery at low, moderate, and high burn severity (BS) levels. To evaluate the spatio-temporal trends of the recovery, the Mann–Kendall statistical test and Theil–Sen’s slope estimator were utilized. The results suggest that post-fire spectral recovery is significantly influenced by the degree of the BS in affected areas. The higher the class of BS, the faster and more extensive the reforestation of the area occurs. About 91% (40,446 ha) of the first 5-year forest recovery after the wildfire belonged to the BS classes of moderate and high severity. A regression model indicated that land surface temperature (LST) plays a more critical role in post-fire recovery compared to precipitation variability (Pr), accounting for approximately 65% of the variance in recovery outcomes. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical distribution of the study area in the western Russian Federation on the Landsat scenes.</p>
Full article ">Figure 2
<p>Flow chart of the methodology.</p>
Full article ">Figure 3
<p>Annual number of the Landsat (TM/ETM+/OLI/OLI2) scene observations over the period May–September 2000–2023.</p>
Full article ">Figure 4
<p>Examples of the forest recovery in three different 2010 wildfire burn severity sites in the Middle Volga region.</p>
Full article ">Figure 5
<p>Schematic of the spectral recovery metrics for the 2010 forest fire burnt area using an example of a fitted NBR trajectory in the LandTrendr program.</p>
Full article ">Figure 6
<p>Box and whisker plot of the time-series recovery metrics by the three burn severity classes: (<b>a</b>) the relative recovery indicator at 5 years after the wildfire (RRI<sub>5</sub>); (<b>b</b>) the average annual recovery (YrYr), and (<b>c</b>) the annual recovery at 5 years after the wildfire (YrYr<sub>5</sub>).</p>
Full article ">Figure 7
<p>Time series of the NBR recovery by burn severity after the wildfires: (<b>a</b>) 2002; (<b>b</b>) 2006; (<b>c</b>) 2010; and (<b>d</b>) 2018. The green horizontal dashed line represents the mean NBR reference for the two years prior to the wildfire, while the red vertical dashed line marks the year of the wildfire.</p>
Full article ">Figure 8
<p>On-site examples of the three recovery stages on the BA of the high BS for the spectral metrics: (<b>a</b>) RRI<sub>5</sub>; (<b>b</b>) Y2R80; and (<b>c</b>) Rec<sub>endTS</sub>.</p>
Full article ">Figure 9
<p>Spatial trends of the post-fire natural forest recovery in the Middle Volga region (2003–2023) based on the Mann–Kendall’s Tau test (<span class="html-italic">p</span> &lt; 0.05) for sites with different burn severities: (<b>a</b>) low; (<b>b</b>) moderate; and (<b>c</b>) high.</p>
Full article ">Figure 10
<p>The spatial trends of post-fire natural forest recovery in the Middle Volga region (2003–2023) based on the Theil–Sen’s slope (<span class="html-italic">p</span> &lt; 0.05) for sites with different burn severities: (<b>a</b>) low; (<b>b</b>) moderate; and (<b>c</b>) high.</p>
Full article ">Figure 11
<p>The Theil–Sen’s slope (magnitude) trend for the forest recovery in the Middle Volga region from 2003–2023 for sites with different burn severities.</p>
Full article ">Figure 12
<p>The fitting results of a multiple linear regression model to describe the relationship between the natural forest recovery (NBR) during the 5 years after a wildfire, the land surface temperature (LST), and the precipitation (Pr) for sites with different burn severities: (<b>a</b>) low; (<b>b</b>) moderate; and (<b>c</b>) high.</p>
Full article ">Figure 13
<p>The fitting results of a multiple linear regression model to describe the relationship between the natural forest recovery (NBR) during the 5 + 10 years after a wildfire, the land surface temperature (LST), and precipitation (PR) for sites with different burn severities: (<b>a</b>) low; (<b>b</b>) moderate; and (<b>c</b>) high.</p>
Full article ">
17 pages, 5901 KiB  
Article
A Cropland Disturbance Monitoring Method Based on Probabilistic Trajectories
by Jiawei Jiang, Juanle Wang, Keming Yang, Denis Fetisov, Kai Li, Meng Liu and Weihao Zou
Remote Sens. 2024, 16(21), 4048; https://doi.org/10.3390/rs16214048 - 30 Oct 2024
Viewed by 367
Abstract
Acquiring the spatiotemporal patterns of cropland disturbance is of great significance for regional sustainable agricultural development and environmental protection. However, effective monitoring of cropland disturbances remains a challenge owing to the complexity of the terrain landscape and the reliability of the training samples. [...] Read more.
Acquiring the spatiotemporal patterns of cropland disturbance is of great significance for regional sustainable agricultural development and environmental protection. However, effective monitoring of cropland disturbances remains a challenge owing to the complexity of the terrain landscape and the reliability of the training samples. This study integrated automatic training sample generation, random forest classification, and the LandTrendr time-series segmentation algorithm to propose an efficient and reliable medium-resolution cropland disturbance monitoring scheme. Taking the Amur state of Russia in the Amur river basin, a transboundary region between Russia and China in east Asia with rich agriculture resources as research area, this approach was conducted on the Google Earth Engine cloud-computing platform using extensive remote-sensing image data. A high-confidence sample dataset was then created and a random forest classification algorithm was applied to generate the cropland classification probabilities. LandTrendr time-series segmentation was performed on the interannual cropland classification probabilities. Finally, the identification, spatial mapping, and analysis of cropland disturbances in Amur state were completed. Further cross-validation comparisons of the accuracy assessment and spatiotemporal distribution details demonstrated the high accuracy of the dataset, and the results indicated the applicability of the method. The study revealed that 2815.52 km2 of cropland was disturbed between 1990 and 2021, primarily focusing on the southern edge of the Amur state. The most significant disturbance occurred in 1991, affecting 1431.48 km2 and accounting for 50.84% of the total disturbed area. On average, 87.98 km2 of croplands are disturbed annually. Additionally, 2495.4 km2 of cropland was identified as having been disturbed at least once during the past 32 years, representing 83% of the total disturbed area. This study introduced a novel approach for identifying cropland disturbance information from long time-series probabilistic images. This methodology can also be extended to monitor the spatial and temporal dynamics of other land disturbances caused by natural and human activities. Full article
Show Figures

Figure 1

Figure 1
<p>Situation of Amur state.</p>
Full article ">Figure 2
<p>Overall technical process.</p>
Full article ">Figure 3
<p>The number of available Landsat images during 1991–2020 in Amur state.</p>
Full article ">Figure 4
<p>LandTrendr time-series segmentation of a cropland disturbance pixel. A: disturbance start cropland probability. B: disturbance stop cropland probability. C: C = A-B, disturbance magnitude of the cropland. E: disturbance start year. F: disturbance end year. D: D = F-E, disturbance time duration.</p>
Full article ">Figure 5
<p>Spatial and temporal distribution of cropland disturbance.</p>
Full article ">Figure 6
<p>Duration of cropland disturbance.</p>
Full article ">Figure 7
<p>Area of the duration of cropland disturbance (unit = year). The legend indicates the duration of the disturbance. The left graph shows the area of cropland disturbance duration, and the right graph shows the percentage of area of cropland disturbance duration.</p>
Full article ">Figure 8
<p>Comparison between cropland disturbance monitoring results with Google satellite images. (<b>A</b>,<b>B</b>) Google satellite images, (<b>C</b>) cropland disturbance monitoring results, and (<b>D</b>) disturbed trajectories of cropland based on LandTrendr algorithm fitting. The color gradient of the disturbance monitoring results represents different disturbance years.</p>
Full article ">
21 pages, 15249 KiB  
Article
Variations of Lake Ice Phenology Derived from MODIS LST Products and the Influencing Factors in Northeast China
by Xiaoguang Shi, Jian Cheng, Qian Yang, Hongxing Li, Xiaohua Hao and Chunxu Wang
Remote Sens. 2024, 16(21), 4025; https://doi.org/10.3390/rs16214025 - 30 Oct 2024
Viewed by 407
Abstract
Lake ice phenology serves as a sensitive indicator of climate change in the lake-rich Northeast China. In this study, the freeze-up date (FUD), break-up date (BUD), and ice cover duration (ICD) of 31 lakes were extracted from a time series of the land [...] Read more.
Lake ice phenology serves as a sensitive indicator of climate change in the lake-rich Northeast China. In this study, the freeze-up date (FUD), break-up date (BUD), and ice cover duration (ICD) of 31 lakes were extracted from a time series of the land water surface temperature (LWST) derived from the combined MOD11A1 and MYD11A1 products for the hydrological years 2001 to 2021. Our analysis showed a high correlation between the ice phenology measures derived by our study and those provided by hydrological records (R2 of 0.89) and public datasets (R2 > 0.7). There was a notable coherence in lake ice phenology in Northeast China, with a trend in later freeze-up (0.21 days/year) and earlier break-up (0.19 days/year) dates, resulting in shorter ice cover duration (0.50 days/year). The lake ice phenology of freshwater lakes exhibited a faster rate of change compared to saltwater lakes during the period from HY2001 to HY2020. We used redundancy analysis and correlation analysis to study the relationships between the LWST and lake ice phenology with various influencing factors, including lake properties, local climate factors, and atmospheric circulation. Solar radiation, latitude, and air temperature were found to be the primary factors. The FUD was more closely related to lake characteristics, while the BUD was linked to local climate factors. The large-scale oscillations were found to influence the changes in lake ice phenology via the coupled influence of air temperature and precipitation. The Antarctic Oscillation and North Atlantic Oscillation correlate more with LWST in winter, and the Arctic Oscillation correlates more with the ICD. Full article
Show Figures

Figure 1

Figure 1
<p>Geographic location of Northeast China and the 31 selected lakes: (<b>a</b>) geographic map of Northeast China; (<b>b</b>) elevation map of Northeast China.</p>
Full article ">Figure 2
<p>Example of lake ice phenology extraction for Chagan Lake.</p>
Full article ">Figure 3
<p>Determination coefficient (R<sup>2</sup>) for two curve fitting methods applied to LWST: (<b>a</b>) Orthogonal Distance Regression (ODR); (<b>b</b>) Levenberg Marquardt (L-M).</p>
Full article ">Figure 4
<p>Comparison of daily lake surface temperatures derived from in situ measurements and remote sensing satellites spanning the period HY2001–HY2020: (<b>a</b>) ground surface temperature vs. lake water surface temperature; (<b>b</b>) air temperature vs lake water surface temperature.</p>
Full article ">Figure 5
<p>Comparison of lake ice phenology with three different sources: (<b>a</b>) FUD between our results and hydrological records; (<b>b</b>) BUD from our results and hydrological records; (<b>c</b>) ICD between our results and hydrological records; (<b>d</b>) FUD between our results and Qiu_2019; (<b>e</b>) BUD between our results and Qiu_2019; (<b>f</b>) ICD between our results and Qiu_2019; (<b>g</b>) FUD between our results and Wang_2021; (<b>h</b>) BUD between our results and Wang_2021; (<b>i</b>) ICD between our results and Wang_2021.</p>
Full article ">Figure 6
<p>Boxplot of lake ice phenology of 31 lakes in Northeast China from 2000 to 2021: (<b>a</b>) freeze-up date; (<b>b</b>) break-up date; (<b>c</b>) ice cover duration.</p>
Full article ">Figure 7
<p>Spatial distribution of lake ice phenology of 31 lakes in Northeast China from 2000 to 2021: (<b>a</b>) freeze-up date (FUD); (<b>b</b>) break-up date (BUD); (<b>c</b>) ice cover duration (ICD); (<b>d</b>) yearly changing rates of the FUD; (<b>e</b>) yearly changing rates of the BUD; (<b>f</b>) yearly changing rates of the ICD; (<b>g</b>) lake ice process of freshwater and saltwater.</p>
Full article ">Figure 8
<p>Inter-annual variations of LWST: (<b>a</b>) daytime LWST; (<b>b</b>) nighttime LWST; (<b>c</b>) mean LWST; (<b>d</b>) DTD.</p>
Full article ">Figure 9
<p>Inter-annual variation of lake ice phenology in freshwater and saltwater lakes in Northeast China from HY2001 to HY2020: (<b>a</b>) freeze-up date (FUD); (<b>b</b>) break-up date (BUD); (<b>c</b>) ice cover duration (ICD).</p>
Full article ">Figure 10
<p>Redundancy analysis between the LWST, lake ice phenology, and impact factors in Northeast China (n = 620): (<b>a</b>) FLWST within the full year; (<b>b</b>) WLWST within the winter; (<b>c</b>) lake ice phenology. The variables considered for the FLWST and WLWST included the daytime LWST, nighttime LWST, and daily mean LWST. As for lake ice phenology, the variables were the FUD, BUD, and ICD. Lake properties included longitude (LON) and latitude (LAT), lake area (AREA), altitude (ALT), and average water depth (AWD). Local climate factors included air temperature (AT), ground surface temperature (GST), precipitation (PR), wind speed (WS), surface net solar radiation (SR), and downward surface solar radiation (DSR). Atmospheric circulation included the North Atlantic Oscillation (NAO), Arctic Oscillation (AO), Antarctic Oscillation (AAO) and Pacific–North American Pattern (PNA).</p>
Full article ">Figure 11
<p>Correlation coefficients between the LWST, lake ice phenology, and impact factors. FLWST<sub>day</sub>, FLWST<sub>night</sub>, and FLWST<sub>mean</sub> represent the daytime, nighttime, and daily mean of LWST for a full year, and WLWST<sub>day</sub>, WLWST<sub>night</sub>, and WLWST<sub>mean</sub> represent the daytime, nighttime and daily mean of LWST for the winter from December to February. The other acronyms are the same as those in <a href="#remotesensing-16-04025-f010" class="html-fig">Figure 10</a>. ** means significant at 99 % level (<span class="html-italic">p</span> &lt; 0.01), and * means significant at 95 % level (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 12
<p>Proportions of null values in different MODIS LST products. MOD<sub>day</sub> and MOD<sub>night</sub> stand for the daytime and nighttime LWST provided by MOD11A1, and MYD<sub>day</sub> and MYD<sub>night</sub> the daytime and nighttime LWST provided by MYD11A1. Day represents the merged production by MOD<sub>day</sub> and MYD<sub>day</sub>, and Night represents the merged production by MOD<sub>day</sub> and MYD<sub>day</sub>.</p>
Full article ">
21 pages, 13748 KiB  
Article
Dynamic Changes in and Driving Factors of Soil Organic Carbon in China from 2001 to 2020
by Fuyan Zou, Min Yan, Liankai Zhang, Jinjiang Yang, Guiren Chen, Keqiang Shan, Chen Zhang, Xiongwei Xu, Zhenhui Wang and Can Xu
Land 2024, 13(11), 1764; https://doi.org/10.3390/land13111764 - 27 Oct 2024
Viewed by 665
Abstract
It remains unclear what changes have occurred in the distribution pattern of and trend in soil organic carbon (SOC) in China against the background of climate and land use change. Clarifying the dynamic changes in SOC and their driving factors in different regions [...] Read more.
It remains unclear what changes have occurred in the distribution pattern of and trend in soil organic carbon (SOC) in China against the background of climate and land use change. Clarifying the dynamic changes in SOC and their driving factors in different regions of China is therefore crucial for assessing the global carbon cycle. In this study, we collected and supplemented a large amount of soil organic carbon density (SOCD) data in China from 2001 to 2020 and extracted data on environmental covariates (ECs) for the corresponding years. A random forest model was used to estimate the SOCD at a depth of 0–20 cm and 0–100 cm in China for the years 2001, 2005, 2010, 2015, and 2020, and we explored the trend of SOCD changes and their key driving factors. The results showed the following: (1) Compared with previous studies, the predictive ability of the 0–100 cm depth model was greatly improved; the coefficient of determination (R2) was 0.61 and Lin’s concordance correlation coefficient (LCCC) was =0.76. (2) From 2001 to 2020, China’s soil organic carbon stocks (SOCS) were 38.11, 39.11, 39.88, 40.16, and 41.12 Pg C for the 0–20 cm depth and 110.49, 112.67, 112.80, 113.06, and 114.96 Pg C for the 0–100 cm depth, respectively. (3) The effects of temperature and precipitation on SOCD in China showed obvious regional variability, and land use changes had mainly positive effects on SOCD in all regions of China, which was related to the large-scale implementation of ecological protection and restoration and the policy of returning farmland to forests and grasslands in China. This study provides strong scientific support for addressing climate change and rationalizing the use of land resources. Full article
(This article belongs to the Section Land Systems and Global Change)
Show Figures

Figure 1

Figure 1
<p>Distribution of sample values of soil carbon density in soil layers at 0–20 cm (<b>a</b>) and 0–100 cm (<b>b</b>) depths. Sampled data refer to the soil data obtained from our project team’s field survey. Collected data refer to the collection of data from the literature and databases.</p>
Full article ">Figure 2
<p>Scatter plots of estimated and predicted SOCD values: depths of 0–20 cm (<b>a</b>) and 0–100 cm (<b>b</b>). R<sup>2</sup>, coefficient of determination; RMSE, root mean squared error; LCCC, Lin’s concordance correlation coefficient.</p>
Full article ">Figure 3
<p>Order of importance of ECs used to predict SOCD in random forest: depths of 0–20 cm (<b>a</b>) and 0–100 cm (<b>b</b>).</p>
Full article ">Figure 4
<p>Spatial distribution of and temporal variation in 0–20 cm SOCD. (<b>a</b>–<b>e</b>) represent SOCD distribution maps; (<b>f</b>) shows the change in SOCD. When slope &gt; 0, the SOCD of the time series shows an increasing trend; when slope &lt; 0, the SOCD of the time series shows a decreasing trend.</p>
Full article ">Figure 5
<p>Spatial distribution of and temporal variation in 0–100 cm SOCD. (<b>a</b>–<b>e</b>) represent SOCD spatial distribution; (<b>f</b>) shows the change in SOCD. When slope &gt; 0, the SOCD of the time series shows an increasing trend; when slope &lt; 0, the SOCD of the time series shows a decreasing trend.</p>
Full article ">Figure 6
<p>Zonal statistics for SOCS, mean SOCD, and mean Theil–Sen median slope: depth of 0–20 cm (<b>a</b>) and 0–100 cm (<b>b</b>).</p>
Full article ">Figure 6 Cont.
<p>Zonal statistics for SOCS, mean SOCD, and mean Theil–Sen median slope: depth of 0–20 cm (<b>a</b>) and 0–100 cm (<b>b</b>).</p>
Full article ">Figure 7
<p>Spatial pattern of partial correlation and correlation coefficients of SOCD at a depth between 0 and 20 cm and influencing factors from 2001 to 2020. These factors include temperature (<b>a</b>), precipitation (<b>b</b>), land use disturbance intensity (<b>c</b>), and (<b>d</b>) the percentage of partial correlation and correlation between SOCD at a 0 to 20 cm depth and the three influencing factors. Upward and downward bars indicate percentages of positive and negative correlation, respectively. Colored areas indicate correlation coefficients greater than 0.5 or less than −0.5.</p>
Full article ">Figure 8
<p>Spatial pattern of partial correlation and correlation coefficients of SOCD at a depth between 0 and 100 cm and influencing factors from 2001 to 2020. These factors include temperature (<b>a</b>), precipitation (<b>b</b>), land use disturbance intensity (<b>c</b>), and (<b>d</b>) the percentage of partial correlation and correlation between SOCD at a 0 to 100 cm depth and the three influencing factors. Upward and downward bars indicate percentages of positive and negative correlation, respectively. Colored areas indicate correlation coefficients greater than 0.5 or less than −0.5.</p>
Full article ">Figure 9
<p>Land use change from 2001 to 2020 (the left side is the area transferred out of different land types, and the right side is the area transferred in from different land types; area unit is km<sup>2</sup>).</p>
Full article ">
23 pages, 12411 KiB  
Article
Does ERA5-Land Effectively Capture Extreme Precipitation in the Yellow River Basin?
by Chunrui Guo, Ning Ning, Hao Guo, Yunfei Tian, Anming Bao and Philippe De Maeyer
Atmosphere 2024, 15(10), 1254; https://doi.org/10.3390/atmos15101254 - 21 Oct 2024
Viewed by 590
Abstract
ERA5-Land is a valuable reanalysis data resource that provides near-real-time, high-resolution, multivariable data for various applications. Using daily precipitation data from 301 meteorological stations in the Yellow River Basin from 2001 to 2013 as benchmark data, this study aims to evaluate ERA5-Land’s capability [...] Read more.
ERA5-Land is a valuable reanalysis data resource that provides near-real-time, high-resolution, multivariable data for various applications. Using daily precipitation data from 301 meteorological stations in the Yellow River Basin from 2001 to 2013 as benchmark data, this study aims to evaluate ERA5-Land’s capability of monitoring extreme precipitation. The evaluation study is conducted from three perspectives: precipitation amount, extreme precipitation indices, and characteristics of extreme precipitation events. The results show that ERA5-Land can effectively capture the spatial distribution patterns and temporal trends in precipitation and extreme precipitation; however, it also exhibits significant overestimation and underestimation errors. ERA5-Land significantly overestimates total precipitation and indices for heavy precipitation and extreme precipitation (R95pTOT and R99pTOT), with errors reaching up to 89%, but underestimates the Simple Daily Intensity Index (SDII). ERA5-Land tends to overestimate the duration of extreme precipitation events but slightly underestimates the total and average precipitation of these events. These findings provide a scientific reference for optimizing the ERA5-Land algorithm and for users in selecting data. Full article
(This article belongs to the Special Issue Advances in Rainfall-Induced Hazard Research)
Show Figures

Figure 1

Figure 1
<p>The map of the Yellow River Basin, including (<b>a</b>) the geographical location of the Yellow River Basin in China, (<b>b</b>) monthly precipitation in the Yellow River Basin and its sub-basins, and (<b>c</b>) the terrain conditions and distribution of the meteorological stations used in this study.</p>
Full article ">Figure 2
<p>The concept diagram of an extreme precipitation event and its characteristics.</p>
Full article ">Figure 3
<p>Spatial distribution comparison of annual precipitation and precipitation days from 2001 to 2013. The upper right corner of each panel features a dual-ring presentation. The outer ring illustrates (<b>a</b>) annual precipitation amount from the ERA5-Land and (<b>b</b>) the count of precipitation days from the ERA5-Land. The inner ring displays (<b>a</b>) annual precipitation amount from station observations and (<b>b</b>) the count of precipitation days from station observations.</p>
Full article ">Figure 4
<p>Spatial distribution of statistical indices between daily precipitation from ERA5-Land and station observations for (<b>a</b>) <span class="html-italic">RB</span>, (<b>b</b>) <span class="html-italic">CC</span>, and (<b>c</b>) <span class="html-italic">RMSE</span> from 2001 to 2013.</p>
Full article ">Figure 5
<p>The spatial pattern of the <span class="html-italic">RB</span> of different extreme precipitation indices including (<b>a</b>) SDII, (<b>b</b>) PRCPTOT, (<b>c</b>) RX1day, (<b>d</b>) RX5day, (<b>e</b>) R95pTOT, (<b>f</b>) R99pTOT, (<b>g</b>) R10mm, (<b>h</b>) R20mm, (<b>i</b>) CDD, and (<b>j</b>) CWD between ERA5-Land and station observations.</p>
Full article ">Figure 6
<p>The spatial pattern of CCs for different extreme precipitation indices including (<b>a</b>) SDII, (<b>b</b>) PRCPTOT, (<b>c</b>) RX1day, (<b>d</b>) RX5day, (<b>e</b>) R95pTOT, (<b>f</b>) R99pTOT, (<b>g</b>) R10mm, (<b>h</b>) R20mm, (<b>i</b>) CDD, and (<b>j</b>) CWD between ERA5-Land and station observations.</p>
Full article ">Figure 7
<p>The spatial pattern of <span class="html-italic">RMSE</span>s for different extreme precipitation indices between ERA5-Land and station observations. The <span class="html-italic">RMSE</span> units are mm/day for (<b>a</b>) SDII; mm for (<b>b</b>–<b>f</b>) PRCPTOT, RX1day, RX5day, R95pTOT, and R99pTOT; and days for (<b>g</b>–<b>j</b>) R10mm, R20mm, CDD, and CWD.</p>
Full article ">Figure 8
<p>The temporal evolution of various extreme precipitation indices based on both ERA5-Land and station observations including (<b>a</b>) SDII, (<b>b</b>) PRCPTOT, (<b>c</b>) RX1day, (<b>d</b>) RX5day, (<b>e</b>) R95PTOT, (<b>f</b>) R99PTOT, (<b>g</b>) R10mm, (<b>h</b>) R20mm, (<b>i</b>) CDD, and (<b>j</b>) CWD. STN indicates station observations. The black lines represent the indices for the entire Yellow River Basin; the orange lines represent those for the upper basin; the yellow lines represent those for the middle basin; the green lines represent those for the lower basin.</p>
Full article ">Figure 9
<p>Spatial pattern of event characteristics for (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) station observations and (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) ERA5-Land, as well as (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) the difference between ERA5-Land and station observations. STN indicates station observations.</p>
Full article ">Figure 10
<p>Annual time distribution of extreme precipitation events for (<b>a</b>) <span class="html-italic">EF</span>, (<b>b</b>) <span class="html-italic">ED</span>, (<b>c</b>) <span class="html-italic">EP</span>, (<b>d</b>) <span class="html-italic">ET</span>, and (<b>e</b>) <span class="html-italic">EM</span>. STN indicates station observations.</p>
Full article ">Figure 11
<p>Monthly frequency of extreme precipitation events spanning from (<b>a</b>) station observations and (<b>b</b>) ERA5-Land from 2001 to 2013. STN indicates station observations.</p>
Full article ">Figure 12
<p>Annual time distribution of extreme precipitation events for (<b>a</b>) <span class="html-italic">EF</span>, (<b>b</b>) <span class="html-italic">ED</span>, (<b>c</b>) <span class="html-italic">EP</span>, (<b>d</b>) <span class="html-italic">ET</span>, and (<b>e</b>) <span class="html-italic">EM</span> in the upper basin, middle basin and lower basin. STN indicates station observations.</p>
Full article ">Figure A1
<p>The theoretical–methodological flowchart.</p>
Full article ">Figure A2
<p>Spatial distribution of annual mean extreme precipitation indices including (<b>a</b>,<b>b</b>) SDII, (<b>c</b>,<b>d</b>) PRCPTOT, (<b>e</b>,<b>f</b>) RX1day, (<b>g</b>,<b>h</b>) RX5day, (<b>i</b>,<b>j</b>) P95pTOT in the Yellow River Basin from 2001 to 2013. STN indicates station observations.</p>
Full article ">Figure A3
<p>Spatial distribution of annual mean extreme precipitation indices including (<b>a</b>,<b>b</b>) R99pTOT, (<b>c</b>,<b>d</b>) R10mm, (<b>e</b>,<b>f</b>) R20mm, (<b>g</b>,<b>h</b>) CDD, (<b>i</b>,<b>j</b>) CWD in the Yellow River Basin from 2001 to 2013. STN indicates station observations.</p>
Full article ">
20 pages, 8446 KiB  
Article
Attribution of Vegetation Dynamics in the Yellow River Water Conservation Area Based on the Deep ConvLSTM Model
by Zhi Liang, Ruochen Sun and Qingyun Duan
Remote Sens. 2024, 16(20), 3875; https://doi.org/10.3390/rs16203875 - 18 Oct 2024
Viewed by 612
Abstract
Climate change and human activities have significantly impacted the long-term growth of vegetation, thereby altering the ecosystem’s response mechanisms. The Yellow River Water Conservation Area (YRWCA) is a critical ecological functional zone in China. Since 1982, the vegetation in the YRWCA has changed [...] Read more.
Climate change and human activities have significantly impacted the long-term growth of vegetation, thereby altering the ecosystem’s response mechanisms. The Yellow River Water Conservation Area (YRWCA) is a critical ecological functional zone in China. Since 1982, the vegetation in the YRWCA has changed significantly, and the primary drivers of vegetation which changed before and after 2000 were identified as climate change and human activities, respectively. However, the extent to which different drivers contribute to the vegetation dynamics of the YRWCA remains uncertain. In this study, we introduced a modified deep Convolutional Long Short-Term Memory (ConvLSTM) model to quantify the contributions of climate change and human activities to vegetation change while considering the spatiotemporal heterogeneity. We identified areas with minimal human activity before 2000 using the residual trend method, and used the regional data from these areas to train the model. Subsequently, we applied the trained deep ConvLSTM model to perform an attribution analysis after 2000. The results show that the deep ConvLSTM effectively captures the impacts of climate change on vegetation growth and outperforms the widely used Random Forest model (RF). Despite the fact that the input data of RF were optimized, ConvLSTM still distinctly outperformed RF, achieving R2, MAE, and RMSE values of 0.99, 0.013, and 0.018, respectively, compared to RF’s corresponding values of 0.94, 0.038, and 0.045. Since 2000, the regional normalized difference vegetation index (NDVI) has shown a broad increasing trend, particularly in dryland, primarily induced by human activities from 2006 to 2015. Furthermore, an analysis of changes in regional land use, particularly in drylands, revealed that the highest magnitude of conversion of farmland back to forest or grass was recorded from 2000 to 2005. However, the most significant contributions from human activities occurred from 2006 to 2015, indicating a time lag in vegetation recovery from these ecological programs. The attribution results provide valuable insights for the implementation of ecological programs, and the introduced deep ConvLSTM proves the suitability of deep learning models that capture spatiotemporal features in vegetation growth simulations, allowing for broader applications. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Spatial location and topographic characteristics of the Yellow River Water Conservation Area (YRWCA).</p>
Full article ">Figure 2
<p>The construction process of the deep Convolutional Long Short-Term Memory (ConvLSTM) model. For a given spatial point, climate factors are cropped into a 5 × 5 dataset, while normalized difference vegetation index (NDVI) is cropped into a 1 × 1 dataset. Subsequently, four consecutive months of climate factors and NDVI data for the last month are, respectively, used as inputs and outputs, with the inclusion of certain terrestrial variables factors as inputs to ensure the model’s robustness. The data are fed into a two-layer ConvLSTM, where the output from the hidden layers undergoes pooling across channels before being passed to a fully connected layer to obtain NDVI data. The predicted NDVI data can then be concatenated to form the predicted NDVI distribution.</p>
Full article ">Figure 3
<p>Vegetation characteristics of the YRWCA.</p>
Full article ">Figure 4
<p>Multi-year trends of various climate factors in growing season. (<b>a</b>–<b>i</b>) denote the trends of the variables: t2m, tmax, tmin, tp, ssrd, d2m, u10, v10, sp, respectively.</p>
Full article ">Figure 5
<p>Changes in human activity intensity during different time periods.</p>
Full article ">Figure 6
<p>Mean time-lag spatial distribution of vegetation by different climate factors. (<b>a</b>–<b>i</b>) denote the mean time-lag of the variables: t2m, tmax, tmin, tp, ssrd, d2m, u10, v10, and sp, respectively.</p>
Full article ">Figure 7
<p>Results of the residual trend method before and after 2000: (<b>a1</b>–<b>a3</b>) Shows the NDVI trends over different years. (<b>b1</b>–<b>b3</b>) Depicts the climate regression trends. (<b>c1</b>–<b>c3</b>) Illustrates residual trends due to human activities.</p>
Full article ">Figure 8
<p>The scatter plot of the NDVI using the trained deep ConvLSTM model during the validation period.</p>
Full article ">Figure 9
<p>Effects of climate change and human activity on NDVI since 2000: (<b>a1</b>–<b>a3</b>) 2000–2015, (<b>b1</b>–<b>b3</b>) 2000–2005, (<b>c1</b>–<b>c3</b>) 2006–2010, and (<b>d1</b>–<b>d3</b>) 2011–2015. NDVI change intensity refers to the difference between observed NDVI values and baseline NDVI predictions. Climate-driven NDVI and human-driven NDVI, respectively, represent the ratios of NDVI change caused by climate change and human activity compared to the baseline period, expressed in percentage.</p>
Full article ">Figure 10
<p>Annual magnitude of climate change and human activity impacts on different vegetation types annually from 2000 to 2015, categorized by the entire region, dryland, forest, shrubland, meadow, grassland, and wetland.</p>
Full article ">
16 pages, 6843 KiB  
Article
Seasonal–Diurnal Distribution of Lightning over Bulgaria and the Black Sea and Its Relationship with Sea Surface Temperature
by Savka Petrova, Rumjana Mitzeva, Vassiliki Kotroni and Elisaveta Peneva
Atmosphere 2024, 15(10), 1233; https://doi.org/10.3390/atmos15101233 - 15 Oct 2024
Viewed by 375
Abstract
A seasonal–diurnal analysis of land-sea contrast in lightning activity over Bulgaria and the Black Sea over 10 years is presented here. The maximum number of flashes over both surface types is registered during the summer (with a peak over Bulgaria in June and [...] Read more.
A seasonal–diurnal analysis of land-sea contrast in lightning activity over Bulgaria and the Black Sea over 10 years is presented here. The maximum number of flashes over both surface types is registered during the summer (with a peak over Bulgaria in June and over the Black Sea in July) and a minimum number in winter (December/February, respectively). During spring, the maximum flash density is observed over Bulgaria (in May), while in autumn, it is over the Black Sea (in September). The results show that only in autumn lightning activity dominates over the Black Sea compared to over land (Bulgaria), while in winter, spring, and summer is vice versa. For this reason, an additional investigation was conducted to determine whether there is a relationship between lightning activity and the sea surface temperature (SST) of the Black Sea in autumn. The analysis reveals that the influence of SST on the formation of thunderstorms over the Black Sea varies depending on the diurnal time interval, with the effect being more significant at night. At nighttime intervals, there is a clear trend of increasing mean flash frequency per case with rising SST (linear correlation coefficients range from R = 0.92 to 0.98), while during the daytime, this trend is not as evident. This indicates that, during the day, other favorable atmospheric processes have a greater influence on the formation of thunderstorms than sea-surface temperature, while in the autumn night hours, the higher SST values probably play a more significant role in thunderstorms formation, in combination with the corresponding orographic conditions. Full article
(This article belongs to the Special Issue Atmospheric Electricity (2nd Edition))
Show Figures

Figure 1

Figure 1
<p>Elevation [m] in the studied area. Positive numbers represent the elevation above the sea level, negative numbers—below the sea level. The data are taken from the GEBCO 1-min global grid (<a href="http://www.gebco.net" target="_blank">www.gebco.net</a>, accessed on 1 May 2022).</p>
Full article ">Figure 2
<p>Spatial distribution of the number of recorded flashes within 0.25° × 0.25° grid boxes over Bulgaria and the Black Sea during winter (<b>a</b>), spring (<b>b</b>), summer (<b>c</b>), and autumn (<b>d</b>) from March 2005 to February 2015. The color scale represents the number of flashes; note that the scale is different in each season.</p>
Full article ">Figure 3
<p>Flash density (flashes/km<sup>2</sup>) over the Black Sea (blue columns) and Bulgaria (white columns) during winter (DJF), spring (MAM), summer (JJA), and autumn (SON) from March 2005 to February 2015.</p>
Full article ">Figure 4
<p>Flash density (flashes/km<sup>2</sup>) for months in winter (<b>a</b>), spring (<b>b</b>), summer (<b>c</b>), and autumn (<b>d</b>) over the Black Sea (blue columns) and Bulgaria (white columns), March 2005–February 2015.</p>
Full article ">Figure 5
<p>Diurnal variation of flash density (flashes/km<sup>2</sup>) at 3-h time intervals during winter (<b>a</b>), spring (<b>b</b>), summer (<b>c</b>), and autumn (<b>d</b>) over the Black Sea (blue columns) and Bulgaria (white columns) from March 2005 to February 2015. Local time is UTC + 2 h or UTC + 3 h.</p>
Full article ">Figure 6
<p>Diurnal spatial distribution of the number of flashes at 3-h time intervals in boxes of 0.25 × 0.25 degrees for each season: (<b>a</b>)—winter (DFJ), (<b>b</b>)—spring (MAM), (<b>c</b>)—summer (JJA), (<b>d</b>)—autumn (SON) of the 10 years (Marth 2005–February 2015). The color scale represents the number of flashes; note that the scale is different in each season. Local time is UTC + 2 h or UTC + 3 h.</p>
Full article ">Figure 7
<p>Spatial distribution of number of recorded flashes in boxes of 0.25° × 0.25° over Bulgaria and the Black Sea for September (2005–2014).</p>
Full article ">Figure 8
<p>Mean values of sea surface temperature (SST) for September (2005–2014).</p>
Full article ">Figure 9
<p>Box and Whisker plot of the sea surface temperature (SST) for the cases with and without flashes for all four investigated time-intervals from September months. (median-blue line; 25th–75th percentile, blue box; 10th–90th percentile). LT = UTC + 3 h.</p>
Full article ">Figure 10
<p>Box and Whisker plot of the flash frequency (number of lightning per case) as a function of sea surface temperature (SST) during nighttime intervals ((1800–2100) UTC; (0000–0300) UTC) and daytime intervals ((0600–0900) UTC; (1200–1500) UTC). (trend line of: mean—red line, median—yellow line; 25th–75th percentile, blue box; 10th–90th percentile, whisker; the value in blue box is number of cloud cases). LT = UTC + 3 h. The analysis includes the September months of the period 2005–2014.</p>
Full article ">Figure 11
<p>Spatial–diurnal distribution of flashes at 3-h time intervals: nighttime (1800–2100) UTC, (0000–0300) UTC and daytime (0600–0900) UTC, (1200–1500) UTC for the September months (2005–2014). LT = UTC + 3 h.</p>
Full article ">
20 pages, 8335 KiB  
Article
Evaluating the Multidimensional Stability of Regional Ecosystems Using the LandTrendr Algorithm
by Lijuan Li, Jiaqiang Du, Jin Wu, Zhilu Sheng, Xiaoqian Zhu, Zebang Song, Guangqing Zhai and Fangfang Chong
Remote Sens. 2024, 16(20), 3762; https://doi.org/10.3390/rs16203762 - 10 Oct 2024
Viewed by 523
Abstract
Stability is a key characteristic for understanding ecosystem processes and evolution. However, research on the stability of complex ecosystems often faces limitations, such as reliance on single parameters and insufficient representation of continuous changes. This study developed a multidimensional stability assessment system for [...] Read more.
Stability is a key characteristic for understanding ecosystem processes and evolution. However, research on the stability of complex ecosystems often faces limitations, such as reliance on single parameters and insufficient representation of continuous changes. This study developed a multidimensional stability assessment system for regional ecosystems based on disturbances. Focusing on the lower reaches of the Yellow River Basin (LR-YRB), we integrated the remote sensing ecological index (RSEI) with texture structural parameters, and applied the Landsat-based detection of trends in disturbance and recovery (LandTrendr) algorithm to analyze the continuous changes in disturbances and recovery from 1986 to 2021, facilitating the quantification and evaluation of resistance, resilience, and temporal stability. The results showed that 72.27% of the pixels experienced 1–9 disturbances, indicating the region’s sensitivity to external factors. The maximum disturbances primarily lasted 2–3 years, with resistance and resilience displaying inverse spatial patterns. Over the 35-year period, 61.01% of the pixels exhibited moderate temporal stability. Approximately 59.83% of the pixels recovered or improved upon returning to pre-disturbance conditions after maximum disturbances, suggesting a strong recovery capability. The correlation among stability dimensions was low and influenced by disturbance intensity, underscoring the necessity for a multidimensional assessment of regional ecosystem stability based on satellite remote sensing. Full article
Show Figures

Figure 1

Figure 1
<p>Study area and spatial distribution of land use/landscape types (2020).</p>
Full article ">Figure 2
<p>Schematic diagram of ecosystem stability parameters under disturbance.</p>
Full article ">Figure 3
<p>Distribution of disturbance frequency.</p>
Full article ">Figure 4
<p>Maximum disturbance occurrence years and duration.</p>
Full article ">Figure 5
<p>Distribution of resistance (avg—the mean value; std—the standard deviation. <a href="#remotesensing-16-03762-f005" class="html-fig">Figure 5</a>, <a href="#remotesensing-16-03762-f006" class="html-fig">Figure 6</a> and <a href="#remotesensing-16-03762-f007" class="html-fig">Figure 7</a> have the same meaning).</p>
Full article ">Figure 6
<p>Distribution of resilience.</p>
Full article ">Figure 7
<p>Distribution of temporal stability.</p>
Full article ">Figure 8
<p>Spatial distribution of regime shift rate of stability.</p>
Full article ">Figure 9
<p>Pairwise scatterplot kernel density maps (<b>a</b>–<b>c</b>) and correlation coefficients (<b>d</b>) between stability dimensions (<span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 10
<p>Disturbance assessment accuracy case: (<b>a</b>) disturbance indicator parameters (top: RSEI; bottom: composite stability parameter index; “+” represents selected points); (<b>b</b>) disturbance and recovery detection results based on LandTrendr (blue line: original values of composite parameters; orange line: LandTrendr fitted curve; yellow vertical line: start year of disturbance; purple vertical line: end year of recovery).</p>
Full article ">
16 pages, 8172 KiB  
Article
Spatiotemporal Variation in Soil Wind Erosion in the Northern Slope of the Tianshan Mountains from 2000 to 2018
by Shiyu Wang and Ximeng Xu
Land 2024, 13(10), 1604; https://doi.org/10.3390/land13101604 - 2 Oct 2024
Viewed by 488
Abstract
The Northern Slope of the Tianshan Mountains (NSTM) is characterized by complex and diverse terrain, which represents a fragile ecological environment. Soil wind erosion is a key factor affecting the natural ecosystem and the social development of the region, but it has not [...] Read more.
The Northern Slope of the Tianshan Mountains (NSTM) is characterized by complex and diverse terrain, which represents a fragile ecological environment. Soil wind erosion is a key factor affecting the natural ecosystem and the social development of the region, but it has not been well understood until now. In this study, the revised wind erosion equation (RWEQ) was employed to display the spatial and temporal characteristics of soil wind erosion in the NSTM from 2000 to 2018. In addition, the main driving factors of wind erosion were analyzed. The results showed that approximately 94.25% of the NSTM experienced soil wind erosion, with a multi-year average actual soil wind erosion modulus of 6556.40 t·km−2·a−1. From 2000 to 2018, the actual soil wind erosion modulus in the NSTM showed a trend of fluctuational increase, with an increase rate of 44.65 t·km−2·a−2, but the area affected by soil wind erosion exhibited a downward trend. The wind erosion rate decreased in 76.38% of the total area, except for some areas such as Hami, with an increasing trend of soil wind erosion. The wind factor in RWEQ showed a significant linear relationship with the soil wind erosion modulus (r = 0.62, p < 0.01). Land use changes also have a critical impact on the soil wind erosion. The results of geographical detectors show that the combined effect of weather factor and vegetation factor can explain more than 60% of the changes in soil wind erosion. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of the Northern Slope of the Tianshan Mountains (NSTM).</p>
Full article ">Figure 2
<p>Spatial distribution of multi-year average actual wind erosion modulus in the NSTM from 2000 to 2018.</p>
Full article ">Figure 3
<p>Actual wind erosion modulus in the NSTM during 2000 to 2018.</p>
Full article ">Figure 4
<p>Interannual variation in wind erosion area with different classifications in the NSTM from 2000 to 2018.</p>
Full article ">Figure 5
<p>The spatial distribution of the actual wind erosion modulus change trend in the NSTM from 2000 to 2018.</p>
Full article ">Figure 6
<p>Impact of wind factor on actual soil wind erosion modulus in the NSTM from 2000 to 2018.</p>
Full article ">Figure 7
<p>Spatial distribution (<b>a</b>) and quantification of land use changes (<b>b</b>) in the NSTM from 2000 to 2018, and the actual soil wind erosion modulus for different land uses in the untransformed land use (<b>c</b>) and the transformed land use (<b>d</b>). The untransformed land use indicates that land use in 2000 and 2018 was the same, while the transformed land use indicates that land use has changed from 2000 to 2018.</p>
Full article ">Figure 8
<p>Geographic detection interaction diagram. The numbers in cells represent the q statistic values. C represents vegetation factor, WF represents weather factor, K’ represents soil roughness factor, SCF represents soil crust factor, and EF represents soil erodibility factor.</p>
Full article ">
15 pages, 10917 KiB  
Article
Geo-Sensing-Based Analysis of Urban Heat Island in the Metropolitan Area of Merida, Mexico
by Francisco A. Sánchez-Sánchez, Marisela Vega-De-Lille, Alejandro A. Castillo-Atoche, José T. López-Maldonado, Mayra Cruz-Fernandez, Enrique Camacho-Pérez and Juvenal Rodríguez-Reséndiz
Sensors 2024, 24(19), 6289; https://doi.org/10.3390/s24196289 - 28 Sep 2024
Viewed by 890
Abstract
Urban Heat Islands are a major environmental and public health concern, causing temperature increase in urban areas. This study used satellite imagery and machine learning to analyze the spatial and temporal patterns of land surface temperature distribution in the Metropolitan Area of Merida [...] Read more.
Urban Heat Islands are a major environmental and public health concern, causing temperature increase in urban areas. This study used satellite imagery and machine learning to analyze the spatial and temporal patterns of land surface temperature distribution in the Metropolitan Area of Merida (MAM), Mexico, from 2001 to 2021. The results show that land surface temperature has increased in the MAM over the study period, while the urban footprint has expanded. The study also found a high correlation (r> 0.8) between changes in land surface temperature and land cover classes (urbanization/deforestation). If the current urbanization trend continues, the difference between the land surface temperature of the MAM and its surroundings is expected to reach 3.12 °C ± 1.11 °C by the year 2030. Hence, the findings of this study suggest that the Urban Heat Island effect is a growing problem in the MAM and highlight the importance of satellite imagery and machine learning for monitoring and developing mitigation strategies. Full article
(This article belongs to the Special Issue Application of Satellite Remote Sensing in Geospatial Monitoring)
Show Figures

Figure 1

Figure 1
<p>Analyzing the interrelationships between keywords and research trends in publications retrieved from Scopus.</p>
Full article ">Figure 2
<p>Study area corresponding to the MAM, Yucatan, Mexico.</p>
Full article ">Figure 3
<p>Flowchart of the urban footprint classification process.</p>
Full article ">Figure 4
<p>LST maps of the MAM.</p>
Full article ">Figure 5
<p>LST distribution in (<b>a</b>) the MAM and (<b>b</b>) the MAM surroundings.</p>
Full article ">Figure 6
<p>Monthly variation in (<b>a</b>) average LST of selected year intervals and (<b>b</b>) their difference with the interval 2001–2005 as reference value.</p>
Full article ">Figure 7
<p>Urban footprint maps of the MAM.</p>
Full article ">Figure 8
<p>Average LST difference between the MAM and its surroundings from 2001 to 2021.</p>
Full article ">
26 pages, 2949 KiB  
Article
Study on Transportation Carbon Emissions in Tibet: Measurement, Prediction Model Development, and Analysis
by Wu Bo, Kunming Zhao, Gang Cheng, Yaping Wang, Jiazhe Zhang, Mingkai Cheng, Can Yang and Wa Da
Sustainability 2024, 16(19), 8419; https://doi.org/10.3390/su16198419 - 27 Sep 2024
Viewed by 648
Abstract
In recent years, the socio-economic development in the Tibet region of China has experienced substantial growth. However, transportation increasingly strains the region’s fragile ecological environment. Most studies overlook the accurate measurement and analysis of factors influencing traffic carbon emissions in Tibet due to [...] Read more.
In recent years, the socio-economic development in the Tibet region of China has experienced substantial growth. However, transportation increasingly strains the region’s fragile ecological environment. Most studies overlook the accurate measurement and analysis of factors influencing traffic carbon emissions in Tibet due to data scarcity. To address this, this paper applies an improved traffic carbon emissions model, using transportation turnover data to estimate emissions in Tibet from 2008 to 2020. Simultaneously, the estimated traffic carbon emissions in Tibet served as the predicted variable, and various machine learning algorithms, including Radial Basis Function Support Vector Machine (RBF-SVM), eXtreme Gradient Boosting (XGBoost), Random Forest, and Gradient Boosting Decision Tree (GBDT) are employed to conduct an initial comparison of the constructed prediction models using three-fold cross-validation and multiple evaluation metrics. The best-performing model undergoes further optimization using Grid Search (GS) and Real-coded Genetic Algorithm (RGA). Finally, the central difference method and Local Interpretable Model-Agnostic Explanation (LIME) algorithm are used for local sensitivity and interpretability analyses on twelve core variables. The results assess each variable’s contribution to the model’s output, enabling a comprehensive analysis of their impact on Tibet’s traffic carbon emissions. The findings demonstrate a significant upward trend in Tibet’s traffic carbon emissions, with road transportation and civil aviation being the main contributors. The RBF-SVM algorithm is most suitable for predicting traffic carbon emissions in this region. After GS optimization, the model’s R2 value exceeded 0.99, indicating high predictive accuracy and stability. Key factors influencing traffic carbon emissions in Tibet include civilian vehicle numbers, transportation land-use area, transportation output value, urban green coverage areas, per capita GDP, and built-up area. This paper provides a systematic framework and empirical support for measuring, predicting, and analyzing factors influencing traffic carbon emissions in Tibet. It employs innovative measurement methods, optimized machine learning models, and detailed sensitivity and interpretability analyses. The results can guide regional carbon reduction targets and promote green sustainable development. Full article
Show Figures

Figure 1

Figure 1
<p>R<sup>2</sup> scores of each model from three-fold cross-validation.</p>
Full article ">Figure 2
<p>Comparison of fitted values and actual values on the training set for the four models.</p>
Full article ">Figure 3
<p>Comparison of predicted values and actual values on the test set for the four models.</p>
Full article ">Figure 4
<p>R<sup>2</sup> scores of the three optimized models across different principal components.</p>
Full article ">Figure 5
<p>Comparison of fitted values and actual values on the training set for Model_gs_rs and Model_rga_rs.</p>
Full article ">Figure 6
<p>Comparison of predicted values and actual values on the test set for Model_gs_rs and Model_rga_rs.</p>
Full article ">Figure 7
<p>Impact of core variables on model output using the central difference method.</p>
Full article ">Figure 8
<p>Feature importance of variables using LIME mean values.</p>
Full article ">
Back to TopTop