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

Search Results (154)

Search Parameters:
Keywords = rugged terrain

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
18 pages, 7848 KiB  
Article
Effects of Climate Change and Human Activities on Streamflow in Arid Alpine Water Source Regions: A Case Study of the Shiyang River, China
by Honghua Xia, Yingqing Su, Linshan Yang, Qi Feng, Wei Liu and Jian Ma
Land 2024, 13(11), 1961; https://doi.org/10.3390/land13111961 - 20 Nov 2024
Viewed by 267
Abstract
Climate change and human activities were identified as the primary drivers of streamflow in arid alpine regions. However, limitations in observational data have resulted in a limited understanding of streamflow changes in these water sources, which hinders efforts to adapt to ongoing climate [...] Read more.
Climate change and human activities were identified as the primary drivers of streamflow in arid alpine regions. However, limitations in observational data have resulted in a limited understanding of streamflow changes in these water sources, which hinders efforts to adapt to ongoing climate change and to formulate effective streamflow management policies. Here, we use the four main tributaries in the upper reach of the Shiyang River in China as a case study to investigate the long-term trends in streamflow within arid alpine water sources, quantifying the individual contributions of climate change and human activities to these changes. The findings revealed that temperatures and precipitation in arid alpine regions have risen over the past 40 years. Although the warming trend has been significant, it has slowed in recent years. Nevertheless, three-quarters of the rivers are experiencing a decline in streamflow. The land types within the watershed remain relatively stable, with land use and cover change (LUCC) primarily occurring in the Gulang River watershed. Climate change has significantly affected streamflow change in high and rugged terrains, with an influence exceeding 70%. For example, Jingta River showed an impact of 118.79%, Zamu River 84.00%, and Huangyang River 71.43%. Human-driven LUCC, such as the expansion of cultivated and urban land, have led to increased water consumption, resulting in reduced streamflow. This effect is particularly pronounced in the low-lying and gently undulating areas of the Gulang River, where LUCC account for 78.68% of the change in streamflow. As human activities intensify and temperatures continue to rise, further declines in streamflow are projected, highlighting the urgent need for effective water resource management. These insights highlight the urgent need for targeted mitigation and adaptation strategies to confront the water scarcity challenges faced by these vulnerable regions. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area. (<b>a</b>) The geographical position of the watershed in China. (<b>b</b>) The environmental background. The abbreviations featured in the figure are listed in <a href="#app1-land-13-01961" class="html-app">Supplementary Material Table S1</a>.</p>
Full article ">Figure 2
<p>Research framework on the effects of climate change and human activities on streamflow.</p>
Full article ">Figure 3
<p>Comparison of streamflow simulated by SWAT model with monthly observation data of hydrologic stations during 1980–2016 in the JTR (<b>a</b>), ZMR (<b>b</b>), HYR (<b>c</b>), and GLR (<b>d</b>).</p>
Full article ">Figure 4
<p>LUCC from 1990 to 2010. (<b>a</b>) and (<b>b</b>) represent the LUCC of the four basins for 1990 and 2010, respectively. (<b>c</b>) indicates the land use dynamic degree of the four basins. (<b>d</b>) refers to the comprehensive land use dynamic degree of the four basins. Note: The abbreviations CL, FL, WB, UrL, UnL, HCG, MCG, and LCG represent cultivated land, forest land, water body, urban land, unutilized land, high-coverage grassland, medium-coverage grassland, and low-coverage grassland, respectively.</p>
Full article ">Figure 5
<p>The 1990–2010 land use transition matrix in the JTR (<b>a</b>), ZMR (<b>b</b>), HYR (<b>c</b>), and GLR (<b>d</b>). The unit of LUCC transfer area is km<sup>2</sup>.</p>
Full article ">Figure 6
<p>Changing trends in temperature, precipitation, and streamflow. Note: β1 and β2 represent Sen’s slope during 1985–2000 and 2001–2016, respectively. “β” in bold denotes the trends for the entire period 1980–2009 (per decade). The single asterisk (“*”) and two asterisks (“**”), represent statistical significance levels of <span class="html-italic">p</span> &lt; 0.1and <span class="html-italic">p</span> &lt; 0.05, respectively.</p>
Full article ">Figure 7
<p>Sen’s slope of temperature (<b>a</b>), precipitation (<b>b</b>), and streamflow (<b>c</b>). Note: The single asterisks (“*”), two asterisks (”**”), and three asterisks (”***”) represent statistical significance levels of <span class="html-italic">p</span> &lt; 0.1, <span class="html-italic">p</span> &lt; 0.05, and <span class="html-italic">p</span> &lt; 0.01, respectively.</p>
Full article ">Figure 8
<p>Changes in monthly streamflow impacted by climate change and LUCC in the JTR (<b>a</b>), ZMR (<b>b</b>), HYR (<b>c</b>), and GLR (<b>d</b>).</p>
Full article ">Figure 9
<p>Streamflow suitability management in arid alpine regions.</p>
Full article ">
28 pages, 31167 KiB  
Article
Optimizing GEDI Canopy Height Estimation and Analyzing Error Impact Factors Under Highly Complex Terrain and High-Density Vegetation Conditions
by Runbo Chen, Xinchuang Wang, Xuejie Liu and Shunzhong Wang
Forests 2024, 15(11), 2024; https://doi.org/10.3390/f15112024 - 17 Nov 2024
Viewed by 555
Abstract
The Global Ecosystem Dynamics Investigation (GEDI) system provides essential data for estimating forest canopy height on a global scale. However, factors such as complex topography and dense canopy can significantly reduce the accuracy of GEDI canopy height estimations. We selected the South Taihang [...] Read more.
The Global Ecosystem Dynamics Investigation (GEDI) system provides essential data for estimating forest canopy height on a global scale. However, factors such as complex topography and dense canopy can significantly reduce the accuracy of GEDI canopy height estimations. We selected the South Taihang region of Henan Province, China, as our study area and proposed an optimization framework to improve GEDI canopy height estimation accuracy. This framework includes correcting geolocation errors in GEDI footprints, screening and analyzing features that affect estimation errors, and combining two regression models with feature selection methods. Our findings reveal a geolocation error of 4 to 6 m in GEDI footprints at the orbital scale, along with an overestimation of GEDI canopy height in the South Taihang region. Relative height (RH), waveform characteristics, topographic features, and canopy cover significantly influenced the estimation error. Some studies have suggested that GEDI canopy height estimates for areas with high canopy cover lead to underestimation, However, our study found that accuracy increased with higher canopy cover in complex terrain and dense vegetation. The model’s performance improved significantly after incorporating the canopy cover parameter into the optimization model. Overall, the R2 of the best-optimized model was improved from 0.06 to 0.61, the RMSE was decreased from 8.73 m to 2.23 m, and the rRMSE decreased from 65% to 17%, resulting in an accuracy improvement of 74.45%. In general, this study reveals the factors affecting the accuracy of GEDI canopy height estimation in areas with complex terrain and dense vegetation cover, on the premise of minimizing GEDI geolocation errors. Employing the proposed optimization framework significantly enhanced the accuracy of GEDI canopy height estimates. This study also highlighted the crucial role of canopy cover in improving the precision of GEDI canopy height estimation, providing an effective approach for forest monitoring in such regions and vegetation conditions. Future studies should further improve the classification of tree species and expand the diversity of sample tree species to test the accuracy of canopy height estimated by GEDI in different forest structures, consider the distortion of optical remote sensing images caused by rugged terrain, and further mine the information in GEDI waveforms so as to enhance the applicability of the optimization framework in more diverse forest environments. Full article
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>In the figure, (<b>a</b>) shows the location of Henan Province in China, (<b>b</b>) illustrates the study area’s location within Henan Province, and (<b>c</b>) presents the DEM of the study area, with each individual area number corresponding to the ALS aerial flight areas.</p>
Full article ">Figure 2
<p>CHM raster maps based on ALS acquisition: bottom images are true color images of Sentinel-2 in May 2023; black dots are GEDI footprints.</p>
Full article ">Figure 3
<p>The distribution of slope and canopy cover within the aerial flight zone after cropping based on ALS slope and canopy cover raster maps. Panel (<b>a</b>) shows the slope distribution following the cropping of the ALS airspace slope raster map using the extent of forested land from the Land Use Survey. Panel (<b>b</b>) shows the slope distribution across ALS airspace. Panels (<b>c</b>,<b>d</b>) show the canopy cover, where the vertical axis represents the number of raster pixels and the horizontal axis indicates the canopy cover (0-1).</p>
Full article ">Figure 4
<p>A square fishing net with a length and width of 5 m was used to calculate statistics of the DEM and slope within each grid, and the mean value, range, standard deviation, and mean slope of the DEM were calculated. Due to the huge amount of data, the data of all grids were not counted, but 5000 grids were randomly selected in each ALS region for statistics. The change from blue to red means the density goes from small to large.</p>
Full article ">Figure 5
<p>The principle of geolocation error correction is illustrated as follows: (<b>a</b>) displays the displacement mode of the footprint, where the red circle in the center represents the original GEDI location, and the cyan spot indicates the position after displacement. The angular step is set at 30°, while the distance step is 2 m. (<b>b</b>) shows the waveform corresponding to the GEDI location. (<b>c</b>) depicts the simulated waveform from ALS, and (<b>d</b>) presents the aligned GEDI waveform and ALS simulated waveform.</p>
Full article ">Figure 6
<p>Overall frame flowchart.</p>
Full article ">Figure 7
<p>The effect of geolocation error correction for a single footprint. Toploc and botloc refer to the start and end positions of the signal, respectively. Panels (<b>a</b>,<b>b</b>) display the original and corrected geolocation waveforms of the complex footprint, while panels (<b>c</b>,<b>d</b>) show the original and corrected geolocation waveforms of the simple footprint.</p>
Full article ">Figure 8
<p>The statistics of all R averages after displacing footprints to the same location for the same acquisition date. Each polar plot represents the average correction effect of geolocation errors for all footprints corresponding to the same acquisition date. The top label of each polar plot indicates the data acquisition date in the format YYYYDDD.</p>
Full article ">Figure 9
<p>R-values between individual features and GEDI canopy height estimation error, All feature parameters in the figure are significantly correlated (<math display="inline"><semantics> <mrow> <mi>p</mi> <mo>≤</mo> <mn>0.05</mn> </mrow> </semantics></math>), with a positive correlation in blue and a negative correlation in orange.</p>
Full article ">Figure 10
<p>Box plots of error distribution in different intervals of each feature with the absolute value of R above 0.3. The left vertical axis is the error (m) and the right vertical axis is the RMSE (m) of RH96 and CHM96.</p>
Full article ">Figure 11
<p>This figure shows the importance of each feature parameter with respect to the residuals: the upper figure shows the top 30 feature parameters in terms of importance, and the lower figure shows the thumbnail of the importance distribution of all feature parameters, where the blue part is the detailed distribution of the importance of the top 30 features in the upper figure.</p>
Full article ">Figure 12
<p>Box plots of error distribution in different intervals of each feature with the absolute value of RF importance above 1%. The left vertical axis is error (m) and the right vertical axis is the RMSE (m) of RH_96 and CHM_96.</p>
Full article ">Figure 13
<p>In the case of selecting different numbers of features, the model effects of various combinations of regression models and feature extraction methods are presented. The results are organized in the vertical coordinates from top to bottom in the order of <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>R</mi> <mi>M</mi> <mi>S</mi> <mi>E</mi> </mrow> </semantics></math> (m), and <math display="inline"><semantics> <mrow> <mi>r</mi> <mi>R</mi> <mi>M</mi> <mi>S</mi> <mi>E</mi> </mrow> </semantics></math> (%). The horizontal coordinates indicate the number of feature parameters used in the model.</p>
Full article ">Figure 14
<p>The left and right panels show the data distribution of RH_96 and RHT_96 with CHM96, respectively.</p>
Full article ">Figure 15
<p>The upper panel is a localized thumbnail of the remote sensing image, the blue part is the non-shadowed area, the white part is the shadowed area, and the lower two panels are the reflectance distributions of the red, green, and blue bands in the shadowed and non-shadow areas, respectively.</p>
Full article ">
18 pages, 1660 KiB  
Article
Evaluating the Soil Properties of Different Land Use Types in the Deviskel Watershed in the Hilly Region of Northeast Türkiye
by Esin Erdoğan Yüksel and Gökhan Yavuz
Sustainability 2024, 16(22), 9732; https://doi.org/10.3390/su16229732 - 8 Nov 2024
Viewed by 551
Abstract
Land use is a remarkable human-induced change that has redesigned the Earth’s surface since the beginning of civilization. Due to the combination of rugged terrain and low-income levels in rural areas, people in watershed regions often resort to overexploiting forests, agricultural land, and [...] Read more.
Land use is a remarkable human-induced change that has redesigned the Earth’s surface since the beginning of civilization. Due to the combination of rugged terrain and low-income levels in rural areas, people in watershed regions often resort to overexploiting forests, agricultural land, and grasslands beyond their capacity. As a result of these spatio-temporal changes in land use, various soil properties undergo changes. This study aims to determine the changes in some physical (texture, bulk weight, particle density, total porosity), hydro-physical (water holding capacity, permeability, field capacity, wilting point), physico-chemical (organic matter, pH, electrical conductivity), and erodibility (dispersion ratio, colloid–moisture equivalent ratio, erosion ratio, clay ratio, aggregate stability and K-factor of Universal Soil Loss Equation-USLE) properties of soil depending on land use in the Deviskel Watershed in the city of Artvin in Türkiye. For this purpose, disturbed (composite) and undisturbed (cylinder) soil samples were taken from a 0 to 20 cm depth at 108 different points in the determined areas (36 from forests, 36 from agricultural areas, and 36 from grassland areas). It was determined that 15 of the 19 soil properties examined showed statistical differences depending on the change in land use. All the examined soil properties, except for clay content, particle density, dispersion ratio, and aggregate stability, were found to be statistically significantly affected by the change in land use, and the reasons behind these changes were discussed. The particle density had the lowest coefficient of variation value (15.26%) while electrical conductivity had the highest coefficient of variation value (91.25%). According to erosion tendencies, all watershed soils were found to be susceptible to erosion. The average aggregate stability was 88.52% in forest soils, 84.84% in agricultural soils, and 85.48% in grassland soils. The average USLE-K factor was determined to be 0.22 for forests, while it was determined to be 0.17 and 0.18 for agriculture and grassland areas, respectively. According to the USLE-K factor, 68.37% of the watershed was dominated by moderately erodible soils, while 31.63% consisted of highly erodible soils. Based on the colloid–moisture equivalent ratio, erosion ratio, and clay ratio, which are statistically different erodibility features, the grassland soils of the research area were found to be more susceptible to erosion than forest and agricultural soils. In terms of aggregate stability, which indicates resistance to water erosion, forest areas had higher values, while agricultural lands were more prone to erosion. Full article
(This article belongs to the Section Soil Conservation and Sustainability)
Show Figures

Figure 1

Figure 1
<p>Location of Deviskel Watershed.</p>
Full article ">Figure 2
<p>The land use map and soil sampling points of the Deviskel Watershed.</p>
Full article ">Figure 3
<p>Erodibility map of the study area.</p>
Full article ">Figure 4
<p>Distribution of the range of K-factors grouped in percentages.</p>
Full article ">
25 pages, 2178 KiB  
Article
User Sentiment Analysis of the Shared Charging Service for China’s G318 Route
by Mei Wang, Siva Shankar Ramasamy, Xi Yu, Mutong Liu, Ahmad Yahya Dawod and Huayue Chen
Electronics 2024, 13(22), 4335; https://doi.org/10.3390/electronics13224335 - 5 Nov 2024
Viewed by 496
Abstract
Shared charging services have gained popularity for their contribution to green travel. Accurately identifying the core factors that influence user experience (UX) not only enhances service quality and optimizes user satisfaction, but also promotes the dissemination of green travel concepts. However, the influencing [...] Read more.
Shared charging services have gained popularity for their contribution to green travel. Accurately identifying the core factors that influence user experience (UX) not only enhances service quality and optimizes user satisfaction, but also promotes the dissemination of green travel concepts. However, the influencing factors and their mechanisms vary significantly across regions, particularly along the Chengdu–Lhasa (G318) route, which features large elevation changes, diverse climatic conditions, rugged terrain, and frequent geological disasters, making the influencing factors particularly complex. This study analyzes comment texts from 38 shared charging stations along the G318 route in the e-Charging APP, totaling 15,214 comments. A comprehensive approach is employed, including high-frequency word analysis, term frequency–inverse document frequency (TF-IDF) comparison, co-occurrence semantic network and co-word matrix feature correlation analysis, Latent Dirichlet Allocation (LDA) topic modeling, and sentiment analysis. This multifaceted analysis explores core themes, user viewpoints, and sentiments in the comments, focusing on users’ perspectives on service quality, usage experience, and environmental impact of the charging stations. The findings indicate that charging speed, service attitude, environment, operational status of hardware and software, and pricing are key factors influencing user sentiment. Users have a high demand for the perfection of supporting facilities of shared charging stations, directly affecting user satisfaction and indirectly influencing the brand image and market competitiveness of enterprises. Full article
(This article belongs to the Special Issue Intelligent Data Analysis and Learning)
Show Figures

Figure 1

Figure 1
<p>Route map of the South Route of G318.</p>
Full article ">Figure 2
<p>Distribution map of 38 shared charging stations.</p>
Full article ">Figure 3
<p>Overall evaluation chart of 38 shared charging stations.</p>
Full article ">Figure 4
<p>Word cloud diagram.</p>
Full article ">Figure 5
<p>Co-occurrence semantic network diagram.</p>
Full article ">Figure 6
<p>The distribution of emotion value and quantity.</p>
Full article ">
19 pages, 4951 KiB  
Article
Combining Remote Sensing Data and Geochemical Properties of Ultramafics to Explore Chromite Ore Deposits in East Oltu Erzurum, Turkey
by Amr Abd El-Raouf, Fikret Doğru, Özgür Bilici, Islam Azab, Sait Taşci, Lincheng Jiang, Kamal Abdelrahman, Mohammed S. Fnais and Omar Amer
Minerals 2024, 14(11), 1116; https://doi.org/10.3390/min14111116 - 2 Nov 2024
Viewed by 530
Abstract
The present research’s main objective was to apply thorough exploration approaches that combine remote sensing data with geochemical sampling and analysis to predict and identify potential chromitite locations in a complex geological site, particularly in rugged mountainous terrain, and differentiate the ultramafic massif [...] Read more.
The present research’s main objective was to apply thorough exploration approaches that combine remote sensing data with geochemical sampling and analysis to predict and identify potential chromitite locations in a complex geological site, particularly in rugged mountainous terrain, and differentiate the ultramafic massif containing chromitite orebodies from other lithologies. The ultramafic massif forming the mantle section of the Kırdağ ophiolite, located within the Erzurum–Kars Ophiolite Zone and emerging in the east of Oltu district (Erzurum, NE Turkey), was selected as the study area. Optimum index factor (OIF), false-color composite (FCC), decorrelation stretch (DS), band rationing (BR), minimum noise fraction (MNF), and principal and independent component analyses (PCA-ICA) were performed to differentiate the lithological features and identify the chromitite host formations. The petrography, mineral chemistry, and whole-rock geochemical properties of the harzburgites, which are the host rocks of chromitites in the research area, were evaluated to verify and confirm the remote sensing results. In addition, detailed petrographic properties of the pyroxenite and chromitite samples are presented. The results support the existence of potential chromitite formations in the mantle section of the Kırdağ ophiolite. Our remote sensing results also demonstrate the successful detection of the spectral anomalies of this ultramafic massif. The mineral and whole-rock geochemical features provide clear evidence of petrological processes, such as partial melting and melt–peridotite interactions during the harzburgite formation. The chromian spinels’ Cr#, Mg#, Fe3+, Al2O3, and TiO2 concentrations indicate that the harzburgite formed in a fore-arc environment. The Al2O3 content and Mg# of the pyroxenes and the whole-rock Al2O3/MgO ratio and V contents of the harzburgite are also compatible with these processes. Consequently, the combined approaches demonstrated clear advantages over conventional chromitite exploration techniques, decreasing the overall costs and supporting the occurrence of chromite production at the site. Full article
(This article belongs to the Section Mineral Exploration Methods and Applications)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Illustration depicting the overall distribution of ophiolite belts and prominent tectonic suture areas across Turkey, with modifications derived from [<a href="#B21-minerals-14-01116" class="html-bibr">21</a>,<a href="#B22-minerals-14-01116" class="html-bibr">22</a>]. (<b>b</b>) Geological map specifically detailing the Kırdağ ophiolite, adapted from [<a href="#B23-minerals-14-01116" class="html-bibr">23</a>].</p>
Full article ">Figure 2
<p>A flow chart showing the methodology of the combined approach applied in the investigated area.</p>
Full article ">Figure 3
<p>(<b>a</b>) False-color composite image in the RGB of ASTER bands (2, 3, 5). Ultra = ultramafic and Gab = gabbro; (<b>b</b>) false-color composite image in the RGB of ASTER bands (8, 3, 1); (<b>c</b>) false-color composite image in the RGB of ASTER bands (1, 2, 3); and (<b>d</b>) false-color composite image in the RGB of the ASTER band ratio (4/8, 4/1, and 3/2 × 4/3).</p>
Full article ">Figure 4
<p>(<b>a</b>) Grayscale image of the ASTER band ratio (3/4), (<b>b</b>) false-color composite image in the RGB of ASTER MNF (1, 2, 3), and (<b>c</b>) false-color composite image in the RGB of ASTER MNF (9, 6, 4).</p>
Full article ">Figure 5
<p>(<b>a</b>) False-color composite image in the RGB of ASTER PCs (1, 2, 3), (<b>b</b>) false-color composite image in the RGB of ASTER ICs (1, 2, 3), (<b>c</b>) false-color composite image in the RGB of ASTER PCs (5, 4, 2), and (<b>d)</b> false-color composite image in the RGB of ASTER (b4, PC1, PC2).</p>
Full article ">Figure 6
<p>(<b>a</b>) Field snapshot exhibiting the juxtaposition between the host harzburgite and adjacent lithologies (including dunite, pyroxenite, and chromitite) within the research locale. (<b>b</b>–<b>h</b>) Detailed close-up images showcasing the characteristics of the harzburgite, dunite, pyroxenite, and chromitite pod within the mantle section of the Kırdağ ophiolite.</p>
Full article ">Figure 7
<p>Thin-section photomicrographs of the harzburgite (<b>a</b>,<b>b</b>), dunite (<b>c</b>,<b>d</b>), pyroxenite (<b>e</b>), and chromitite (<b>f</b>). Microphotos (<b>a</b>,<b>c</b>,<b>e</b>) were taken under cross-polarized light, (<b>b</b>,<b>d</b>) were taken under plane-polarized light, and (<b>f</b>) was taken using a reflecting microscope for a chromitite ore sample. In the figure, spnl = chromian spinel, ol = olivine, opx = orthopyroxene, cpx = clinopyroxene, and srpn = serpentine.</p>
Full article ">Figure 8
<p>Graphical representations (<b>a</b>,<b>b</b>) illustrating the relationship between the Cr-number, Mg-number, and Cr-number and TiO<sub>2</sub> for chromian spinels within the harzburgite. Abyssal, fore-arc peridotite, and boninite fields are derived from [<a href="#B68-minerals-14-01116" class="html-bibr">68</a>], whereas the reaction fields and partial melting trend are referenced from [<a href="#B66-minerals-14-01116" class="html-bibr">66</a>,<a href="#B69-minerals-14-01116" class="html-bibr">69</a>], respectively. The diagram (<b>c</b>) presents the relationship between TiO<sub>2</sub> and Fe<sup>3+</sup>-number for chromian spinels. The fields representing the Mid-Ocean Ridge (MOR) and Supra-Subduction Zone (SSZ) contexts are based on data from [<a href="#B64-minerals-14-01116" class="html-bibr">64</a>]. Diagram (<b>d</b>) illustrates the correlation between Al<sub>2</sub>O<sub>3</sub> and Mg-number for orthopyroxene. In addition, (<b>e</b>) depicts the exact correlation for clinopyroxene. The fields representing abyssal and SSZ peridotites are based on data from [<a href="#B67-minerals-14-01116" class="html-bibr">67</a>]. The diagram (<b>f</b>) illustrates the relationship between V (ppm) and Al<sub>2</sub>O<sub>3</sub>/MgO for the harzburgites from the Kırdağ ophiolite. The fields representing the fore-arc and abyssal peridotites are derived from data compiled by [<a href="#B65-minerals-14-01116" class="html-bibr">65</a>,<a href="#B67-minerals-14-01116" class="html-bibr">67</a>], respectively.</p>
Full article ">Figure 9
<p>Sentinel-2-L2A True Color B4, B3, B2 showing the location of high-potential chromite-bearing mineralized zones based on integrating remote sensing and geochemical results with field validation.</p>
Full article ">
17 pages, 11423 KiB  
Article
Spatiotemporal Variability of Soil Erosion in the Pisha Sandstone Region: Influences of Precipitation and Vegetation
by Zhenqi Yang, Jianying Guo, Fucang Qin, Yan Li, Xin Wang, Long Li and Xinyu Liu
Sustainability 2024, 16(21), 9313; https://doi.org/10.3390/su16219313 - 26 Oct 2024
Viewed by 691
Abstract
The Pisha sandstone area, situated in the upper and middle reaches of the Yellow River in China, is characterized by severe soil and water erosion, making it one of the most critical regions on the Loess Plateau. The rugged terrain and exposed bedrock [...] Read more.
The Pisha sandstone area, situated in the upper and middle reaches of the Yellow River in China, is characterized by severe soil and water erosion, making it one of the most critical regions on the Loess Plateau. The rugged terrain and exposed bedrock complicate management efforts for this area, posing challenges for accurate forecasting using soil erosion models. Through an analysis of terrain, vegetation, and precipitation impacts on soil erosion, this study offers theoretical support for predicting soil erosion within the exposed Pisha sandstone area of the Loess Plateau. This has substantial implications for guiding water and soil conservation measures in this region. Focusing on China’s exposed sandstone area within the Geqiugou watershed, temporal and spatial changes in vegetation cover and land use from 1990 to 2020 were analyzed. The result shows that, from 1990 to 2020, the grassland area has exhibited a consistent downward trend, with successive reductions of 64.86% to 59.46%. The area of low vegetation cover witnessed a significant decline of 59.29% in 2020 compared to that in 1990. The moderate erosion area decreased from 84.52 to 57.17 km2. The significant reduction in soil and water loss can be attributed to the expansion of forest and grassland areas, with the implementation of the Grain for Green project serving as a key policy driver for facilitating this expansion. This study provided a good example of combining rainfall with vegetation coverage to fast estimation soil erosion. A mathematical relationship between the vegetation rainfall coupling index (RV) and soil erosion was established with strong fitting effects, enabling estimation of the soil erosion volume under varying slope conditions within Pisha sandstone areas. The main focus of future soil and water conservation in the Pisha sandstone area should be on effectively managing the channel slope and minimizing exposed bedrock areas through a combination of slope cutting, the application of anticorrosive materials, and the implementation of artificial vegetation planting. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of the study site.</p>
Full article ">Figure 2
<p>The field for monitoring soil and water loss in Ordos City.</p>
Full article ">Figure 3
<p>Interannual change characteristics of ground utilization form in the study area from 1990 to 2020.</p>
Full article ">Figure 4
<p>Spatial distribution map of <span class="html-italic">VC</span> in the watershed.</p>
Full article ">Figure 5
<p>Spatial distribution map of the hydraulic erosion modulus.</p>
Full article ">Figure 6
<p>Relationship between vegetation coverage and soil erosion under different slopes. (<b>a</b>) The vegetation coverage of various slopes from 1990 to 2020. (<b>b</b>) The soil erosion amount of various slopes from 1990 to 2020. (<b>c</b>) The correlation between vegetation coverage and soil erosion under varying slope conditions.</p>
Full article ">Figure 6 Cont.
<p>Relationship between vegetation coverage and soil erosion under different slopes. (<b>a</b>) The vegetation coverage of various slopes from 1990 to 2020. (<b>b</b>) The soil erosion amount of various slopes from 1990 to 2020. (<b>c</b>) The correlation between vegetation coverage and soil erosion under varying slope conditions.</p>
Full article ">Figure 7
<p>Relationship between annual rainfall and soil erosion at different gradients.</p>
Full article ">Figure 8
<p>Relationship between RV and soil erosion at different gradients.</p>
Full article ">
15 pages, 3362 KiB  
Article
Assessing Atlantic Kelp Forest Restoration Efforts in Southern Europe
by Alexandre F. S. Marques, Álvaro Sanchéz-Gallego, Rodrigo R. Correia, Isabel Sousa-Pinto, Silvia Chemello, Inês Louro, Marco F. L. Lemos and João N. Franco
Sustainability 2024, 16(21), 9176; https://doi.org/10.3390/su16219176 - 23 Oct 2024
Viewed by 948
Abstract
Kelp forests are essential marine ecosystems increasingly compromised by human activities. Effective reforestation strategies are urgently needed, and the “green gravel” method is a viable tool already used in some European regions. This study aimed to assess the success of this method using [...] Read more.
Kelp forests are essential marine ecosystems increasingly compromised by human activities. Effective reforestation strategies are urgently needed, and the “green gravel” method is a viable tool already used in some European regions. This study aimed to assess the success of this method using the native Kelp species Laminaria ochroleuca on the Portuguese coastline. Cultures of green gravel were reared until the specimens reached a size of approximately 3 cm. The gravel was then deployed at selected sites in Peniche, Berlengas, and Cascais. Over an eight-month period, scientific scuba divers monitored the integration of Kelp, along with associated fish, invertebrate, and algae communities. Nutrient availability, temperature, water movement, substrate type, and Rugosity Index (RI) were also measured. The highest success rate was 12% in Consolação, with Elefante and Galos (Berlengas) reaching 7% and 4%, respectively. By the end of the monitoring period, Cascais had no remaining Kelp on green gravel. Present data suggest that higher success is dependent on less rugged and higher RI topography. Higher grazing pressure, rougher terrain, and unexpected sedimentation appear to be the main obstacles to deployment success. Solid knowledge (biologic and topographic) on the restoration site, starting restoration actions near already established Kelp forests, and significantly scaling up restoration efforts could substantially improve the success of the green gravel method in future reforestation campaigns. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>A</b>) Section of the western continental Portuguese coast where the deployments took place. Three distinct zones were selected: (<b>B</b>) The coast of Peniche, with Marques-Neves as the reference Kelp forest (I) and Consolação (II); (<b>C</b>) In the Berlengas Islands at the sites Elefante (III) and Galos (IV); (<b>D</b>) and in the Cascais area at Boca do Inferno (V).</p>
Full article ">Figure 2
<p>(<b>A</b>) Technique used to deploy the gravel from the surface; (<b>B</b>) Scientific divers deployed buoys marking the area for deployment (100 m<sup>2</sup>) in Berlengas at the Galos site; (<b>C</b>) close-up of the deployed green gravel.</p>
Full article ">Figure 3
<p>Percentage of each bottom substratum category (reef plateau, boulders, and sand) measured in 25 m transects at the different studied sites (<span class="html-italic">n</span> = 5). Numbers within each bar indicate the site’s Rugosity Index (RI) achieved (<span class="html-italic">n</span> = 5).</p>
Full article ">Figure 4
<p>(<b>A</b>) Average seawater temperature (measured in situ); (<b>B</b>) sea surface temperature—SST; (<b>C</b>) near surface chlorophyll <span class="html-italic">a</span> level—Chla <span class="html-italic">a</span>, and (<b>D</b>) absolute average water movement in all monitored sites. Statistically significant different groups (a, b, and c) are shown (Tukey’s HSD; <span class="html-italic">p</span> ≤ 0.05). Error bars are ±2 sd of means (n.a. = not assessed).</p>
Full article ">Figure 5
<p>(<b>A</b>) Concentration of phosphates, (<b>B</b>) nitrates, and (<b>C</b>) nitrites at depth at all monitored sites. Statistically significant different groups (a and b) are shown (Tukey’s HSD; <span class="html-italic">p</span> ≤ 0.05). Error bars are ±2 sd of means. (<span class="html-italic">n</span> = 9 per site).</p>
Full article ">Figure 6
<p>(<b>A</b>) Gravel retention and (<b>B</b>) deployment success at the four deployment sites in the three monitoring moments after 3, 6, and 8 months of deployment.</p>
Full article ">Figure 7
<p>Coefficient plot of the first component of the PLS regression model. The model included data from all the deployment sites to explain the green gravel success index. Importance was deemed as high if the coefficient was higher than 0.15 (<span class="html-italic">n</span> = 80).</p>
Full article ">
22 pages, 25663 KiB  
Article
Trade-Off and Coordination between Development and Ecological Protection of Urban Agglomerations along Rivers: A Case Study of Urban Agglomerations in the Shandong Section of the Lower Yellow River
by Anbei Liu, Tingting Yan, Shengxiang Shi, Weijun Zhao, Sihang Ke and Fangshu Zhang
Land 2024, 13(9), 1368; https://doi.org/10.3390/land13091368 - 26 Aug 2024
Viewed by 596
Abstract
Urban development of clusters situated along rivers significantly affects the health of the river ecosystems, the quality of urban environments, and the overall well-being of local communities. Ecosystem service supply value (ESSV) measures the delivery of ecosystem goods and services within a specific [...] Read more.
Urban development of clusters situated along rivers significantly affects the health of the river ecosystems, the quality of urban environments, and the overall well-being of local communities. Ecosystem service supply value (ESSV) measures the delivery of ecosystem goods and services within a specific timeframe in a particular area. Using the lower Yellow River urban agglomeration (Shandong section) as a case, we comprehensively applied land use structure and intensity change analysis, quantitative calculation of ESS, and geographical probe methods to unveil ESS and its mechanism of response to the spatio-temporal evolution of the intensity of land use in urban agglomeration along the river. The key results were as follows: (1) Over the past two decades, farmland and construction land areas have continued to decrease and increase, respectively, with the intensity of land use change being highest from 2005 to 2010. (2) ESS has continued to rise over the past 20 years, with the income in 2020 being 11.142 billion yuan, an increase of 31.13%. The “low-value areas” are mainly concentrated in Liaocheng City, Dezhou City, and Tai’an City, which are characterized by predominantly flat terrains where farmland constitutes the principal land use type. Conversely, “high-value areas” are largely in the counties bordering the Yellow River, including the upper estuary in the north and the rugged, southeastern terrains. (3) Areas with concentrated ESSV were primarily localized in the northern estuary area and along the Yellow River in a scattered point-like pattern. The spatial distribution of hotspots has become increasingly concentrated, transitioning from points to planes. Conversely, cold spots initially increased in number before subsequently decreasing. Waterbody was the most sensitive ESSV-determining factor. (4) The spatial heterogeneity of ESSV emerges as a consequence of the interaction of multiple factors, and among these interactions, those involving NDVI and POP contain the greatest explanatory power. Our findings are expected to offer a scientific foundation for optimizing land spatial patterns and enhancing ecological management in the lower Yellow River region. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>Land use in the study area in 2000–2020.</p>
Full article ">Figure 2 Cont.
<p>Land use in the study area in 2000–2020.</p>
Full article ">Figure 3
<p>Sankey map depicting the change in land use from 2000 to 2020.</p>
Full article ">Figure 4
<p>Time-intensity analysis of four time intervals: 2000–2005, 2005–2010, 2010–2015, and 2015–2020.</p>
Full article ">Figure 5
<p>The category intensities in 2000–2005, 2005–2010, 2010–2015, and 2015–2020 are shown in figures (<b>a</b>–<b>d</b>), respectively.</p>
Full article ">Figure 6
<p>Transition intensity of the given category gains during the four time intervals. The green lines and orange lines in the figures represent the intensity of the transition from M to other categories, and from other categories to M, respectively.</p>
Full article ">Figure 6 Cont.
<p>Transition intensity of the given category gains during the four time intervals. The green lines and orange lines in the figures represent the intensity of the transition from M to other categories, and from other categories to M, respectively.</p>
Full article ">Figure 7
<p>Characteristics of the spatial distribution of ESSV from 2000 to 2020. The blue and red colors indicate the low-value and high-value areas, respectively.</p>
Full article ">Figure 8
<p>Land use type transfer and distribution of ESSV hot spots in the study area.</p>
Full article ">Figure 9
<p>Results of the spatial differentiation of factors driving ESSV in the study area.</p>
Full article ">Figure 10
<p>Interactive effect of driving factors from 2000 to 2020. Note: X<sub>1</sub>: population; X2: per capita GDP; X3: precipitation; X4: temperature; X5: normalized difference vegetation index; X6: night light index; X7: distance from the Yellow River.</p>
Full article ">
27 pages, 22313 KiB  
Article
Landslide Risk Assessments through Multicriteria Analysis
by Fatma Zohra Chaabane, Salim Lamine, Mohamed Said Guettouche, Nour El Islam Bachari and Nassim Hallal
ISPRS Int. J. Geo-Inf. 2024, 13(9), 303; https://doi.org/10.3390/ijgi13090303 - 25 Aug 2024
Viewed by 1703
Abstract
Natural risks comprise a whole range of disasters and dangers, requiring comprehensive management through advanced assessment, forecasting, and warning systems. Our specific focus is on landslides in difficult terrains. The evaluation of landslide risks employs sophisticated multicriteria models, such as the weighted sum [...] Read more.
Natural risks comprise a whole range of disasters and dangers, requiring comprehensive management through advanced assessment, forecasting, and warning systems. Our specific focus is on landslides in difficult terrains. The evaluation of landslide risks employs sophisticated multicriteria models, such as the weighted sum GIS approach, which integrates qualitative parameters. Despite the challenges posed by the rugged terrain in Northern Algeria, it is paradoxically home to a dense population attracted by valuable hydro-agricultural resources. The goal of our research is to study landslide risks in these areas, particularly in the Mila region, with the aim of constructing a mathematical model that integrates both hazard and vulnerability considerations. This complex process identifies threats and their determining factors, including geomorphology and socio-economic conditions. We developed two algorithms, the analytic hierarchy process (AHP) and the fuzzy analytic hierarchy process (FAHP), to prioritize criteria and sub-criteria by assigning weights to them, aiming to find the optimal solution. By integrating multi-source data, including satellite images and in situ measurements, into a GIS and applying the two algorithms, we successfully generated landslide susceptibility maps. The FAHP method demonstrated a higher capacity to manage uncertainty and specialist assessment errors. Finally, a comparison between the developed risk map and the observed risk inventory map revealed a strong correlation between the thematic datasets. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of Mila province in the northeast of Algeria.</p>
Full article ">Figure 2
<p>Hierarchical structure of the AHP method.</p>
Full article ">Figure 3
<p>Landsat 8 satellite imagery (<b>left</b>) and the DEM of Mila province (<b>right</b>).</p>
Full article ">Figure 4
<p>Rainfall map.</p>
Full article ">Figure 5
<p>Altitude classes map.</p>
Full article ">Figure 6
<p>Hydrographic distance map.</p>
Full article ">Figure 7
<p>Aspect classes map.</p>
Full article ">Figure 8
<p>NDMI map.</p>
Full article ">Figure 9
<p>Lithology map.</p>
Full article ">Figure 10
<p>NDVI map.</p>
Full article ">Figure 11
<p>Drainage density map.</p>
Full article ">Figure 12
<p>Slope map.</p>
Full article ">Figure 13
<p>Land use map.</p>
Full article ">Figure 14
<p>Slide inventory map of Mila province.</p>
Full article ">Figure 15
<p>Landslide susceptibility map (<b>AHP</b>).</p>
Full article ">Figure 16
<p>Landslide susceptibility map (<b>FAHP</b>).</p>
Full article ">
22 pages, 5097 KiB  
Article
Disentangling the Response of Vegetation Dynamics to Natural and Anthropogenic Drivers over the Minjiang River Basin Using Dimensionality Reduction and a Structural Equation Model
by Yujie Kang, Ziqin Wang, Binni Xu, Wenjie Shen, Yu Chen, Xiaohui Zhou, Yanguo Liu, Tingbin Zhang, Guoyan Wang, Yuling Jia and Jingji Li
Forests 2024, 15(8), 1438; https://doi.org/10.3390/f15081438 - 15 Aug 2024
Viewed by 837
Abstract
Located at an average elevation of approximately 2000 m, the Minjiang River Basin (MJB), a key tributary of the Upper Yangtze River, straddles the Western Sichuan Plateau and the Sichuan Basin. Vegetation here is crucial for human life, providing oxygen and energy. However, [...] Read more.
Located at an average elevation of approximately 2000 m, the Minjiang River Basin (MJB), a key tributary of the Upper Yangtze River, straddles the Western Sichuan Plateau and the Sichuan Basin. Vegetation here is crucial for human life, providing oxygen and energy. However, the influence of climatic variables, human activities, and rugged terrain on vegetation vitality is still debated. This study mainly leverages data from the Normalized Difference Vegetation Index (NDVI), meteorological stations data, and land use data. Analytical techniques include trend analysis, partial correlation coefficient analysis (PCC), principal component analysis (PCA), and partial least squares structural equation modeling (PLS-SEM). Results indicate a stable upward trend in vegetation growth with minimal fluctuations, with a growth rate of 0.95 × 10−3/a (p < 0.01). PCC analysis shows a positive correlation between NDVI and key climatic elements in over 60% of the area. The areas with significant vegetation growth had the highest average PCC. PCA and PLS-SEM identify temperature and precipitation as primary growth drivers, while elevation and land use intensity hinder growth. The MJB landscape reveals thresholds and tipping points, with specific temperature and precipitation benchmarks varying by elevation, delineating the boundary between flourishing vegetation and growth inhibition. Full article
Show Figures

Figure 1

Figure 1
<p>Location and elevation of MJB.</p>
Full article ">Figure 2
<p>Spatial distribution of datasets: (<b>a</b>). The land use types in 2020; (<b>b</b>). The average annual precipitation from 2000 to 2020 (mm); (<b>c</b>). The average annual GDP from 2000 to 2020 (10,000 yuan/km<sup>2</sup>); (<b>d</b>). The average annual NDVI from 2000 to 2020; (<b>e</b>). The average annual temperature from 2000 to 2020 (°C); (<b>f</b>). The average annual POP from 2000 to 2020 (person/km<sup>2</sup>).</p>
Full article ">Figure 3
<p>The flowchart of the study.</p>
Full article ">Figure 4
<p>Spatial distribution patterns of LUI: (<b>a</b>) 2000; (<b>b</b>) 2005; (<b>c</b>) 2010; (<b>d</b>) 2015; (<b>e</b>) 2020; and (<b>f</b>) Annual average.</p>
Full article ">Figure 5
<p>(<b>a</b>). Spatial distribution of CV; (<b>b</b>). Annual NDVI changes from 2000 to 2020; (<b>c</b>). The spatial distribution of trends; (<b>d</b>). The spatial distribution of trend and Hurst exponent coupling.</p>
Full article ">Figure 6
<p>(<b>a</b>). The PCC between NDVI and precipitation; (<b>b</b>). The significance of PCC/precipitation; (<b>c</b>). The correlation between NDVI and precipitation; (<b>d</b>). The PCC between NDVI and temperature; (<b>e</b>). the significance of PCC/temperature; (<b>f</b>). The correlation between NDVI and mean temperature.</p>
Full article ">Figure 7
<p>The biplot of PCA and the percentage of variance are explained by the principal components: (<b>a</b>). The graph displays the percentage of variance explained by PC1 and PC2 in a PCA. (<b>b</b>). The <span class="html-italic">x</span>-axis represents the number of principal components, while the <span class="html-italic">y</span>-axis indicates the proportion of variance explained by each principal component.</p>
Full article ">Figure 8
<p>PLS-SEM analysis of vegetation NDVI responses to climate drivers (precipitation and temperature), topographic factors (elevation and slope), and human activities (LUI, GDP, and POP) is presented. The thickness of the lines represents the absolute value of the path coefficients: thicker lines indicate larger absolute values.</p>
Full article ">Figure 9
<p>The interrelation between NDVI/PCC and climatic factors across different elevations.</p>
Full article ">
17 pages, 5892 KiB  
Article
Improved A* Algorithm for Mobile Robots under Rough Terrain Based on Ground Trafficability Model and Ground Ruggedness Model
by Zhiguang Liu, Song Guo, Fei Yu, Jianhong Hao and Peng Zhang
Sensors 2024, 24(15), 4884; https://doi.org/10.3390/s24154884 - 27 Jul 2024
Cited by 1 | Viewed by 912
Abstract
Considering that the existing path planning algorithms for mobile robots under rugged terrain do not consider the ground flatness and the lack of optimality, which leads to the instability of the center of mass of the mobile robot, this paper proposes an improved [...] Read more.
Considering that the existing path planning algorithms for mobile robots under rugged terrain do not consider the ground flatness and the lack of optimality, which leads to the instability of the center of mass of the mobile robot, this paper proposes an improved A* algorithm for mobile robots under rugged terrain based on the ground accessibility model and the ground ruggedness model. Firstly, the ground accessibility and ruggedness models are established based on the elevation map, expressing the ground flatness. Secondly, the elevation cost function that can obtain the optimal path is designed based on the two types of models combined with the characteristics of the A* algorithm, and the continuous cost function is established by connecting with the original distance cost function, which avoids the center-of-mass instability caused by the non-optimal path. Finally, the effectiveness of the improved algorithm is verified by simulation and experiment. The results show that compared with the existing commonly used path planning algorithms under rugged terrain, the enhanced algorithm improves the smoothness of paths and the optimization degree of paths in the path planning process under rough terrain. Full article
(This article belongs to the Topic Advances in Mobile Robotics Navigation, 2nd Volume)
Show Figures

Figure 1

Figure 1
<p>The autonomous navigation system framework. The mobile robot establishes an elevation map and a target node based on the goal pose and point cloud information, after which a new path is generated by the path planner module and tracked by the controller until the mobile robot reaches the node where the goal pose is located.</p>
Full article ">Figure 2
<p>Ground trafficability model. The elevation difference model represents the difference between the child and parent nodes, where (<span class="html-italic">x</span><sub>1</sub>, <span class="html-italic">y</span><sub>1</sub>, <span class="html-italic">z</span><sub>1</sub>) and (<span class="html-italic">x</span><sub>2</sub>, <span class="html-italic">y</span><sub>2</sub>, <span class="html-italic">z</span><sub>2</sub>) represent the parent and child node coordinates, respectively. The slope model represents the slope between the child node and the parent node, and (<span class="html-italic">x</span><sub>1</sub>, <span class="html-italic">y</span><sub>1</sub>, <span class="html-italic">z</span><sub>1</sub>) and (<span class="html-italic">x</span><sub>3</sub>, <span class="html-italic">y</span><sub>3</sub>, <span class="html-italic">z</span><sub>3</sub>) represent the parent and child node coordinates, respectively.</p>
Full article ">Figure 3
<p>Ground ruggedness model. The ground ruggedness model describes the degree of fluctuation of the search ground. The child node is (<span class="html-italic">x<sub>i</sub></span>, <span class="html-italic">y<sub>j</sub></span>, <span class="html-italic">z<sub>i</sub></span><sub>,<span class="html-italic">j</span></sub>).</p>
Full article ">Figure 4
<p>Ackermann steering mobile robot platforms. The robot is built with 3D LiDAR for sensing and is processed and controlled by Nvidia Jetson TX2 running Ubuntu 18.04 and ROS 1.12.17.</p>
Full article ">Figure 5
<p>MATLAB rough terrain simulation map.</p>
Full article ">Figure 6
<p>Paths formed by the mobile robot path planning in MATLAB rough terrain simulation map: (<b>a</b>) shows the path formed by the OA* algorithm; (<b>b</b>) shows the path formed by the AOEA* algorithm; and (<b>c</b>) shows the path formed by the IA* algorithm.</p>
Full article ">Figure 7
<p>Paths are formed by the mobile robot path planning in the Gazebo rough terrain simulation, in which the red path is the path created by the AOEA* algorithm, and the yellow path is the path made by the IA* algorithm. Arrows and dots represent the start and end points of paths.</p>
Full article ">Figure 8
<p>The <span class="html-italic">x</span>-, <span class="html-italic">y</span>-, and <span class="html-italic">z</span>-axis angles of each node of the mobile robot to the ground coordinate system for the paths planned by the two algorithms in the Gazebo rough terrain simulation: (<b>a</b>) shows the <span class="html-italic">x</span> angles of each node of the mobile robot to the ground coordinate system for the paths planned by the two algorithms in the Gazebo rough terrain simulation; (<b>b</b>) shows the <span class="html-italic">y</span> angles of each node of the mobile robot to the ground coordinate system for the paths planned by the two algorithms in the Gazebo rough terrain simulation; and (<b>c</b>) shows the <span class="html-italic">z</span> angles of each node of the mobile robot to the ground coordinate system for the paths planned by the two algorithms in the Gazebo rough terrain simulation.</p>
Full article ">Figure 9
<p>Paths are formed by the mobile robot path planning in the real-life environment, in which the yellow path is the path created by the AOEA* algorithm, and the blue path is the path made by the IA* algorithm. Arrows and dots represent the start and end points of paths.</p>
Full article ">Figure 10
<p>The <span class="html-italic">x</span>-, <span class="html-italic">y</span>-, and <span class="html-italic">z</span>-axis angles of each node of the mobile robot to the ground coordinate system for the paths planned by the two algorithms in the Gazebo rough terrain simulation: (<b>a</b>) shows the <span class="html-italic">x</span> angles of each node of the mobile robot to the ground coordinate system for the paths planned by the two algorithms in the Gazebo rough terrain simulation; (<b>b</b>) shows the <span class="html-italic">y</span> angles of each node of the mobile robot to the ground coordinate system for the paths planned by the two algorithms in the Gazebo rough terrain simulation; and (<b>c</b>) shows the <span class="html-italic">z</span> angles of each node of the mobile robot to the ground coordinate system for the paths planned by the two algorithms in the Gazebo rough terrain simulation.</p>
Full article ">
23 pages, 19174 KiB  
Article
Unmanned Aerial Vehicle Landing on Rugged Terrain by On-Board LIDAR–Camera Positioning System
by Cheng Zou, Yezhen Sun and Linghua Kong
Appl. Sci. 2024, 14(14), 6079; https://doi.org/10.3390/app14146079 - 12 Jul 2024
Viewed by 964
Abstract
Safely landing unmanned aerial vehicles (UAVs) in unknown environments that are denied by GPS is challenging but crucial. In most cases, traditional landing methods are not suitable, especially under complex terrain conditions with insufficient map information. This report proposes an innovative multi-stage UAV [...] Read more.
Safely landing unmanned aerial vehicles (UAVs) in unknown environments that are denied by GPS is challenging but crucial. In most cases, traditional landing methods are not suitable, especially under complex terrain conditions with insufficient map information. This report proposes an innovative multi-stage UAV landing framework involving (i) point cloud and image fusion positioning, (ii) terrain analysis, and (iii) neural network semantic recognition to optimize landing site selection. In the first step, 3D point cloud and image data are fused to attain a comprehensive perception of the environment. In the second step, an energy cost function considering texture and flatness is employed to identify potential landing sites based on energy scores. To navigate the complexities of classification for precise landings, the results are stratified by the difficulty of various UAV landing scenarios. In the third step, a network model is applied to analyze UAV landing site options by integrating the ResNet50 network with a convolutional block attention module. Experimental results indicate a reduction in computational load and improved landing site identification accuracy. The developed framework fuses multi-modal data to enhance the safety and feasibility of UAV landings in complex environments. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Landing framework involving a LIDAR–camera positioning system.</p>
Full article ">Figure 2
<p>Fusion of LIDAR point cloud and image.</p>
Full article ">Figure 3
<p>LIDAR and camera positioning framework.</p>
Full article ">Figure 4
<p>Example showing that the local smoothness of LIDAR points is estimated on the image plane. Left: the original image is captured from the UAV, and the SED image is applied to evaluate the edge energy. Two LIDAR points (A and B) are selected to show the computation of <math display="inline"><semantics> <msub> <mi>f</mi> <mi>c</mi> </msub> </semantics></math>. Right: Plane and edge points are evaluated in terms of their local smoothness by summing the SED energy along the minimum route to their neighbors.</p>
Full article ">Figure 5
<p>The <math display="inline"><semantics> <msub> <mi>D</mi> <mrow> <mi>S</mi> <mi>c</mi> <mi>e</mi> <mi>n</mi> <mi>e</mi> </mrow> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>D</mi> <mrow> <mi>L</mi> <mi>a</mi> <mi>n</mi> <mi>d</mi> </mrow> </msub> </semantics></math> datasets were derived from the UC Merced Land Use, AID, and Places365 datasets.</p>
Full article ">Figure 6
<p>ResNet50 neural network residual function with convolutional block attention module network architecture, (<b>a</b>) Internal structure diagram of CAM module, (<b>b</b>) Internal structure diagram of SAM module.</p>
Full article ">Figure 7
<p>Schematic diagram of the UAV experimental platform.</p>
Full article ">Figure 8
<p>Simulation of UAVs and constructed scenes. (<b>a</b>) System configuration of the UAV in the simulation including Livox Horizon and RealSense with scanning ranges parallel to the ground; the images in (<b>b</b>–<b>e</b>) are virtual environments based on four common scenarios.</p>
Full article ">Figure 9
<p>Running trajectories of the UAV: (<b>a</b>) terrain scanning trajectory of the UAV; (<b>b</b>) landing designated trajectory of the UAV.</p>
Full article ">Figure 10
<p>Comparative analysis of the absolute error of the UAV trajectory using four positioning methods during algorithm execution: (<b>a</b>) comparing each positioning technology’s trajectory with the expected value; (<b>b</b>–<b>e</b>) evaluating the mean absolute error of the UAV trajectory for each method over 10 simulations.</p>
Full article ">Figure 11
<p>Based on the image and point cloud information acquired during the experiment and the confirmation of the landing site, (<b>i</b>–<b>iii</b>) illustrate the steps of the experimental process, and (<b>a</b>–<b>e</b>) show the types of experimental sites selected for the outdoor tests (part of the point cloud map is rotated). The blue points marked in the lower column show the approximate landing site selection, and the red points show the final landing site.</p>
Full article ">Figure 12
<p>Comparison of UAV landing positions with different methods.</p>
Full article ">
16 pages, 7662 KiB  
Article
Exploring the Influence of Terrain Blockage on Spatiotemporal Variations in Land Surface Temperature from the Perspective of Heat Energy Redistribution
by Hong Gao, Yong Dong, Liang Zhou and Xi Wang
ISPRS Int. J. Geo-Inf. 2024, 13(6), 200; https://doi.org/10.3390/ijgi13060200 - 14 Jun 2024
Viewed by 828
Abstract
Land surface temperature (LST) is a critical indicator of the earth’s surface environment, which has significant implications for research on the ecological environment and climate change. The influence of terrain on LST is complex due to its rugged and varied surface topography. The [...] Read more.
Land surface temperature (LST) is a critical indicator of the earth’s surface environment, which has significant implications for research on the ecological environment and climate change. The influence of terrain on LST is complex due to its rugged and varied surface topography. The relationship between traditional terrain features and LST has been comprehensively discussed in the literature; however, terrain blockage has received less attention and could influence LST by hindering the redistribution of heat energy in mountain regions. Here, we investigate the influence of terrain blockage on the spatiotemporal variation in LST in mountain regions. We first propose a terrain feature framework to characterize the effect of terrain blockage from the perspective of heat energy redistribution and then adopt a random forest model to analyze the relationship between terrain blockage features and LST over a whole year. The results show that terrain blockage significantly influences the spatial heterogeneity of LST, which can be effectively simulated based on terrain blockage features, with a mean deviation of less than 0.15 K. Terrain blockage has a more pronounced influence on LST during the four months from June to September. This influence is also more evident during nighttime than daytime. Regarding LST in mountain regions, local terrain blockage features have a greater influence than global terrain blockage features. In spatial terms, the influence of terrain blockage on LST is uniform. Moreover, the diurnal variation in LST can also be effectively simulated based on terrain blockage. The contribution of this study lies in the finding that terrain blockage can influence the spatiotemporal variation in LST through the process of heat energy redistribution. The terrain blockage features proposed in this study may be useful for other studies of the ecological environment in mountain regions. Full article
Show Figures

Figure 1

Figure 1
<p>Location of study area.</p>
Full article ">Figure 2
<p>The flow diagram of this study.</p>
Full article ">Figure 3
<p>Illustration for extracting section elevation series in four directions. (<b>a</b>) Determination for the four main directions according to LST data. (<b>b</b>) Results of section elevation series for the four main directions in the spatial grid named G1000.</p>
Full article ">Figure 4
<p>Coefficients of correlation between different TBFs and LST.</p>
Full article ">Figure 5
<p>Evaluation indices for LST simulation in different months.</p>
Full article ">Figure 6
<p>Importance of different types of TBFs.</p>
Full article ">Figure 7
<p>Importance of different directions of TBFs.</p>
Full article ">Figure 8
<p>Scatter plot of simulation for single days and nights (the dotted line is the 1:1 line, where the simulated value equals the original value).</p>
Full article ">Figure 9
<p>Spatial distribution of the origin and simulation and the associated errors.</p>
Full article ">Figure 10
<p>Simulation accuracy for LST diurnal variation.</p>
Full article ">
30 pages, 6730 KiB  
Article
An Enhanced Multiple Unmanned Aerial Vehicle Swarm Formation Control Using a Novel Fractional Swarming Strategy Approach
by Abdul Wadood, Al-Fahad Yousaf and Aadel Mohammed Alatwi
Fractal Fract. 2024, 8(6), 334; https://doi.org/10.3390/fractalfract8060334 - 3 Jun 2024
Cited by 1 | Viewed by 879
Abstract
This paper addresses the enhancement of multiple Unmanned Aerial Vehicle (UAV) swarm formation control in challenging terrains through the novel fractional memetic computing approach known as fractional-order velocity-pausing particle swarm optimization (FO-VPPSO). Existing particle swarm optimization (PSO) algorithms often suffer from premature convergence [...] Read more.
This paper addresses the enhancement of multiple Unmanned Aerial Vehicle (UAV) swarm formation control in challenging terrains through the novel fractional memetic computing approach known as fractional-order velocity-pausing particle swarm optimization (FO-VPPSO). Existing particle swarm optimization (PSO) algorithms often suffer from premature convergence and an imbalanced exploration–exploitation trade-off, which limits their effectiveness in complex optimization problems such as UAV swarm control in rugged terrains. To overcome these limitations, FO-VPPSO introduces an adaptive fractional order β and a velocity pausing mechanism, which collectively enhance the algorithm’s adaptability and robustness. This study leverages the advantages of a meta-heuristic computing approach; specifically, fractional-order velocity-pausing particle swarm optimization is utilized to optimize the flying path length, mitigate the mountain terrain costs, and prevent collisions within the UAV swarm. Leveraging fractional-order dynamics, the proposed hybrid algorithm exhibits accelerated convergence rates and improved solution optimality compared to traditional PSO methods. The methodology involves integrating terrain considerations and diverse UAV control parameters. Simulations under varying conditions, including complex terrains and dynamic threats, substantiate the effectiveness of the approach, resulting in superior fitness functions for multi-UAV swarms. To validate the performance and efficiency of the proposed optimizer, it was also applied to 13 benchmark functions, including uni- and multimodal functions in terms of the mean average fitness value over 100 independent trials, and furthermore, an improvement at percentages of 29.05% and 2.26% is also obtained against PSO and VPPSO in the case of the minimum flight length, as well as 16.46% and 1.60% in mountain terrain costs and 55.88% and 31.63% in collision avoidance. This study contributes valuable insights to the optimization challenges in UAV swarm-formation control, particularly in demanding terrains. The FO-VPPSO algorithm showcases potential advancements in swarm intelligence for real-world applications. Full article
Show Figures

Figure 1

Figure 1
<p>Mission area comprised of mountains.</p>
Full article ">Figure 2
<p>Mountain terrain cost.</p>
Full article ">Figure 3
<p>Leader–follower formation mode.</p>
Full article ">Figure 4
<p>Flowchart of FO_VPPSO.</p>
Full article ">Figure 5
<p>Adaptive beta factor plot.</p>
Full article ">Figure 6
<p>Graphical demonstration of process.</p>
Full article ">Figure 7
<p>XYZ-axis plot of trajectories.</p>
Full article ">Figure 8
<p>XZ-axes of trajectories.</p>
Full article ">Figure 9
<p>YZ-axes plot of trajectories.</p>
Full article ">Figure 10
<p>Fitness graph of all UAVs (FOVPPSO).</p>
Full article ">Figure 11
<p>Fitness graph of all UAVs (VPPSO).</p>
Full article ">Figure 12
<p>Dynamic obstacle movements (trajectories).</p>
Full article ">Figure 13
<p>Path deviation concept.</p>
Full article ">Figure 14
<p>Path deviation of swarm.</p>
Full article ">Figure 15
<p>Path deviation for leader UAV (dynamic obstacle).</p>
Full article ">Figure 16
<p>Convergence behaviors of all UAVS (FOVPPSO).</p>
Full article ">Figure 17
<p>Convergence behaviors of all UAVs (VPPSO).</p>
Full article ">
15 pages, 9622 KiB  
Technical Note
Estimation of Antarctic Ice Sheet Thickness Based on 3D Density Interface Inversion Considering Terrain and Undulating Observation Surface Simultaneously
by Yandong Liu, Jun Wang, Fang Li and Xiaohong Meng
Remote Sens. 2024, 16(11), 1905; https://doi.org/10.3390/rs16111905 - 25 May 2024
Viewed by 873
Abstract
The thickness of the Antarctic ice sheet is a crucial parameter for inferring glacier mass and its evolution process. In the literature, the gravity method has been proven to be one of the effective means for estimating ice sheet thickness. And it is [...] Read more.
The thickness of the Antarctic ice sheet is a crucial parameter for inferring glacier mass and its evolution process. In the literature, the gravity method has been proven to be one of the effective means for estimating ice sheet thickness. And it is a preferred approach when direct measurements are not available. However, few gravity inversion methods are valid in rugged terrain areas with undulating observation surfaces (UOSs). To solve this problem, this paper proposes an improved high-precision 3D density interface inversion method considering terrain and UOSs simultaneously. The proposed method utilizes airborne gravity data at their flight altitudes, instead of the continued data yield from the unstable downward continuation procedure. In addition, based on the undulating right rectangular prism model, the large reliefs of the terrain are included in the iterative inversion. The proposed method is verified on two synthetic examples and is successfully applied to real data in East Antarctica. Full article
Show Figures

Figure 1

Figure 1
<p>Interpretation model of an ice sheet, considering terrain and UOSs. The UOSs are drawn with green lines, and the solid black dots represent the observed points. The ice sheet is divided into a series of juxtaposed vertical prisms in blue, and the parameter to be calculated is the elevation of the prisms’ bottom.</p>
Full article ">Figure 2
<p>A flow chart of the interface inversion considering terrain and UOSs simultaneously.</p>
Full article ">Figure 3
<p>The models and gravity anomaly in synthetic example 1: (<b>a</b>) The elevation of observation surface. (<b>b</b>) The elevation of terrain. (<b>c</b>) The elevation of ice–rock interface. (<b>d</b>) The ice sheet thickness. (<b>e</b>) The theoretical gravity anomaly due to these models.</p>
Full article ">Figure 4
<p>The results of the simple models obtained by the conventional method: (<b>a</b>) The calculated gravity anomaly. (<b>b</b>) The residuals of data fitting. (<b>c</b>) The estimated ice thickness. (<b>d</b>) The residuals between the inverted thickness and the true value.</p>
Full article ">Figure 5
<p>The results of the simple models obtained by the proposed method: (<b>a</b>) The calculated gravity anomaly. (<b>b</b>) The residuals of data fitting. (<b>c</b>) The estimated ice thickness. (<b>d</b>) The residuals between the inverted thickness and the true value.</p>
Full article ">Figure 6
<p>The models and gravity anomaly in synthetic example 2: (<b>a</b>) The elevation of observation surface. (<b>b</b>) The elevation of terrain. (<b>c</b>) The elevation of ice–rock interface. (<b>d</b>) The ice sheet thickness. (<b>e</b>) The theoretical gravity anomaly due to these models.</p>
Full article ">Figure 7
<p>The results of the complex models obtained by the conventional method: (<b>a</b>) The calculated gravity anomaly. (<b>b</b>) The residuals of data fitting. (<b>c</b>) The estimated ice thickness. (<b>d</b>) The residuals between the inverted ice thickness and the true value.</p>
Full article ">Figure 8
<p>The results of the complex models obtained by the proposed method: (<b>a</b>) The calculated gravity anomaly. (<b>b</b>) The residuals of data fitting. (<b>c</b>) The estimated ice thickness. (<b>d</b>) The residuals between the inverted thickness and the true value.</p>
Full article ">Figure 9
<p>Location of the study area in Antarctica. The base map is the elevation of the Antarctic topography. The study area is bounded by a green rectangle.</p>
Full article ">Figure 10
<p>The airborne geophysical data from the AGAP Project. Aerogravity data: (<b>a</b>) the residual gravity anomaly caused by the ice–rock interface; (<b>b</b>) the aircraft altitude. Radio-echo sounding data: (<b>c</b>) ice surface elevation, (<b>d</b>) ice bed elevation, and (<b>e</b>) ice thickness.</p>
Full article ">Figure 11
<p>The results obtained by the conventional method in the survey area: (<b>a</b>) The calculated gravity anomaly. (<b>b</b>) The residuals of data fitting. (<b>c</b>) The estimated ice thickness. (<b>d</b>) The residuals between the estimated ice thickness and the radar-derived ice thickness.</p>
Full article ">Figure 12
<p>The results obtained by the proposed method in the survey area: (<b>a</b>) The calculated gravity anomaly. (<b>b</b>) The residuals of data fitting. (<b>c</b>) The estimated ice thickness. (<b>d</b>) The residuals between the inverted thickness and the radar-derived value.</p>
Full article ">
Back to TopTop