[go: up one dir, main page]

 
 
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

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
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

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
remove_circle_outline

Search Results (3,463)

Search Parameters:
Keywords = DEM

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
22 pages, 11550 KiB  
Article
Numerical Study of Cone Penetration Tests in Lunar Regolith for Strength Index
by Xueliang Zhao, Zixiong Liu, Yu Li, Hao Wang and Zhaodong Xu
Appl. Sci. 2024, 14(22), 10645; https://doi.org/10.3390/app142210645 (registering DOI) - 18 Nov 2024
Abstract
The cohesive properties of lunar regolith, combined with a low-gravity environment, result in it having a distinct mechanical behavior from sandy soil on Earth. Consequently, empirical formulas derived from cone penetration tests (CPTs) for calculating the shear strength parameters of Earth’s sand cannot [...] Read more.
The cohesive properties of lunar regolith, combined with a low-gravity environment, result in it having a distinct mechanical behavior from sandy soil on Earth. Consequently, empirical formulas derived from cone penetration tests (CPTs) for calculating the shear strength parameters of Earth’s sand cannot be directly applied to lunar regolith. This study utilized the three-dimensional discrete element method (DEM) to numerically simulate triaxial shear tests and cone penetration tests in a lunar environment. The particle contact model for lunar regolith in the discrete element method (DEM) simulation incorporated the hysteresis effect of van der Waals forces, thereby simulating the cohesive properties of lunar regolith in a lunar environment. We proposed a relationship for calculating the shear strength index of lunar regolith based on normalized cone tip resistance using the results from triaxial and CPT simulations and referencing empirical formulas derived from ground-based CPT data. The results of this study provide a valuable reference for future lunar CPTs. Full article
(This article belongs to the Section Civil Engineering)
Show Figures

Figure 1

Figure 1
<p>A contact mechanics model of lunar regolith particles.</p>
Full article ">Figure 2
<p>Contact model mesoscopic parameter analysis results.</p>
Full article ">Figure 3
<p>Particle gradation curve.</p>
Full article ">Figure 4
<p>A comparison of stress–strain curves after triaxial simulation calibration and triaxial test results.</p>
Full article ">Figure 5
<p>Simulated foundation diagram.</p>
Full article ">Figure 6
<p>Schematic diagram of particle size distribution of each section.</p>
Full article ">Figure 7
<p>Probe size diagram.</p>
Full article ">Figure 8
<p>A three-dimensional simulation model of the static penetration test.</p>
Full article ">Figure 9
<p>A triaxial shear simulation model.</p>
Full article ">Figure 10
<p>Stress–strain curves from triaxial simulations with different relative densities.</p>
Full article ">Figure 11
<p>The peak internal friction angle <span class="html-italic">φ</span> and cohesion <span class="html-italic">c</span>.</p>
Full article ">Figure 12
<p>Cone tip resistance–penetration depth curves for different relative densities.</p>
Full article ">Figure 13
<p>Activation of cone tip resistance in the lunar surface–depth curve.</p>
Full article ">Figure 14
<p>A comparison of normalized cone tip resistance–depth curves in a lunar environment within the range of <span class="html-italic">Z</span> = 4–8.</p>
Full article ">Figure 15
<p>The fitting formula developed in this study compared with other empirical formulas, Jamiolkowski et al. [<a href="#B13-applsci-14-10645" class="html-bibr">13</a>].</p>
Full article ">Figure 16
<p>A fitted curve and empirical formula [<a href="#B25-applsci-14-10645" class="html-bibr">25</a>] comparison.</p>
Full article ">Figure 17
<p>Direct fitting with indirect comparison.</p>
Full article ">Figure 18
<p>A comparison of fitting formulas with the proposed empirical formula, Marco et al. [<a href="#B10-applsci-14-10645" class="html-bibr">10</a>].</p>
Full article ">Figure 19
<p>The relationship of relative density with adhesion.</p>
Full article ">Figure 20
<p>The relationship between cohesion <span class="html-italic">c</span> and stable normalized cone tip resistance.</p>
Full article ">
16 pages, 5187 KiB  
Article
Effect of High-Stress Levels on the Shear Behavior of Geosynthetic-Reinforced Marine Coral Sands
by Lixing Liu, Zhixiong Chen, Xuanming Ding and Qiang Ou
J. Mar. Sci. Eng. 2024, 12(11), 2081; https://doi.org/10.3390/jmse12112081 (registering DOI) - 18 Nov 2024
Abstract
As an important construction material, the mechanical and deformation properties of marine coral sand determine the safety and stability of related island and coastal engineering construction. The porous and easily broken characteristics of coral sand often make it difficult to meet engineering construction [...] Read more.
As an important construction material, the mechanical and deformation properties of marine coral sand determine the safety and stability of related island and coastal engineering construction. The porous and easily broken characteristics of coral sand often make it difficult to meet engineering construction needs. In particular, coral sand undergoes a large amount of particle breakage under high-stress conditions, which in turn negatively affects its mechanical and deformation properties. In this study, the macro- and micro-mechanical behavior of geosynthetic-reinforced coral sand under high confining pressure was investigated and compared with unreinforced cases using the three-dimensional discrete element method (DEM), which was verified by indoor triaxial tests. The results showed that the stress–strain responses of unreinforced and reinforced coral sand under high confining pressure showed completely different trends, i.e., the hardening tendency shown in the reinforced case. Geosynthetic reinforcement can significantly inhibit the stress–strain softening and bulging deformation of coral sand under high confining pressure, thus improving the shear mechanical performance of the reinforced sample. At the microscopic scale, high confining pressure and reinforcement affected the contact force distribution pattern and stress level between particles, determining the macroscopic mechanical and deformation performance. In addition, the breakage of particles under high confining pressure was mainly affected by shear strain and reinforcement. The particle fragment distribution, particle gradation, and relative breakage index exhibited different trends at different confining pressure levels. These breakage characteristics were closely related to the deformation and stress levels of unreinforced and reinforced samples. Full article
(This article belongs to the Section Ocean Engineering)
Show Figures

Figure 1

Figure 1
<p>Shape characteristics of coral sand particles: (<b>a</b>) coral sand from the South China Sea, and (<b>b</b>) angularity and porosity characteristics.</p>
Full article ">Figure 2
<p>3D coral sand particle construction and microscopic contact models based on DEM: (<b>a</b>) 3D DEM particle, (<b>b</b>) linear contact model, and (<b>c</b>) linear parallel bonding model.</p>
Full article ">Figure 3
<p>Specimen preparation, load application, and failure modes in triaxial testing and DEM simulation.</p>
Full article ">Figure 4
<p>Effect of high confining pressure on failure patterns: (<b>a</b>) unreinforced case and (<b>b</b>) reinforcement.</p>
Full article ">Figure 5
<p>Effect of high confining pressure on deviatoric stress–axial strain curves: (<b>a</b>) unreinforced case and (<b>b</b>) reinforcement.</p>
Full article ">Figure 6
<p>Effect of high confining pressure on the shear strength: (<b>a</b>) unreinforced case and (<b>b</b>) reinforcement.</p>
Full article ">Figure 7
<p>Effect of high confining pressure on the shear mechanical parameters compared with the low pressure condition [<a href="#B12-jmse-12-02081" class="html-bibr">12</a>]: (<b>a</b>) internal friction angle and (<b>b</b>) pseudo-cohesion.</p>
Full article ">Figure 8
<p>Effect of high confining pressure on microscopic contact force: (<b>a</b>) unreinforced case and (<b>b</b>) reinforcement.</p>
Full article ">Figure 9
<p>Tensile force of geogrids under varied confining pressures and shear strains: (<b>a</b>) 6% strain, (<b>b</b>) 14% strain, and (<b>c</b>) 20% strain.</p>
Full article ">Figure 10
<p>3D distribution of particle fragments under 1600 kPa confining pressure and different axial strains: (<b>a</b>) unreinforced case and (<b>b</b>) reinforcement.</p>
Full article ">Figure 11
<p>3D distribution of particle fragments under 20% shear strain and different confining pressures: (<b>a</b>) unreinforced case and (<b>b</b>) reinforcement.</p>
Full article ">Figure 12
<p>Effects of high confining pressure and reinforcement on the particle size distribution: (<b>a</b>) unreinforced case, <b>c</b> = 800 kPa, (<b>b</b>) unreinforced case, <b>c</b> = 1600 kPa, (<b>c</b>) reinforced case, <b>c</b> = 800 kPa, and (<b>d</b>) reinforced case, <b>c</b> = 1600 kPa.</p>
Full article ">Figure 13
<p>Evolution of the relative breakage index <span class="html-italic">B<sub>r</sub></span>: (<b>a</b>,<b>b</b>) evolution of the relative breakage index with axial strain in unreinforced and reinforced cases, and (<b>c</b>,<b>d</b>) final value of the relative breakage index in unreinforced and reinforced cases.</p>
Full article ">
21 pages, 8292 KiB  
Article
Modeling and Dynamic Characteristics of Tracked Vehicle Equipped with Symmetrical Suspensions Based on Multi-Body Dynamics and Discrete Element Coupling Method
by Jing Tao, Zhiyuan Deng, Xiuquan Cao, Guangzhong Hu and Ping Wang
Appl. Sci. 2024, 14(22), 10618; https://doi.org/10.3390/app142210618 - 18 Nov 2024
Viewed by 118
Abstract
For improving the adaptability of a tracked vehicle equipped with a symmetrical suspension system on complex hilly terrain, based on the coupling method of multi-body dynamics (MBD) and discrete element method (DEM), an MBD-DEM coupling model was built and verified to explore its [...] Read more.
For improving the adaptability of a tracked vehicle equipped with a symmetrical suspension system on complex hilly terrain, based on the coupling method of multi-body dynamics (MBD) and discrete element method (DEM), an MBD-DEM coupling model was built and verified to explore its dynamic behaviors on soil. Firstly, according to the basic parameters of the tracked vehicle equipped with a symmetrical suspension system, a corresponding MBD model was built in Recurdyn V9R4 software. Based on the Euler–Lagrange method, the mathematical structure of the symmetrical suspension system was analyzed, and a corresponding mathematical simulation model was built in Matlab2016/Simulink to verify the MBD model. Secondly, based on the DEM theory and the parameter of the soil in a hilly area, a granular pavement model was built. Then, based on the coupling method of MBD and DEM, the corresponding MBD-DEM coupling model was built. Finally, using the MBD-DEM coupling model, the dynamic behaviors of the tracked vehicle equipped with a symmetrical suspension system under horizontal condition, the climbing condition and the obstacle crossing condition were obtained and discussed. The study results show that the proposed MBD-DEM coupling model could be used effectively to analyze the dynamic characteristics of the tracked vehicle. In addition, according to the analysis of the dynamic characteristics of the proposed tracked vehicle, the tracked vehicle equipped with a symmetrical suspension presents good adaptabilities under various working conditions. Full article
Show Figures

Figure 1

Figure 1
<p>Diagram of tracked vehicle.</p>
Full article ">Figure 2
<p>Diagram of track subsystem.</p>
Full article ">Figure 3
<p>Track cross-section.</p>
Full article ">Figure 4
<p>Schematic diagram of symmetrical suspension.</p>
Full article ">Figure 5
<p>MBD model of tracked vehicle.</p>
Full article ">Figure 6
<p>An MBD model of the tracked vehicle and the road surface.</p>
Full article ">Figure 7
<p>Simplification of symmetrical suspension.</p>
Full article ">Figure 8
<p>Road excitation model in Matlab software.</p>
Full article ">Figure 9
<p>The vertical accelerations of the jointed point C obtained by the MBD model and the mathematical model.</p>
Full article ">Figure 10
<p>Schematic of contact between track and particles.</p>
Full article ">Figure 11
<p>MBD-DEM coupling model.</p>
Full article ">Figure 12
<p>The vertical acceleration of the vehicle body at the mass center and the symmetrical suspensions at the hinge point under various velocities.</p>
Full article ">Figure 13
<p>Acceleration amplitude evolution with various velocities.</p>
Full article ">Figure 14
<p>The vertical acceleration of the vehicle body at the mass center and the symmetrical suspensions at the hinge point under the climbing conditions.</p>
Full article ">Figure 15
<p>The pitch angles of the tracked vehicle under the climbing conditions.</p>
Full article ">Figure 16
<p>The maximum pitch angles of the tracked vehicle evolution with the climbing slope.</p>
Full article ">Figure 17
<p>The vertical acceleration of the vehicle body at the mass center and the symmetrical suspensions at the hinge point under the obstacle crossing conditions.</p>
Full article ">Figure 18
<p>The dynamic deflection of the symmetrical suspension under the obstacle crossing condition.</p>
Full article ">Figure 19
<p>The suspension deflection of valley evolution with the height of the obstacle.</p>
Full article ">Figure 20
<p>Vertical displacement of vehicle at mass center.</p>
Full article ">
22 pages, 24817 KiB  
Article
Construction of Mining Subsidence Basin and Inversion of Predicted Subsidence Parameters Based on UAV Photogrammetry Products Considering Horizontal Displacement
by Jinqi Zhao, Yufen Niu, Zhengpei Zhou, Zhong Lu, Zhimou Wang, Zhaojiang Zhang, Yiyao Li and Ziheng Ju
Remote Sens. 2024, 16(22), 4283; https://doi.org/10.3390/rs16224283 (registering DOI) - 17 Nov 2024
Viewed by 195
Abstract
Constructing high-precision subsidence basins is of paramount importance for mining subsidence monitoring. Traditional unmanned aerial vehicle (UAV) photogrammetry techniques typically construct subsidence basins by directly differencing digital elevation models (DEMs) from different monitoring periods. However, this method often neglects the influence of horizontal [...] Read more.
Constructing high-precision subsidence basins is of paramount importance for mining subsidence monitoring. Traditional unmanned aerial vehicle (UAV) photogrammetry techniques typically construct subsidence basins by directly differencing digital elevation models (DEMs) from different monitoring periods. However, this method often neglects the influence of horizontal displacement on the accuracy of the subsidence basin. Taking a mining area in Ordos, Inner Mongolia, as an example, this study employed the normalized cross-correlation (NCC) matching algorithm to extract horizontal displacement information between two epochs of a digital orthophoto map (DOM) and subsequently corrected the horizontal position of the second-epoch DEM. This ensured that the planar positions of ground feature points remained consistent in the DEM before and after subsidence. Based on this, the vertical displacement in the subsidence area (the subsidence basin) was obtained via DEM differencing, and the parameters of the post-correction subsidence basin were inverted using the probability integral method (PIM). The experimental results indicate that (1) the horizontal displacement was influenced by the gully topography, causing the displacement within the working face to be segmented on both sides of the gully; (2) the influence of the terrain on the subsidence basin was significantly reduced after correction; (3) the post-correction surface subsidence curve was smoother than the pre-correction curve, with abrupt error effects markedly diminished; (4) the accuracy of the post-correction subsidence basin increased by 43.12% compared with the total station data; and (5) comparing the measured horizontal displacement curve with that derived using the probability integral method revealed that the horizontal displacement on the side of an old goaf adjacent to the newly excavated working face shifted toward the advancing direction of the new working face as mining progressed. This study provides a novel approach and insights for using low-cost UAVs to construct high-precision subsidence basins. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of the study area location. (<b>a</b>) Map of China; (<b>b</b>) DEM of Ordos; (<b>c</b>) study area.</p>
Full article ">Figure 2
<p>Technical flow chart of this research.</p>
Full article ">Figure 3
<p>Schematic diagram of the DEM correction process.</p>
Full article ">Figure 4
<p>(<b>a</b>) East–west displacement; (<b>b</b>) north–south displacement.</p>
Full article ">Figure 5
<p>Illustration of the relationship between horizontal displacement and topography. (<b>a</b>,<b>b</b>) are cross-sectional views of profile A-A′; (<b>c</b>,<b>d</b>) are cross-sectional views of profile B-B′.</p>
Full article ">Figure 6
<p>Horizontal displacement in gully topography. (<b>a</b>) A-A′ cross-section; (<b>b</b>) local displacement field.</p>
Full article ">Figure 7
<p>Subsidence basin. (<b>a</b>) Pre-correction subsidence basin; (<b>b</b>) post-correction subsidence basin.</p>
Full article ">Figure 8
<p>Local maps of areas I and II. (<b>a</b>) Magnified view of area I pre-correction; (<b>b</b>) magnified view of area I post-correction; (<b>c</b>) magnified view of area II pre-correction; (<b>d</b>) magnified view of area II post-correction; (<b>e</b>) 1-1′ cross-section; (<b>f</b>) 2-2′ cross-section.</p>
Full article ">Figure 9
<p>Subsidence curves of pre-correction and post-correction. (<b>a</b>) A-A′ cross-section; (<b>b</b>) C-C′ cross-section.</p>
Full article ">Figure 10
<p>Inverted subsidence basin.</p>
Full article ">Figure 11
<p>Measured subsidence basin.</p>
Full article ">Figure 12
<p>(<b>a</b>) Strike main profile; (<b>b</b>) dip main profile.</p>
Full article ">Figure 13
<p>Horizontal displacement of strike main profile. (<b>a</b>) Strike main profile; (<b>b</b>) partial enlarged detail.</p>
Full article ">Figure 14
<p>Horizontal displacement of dip main profile. (<b>a</b>) Dip main profile; (<b>b</b>) partial enlarged detail.</p>
Full article ">Figure 15
<p>Horizontal displacement error.</p>
Full article ">Figure 16
<p>Statistical chart of residuals for subsidence basin.</p>
Full article ">Figure 17
<p>Statistical chart of strike residuals.</p>
Full article ">Figure 18
<p>Statistical chart of dip residuals.</p>
Full article ">Figure 19
<p>Statistical analysis of errors in subsidence basin.</p>
Full article ">
22 pages, 42906 KiB  
Article
Monitoring the Subsidence in Wan’an Town of Deyang Based on PS-InSAR Technology (Sichuan, China)
by Hongyi Guo, Antonio Miguel Martínez-Graña and José Angel González-Delgado
Sustainability 2024, 16(22), 10010; https://doi.org/10.3390/su162210010 - 16 Nov 2024
Viewed by 526
Abstract
In recent years, land subsidence has become a crucial factor affecting urban safety and sustainable development, especially in Wan’an Town. To accurately monitor and analyze the land subsidence in Wan’an Town, this study uses the PS-InSAR technique combined with an improved DEM for [...] Read more.
In recent years, land subsidence has become a crucial factor affecting urban safety and sustainable development, especially in Wan’an Town. To accurately monitor and analyze the land subsidence in Wan’an Town, this study uses the PS-InSAR technique combined with an improved DEM for detailed research on land subsidence in Wan’an Town. PS-InSAR, or Permanent Scatterer Interferometric SAR, is suitable for high-precision monitoring of surface deformation. The natural neighbor interpolation method optimizes DEM data, improving its spatial resolution and accuracy. In this study, multiple periods of SAR imagery data of Wan’an Town were collected and preprocessed through radiometric calibration, phase unwrapping, and other steps. Using the PS-InSAR technique, the phase information of permanent scatterers (PS points) on the surface was extracted to establish a deformation model and preliminarily analyze the land subsidence in Wan’an Town. Concurrently, the DEM data were optimized using the natural neighbor interpolation method to enhance its accuracy. Finally, the optimized DEM data were combined with the surface deformation information extracted through the PS-InSAR technique for a detailed analysis of the land subsidence in Wan’an Town. The research results indicate that the DEM data optimized by the natural neighbor interpolation method have higher accuracy and spatial resolution, providing a more accurate reflection of the topographical features of Wan’an Town. The research found that the optimized DEM provided a more accurate reflection of Wan’an Town’s topographical features. By combining PS-InSAR data, subsidence information from 2016 to 2024 was calculated. The study area showed varying degrees of subsidence, with rates ranging from 6 mm/year to 10 mm/year. Four characteristic deformation areas were analyzed for causes and influencing factors. The findings contribute to understanding urban land subsidence, guiding urban planning, and providing data support for geological disaster warning and prevention. Full article
Show Figures

Figure 1

Figure 1
<p>Digital elevation model of the study area.</p>
Full article ">Figure 2
<p>Topography of the study area and radar image coverage area.</p>
Full article ">Figure 3
<p>Geology map of the study area.</p>
Full article ">Figure 4
<p>Elevation contrast chart.</p>
Full article ">Figure 5
<p>Workflow of PS processing.</p>
Full article ">Figure 6
<p>Spatial and temporal baseline distribution map.</p>
Full article ">Figure 7
<p>Differential interferogram.</p>
Full article ">Figure 8
<p>Spatial distribution of the average subsidence rate in the study area from 2016 to 2024.</p>
Full article ">Figure 9
<p>Settlement comparison diagram.</p>
Full article ">Figure 10
<p>Total subsidence in the study area from 2016 to 2024.</p>
Full article ">Figure 11
<p>Natural neighbor interpolation.</p>
Full article ">Figure 12
<p>Time-series deformation map of the study area from 2016 to 2024.</p>
Full article ">Figure 13
<p>GPS survey map. (<b>A</b>) Field survey map of target A, (<b>B</b>) field survey map of target B, (<b>C</b>) field survey map of target C, (<b>D</b>) field survey map of target D.</p>
Full article ">Figure 14
<p>(<b>A</b>) Deformation analysis of Deformation Target Area A. (<b>B</b>) Deformation analysis of Deformation Target Area B. (<b>C</b>) Deformation analysis of Deformation Target Area C. (<b>D</b>) Deformation analysis of Deformation Target Area D.</p>
Full article ">Figure 14 Cont.
<p>(<b>A</b>) Deformation analysis of Deformation Target Area A. (<b>B</b>) Deformation analysis of Deformation Target Area B. (<b>C</b>) Deformation analysis of Deformation Target Area C. (<b>D</b>) Deformation analysis of Deformation Target Area D.</p>
Full article ">
23 pages, 19709 KiB  
Article
Analysis of Gravel Migration Patterns During Vibration Rolling and Their Impact on GCL Performance Based on DEM
by Hu Huang, Feihao Chen, Qingming Qiu, Ruihang Li and Lixia Guo
Buildings 2024, 14(11), 3640; https://doi.org/10.3390/buildings14113640 (registering DOI) - 15 Nov 2024
Viewed by 248
Abstract
In this study, a multilayer composite rolling model consisting of a rolling wheel, a protective layer, a GCL, and a support layer was constructed by the discrete element method (DEM). Soil compaction and gravel migration, and their effects on the GCL, were analyzed [...] Read more.
In this study, a multilayer composite rolling model consisting of a rolling wheel, a protective layer, a GCL, and a support layer was constructed by the discrete element method (DEM). Soil compaction and gravel migration, and their effects on the GCL, were analyzed from a fine viewpoint, and three key indexes for the safety assessment of the GCL were proposed: local elongation, gravel embedment value, and bentonite allotment number. The results show that the soil porosity and cumulative settlement do not decrease all the time with the number of rolling passes, and there exists an optimal number of rolling passes during the rolling process; the protective layer of gravel soil moves more frequently than the support layer; and the nearly rectangular and nearly elliptical gravels are more likely to rotate. The maximum local elongation of the GCL was 3.79% during the lapping process, and all gravels in contact with the upper boundary of the GCL extruded the GCL to varying degrees during the lapping process. The distribution of bentonite particles is closely related to the contact mode between gravel and GCL. Full article
(This article belongs to the Section Building Structures)
Show Figures

Figure 1

Figure 1
<p>The proportions of test material particle sizes.</p>
Full article ">Figure 2
<p>Field soil grading curve.</p>
Full article ">Figure 3
<p>Triaxial compression test sample preparation process: (<b>a</b>) sample preparation steps; (<b>b</b>) sample preparation process.</p>
Full article ">Figure 4
<p>Test equipment.</p>
Full article ">Figure 5
<p>Stress-strain curves of lime-mixed gravel soil under different confining pressures in the triaxial compression test.</p>
Full article ">Figure 6
<p>Simulation method of gravel soil.</p>
Full article ">Figure 7
<p>The process of simulating gravels.</p>
Full article ">Figure 8
<p>Vibration rolling model schematic diagram with site construction drawings: (<b>a</b>) construction diagram; (<b>b</b>) modeling process; (<b>c</b>) backfill construction drawings.</p>
Full article ">Figure 8 Cont.
<p>Vibration rolling model schematic diagram with site construction drawings: (<b>a</b>) construction diagram; (<b>b</b>) modeling process; (<b>c</b>) backfill construction drawings.</p>
Full article ">Figure 9
<p>Simulation model gradation curve.</p>
Full article ">Figure 10
<p>Triaxial compression test: (<b>a</b>) numerical simulation model; (<b>b</b>) numerical simulation results and test results comparison diagram.</p>
Full article ">Figure 11
<p>Calibration process of the GCL: (<b>a</b>) GCL tensile test model; (<b>b</b>) test results.</p>
Full article ">Figure 12
<p>Porosity and cumulative settlement with increasing number of rolling passes.</p>
Full article ">Figure 13
<p>Vector displacement field of lime-mixed gravel soil and GCL under different numbers of rolling passes: (<b>a</b>) the fourth time; (<b>b</b>) the sixth time; (<b>c</b>) the ninth time; (<b>d</b>) the eleventh time; (<b>e</b>) the fifteenth time.</p>
Full article ">Figure 14
<p>Typical gravel distribution and shape. (<b>a</b>) Typical gravel distribution location map; (<b>b</b>) Typical gravel shape.</p>
Full article ">Figure 15
<p>The movements and trajectories of gravels in the PL group: (<b>a</b>) PL1-1; (<b>b</b>) PL1-2; (<b>c</b>) PL1-3; (<b>d</b>) PL2-1; (<b>e</b>) PL2-2; (<b>f</b>) PL2-3.</p>
Full article ">Figure 15 Cont.
<p>The movements and trajectories of gravels in the PL group: (<b>a</b>) PL1-1; (<b>b</b>) PL1-2; (<b>c</b>) PL1-3; (<b>d</b>) PL2-1; (<b>e</b>) PL2-2; (<b>f</b>) PL2-3.</p>
Full article ">Figure 16
<p>The movements and trajectories of gravels in the GL group: (<b>a</b>) GL1; (<b>b</b>) GL2; (<b>c</b>) GL3.</p>
Full article ">Figure 17
<p>The movements and trajectories of gravels in the SL group: (<b>a</b>) SL1; (<b>b</b>) SL2; (<b>c</b>) SL3.</p>
Full article ">Figure 18
<p>The maximum displacement and maximum rotation speed curves of gravels in each group: (<b>a</b>) maximum displacement; (<b>b</b><span class="html-italic">)</span> maximum rotation speed.</p>
Full article ">Figure 19
<p>Displacement diagram of the contact points between typical gravels and the GCL during the rolling process.</p>
Full article ">Figure 20
<p>Local length diagram of the GCL.</p>
Full article ">Figure 21
<p>Local elongation of the GCL under actual working conditions.</p>
Full article ">Figure 22
<p>Definition of gravel embedding value.</p>
Full article ">Figure 23
<p>The change in gravel embedding value with the number of rolling passes.</p>
Full article ">Figure 24
<p>Observation points of the distribution of coordination number of bentonite.</p>
Full article ">Figure 25
<p>The change in the coordination number at the typical position of the gravel–GCL contact during rolling: (<b>a</b>) A contact position; (<b>b</b>) B contact position; (<b>c</b>) C contact position.</p>
Full article ">Figure 26
<p>Schematic diagram of GCL changes. (<b>a</b>) Effect of gravel surface contact on bentonite particles; (<b>b</b>) Effect of gravel point contact on bentonite particles.</p>
Full article ">
19 pages, 7362 KiB  
Article
Highly Efficient JR Optimization Technique for Solving Prediction Problem of Soil Organic Carbon on Large Scale
by Harsh Vazirani, Xiaofeng Wu, Anurag Srivastava, Debajyoti Dhar and Divyansh Pathak
Sensors 2024, 24(22), 7317; https://doi.org/10.3390/s24227317 (registering DOI) - 15 Nov 2024
Viewed by 245
Abstract
We utilized remote sensing and ground cover data to predict soil organic carbon (SOC) content across a vast geographic region. Employing a combination of machine learning and deep learning techniques, we developed a novel data fusion approach that integrated Digital Elevation Model (DEM) [...] Read more.
We utilized remote sensing and ground cover data to predict soil organic carbon (SOC) content across a vast geographic region. Employing a combination of machine learning and deep learning techniques, we developed a novel data fusion approach that integrated Digital Elevation Model (DEM) data, MODIS satellite imagery, WOSIS soil profile data, and CHELSA environmental data. This combined dataset, named GeoBlendMDWC, was specifically designed for SOC prediction. The primary aim of this research is to develop and evaluate a novel optimization algorithm for accurate SOC prediction by leveraging multi-source environmental data. Specifically, this study aims to (1) create an integrated dataset combining remote sensing and ground data for comprehensive SOC analysis, (2) develop a new optimization technique that enhances both machine learning and deep learning model performance, and (3) evaluate the algorithm’s efficiency and accuracy against established optimization methods like Jaya and GridSearchCV. This study focused on India, Australia, and South Africa, countries known for their significant agricultural activities. We introduced a novel optimization technique for both machine learning and deep neural networks, comparing its performance to established methods like the Jaya optimization technique and GridSearchCV. The models evaluated included XGBoost Regression, LightGBM, Gradient Boosting Regression (GBR), Random Forest Regression, Decision Tree Regression, and a Multilayer Perceptron (MLP) model. Our research demonstrated that the proposed optimization algorithm consistently outperformed existing methods in terms of execution time and performance. It achieved results comparable to GridSearchCV, reaching an R2 of 90.16, which was a significant improvement over the base XGBoost model’s R2 of 79.08. In deep learning optimization, it significantly outperformed the Jaya algorithm, achieving an R2 of 61.34 compared to Jaya’s 30.04. Moreover, it was 20–30 times faster than GridSearchCV. Given its speed and accuracy, this algorithm can be applied to real-time data processing in remote sensing satellites. This advanced methodology will greatly benefit the agriculture and farming sectors by providing precise SOC predictions. Full article
(This article belongs to the Section Remote Sensors)
Show Figures

Figure 1

Figure 1
<p>Soil organic carbon content by region in tons-per-hectare (ton/ha) in India [<a href="#B34-sensors-24-07317" class="html-bibr">34</a>], Australia [<a href="#B35-sensors-24-07317" class="html-bibr">35</a>], and Africa [<a href="#B36-sensors-24-07317" class="html-bibr">36</a>] respectively.</p>
Full article ">Figure 2
<p>Research workflow.</p>
Full article ">Figure 3
<p>Flowchart of the optimization algorithm.</p>
Full article ">Figure 4
<p>Correlation matrix for the Indian–Australian–African combined dataset. Pearson correlation methodology is used to calculate the correlation values.</p>
Full article ">Figure 5
<p>Average execution time comparison in milliseconds between the machine learning models when using different optimization techniques.</p>
Full article ">
27 pages, 13332 KiB  
Article
Numerical Prediction of Solid Particle Erosion in Jet Pumps Based on a Calibrated Model
by Xuanchen Wan, Mengxue Dong, Maosen Xu, Chuanhao Fan, Jiegang Mou and Shuai Han
Energies 2024, 17(22), 5720; https://doi.org/10.3390/en17225720 - 15 Nov 2024
Viewed by 243
Abstract
Jet pumps are widely used in petrochemical processes, nuclear cooling, and wastewater treatment due to their simple structure, high reliability, and stable performance under extreme conditions. However, when transporting solid-laden two-phase flows, they face severe erosion problems, leading to reduced efficiency, malfunctions, or [...] Read more.
Jet pumps are widely used in petrochemical processes, nuclear cooling, and wastewater treatment due to their simple structure, high reliability, and stable performance under extreme conditions. However, when transporting solid-laden two-phase flows, they face severe erosion problems, leading to reduced efficiency, malfunctions, or even failure. Therefore, optimizing jet pump performance and extending its service life is crucial. In this study, an experimental platform was established to conduct experiments on wall erosion in jet pumps. The CFD-DEM method was used to simulate the solid–liquid two-phase flow in the jet pump, comparing six erosion models for predicting erosion rates. The Grey Wolf Optimization algorithm was applied to calibrate model coefficients. The results indicate that the Neilson erosion model shows the best consistency with the experimental results. The inlet flow rate significantly influenced the erosion rates, while the flow rate ratio had a smaller effect. The particle concentration exhibited a nonlinear relationship with erosion, with diminishing impact beyond a certain threshold. As the factors varied, the erosion distribution tended to be uniform, but high erosion areas remained locally concentrated, indicating intensified localized erosion. Full article
Show Figures

Figure 1

Figure 1
<p>Structure and size of jet pump.</p>
Full article ">Figure 2
<p>Slot hole location and monitoring surface naming: (<b>a</b>) slot hole location and cross section; (<b>b</b>) monitoring surface naming.</p>
Full article ">Figure 3
<p>Experimental platform. 1—transparent jet pump; 2—power pump; 3—valve; 4—clean water tank; 5—sand-filled water tank; 6—stirrer; 7—inlet flow meter; 8—outlet flow meter; 9—inlet pressure sensor; 10—suction pressure sensor; 11—outlet pressure sensor; 12—high-speed camera.</p>
Full article ">Figure 4
<p>Center pressure coefficient curve.</p>
Full article ">Figure 5
<p>Computational domain and mesh details: (<b>a</b>) computational domain mesh; (<b>b</b>) mesh details.</p>
Full article ">Figure 6
<p>Jet pump external characteristic curve.</p>
Full article ">Figure 7
<p>Percentage of wear and tear on monitored surfaces.</p>
Full article ">Figure 8
<p>Erosion rate variation with inlet flow rate.</p>
Full article ">Figure 9
<p>Erosion rate variation with flow rate ratio.</p>
Full article ">Figure 10
<p>Erosion rate variation with particle size ratio.</p>
Full article ">Figure 11
<p>Erosion rate variation with particle volume concentration.</p>
Full article ">Figure A1
<p>Erosion rate versus time curve: (<b>a</b>) Archard erosion model; (<b>b</b>) DNV erosion model; (<b>c</b>) Zhang erosion model; (<b>d</b>) Oka erosion model; (<b>e</b>) Neilson erosion model; (<b>f</b>) Ahlert erosion model.</p>
Full article ">Figure A2
<p>Spatial division of throat with different flow directions and gravity effects: (<b>a</b>) X-axis direction; (<b>b</b>) Y-axis direction; (<b>c</b>) Z-axis direction.</p>
Full article ">Figure A3
<p>Erosion contour plots for different inlet flow rates: (<b>a</b>) <span class="html-italic">Q<sub>p</sub></span> = 1.03 m<sup>3</sup>/h; (<b>b</b>) <span class="html-italic">Q<sub>p</sub></span> = 2.07 m<sup>3</sup>/h; (<b>c</b>) <span class="html-italic">Q<sub>p</sub></span> = 4.14 m<sup>3</sup>/h; (<b>d</b>) <span class="html-italic">Q<sub>p</sub></span> = 8.28 m<sup>3</sup>/h; (<b>e</b>) <span class="html-italic">Q<sub>p</sub></span> = 16.56 m<sup>3</sup>/h.</p>
Full article ">Figure A4
<p>Erosion contour plots for different flow rate ratios: (<b>a</b>) <span class="html-italic">q</span> = 0.6; (<b>b</b>) <span class="html-italic">q</span> = 0.8; (<b>c</b>) <span class="html-italic">q</span> = 1.0; (<b>d</b>) <span class="html-italic">q</span> = 1.2; (<b>e</b>) <span class="html-italic">q</span> = 1.4.</p>
Full article ">Figure A5
<p>Erosion contour plots for different particle sizes: (<b>a</b>) <span class="html-italic">D</span> = 0.5 mm; (<b>b</b>) <span class="html-italic">D</span> = 0.75 mm; (<b>c</b>) <span class="html-italic">D</span> = 1.0 mm; (<b>d</b>) <span class="html-italic">D</span> = 1.25 mm; (<b>e</b>) <span class="html-italic">D</span> = 1.50 mm.</p>
Full article ">Figure A6
<p>Erosion contour plots for different particle concentrations: (<b>a</b>) <span class="html-italic">C<sub>v</sub></span> = 0.1%; (<b>b</b>) <span class="html-italic">C<sub>v</sub></span> = 0.5%; (<b>c</b>) <span class="html-italic">C<sub>v</sub></span> = 1.0%; (<b>d</b>) <span class="html-italic">C<sub>v</sub></span> = 2.0%; (<b>e</b>) <span class="html-italic">C<sub>v</sub></span> = 5.0%.</p>
Full article ">
17 pages, 1631 KiB  
Review
Three-Dimensional Printing of Metallic Parts by Means of Fused Filament Fabrication (FFF)
by Irene Buj-Corral, Felip Fenollosa-Artés and Joaquim Minguella-Canela
Metals 2024, 14(11), 1291; https://doi.org/10.3390/met14111291 - 14 Nov 2024
Viewed by 363
Abstract
Obtaining metallic parts via Additive Manufacturing can yield several advantages over using other traditional manufacturing methods such as machining. Material extrusion (MEX) can handle complex shapes with porous structures and, at the present time, much low-end and desktop equipment is available. In the [...] Read more.
Obtaining metallic parts via Additive Manufacturing can yield several advantages over using other traditional manufacturing methods such as machining. Material extrusion (MEX) can handle complex shapes with porous structures and, at the present time, much low-end and desktop equipment is available. In the present work, different industrial and medical applications of metallic Fused Filament Fabrication (FFF) parts are presented. First, an overview of the process, equipment, and of the metal-filled filaments currently available is provided. Then, the properties of parts and different applications are shown. For example, metal-filled filaments with a low metal content that can be used to obtain plastic parts with metallic appearance (with either steel, copper, or bronze), and filaments with a high metallic content allow obtaining metallic parts with high mechanical strength after a sintering operation. The present contribution aims to be an up-to-date panorama for current industrial and medical results and lessons learnt from the application of FFF to obtain metallic parts. Full article
(This article belongs to the Section Additive Manufacturing)
Show Figures

Figure 1

Figure 1
<p>Material extrusion equipment: (<b>a</b>) mono-extruder open-frame equipment, (<b>b</b>) multiple-extruder closed-frame equipment, (<b>c</b>) equipment for metal compound AM, and (<b>d</b>) oven post-processing of metal parts. Source: Laboratory of Manufacturing Technologies of ETSEIB and Fundació Privada Centre CIM.</p>
Full article ">Figure 2
<p>Metallic material extrusion part: (<b>a</b>) while being printed and (<b>b</b>) after being processed via debinding and sintering. Reprinted with permission from Ref. [<a href="#B56-metals-14-01291" class="html-bibr">56</a>]. Diseño y Fabricación de Maquinaria para Alimentación S.L. (DIFMAC ROURE), Santa Perpètua de Mogoda, Spain, is the owner of the parts.</p>
Full article ">Figure 3
<p>Stiffening strategy for FDMS parts to avoid deformation: (<b>a</b>) printed block of part and supports with interface material and (<b>b</b>) final part after sintering. Reprinted with permission from Ref. [<a href="#B56-metals-14-01291" class="html-bibr">56</a>]. Diseño y Fabricación de Maquinaria para Alimentación S.L. (DIFMAC ROURE), Santa Perpètua de Mogoda, Spain, is the owner of the parts.</p>
Full article ">Figure 4
<p>Dimensional measurement of a steel-filled PLA dogbone specimen with a Mitutoyo micrometer.</p>
Full article ">Figure 5
<p>Three-dimensional-printed disks in (<b>a</b>) bronze, (<b>b</b>) copper, and (<b>c</b>) stainless-steel. Reprinted from Ref. [<a href="#B77-metals-14-01291" class="html-bibr">77</a>].</p>
Full article ">
34 pages, 41034 KiB  
Article
The Dynamics of Air Pollution in the Southwestern Part of the Caspian Sea Basin (Based on the Analysis of Sentinel-5 Satellite Data Utilizing the Google Earth Engine Cloud-Computing Platform)
by Vladimir Tabunshchik, Aleksandra Nikiforova, Nastasia Lineva, Polina Drygval, Roman Gorbunov, Tatiana Gorbunova, Ibragim Kerimov, Cam Nhung Pham, Nikolai Bratanov and Mariia Kiseleva
Atmosphere 2024, 15(11), 1371; https://doi.org/10.3390/atmos15111371 - 14 Nov 2024
Viewed by 271
Abstract
The Caspian region represents a complex and unique system of terrestrial, coastal, and aquatic environments, marked by an exceptional landscape and biological diversity. This diversity, however, is increasingly threatened by substantial anthropogenic pressures. One notable impact of this human influence is the rising [...] Read more.
The Caspian region represents a complex and unique system of terrestrial, coastal, and aquatic environments, marked by an exceptional landscape and biological diversity. This diversity, however, is increasingly threatened by substantial anthropogenic pressures. One notable impact of this human influence is the rising concentration of pollutants atypical for the atmosphere. Advances in science and technology now make it possible to detect certain atmospheric pollutants using remote Earth observation techniques, specifically through data from the Sentinel-5 satellite, which provides continuous insights into atmospheric contamination. This article investigates the dynamics of atmospheric pollution in the southwestern part of the Caspian Sea basin using Sentinel-5P satellite data and the cloud-computing capabilities of the Google Earth Engine (GEE) platform. The study encompasses an analysis of concentrations of seven key pollutants: nitrogen dioxide (NO2), formaldehyde (HCHO), carbon monoxide (CO), ozone (O3), sulfur dioxide (SO2), methane (CH4), and the Aerosol Index (AI). Spatial and temporal variations in pollution fields were examined for the Caspian region and the basins of the seven rivers (key areas) flowing into the Caspian Sea: Sunzha, Sulak, Ulluchay, Karachay, Atachay, Haraz, and Gorgan. The research methodology is based on the use of data from the Sentinel-5 satellite, SRTM DEM data on absolute elevations, surface temperature data, and population density data. Data processing is performed using the Google Earth Engine cloud-computing platform and the ArcGIS software suite. The main aim of this study is to evaluate the spatiotemporal variability of pollutant concentration fields in these regions from 2018 to 2023 and to identify the primary factors influencing pollution distribution. The study’s findings reveal that the Heraz and Gorgan River basins have the highest concentrations of nitrogen dioxide and Aerosol Index levels, marking these basins as the most vulnerable to atmospheric pollution among those assessed. Additionally, the Gorgan basin exhibited elevated carbon monoxide levels, while the highest ozone concentrations were detected in the Sunzha basin. Our temporal analysis demonstrated a substantial influence of the COVID-19 pandemic on pollutant dispersion patterns. Our correlation analysis identified absolute elevation as a key factor affecting pollutant distribution, particularly for carbon monoxide, ozone, and aerosol indices. Population density showed the strongest correlation with nitrogen dioxide distribution. Other pollutants exhibited more complex distribution patterns, influenced by diverse mechanisms associated with local emission sources and atmospheric dynamics. Full article
(This article belongs to the Special Issue Study of Air Pollution Based on Remote Sensing (2nd Edition))
Show Figures

Figure 1

Figure 1
<p>Geographic location of the study area.</p>
Full article ">Figure 2
<p>Overall methodology for air pollution assessment and the evaluation of its relationship with geographical factors.</p>
Full article ">Figure 3
<p>Distribution of Aerosol Index: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 3 Cont.
<p>Distribution of Aerosol Index: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 4
<p>Distribution of nitrogen dioxide concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 5
<p>Distribution of ozone concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 5 Cont.
<p>Distribution of ozone concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 6
<p>Distribution of carbon monoxide concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 6 Cont.
<p>Distribution of carbon monoxide concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 7
<p>Distribution of formaldehyde concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 8
<p>Distribution of sulfur dioxide concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 9
<p>Distribution of methane concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 9 Cont.
<p>Distribution of methane concentrations: (<b>a</b>) 2018, (<b>b</b>) 2019, (<b>c</b>) 2020, (<b>d</b>) 2021, (<b>e</b>) 2022, and (<b>f</b>) 2023.</p>
Full article ">Figure 10
<p>Intra-annual dynamics of average Aerosol Index concentrations in the river basins of Sunzha, Sulak, Ulluchay, Karachay, Atachay, Haraz, and Gorgan.</p>
Full article ">Figure 11
<p>Intra-annual dynamics of average carbon monoxide concentrations in the river basins of Sunzha, Sulak, Ulluchay, Karachay, Atachay, Haraz, and Gorgan.</p>
Full article ">Figure 12
<p>Intra-annual dynamics of average formaldehyde concentrations in the river basins of Sunzha, Sulak, Ulluchay, Karachay, Atachay, Haraz, and Gorgan.</p>
Full article ">Figure 13
<p>Intra-annual dynamics of average carbon ozone concentrations in the river basins of Sunzha, Sulak, Ulluchay, Karachay, Atachay, Haraz, and Gorgan.</p>
Full article ">Figure 14
<p>Intra-annual dynamics of average nitrogen dioxide concentrations in the river basins of Sunzha, Sulak, Ulluchay, Karachay, Atachay, Haraz, and Gorgan.</p>
Full article ">Figure 15
<p>Intra-annual dynamics of average sulfur dioxide concentrations in the river basins of Sunzha, Sulak, Ulluchay, Karachay, Atachay, Haraz, and Gorgan.</p>
Full article ">Figure 16
<p>Correlation diagrams between pollution indicators and geographic factors in the Caspian region: (<b>a</b>) 2019, (<b>b</b>) 2020, (<b>c</b>) 2021, and (<b>d</b>) 2022.</p>
Full article ">Figure 17
<p>Correlation diagrams between pollution indicators and geographic factors in the Sunzha River basin: (<b>a</b>) 2019, (<b>b</b>) 2020, (<b>c</b>) 2021, and (<b>d</b>) 2022.</p>
Full article ">Figure 18
<p>Correlation diagrams between pollution indicators and geographic factors in the Sulak River basin: (<b>a</b>) 2019, (<b>b</b>) 2020, (<b>c</b>) 2021, and (<b>d</b>) 2022.</p>
Full article ">Figure 19
<p>Correlation diagrams between pollution indicators and geographic factors in the Ulluchay River basin: (<b>a</b>) 2019, (<b>b</b>) 2020, (<b>c</b>) 2021, and (<b>d</b>) 2022.</p>
Full article ">Figure 20
<p>Correlation diagrams between pollution indicators and geographic factors in the Karachay River basin: (<b>a</b>) 2019, (<b>b</b>) 2020, (<b>c</b>) 2021, and (<b>d</b>) 2022.</p>
Full article ">Figure 21
<p>Correlation diagrams between pollution indicators and geographic factors in the Atachay River basin: (<b>a</b>) 2019, (<b>b</b>) 2020, (<b>c</b>) 2021, and (<b>d</b>) 2022.</p>
Full article ">Figure 22
<p>Correlation diagrams between pollution indicators and geographic factors in the Haraz River basin: (<b>a</b>) 2019, (<b>b</b>) 2020, (<b>c</b>) 2021, and (<b>d</b>) 2022.</p>
Full article ">Figure 23
<p>Correlation diagrams between pollution indicators and geographic factors in the Gorgan River basin: (<b>a</b>) 2019, (<b>b</b>) 2020, (<b>c</b>) 2021, and (<b>d</b>) 2022.</p>
Full article ">Figure A1
<p>Average Aerosol Index values from 2018 to 2023 in the basins of small and medium rivers of the Caspian Sea.</p>
Full article ">Figure A2
<p>Average nitrogen dioxide concentrations from 2018 to 2023 in the basins of small- and medium-sized rivers of the Caspian Sea.</p>
Full article ">Figure A3
<p>Average tropospheric ozone concentrations from 2018 to 2023 in the basins of small- and medium-sized rivers of the Caspian Sea.</p>
Full article ">Figure A4
<p>Average carbon monoxide concentrations in the basins of small- and medium-sized rivers in the Caspian Sea from 2018 to 2023.</p>
Full article ">Figure A5
<p>Average formaldehyde concentrations from 2018 to 2023 in the basins of small- and medium-sized rivers in the Caspian Sea region.</p>
Full article ">Figure A6
<p>Average sulfur dioxide concentrations in the basins of small- and medium-sized rivers in the Caspian Sea from 2018 to 2023.</p>
Full article ">
15 pages, 14788 KiB  
Article
The DEM Registration Method Without Ground Control Points for Landslide Deformation Monitoring
by Yunchuan Wang, Jia Li, Ping Duan, Rui Wang and Xinrui Yu
Remote Sens. 2024, 16(22), 4236; https://doi.org/10.3390/rs16224236 - 14 Nov 2024
Viewed by 279
Abstract
Landslides are geological disasters that are harmful to both humans and society. Digital elevation model (DEM) time series data are usually used to monitor dynamic changes or surface damage. To solve the problem of landslide deformation monitoring without ground control points (GCPs), a [...] Read more.
Landslides are geological disasters that are harmful to both humans and society. Digital elevation model (DEM) time series data are usually used to monitor dynamic changes or surface damage. To solve the problem of landslide deformation monitoring without ground control points (GCPs), a multidimensional feature-based coregistration method (MFBR) was studied to achieve accurate registration of multitemporal DEMs without GCPs and obtain landslide deformation information. The method first derives the elevation information of the DEM into image pixel information, and the feature points are extracted on the basis of the image. The initial plane position registration of the DEM is implemented. Therefore, the expected maximum algorithm is applied to calculate the stable regions that have not changed between multitemporal DEMs and to perform accurate registrations. Finally, the shape variables are calculated by constructing a DEM differential model. The method was evaluated using simulated data and data from two real landslide cases, and the experimental results revealed that the registration accuracies of the three datasets were 0.963 m, 0.368 m, and 2.459 m, which are 92%, 50%, and 24% better than the 12.189 m, 0.745 m, and 3.258 m accuracies of the iterative closest-point algorithm, respectively. Compared with the GCP-based method, the MFBR method can achieve 70% deformation acquisition capability, which indicates that the MFBR method has better applicability in the field of landslide monitoring. This study provides an idea for landslide deformation monitoring without GCPs and is helpful for further understanding the state and behavior of landslides. Full article
(This article belongs to the Special Issue Advances in GIS and Remote Sensing Applications in Natural Hazards)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Workflow of the MFBR method.</p>
Full article ">Figure 2
<p>Image-based DEM position registration method.</p>
Full article ">Figure 3
<p>3D spatial-feature-based precision registration method.</p>
Full article ">Figure 4
<p>Simulation dataset.</p>
Full article ">Figure 5
<p>Luchun landslide and data collection.</p>
Full article ">Figure 6
<p>The Gongshan landslide.</p>
Full article ">Figure 7
<p>Deformation detection results of the simulation data. (<b>a</b>) Raw deformation. (<b>b</b>) Based on the MFBR method.</p>
Full article ">Figure 8
<p>Deformation detection results for the Luchun landslide. (<b>a</b>) Based on GCPs. (<b>b</b>) Based on MFBR. (<b>c</b>) ICP-based method.</p>
Full article ">Figure 9
<p>Deformation detection results for the Gongshan landslide. (<b>a</b>) Based on GCPs. (<b>b</b>) Based on MFBR. (<b>c</b>) ICP-based method.</p>
Full article ">Figure 10
<p>Stable region extraction results.</p>
Full article ">
29 pages, 6797 KiB  
Article
Integrating Physiology, Transcriptome, and Metabolome Analyses Reveals the Drought Response in Two Quinoa Cultivars with Contrasting Drought Tolerance
by Yang Wang, Yang Wu, Qinghan Bao, Huimin Shi and Yongping Zhang
Int. J. Mol. Sci. 2024, 25(22), 12188; https://doi.org/10.3390/ijms252212188 - 13 Nov 2024
Viewed by 325
Abstract
Quinoa (Chenopodium quinoa Willd.) is an annual broadleaf plant belonging to the Amaranthaceae family. It is a nutritious food crop and is considered to be drought-tolerant, but drought is still one of the most important abiotic stress factors limiting its yield. Quinoa [...] Read more.
Quinoa (Chenopodium quinoa Willd.) is an annual broadleaf plant belonging to the Amaranthaceae family. It is a nutritious food crop and is considered to be drought-tolerant, but drought is still one of the most important abiotic stress factors limiting its yield. Quinoa responses to drought are related to drought intensity and genotype. This study used two different drought-responsive quinoa cultivars, LL1 (drought-tolerant) and ZK1 (drought-sensitive), to reveal the important mechanisms of drought response in quinoa by combining physiological, transcriptomic, and metabolomic analyses. The physiological analysis indicated that Chla/Chlb might be important for drought tolerance in quinoa. A total of 1756 and 764 differentially expressed genes (DEGs) were identified in LL1 and ZK1, respectively. GO (Gene Ontology) enrichment analysis identified 52 common GO terms, but response to abscisic acid (GO:0009737) and response to osmotic stress (GO:0006970) were only enriched in LL1. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis revealed that glycerophospholipid metabolism (ko00564) and cysteine and methionine metabolism (ko00270) ranked at the top of the list in both cultivars. A total of 1844 metabolites were identified by metabolomic analysis. “Lipids and lipid-like” molecules had the highest proportions. The DEMs in LL1 and ZK1 were mainly categorized 6 and 4 Human Metabolome Database (HMDB) superclasses, respectively. KEGG analysis revealed that the ‘α-linolenic acid metabolism’ was enriched in both LL1 and ZK1. Joint KEGG analysis also revealed that the ‘α-linolenic acid metabolism’ pathway was enriched by both the DEGs and DEMs of LL1. There were 17 DEGs and 8 DEMs enriched in this pathway, and methyl jasmonate (MeJA) may play an important role in the drought response of quinoa. This study will provide information for the identification of drought resistance in quinoa, research on the molecular mechanism of drought resistance, and genetic breeding for drought resistance in quinoa. Full article
(This article belongs to the Special Issue Advance in Plant Abiotic Stress: 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>The effects of different drought stress intensities on quinoa seedling height and leaf area. (<b>A</b>) Plant height. (<b>B</b>) Leaf surface area. Note: W1 represents control group, W2, W3, and W4 represent mild, moderate, and severe drought stress, respectively. The same as below. Vertical bars indicate the mean value ± SD (<span class="html-italic">n</span> = 3). The different lowercase letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 2
<p>The effects of different drought stress intensities on quinoa biomass and root–shoot ratio. (<b>A</b>) The above-ground fresh weight (AGFW). (<b>B</b>) The under-ground fresh weight (UGFW). (<b>C</b>) The above-ground dry weight (AGDW). (<b>D</b>) The under-ground dry weight (UGDW). (<b>E</b>) The root–shoot ratio (RSR). Note: vertical bars indicate the mean value ± SD (<span class="html-italic">n</span> = 3). The different lowercase letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 3
<p>The effects of different drought stress intensities on quinoa relative water content. Note: vertical bars indicate the mean value ± SD (<span class="html-italic">n</span> = 3). The different lowercase letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 4
<p>The effects of different drought stress intensities on quinoa root vigor. Note: vertical bars indicate the mean value ± SD (<span class="html-italic">n</span> = 3). The different lowercase letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 5
<p>The effects of different drought stress intensities on quinoa soluble sugars and antioxidase. (<b>A</b>) The soluble sugar content. (<b>B</b>) Enzymatic activity of CAT. (<b>C</b>) Enzymatic activity of SOD. (<b>D</b>) Enzymatic activity of POD. Note: vertical bars indicate the mean value ± SD (<span class="html-italic">n</span> = 3). The different lowercase letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 6
<p>The effects of different drought stress intensities on quinoa chlorophyll. (<b>A</b>) Chlorophyll-a content. (<b>B</b>) Chlorophyll-b content. (<b>C</b>) Total chlorophyll. (<b>D</b>) Chlorophyll-a/chlorophyll-b. Note: vertical bars indicate the mean value ± SD (<span class="html-italic">n</span> = 3). The different lowercase letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 7
<p>The effects of different drought stress intensities on quinoa photosynthetic properties. (<b>A</b>) The assimilation rate. (<b>B</b>) The transpiration rate. (<b>C</b>) The internal CO<sub>2</sub> content. (<b>D</b>) Stomatal conductance. (<b>E</b>) Water use efficiency. Note: vertical bars indicate the mean value ± SD (<span class="html-italic">n</span> = 3). The different lowercase letters indicate a significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 8
<p>Correlation analysis of physiological and biochemical indicators.</p>
Full article ">Figure 9
<p>Gene expression clustering, principal components analysis (PCA), and correlation analysis. (<b>A</b>) Cluster analysis of gene expression. (<b>B</b>) Principle components analysis. (<b>C</b>) Correlation analysis. (<b>D</b>) LL1 gene expression volcano plot. (<b>E</b>) ZK1 gene expression volcano plot.</p>
Full article ">Figure 9 Cont.
<p>Gene expression clustering, principal components analysis (PCA), and correlation analysis. (<b>A</b>) Cluster analysis of gene expression. (<b>B</b>) Principle components analysis. (<b>C</b>) Correlation analysis. (<b>D</b>) LL1 gene expression volcano plot. (<b>E</b>) ZK1 gene expression volcano plot.</p>
Full article ">Figure 10
<p>Bar plot for GO enrichment of DEGs. (<b>A</b>) GO enrichment bar plot of the LL1 DEGs. (<b>B</b>) GO enrichment bar plot of the ZK1 DEGs.</p>
Full article ">Figure 11
<p>Bar plot for KEGG enrichment of DEGs. (<b>A</b>) KEGG enrichment plot of the LL1 DEGs. (<b>B</b>) KEGG enrichment plot of the ZK1 DEGs. Note: bar length represents the number of genes.</p>
Full article ">Figure 12
<p>DEGs in the common KEGG pathway.</p>
Full article ">Figure 13
<p>HMDB annotation results. (<b>A</b>) HMDB annotation results of LL1. (<b>B</b>) HMDB annotation results of ZK1.</p>
Full article ">Figure 14
<p>DEMs KEGG enrichment analysis. (<b>A</b>) DEMs KEGG enrichment analysis of LL1. (<b>B</b>) DEMs KEGG enrichment analysis of ZK1.</p>
Full article ">Figure 15
<p>Expressions of DEGs and DEMs in α-linolenic acid metabolism pathway. Heatmap using gene FPKM values and metabolite abundance. Note: Red triangles represent downregulated expression of DEGs and DEMs, green triangles represent up-regulated.</p>
Full article ">
17 pages, 23726 KiB  
Article
Construction and Analysis of miRNA–mRNA Interaction Network in Ovarian Tissue of Wanxi White Geese Across Different Breeding Stages
by Ruidong Li, Yuhua Wang, Fei Xie, Xinwei Tong, Xiaojin Li, Man Ren, Qianqian Hu and Shenghe Li
Animals 2024, 14(22), 3258; https://doi.org/10.3390/ani14223258 - 13 Nov 2024
Viewed by 315
Abstract
Ovarian development significantly influences the laying performance of geese. In this study, the transcriptome analysis was conducted on the ovarian tissues of Wanxi White Geese during the pre-laying (KL), laying (CL), and ceased-laying period (XL). Short Time-series Expression Miner (STEM) analysis and miRNA–mRNA [...] Read more.
Ovarian development significantly influences the laying performance of geese. In this study, the transcriptome analysis was conducted on the ovarian tissues of Wanxi White Geese during the pre-laying (KL), laying (CL), and ceased-laying period (XL). Short Time-series Expression Miner (STEM) analysis and miRNA–mRNA regulatory network construction were performed to identify the key genes and miRNAs regulating laying traits. Comparative analysis of KL vs. CL, CL vs. XL, and XL vs. KL groups resulted in the identification of 337, 136, and 525 differentially expressed genes (DEGs), and 258, 1131, and 909 differentially expressed miRNAs (DEMs), respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis (p < 0.05) revealed that the main enrichment pathways of DEGs and DEMs at different breeding periods were Neuroactive ligand–receptor interaction, GnRH signaling pathway and Wnt signaling pathway, all associated with ovarian development. According to the three groups of common pathways, four DEGs were screened out, including INHBB, BMP5, PRL, and CGA, along with five DEMs, including let-7-x, miR-124-y, miR-1-y, and miR-10926-z, all of them may affect ovarian development. A miRNA–mRNA regulatory network was constructed through integrated analysis of DEGs and DEMs, revealing nine miRNAs highly associated with ovarian development: miR-101-y, let-7-x, miR-1-x, miR-17-y, miR-103-z, miR-204-x, miR-101-x, miR-301-y, and miR-151-x. The dual-luciferase reporter gene verified the target relationship between WIF1 and miR-204-x, suggesting that these miRNAs may influence ovarian development in Wanxi White Goose by regulating the expression levels of their target genes within ovarian tissue. This study provides a theoretical foundation for analyzing the mechanisms of ovarian development across different breeding periods and accelerating the cultivation of new breeds through post-transcriptional regulation levels. Full article
(This article belongs to the Section Animal Genetics and Genomics)
Show Figures

Figure 1

Figure 1
<p>Ovarian histological analysis of Wanxi White Geese across different egg-laying periods: (<b>A</b>) pre-laying period; (<b>B</b>) laying period; (<b>C</b>) eased period.</p>
Full article ">Figure 2
<p>The DEGs in ovary tissues of Wanxi White Geese in different periods: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL.</p>
Full article ">Figure 3
<p>Functional analysis of DEGs GO in different periods: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL.</p>
Full article ">Figure 4
<p>Functional enrichment analysis of DEGs KEGG in different egg-laying periods: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL.</p>
Full article ">Figure 5
<p>The DEMs in ovary tissues of Wanxi White Geese in different periods: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL.</p>
Full article ">Figure 6
<p>GO function analysis of DEMs in different periods: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL.</p>
Full article ">Figure 7
<p>KEGG function analysis of DEMs in different reproductive periods: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL.</p>
Full article ">Figure 8
<p>mRNA STEM analysis diagram.</p>
Full article ">Figure 9
<p>KEGG functional enrichment analysis of dynamic DEGs: (<b>A</b>) Profile 2; (<b>B</b>) Profile 5.</p>
Full article ">Figure 10
<p>miRNA STEM analysis.</p>
Full article ">Figure 11
<p>Profile 5 KEGG functional enrichment analysis.</p>
Full article ">Figure 12
<p>The intersection of differentially expressed genes and miRNA target genes at different reproductive stages Wayne diagram: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL.</p>
Full article ">Figure 13
<p>KL vs. CL miRNA–mRNA interaction network analysis diagram: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL. Red is up-regulated miRNA, yellow is down-regulated miRNA, and blue is mRNA.</p>
Full article ">Figure 14
<p>KEGG enrichment analysis of interaction network: (<b>A</b>) KL vs. CL; (<b>B</b>) CL vs. XL; (<b>C</b>) KL vs. XL.</p>
Full article ">Figure 15
<p>RT-qPCR validation of RNA-seq results.</p>
Full article ">Figure 16
<p>Double luciferase data analysis diagram.</p>
Full article ">
14 pages, 3818 KiB  
Article
Compaction Evolution and Mechanisms of Granular Materials Due to Gyratory Shearing
by Teng Man
Materials 2024, 17(22), 5525; https://doi.org/10.3390/ma17225525 - 12 Nov 2024
Viewed by 352
Abstract
Granular systems, no matter whether they are dry or saturated, are commonly encountered in both natural scenarios and engineering applications. In this work, we tackle the compaction problem of both dry and saturated granular systems under gyratory shearing compaction, where particles are subjected [...] Read more.
Granular systems, no matter whether they are dry or saturated, are commonly encountered in both natural scenarios and engineering applications. In this work, we tackle the compaction problem of both dry and saturated granular systems under gyratory shearing compaction, where particles are subjected to constant pressure and continuous shear rate, which is quite different from the traditional cyclic shearing compaction. Such phenomena are crucial to the compaction of asphalt mixtures or soils in civil engineering and can be extended to other areas, such as powder processing and pharmaceutical engineering. In this study, we investigated the behavior of both dry and fully saturated mono-dispersed granular materials under gyratory shearing compaction using the discrete element method (DEM) and found that the gyratory speed or interstitial fluid viscosity has almost no impact on the compaction behavior, while the pressure and the particle size play more important roles. Additionally, it is the inertial time scale which dictates the compaction behavior under gyratory shearing in most cases; meanwhile, the viscous time scale can also have influence in some conditions. This work determines the similarity and unity between the granular gyratory compaction and the rheology of granular systems, which has direct relevance to various natural and engineering systems. Full article
Show Figures

Figure 1

Figure 1
<p>Configuration of the gyratory compactor and the model setup in the DEM simulation with four steps: (<b>I</b>,<b>II</b>) dropping particles, (<b>III</b>) adding the top plate and tilt the cylindrical ring, and (<b>IV</b>) starting the gyratory compaction.</p>
Full article ">Figure 2
<p>Compaction curves of granular materials with different conditions. In (<b>a</b>), we plotted the results of simulations with different gyratory speeds while keeping the pressure and the interstitial fluid viscosity constant at 100 kPa and 100 cP, respectively. We plotted the simulations with different interstitial fluid viscosity in (<b>b</b>,<b>c</b>). However, the pressure is different in these two sets of simulations. In (<b>d</b>), we have the relationship between solid fraction and dimensionless time, <math display="inline"><semantics> <mover accent="true"> <mi>T</mi> <mo stretchy="false">^</mo> </mover> </semantics></math>, of simulations with different pressures.</p>
Full article ">Figure 3
<p>Schematic plot of both the shearing dilation and the compressive relaxation. (<b>a</b>) shows the start of a shearing motion and (<b>b</b>) shows the final state after a shear step.</p>
Full article ">Figure 4
<p>Histogram of the contact forces between particles: (<b>a</b>) the histogram in a linear-log coordinate system and (<b>b</b>) the histogram in a log-log coordinate system. Markers <span class="html-fig-inline" id="materials-17-05525-i001"><img alt="Materials 17 05525 i001" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i001.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i002"><img alt="Materials 17 05525 i002" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i002.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i003"><img alt="Materials 17 05525 i003" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i003.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i004"><img alt="Materials 17 05525 i004" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i004.png"/></span> represent simulation results of different gyratory speeds. <span class="html-fig-inline" id="materials-17-05525-i005"><img alt="Materials 17 05525 i005" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i005.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i006"><img alt="Materials 17 05525 i006" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i006.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i007"><img alt="Materials 17 05525 i007" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i007.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i008"><img alt="Materials 17 05525 i008" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i008.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i009"><img alt="Materials 17 05525 i009" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i009.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i010"><img alt="Materials 17 05525 i010" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i010.png"/></span> represent simulation results of different viscosity when <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <mi>n</mi> </msub> <mo>=</mo> <mn>100</mn> </mrow> </semantics></math> kPa, while <span class="html-fig-inline" id="materials-17-05525-i011"><img alt="Materials 17 05525 i011" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i011.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i012"><img alt="Materials 17 05525 i012" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i012.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i013"><img alt="Materials 17 05525 i013" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i013.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i014"><img alt="Materials 17 05525 i014" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i014.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i015"><img alt="Materials 17 05525 i015" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i015.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i016"><img alt="Materials 17 05525 i016" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i016.png"/></span> represent simulation results of different viscosity when <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <mi>n</mi> </msub> <mo>=</mo> <mn>600</mn> </mrow> </semantics></math> kPa. <span class="html-fig-inline" id="materials-17-05525-i017"><img alt="Materials 17 05525 i017" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i017.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i018"><img alt="Materials 17 05525 i018" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i018.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i019"><img alt="Materials 17 05525 i019" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i019.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i020"><img alt="Materials 17 05525 i020" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i020.png"/></span>, <span class="html-fig-inline" id="materials-17-05525-i021"><img alt="Materials 17 05525 i021" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i021.png"/></span> represent simulation with different pressures.</p>
Full article ">Figure 5
<p>(<b>a</b>) Relationship between the solid fraction, <math display="inline"><semantics> <mi>ϕ</mi> </semantics></math>, and the dimensionless time, <math display="inline"><semantics> <mrow> <mover accent="true"> <mi>T</mi> <mo stretchy="false">˜</mo> </mover> <mo>=</mo> <mi>t</mi> <msqrt> <mrow> <msub> <mi>σ</mi> <mi>n</mi> </msub> <mo>/</mo> <msub> <mi>ρ</mi> <mi>p</mi> </msub> </mrow> </msqrt> <mo>/</mo> <mover accent="true"> <mi>d</mi> <mo stretchy="false">¯</mo> </mover> </mrow> </semantics></math>. (<b>b</b>) Relationship between the “visco-plastic” solid fraction, <math display="inline"><semantics> <mrow> <mi>ϕ</mi> <mo>−</mo> <mo>Δ</mo> <msub> <mi>ϕ</mi> <mi>e</mi> </msub> </mrow> </semantics></math>, and the dimensionless time, <math display="inline"><semantics> <mover accent="true"> <mi>T</mi> <mo stretchy="false">˜</mo> </mover> </semantics></math>.</p>
Full article ">Figure 6
<p>Relationship between the change in solid fraction, <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>ϕ</mi> <mi>e</mi> </msub> </mrow> </semantics></math>, and the relative deformation between adjacent particles, <math display="inline"><semantics> <mrow> <msub> <mi>δ</mi> <mi>n</mi> </msub> <mo>/</mo> <msub> <mi>r</mi> <mi>p</mi> </msub> </mrow> </semantics></math>, where <math display="inline"><semantics> <msub> <mi>r</mi> <mi>p</mi> </msub> </semantics></math> is the particle radius. The inset shows the basic structure of closely packed particles used for calculating the change in solid fractions.</p>
Full article ">Figure 7
<p>(<b>a</b>) The relationship between the solid fraction and time for systems with different particle sizes. (<b>b</b>) We vary the <span class="html-italic">x</span>-axis of Figure (<b>a</b>) to the dimensionless time <math display="inline"><semantics> <mrow> <mover accent="true"> <mi>T</mi> <mo stretchy="false">˜</mo> </mover> <mo>=</mo> <mi>t</mi> <mo>/</mo> <mrow> <mo>(</mo> <mi>d</mi> <mo>/</mo> <msqrt> <mrow> <msub> <mi>σ</mi> <mi>n</mi> </msub> <mo>/</mo> <msub> <mi>ρ</mi> <mi>p</mi> </msub> </mrow> </msqrt> <mo>)</mo> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 8
<p>(<b>a</b>) Histogram of the ratio between the tangential and the normal contact forces, <math display="inline"><semantics> <mrow> <msub> <mi>f</mi> <mi>t</mi> </msub> <mo>/</mo> <msub> <mi>f</mi> <mi>n</mi> </msub> </mrow> </semantics></math>, when three systems have the same compacted solid fraction. (<b>b</b>) Histogram of <math display="inline"><semantics> <mrow> <msub> <mi>f</mi> <mi>t</mi> </msub> <mo>/</mo> <msub> <mi>f</mi> <mi>n</mi> </msub> </mrow> </semantics></math>, when three systems are the the same dimensionless time, <math display="inline"><semantics> <mover accent="true"> <mi>T</mi> <mo stretchy="false">˜</mo> </mover> </semantics></math>.</p>
Full article ">Figure 9
<p>(<b>a</b>) shows the loading condition, where the compaction pressure varies between 20 kPa and 5120 kPa. (<b>b</b>) shows the relationship between compaction pressure and the steady-state solid fractions. Marker (<span class="html-fig-inline" id="materials-17-05525-i025"><img alt="Materials 17 05525 i025" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i025.png"/></span>) shows the cases of the change in solid fractions when we increase the pressure from 20 kPa to around 5120 kPa (the first loading process). Marker (<span class="html-fig-inline" id="materials-17-05525-i026"><img alt="Materials 17 05525 i026" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i026.png"/></span>) shows the results when we decrease the pressure from 5120 kPa to 20 kPa (the unloading process). Marker (<span class="html-fig-inline" id="materials-17-05525-i027"><img alt="Materials 17 05525 i027" src="/materials/materials-17-05525/article_deploy/html/images/materials-17-05525-i027.png"/></span>) shows the results when we increase the pressure from 20 kPa to 5120 kPa again (the second loading process).</p>
Full article ">
20 pages, 6927 KiB  
Article
High-Resolution Spaceborne SAR Geolocation Accuracy Analysis and Error Correction
by Facheng Li and Qiming Zeng
Remote Sens. 2024, 16(22), 4210; https://doi.org/10.3390/rs16224210 - 12 Nov 2024
Viewed by 437
Abstract
High-accuracy geolocation is crucial for high-resolution spaceborne SAR images. Most advanced SAR satellites have a theoretical geolocation accuracy better than 1 m, but this may be unrealizable with less accurate external data, such as atmospheric parameters and ground elevations. To investigate the actual [...] Read more.
High-accuracy geolocation is crucial for high-resolution spaceborne SAR images. Most advanced SAR satellites have a theoretical geolocation accuracy better than 1 m, but this may be unrealizable with less accurate external data, such as atmospheric parameters and ground elevations. To investigate the actual SAR geolocation accuracy in common applications, we analyze the properties of different geolocation errors, propose a geolocation procedure, and conduct experiments on TerraSAR-X images and a pair of Tianhui-2 images. The results show that based on GNSS elevations, the geolocation accuracy is better than 1 m for TerraSAR-X and 2 m/4 m for the Tianhui-2 reference/secondary satellites. Based on the WorldDEM and the SRTM, additional geolocation errors of 2 m and 4 m are introduced, respectively. By comparing the effectiveness of different tropospheric correction methods, we find that the GACOS mapping method has advantages in terms of resolution and computational efficiency. We conclude that tropospheric errors and ground elevation errors are the primary factors influencing geolocation accuracy, and the key to improving accuracy is to use higher-accuracy DEMs. Additionally, we propose and validate a geolocation model for the Tianhui-2 secondary satellite. Full article
(This article belongs to the Special Issue SAR Images Processing and Analysis (2nd Edition))
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The three kinds of SAR geolocation errors in the range: (1) <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>r</mi> </mrow> <mrow> <mi>r</mi> <mi>g</mi> </mrow> </msub> </mrow> </semantics></math> is related to the satellite orbit error; <math display="inline"><semantics> <mrow> <mi>S</mi> </mrow> </semantics></math> is the correct position, while <math display="inline"><semantics> <mrow> <msup> <mi>S</mi> <mo>′</mo> </msup> </mrow> </semantics></math> is the incorrect one; <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>r</mi> </mrow> <mrow> <mi>r</mi> <mi>g</mi> </mrow> </msub> <mo>=</mo> <mfenced open="|" close="|" separators="|"> <mrow> <mi>S</mi> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> </mfenced> <mo>−</mo> <mfenced open="|" close="|" separators="|"> <mrow> <msup> <mi>S</mi> <mo>′</mo> </msup> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> </mfenced> </mrow> </semantics></math>. (2) <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>a</mi> <mi>t</mi> <mi>m</mi> </mrow> </msub> </mrow> </semantics></math> is associated with the atmospheric error. The signal’s velocity in the atmosphere, <math display="inline"><semantics> <mrow> <mi>v</mi> </mrow> </semantics></math>, is smaller than in a vacuum, <math display="inline"><semantics> <mrow> <mi>c</mi> </mrow> </semantics></math>; <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>a</mi> <mi>t</mi> <mi>m</mi> </mrow> </msub> <mo>=</mo> <mstyle scriptlevel="0" displaystyle="true"> <mfrac> <mrow> <mi>c</mi> <mo>−</mo> <mi>v</mi> </mrow> <mrow> <mi>c</mi> </mrow> </mfrac> </mstyle> <mfenced open="|" close="|" separators="|"> <mrow> <mi>S</mi> <msub> <mrow> <mi>T</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> </mrow> </mfenced> </mrow> </semantics></math>. (3) <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> </mrow> </msub> </mrow> </semantics></math> is related to the elevation error, <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <mi>h</mi> </mrow> </semantics></math>. <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> </mrow> </msub> <mo>=</mo> <mi mathvariant="sans-serif">Δ</mi> <mi>h</mi> <mrow> <mrow> <mi mathvariant="italic">cos</mi> </mrow> <mo>⁡</mo> <mrow> <mi mathvariant="sans-serif">θ</mi> </mrow> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>The signal propagation model upon which the RT method is based. The propagation path in the troposphere is divided into multiple segments.</p>
Full article ">Figure 3
<p>The geolocation procedure.</p>
Full article ">Figure 4
<p>The geometry for solving the imaging time, <math display="inline"><semantics> <mrow> <mi>t</mi> </mrow> </semantics></math>: (<b>a</b>) single satellite model; (<b>b</b>) Tianhui-2 secondary satellite model.</p>
Full article ">Figure 5
<p>The relationship between the elevation-related geolocation error in range and on the ground. The range error is <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> </mrow> </msub> <mo>=</mo> <mi mathvariant="sans-serif">Δ</mi> <mi>h</mi> <mrow> <mrow> <mi mathvariant="italic">cos</mi> </mrow> <mo>⁡</mo> <mrow> <mi>θ</mi> </mrow> </mrow> </mrow> </semantics></math>, while on the ground it is <math display="inline"><semantics> <mrow> <mo stretchy="false">|</mo> <mi>O</mi> <mi>B</mi> <mo stretchy="false">|</mo> <mo>=</mo> <mi mathvariant="sans-serif">Δ</mi> <mi>h</mi> <mrow> <mrow> <mi mathvariant="italic">cot</mi> </mrow> <mo>⁡</mo> <mrow> <mi>θ</mi> </mrow> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>The geographical distribution of the measurement points.</p>
Full article ">Figure 7
<p>Geolocation results for two representative measurement points: (<b>a</b>) measurement points depicted in a Google Earth image; (<b>b</b>) TerraSAR-X image results without atmospheric correction; (<b>c</b>) TerraSAR-X image results with atmospheric correction. The red boxes indicate the measurement points, and the arrows denote the range direction.</p>
Full article ">Figure 8
<p>Offsets between the reference image and other images in (<b>a</b>) azimuth and (<b>b</b>) range. The <span class="html-italic">x</span>-axis represents the imaging time.</p>
Full article ">Figure 9
<p>Tianhui-2 image geolocation results corresponding to <a href="#remotesensing-16-04210-f007" class="html-fig">Figure 7</a>a with atmospheric correction: (<b>a</b>) reference image; (<b>b</b>) reference image with the overall shift; (<b>c</b>) secondary image; (<b>d</b>) secondary image with the overall shift.</p>
Full article ">Figure 10
<p>Comparison between the zenith delay calculated using the ERA5 and ZPD data: (<b>a</b>) difference changing with the day of the year; (<b>b</b>) statistical results of the original data.</p>
Full article ">Figure 11
<p>(<b>a</b>) DEM of the study area; (<b>b</b>) tropospheric delay on 17 January 2019, derived using the RT method; (<b>c</b>) difference between the ZDM result and the RT result; (<b>d</b>) relationship between the radio refractivity, <math display="inline"><semantics> <mrow> <mi>N</mi> </mrow> </semantics></math>, and the elevation in the study area. <math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mo>−</mo> <mn>3.1</mn> <mo>×</mo> <msup> <mrow> <mn>10</mn> </mrow> <mrow> <mo>−</mo> <mn>2</mn> </mrow> </msup> <mi>h</mi> <mo>+</mo> <mn>293</mn> </mrow> </semantics></math>, when the elevation is greater than 500 m.</p>
Full article ">Figure 12
<p>(<b>a</b>) Tropospheric delay calculated using the GACOS-based ZDM method; (<b>b</b>) GACOS-based ZDM result subtracting the ERA5-based ZDM result.</p>
Full article ">Figure 13
<p>(<b>a</b>) Relationship between the ground elevation and tropospheric delay; (<b>b</b>) kernel density estimation (KDE) curve of the tropospheric delay residual estimated using the regression equation.</p>
Full article ">Figure 14
<p>(<b>a</b>) Variation in the VTEC at 12 PM and 6 PM local time in the study area during 2019; (<b>b</b>) average VTEC at 6 PM local time in the study area in 2019. The red dot indicates the location of the study area.</p>
Full article ">Figure 15
<p>(<b>a</b>) KDE curves of the geolocation error, <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> </mrow> </msub> </mrow> </semantics></math>, based on the WorldDEM and SRTM; (<b>b</b>) cumulative percentage curves of the geolocation error, <math display="inline"><semantics> <mrow> <mfenced open="|" close="|" separators="|"> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> </mrow> </msub> </mrow> </mfenced> </mrow> </semantics></math>, when based on the two DEMs.</p>
Full article ">Figure 16
<p>(<b>a</b>) KDE curves of the geolocation error, <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> </mrow> </msub> </mrow> </semantics></math>, based on the WorldDEM. Measurement points were divided into four groups based on the slope, using threshold values of 0, 25%, 50%, and 75%; (<b>b</b>) KDE curves of <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <msub> <mrow> <mi>s</mi> </mrow> <mrow> <mi>e</mi> <mi>l</mi> <mi>e</mi> </mrow> </msub> </mrow> </semantics></math> based on the SRTM.</p>
Full article ">Figure 17
<p>Registration procedure using the ISCE. The DEM is projected to the radar center coordinate system as z data and is used to generate a simulated amplitude image. The orbit can be corrected using the rgshift calculated by registering the SAR image and the simulated image.</p>
Full article ">Figure 18
<p>The registration’s SNR: (<b>a</b>) TerraSAR-X image; (<b>b</b>) Tianhui-2 reference image.</p>
Full article ">
Back to TopTop