With the renewed interest for lunar surface exploration, the European Space Agency envisions to stimulate the creation of lunar communications and navigation services (LCNS) to enable, among others, autonomous navigation capabilities for lunar rovers. As the number of satellites foreseen in such a service is much smaller compared to Earth based global navigation satellite systems (GNSS), different complementary technologies are pursued to improve the attainable navigation accuracy for lunar rovers. One way to improve the position accuracy provided by the LCNS satellites is to constrain the vertical position using a high resolution digital elevation model (DEM). This article presents the results of a variance covariance analysis of an extended Kalman filter (EKF) implementation in which the LCNS ranging measurements are used together with DEM from the LRO LOLA instrument. Assuming a realistic orbit determination and time synchronization (ODTS) accuracy of the LCNS satellites, the usage of a navigation grade IMU and an oven controlled crystal oscillator (OCXO), a 3-sigma position accuracy of less than 10 meters can be obtained. Furthermore, the availability is substantially improved as the DEM aided solution enables a position solution in case of only 3 visible satellites.

The digital elevation model (DEM) is one of the key geospatial datasets used in many fields of engineering and science for countless applications. In this contribution, we assess the vertical accuracy of the Advanced Land Observing Satellite (ALOS) World 3D-30m (AW3D30) DEM using the runway method (RWYM). The RWYM utilizes the longitudinal profiles of runways which are reliable and ubiquitous reference data. A reference dataset used in this project consists of 36 runways located at various points throughout the world. The same dataset was previously used to test the accuracy of WorldDEM^TM. Our study indicates that AW3D30 has a remarkably high RMSE of 1.78 m (one σ). However, while analyzing the results, it has become apparent that it also contains a widespread elevation anomaly. We conclude that this anomaly is the result of uncompensated sensor noise and the data processing algorithm (downsampling of the higher resolution data). We believe that this issue should be communicated to the user community. Also, we would like to note that the traditional accuracy assessment of a DEM, e.g., statistical assessment of the elevation differences = model – reference, does not allow for identification of these type of anomalies in a DEM.

Marine vessels operating on the Arctic Sea route are constantly prone to collisions and friction with ice. This study discusses the wear of the hull plate caused by the collision of ice against vessels operating in Arctic Sea routes. The abrasive wear of the hull due to ice impact was numerically assessed based on both the incident behavior of ice particles interacting with the flow around the hull and the wear loss of the hull surface caused by the contact force of ice particles. A multi-phase approach was adopted to account for the behavior of ice particles continuously affected by the fluid force around the hull. The fluid force acting on the ice floe was evaluated using computational fluid dynamics (CFD) and the dynamic motion of the drift ice was evaluated using the discrete element method (DEM). The motion of the floating ice particles was updated in real time by iteratively coupling the fluid force and the motion of the ice floe at each time step of the numerical simulation. The results of the wear simulation models were presented in terms of the shape change of the hull surface due to wear. At first, the wear was evaluated for cases in which only the surface paint of the hull was damaged. Thereafter, a computation model considering the shape change of the hull surface experiencing long-term friction of ice particles was introduced. Finally, the numerical procedures to predict the abrasive wear of the hull surface by ice impact were discussed.

Digital Elevation Model (DEM) is a representation of elevation data that represent terrain with or without overlaying objects of the Earth. It is the ideal and most widely used method for determining topography. DEMs are generated from various techniques such as traditional Surveying, Photogrammetry, InSAR, LiDAR, Clinometry and radargrammetry. It has been observed that mostly LiDAR-generated DEMs provide the best accuracy. The unavailability of LiDAR data in most of the region restricts global researchers from high-resolution and accurate DEMs. The recent launch of ICESat-2 with a 13m beam footprint and 0.7m pulse interval, promises elevations at high orbital precision. Its accuracy is of the order of few centimeters in complex topography, because of this ICESat-2 proves to be a good source to generate high-accuracy DEMs. ICESat-2 provides discrete photon data with elevations of points on the Earth’s surface. Traditional interpolation techniques tend to over-smooth the estimated space and still are unable to justify the complicated continuity in the topographical data. Machine learning algorithms are widely being used to extract patterns and spatial extent in geographic data. To estimate a DEM from ICESat-2 LiDAR point data, machine learning regression algorithms are implemented in this study. The present study has been performed for a plain region of Ghaziabad, Uttar Pradesh, India. Studies have shown that Cartosat-1 DEM V3 R1 product provides an accuracy of the order of 2m in predominantly plain regions, hence taken for this region. Current work focuses on comparing various regression-based machine-learning techniques to interpolate DEM from ICESat-2 data The RMSE of the interpolated DEM resulted from the Gradient Boosting Regressor, Random Forest Regressor, Decision Tree Regressor, and Multi-Layer Perceptron Regressor was 7.13m, 7.01m, 7.15m, and 3.76m, respectively when evaluated against the TANDEM-X DEM of the same region. The MLP Regressor is found to perform the best among the four algorithms tested.

Interferometry Synthetic Aperture Radar (InSAR) is an advanced remote sensing technique for studying the earth's surface topography and deformations. It is used to generate high-quality Digital Elevation Models (DEMs). DEMs are a crucial and primary input to various topographical quantification and modelling applications. The quality of input DEMs can be further improved using fusion methods, which combine multi-sensor or multi-temporal datasets intelligently to retrieve the best information amongst the input data. This research study is based on developing a Neural Network based fusion approach for improving InSAR based DEMs in plain and hilly terrains. The study areas comprise of relatively plain terrain from Ghaziabad and hilly terrain of Dehradun and their surrounding regions. The training dataset consists of DEM elevations and derived topographic attributes like slope, aspect, topographic position index (TPI), terrain ruggedness index (TRI), and vector roughness measure (VRM) in different land use land cover classes of the study areas. The spaceborne altimetry ICESat-2 ATL08 photon data is used as a reference elevation. A Feed Forward Neural Network with backpropagation algorithm is trained based on the prepared training samples. The trained model produces fused DEMs by learning the relationship between the input and target samples. This is used to predict elevations in the test areas. The accuracy of results from the models are assessed with TanDEM-X 90 m DEM. The fused DEMs show significant improvement in terms of RMSE over the input DEMs with improvement factor of 94.65 % in plain area and 82.62 % in hilly area. The study concludes that the ANN with its universal approximation property is able to significantly improve the fused DEM.

The discrete element method (DEM) coupled with computational dynamics (CFD) has been considered one of the most sensitive ways of studying fluidized beds. This paper presents gas-solid fluidized bed simulations by use of CFD-DEM. The two-phase interaction is calculated in a way of decoupling. Within the effective neighborhood of the particles, the kernel approximation method is used to calculate the local porosity, thereby calculating the drag force on the particles. At the grid scale, the momentum exchange coefficient is calculated using a two fluid model based method, thereby calculating the source term of the phase interaction in the fluid control equation. The simulated form, size and motion process of the big bubble are all consistent with experimental observations. The simulated solid volume fraction and relative pressure averaged in a horizontal plane at 45 mm above the bottom are in reasonable agreement with the experimental data. The fluctuation of the simulated bed height over time is similar to the experimental measurement results. The reproduction of these qualitative and quantitative results indicates that the proposed method effectively improves simulation performance and accuracy.

This paper presents the results of simulations of particle-laden air-solid jet flow in long straight tubes using CFD-DEM, along with an analysis of coarse-graining. Although previous studies have used CFD-DEM for similar flows, these have typically been in a dilute flow regime where uncoupled simulations are accurate. However, fully coupled simulations can introduce issues, necessitating validation studies to ensure that all coupling parameters are effectively used and that the physics is accurately represented. This paper validates the simulations against two different experimental studies, with fluid Reynolds numbers between 10,000 and 40,000 and Stokes numbers between 5.6 and 50. Interestingly, the profiles of mean particle velocity exhibit fewer discrepancies as the Stokes number increases but more discrepancies for the root mean square velocity compared to the experiments. Particle number flux is consistent with the experiments after the nozzle exit. Coarse-graining is also applied to the same simulations, achieving relatively accurate results. However, as expected, scaling of contact collision frequencies, forces, and stresses cannot be achieved, meaning that coarse-graining may be useful for comparing designs or operating parameters on an industrial scale but falls short when measuring the total energy dissipation of one experiment.

This paper presents an analysis of linear viscous stress Favre-Averaged turbulence models for computational fluid dynamics (CFD) of fully turbulent round jets with a long straight tube geometry in the near field. Although similar work has been performed in the past with very relevant solutions, considerations were not given for the issues and limitations involved with coupling between an Eulerian and Lagrangian phase, such as in fully two-way coupled CFD-DEM. These issues include limitations on solution domain, mesh cell size, wall modelling, and momentum coupling between the two phases in relation to the particles size. Therefore, within these considerations, solutions are provided to the Navier-Stokes equations with various turbulence models using a three-dimensional wedge long straight tube geometry for fully developed turbulence flow. Simulations are performed with a Reynolds number of 15000 and 50000 using two different tube diameters. It is found that a modified k−ε turbulence model achieved the most agreeable results for both the velocity and turbulent flow fields between these two flow regimes, while a modified k−ω SST/BSL also provided suitable results.

Cloudburst is one of the most devastating and frequently occurring natural hazardous events in Indian Himalayan region. Localized deep cumulus convective clouds have a capability of giving enormous amount of rainfall over a limited horizontal area, within a short span of time. Whenever, such events occur, lead to flash floods causing landslides, house collapses, dislocation of traffic, and human casualties on a large scale. Therefore, it is necessary to predict the cloudburst inundation zones accurately to avoid damage associated with them. For this, high resolution Digital Elevation Model generated from CartoSat-1 (Stereo pair) were integrated in MIKE 11 Hydrodynamic 1D model to generate longitudinal profile of the study area and to find water level, peak discharge, flow velocity, flow width at different reaches along the Asi ganga and Bhagirathi river, to know the Cloudburst flood inundation scenario. On 3^rd August 2012 one of the major Cloudburst event occurred in Asi Ganga Valley in Indian Himalayan region which was considered for simulation of hydrodynamic model. For a Cloudburst event, 100 mm/hr rainfall was considered for the simulation of the hydrodynamic model. It is observed that the discharge rise from 50 m³/s to 549.164 m³/s (an abrupt increase of about 10 times) within 1 hr at Sangamchetty in Asiganga river and at Joshiyara area rise from 600 m³/s to 3378.69 m³/s (an abrupt increase of about 5 times) within 4 hr in Bhagirathi river. Similarly the water level rises around 3 m and 6m in Asi Ganga and Bhagirathi rivers respectively. Flash Flood inundation areas due to Cloudburst on 3^rd August 2012 were demarcated from the simulation results in GIS environment.

The TanDEM-X mission is acquiring a new dataset to provide a temporally independent DEM, called "TanDEM-X Change DEM". This set of acquisitions taken between 2017 and 2020 has a clear temporal separation to the TanDEM-X global DEM data which were acquired between 2010 and 2015. Therefore, this new DEM aims to enable the characterization of terrain changes. Improvements in the acquisition planning and the data processing were necessary to generate this Change DEM with fewer acquisitions but still very high accuracy. For this, the use of an edited TanDEM-X DEM as a "starting point" for the interferometric processing is mandatory.

Flooding, exacerbated by climate change, poses a significant threat to certain areas, increasing in frequency and severity. In response, the construction of supplementary dams has emerged as a reliable solution for flood management. This study employs a Remote Sensing (RS) approach integrated with GIS to assess the feasibility of constructing a supplementary dam near Linville, Brisbane, Australia, with the aim of mit-igating floods and preventing overtopping failure at the Wivenhoe Dam. Using QGIS software and a 25-meter resolution DEM from the Queensland Spatial Catalogue ‘QSpatial’ website, four potential dam sites were analysed, considering cross-sections, watershed characteristics, and water volume calculations. Systematic selection criteria were applied to identify the most optimal option based on dam wall di-mensions, volume-to-area, and volume-to-cost ratios. The selected option was further assessed against predefined criteria, yielding one optimal choice. The study provides insights into the feasibility and ef-fectiveness of supplementary dam construction for flood mitigation in the region, with recommendations for future research and implementation.

Tablets are the most used solid dosage form in medical treatments. The realistic modeling 1 of such specific object shapes in computer simulations has been demonstrated already by the author 2 for bi-convex round tablets but needs to be extended for a wider range of tablet shapes. Moreover, 3 the number of tablets must be large enough to simulate industrial drum coaters, which are typically 4 used for film coating processes. Whilst this long lasting process where tablets flow continuously, 5 tablets are loaded by frequent changes of individual contacts. Using a tailored simulation tool named 6 pharmAssist , written for GPUs (modern massively parallel processors in graphics cards and HPC 7 accelerators), we implemented the widely used oblong tablet shape. We are able to simulate flows 8 involving industrial scale number of such tablets. Data mining allows for statistical analyzis of 9 pair-wise contacts of different geometric features. Our simulation results serves as an observation 10 space, which cannot be addressed by experimental analyzis and will help to avoid degradation of 11 tablets like chipping or other defect phenomena. This paper reflects current development work in 12 order to extend the DEM with a new shape class, which comprises the modelling of compacted or 13 extruded objects, which have a smooth exterior, a flat band and sharp contour edges

Snow cover plays a critical role in the global energy and water cycles. Snow cover on the Tibetan Plateau (TP) provides vital water sources in western China and Himalayan regions in addition to its weather and climate significance. The massive high mountain topography of the TP is main conditions for the presence and persistence of snow cover on the plateau at the mid-low latitudes of the Northern Hemisphere (NH). However, how mountain topography controls snow cover distribution on the TP is largely remain unclear and the relationship is not well quantified. Here, the spatial distribution of snow cover and topographic controls on snow cover on the TP are examined based on snow cover frequency (SCF) derived from MODIS snow cover product (MOD10A2 v5) and digital elevation model (DEM). The results show that snow cover on the TP is spatially unevenly distributed and is characterized by rich snow and high SCF on the interior and surrounding high mountain ranges, and less snow and low SCF in inland basins and river valleys. Snow cover on the TP presents elevation dependence with the higher the altitude, the higher the SCF, the longer the snow cover duration and the more stable the intra-annual variation. Annual mean SCF below 3000 m above sea level (masl) is less than 4% and it reaches 77% above 6000 masl. The intra-annual snow cover variation below 4000 masl features a unimodal distribution, while above 4000 masl it presents a bimodal distribution. The mean minimum SCF below 6000 masl occurs in summer, while above 6000 masl it occurs in winter. Because of differences in solar radiation and water vapor sources, mean SCF generally increases with mountain slopes and it is the highest on the north-facing aspect, whereas the lowest is observed on the south-facing aspect.

Digital Elevation Models (DEMs) are widely used in geographic and environmental studies. In the current work, the fusion of multi-source DEMs is investigated to improve the overall accuracy of public domain DEMs. Multi-scale decomposition is an important analytical method in data fusion. Three multi-scale decomposition methods – the wavelet transform (WT), bidimensional empirical mode decomposition (BEMD) and nonlinear adaptive multi-scale decomposition (N-AMD) - are applied to the 1-arc-second Shuttle Radar Topography Mission Global digital elevation model (SRTM-1 DEM) and the Advanced Land Observing Satellite World 3D – 30 m digital surface model (AW3D30 DSM) in China. Of these, the WT and BEMD are popular image fusion methods. A new approach for DEM fusion is developed using N-AMD (which is originally invented to remove the cycle from sunspots). Subsequently, a window-based rule is proposed for the fusion of corresponding frequency components obtained by these methods. Quantitative results show that N-AMD is more suitable for multi-scale fusion of multi-source DEMs, taking the ice cloud and land elevation satellite (ICESat) global land surface altimetry data as a reference. The vertical accuracy of the fused DEM shows significant improvements of 29.6% and 19.3% in a mountainous region and 27.4% and 15.5% in a low-relief region, compared to the SRTM-1 and AW3D30 respectively. Furthermore, a slope position-based linear regression method is developed to calibrate the fused DEM for different slope position classes, by investigating the distribution of the fused DEM error with topography. The results indicate that the accuracy of the DEM calibrated by this method is improved by 16% and 13.6%, compared to the fused DEM in the mountainous region and low-relief region respectively, proving that it is a practical and simple means of further increasing the accuracy of the fused DEM.

Accurate seasonal weather forecasts are vital across a broad spectrum of applications, bearing significant environmental and socioeconomic implications. This importance renders the subject a matter of primary interest to a wide range of stakeholders, including the general public, agricultural sector, emergency responders, financial institutions, and policy strategists. The need for precision in long-term predictive models necessitates the development of innovative methodologies. These methodologies should be capable of deciphering atmospheric patterns and mechanisms with detail, especially at a local level, without the resource constraints that dynamical downscaling imposes. In response to this expanding demand, this study presents a novel solution, combining a stochastic deep learning methodology, specifically Generative Adversarial Networks (GANs), with a Digital Elevation Model (DEM). The cornerstone of the proposed model is the enhancement of gridded temperature fields from seasonal forecasts. The area of focus was the Hellenic region, wherein the spatial resolution is amplified from a coarse 1° x 1° grid to an impressively detailed 0.1° x 0.1° grid. This offers a transformative perspective for interpreting and employing this crucial meteorological data. The results suggest that the downscaled fields adequately approximate the actual spatial distribution, although the predicted values tend to slightly overestimate and in some cases underestimate the original ones. This study underlines the potential of this approach to significantly enhance the resolution and utility of weather forecasts, thereby contributing to a variety of sectors dependent on reliable meteorological data.

The fluid flow field at the upper airways is highly complex due to the complex structure of the airway. The inhaled particle flow, the air streamline and the interaction of the continuum and discrete phase could significantly affect the transport behaviour of the inhaled particles. A range of analytical, mathematical and computational fluid dynamics (CFD) models analyzed the airflow and particle transport in different idealized and asymmetric airway models. A precise understanding of the continuum and discrete phase interaction in realistic human airways is missing, and this study aims to develop a CFD-DEM model for particle transport in realistic airways. This study uses the CFD model for the continuum phase and the discrete element method (DEM) for the discrete phase. A soft sphere approach is used for the interaction of the discrete phase. Proper validation is performed for particle transport efficiency. The CFD-DEM model analyzed the particle transport in an idealized and realistic airway model, and different methods are used to analyze the transport behaviour. During the particle-particle interaction, a stagnation point and a high-pressure zone are observed at the airway model's carinal angle. The numerical results report higher deposition efficiency (DE) for particle-particle interaction than without interaction. The flow field becomes highly complex with the spring constant values, and higher DE is found for high spring constant values. The spring dashpot friction-dshf method shows higher deposition at the upper part of the airways than other interaction methods. The findings of this study and more case-specific analysis would improve the knowledge of aerosol transport in airways and the health risk assessment of the patient.

In this case study, an active runway of a civilian airport in Zonguldak, Turkey, is used to assess the suitability of spaceborne digital elevation models (DEMs) to model an anthropogenic structure. The tested DEMs include the ASTER, the AW3D30 m, the SRTM-1”, the SRTM-3”, the SRTM-X, the TanDEM-3”, and the WorldDEMTM. A photogrammetric high accuracy DEM was also available for the tests. As a reference dataset, a line leveling survey of the runway using a Leica Sprinter 150/150M instrument was performed. The selection of a runway as a testbed for this type of investigation is justified by its unique characteristics, including its flat surface, homogenous surface material, and availability for a ground survey. These characteristics are significant because DEMs over similar structures are free from environment- and target-induced error sources. The most accurate DEM is the WorldDEMTM followed by the SRTM-3” and TanDEM-3”, with vertical errors (LE90) equal to 1.291 m, 1.542 m, and 1.56 m, respectively. This investigation uses a method for identifying the vertical errors in DEMs that is known as the runway method.

Digital Elevation Models (DEMs) are essential in watershed delineation, but the sensitivity of simulated runoff to DEM resolution is poorly understood. This study investigates the impact of DEM resolution on topological attributes and simulated runoff in the Mahabad Dam watershed, Iran. To delineate the watershed, DEMs with 12.5 m, 30 m, and 90 m resolutions were acquired from the ALOS PALSAR, Space Shuttle Radar Topography Mission (SRTM), and ASTER global DEM data source, respectively. Watershed and streamlines were delineated in ArcGIS, with hydrologic analyses performed using the Soil and Water Assessment Tool (SWAT). Sensitivity analysis on parameters contributing to runoff was done using Sequential Uncertainties Fitting Ver-2 (SUFI-2) Algorithm, in SWAT Calibration and Uncertainty Procedures (SWAT-CUP) software. Results showed the watershed area, reach lengths, and elevations in the watershed varied due to DEM resolutions. Higher amounts of runoff were generated when DEMs with finer resolutions were implemented. The 12.5 m DEM generated 3.48% and 0.42% more runoff compared with 90 m and 30 m DEMs, respectively. SWAT-CUP results showed the sensitivity of parameters contributing to runoff changes under different DEM resolutions. Regardless of DEM resolution, surface properties, available water capacity, and moisture levels in the soil are the most sensitive parameters. As the distribution of slope changes in different DEM resolutions, surface parameters are most affected. The findings indicate to reduce computation time and speed up computation procedures, researchers may use DEMs with coarser resolutions at the expense of minor decreases in accuracy.

In recent years discrete element method (DEM) has gradually been applied to the traditional fluidization simulation of fine particles in micro fluidized bed (MFB). The application of DEM in simulating fast fluidization of fine particles in MFB has not yet received attention. This article presents a drag model that relies on the surrounding environment of particles, namely particle circumstance-dependent drag modle or PCDD model. Fast fluidization in MFB of fine particles is simulated using DEM based on the PCDD model. It is found that the gas-solid two-phase flow in MFB meets the common law of fast fluidization in circulating fluidized bed (CFB) to a certain extent, and also presents a special law that is significantly different from that in CFB. Simulations indicate that the local structure in MFB exhibits particle aggregation which is a natural property of fast fluidized, forming a structure where continuous dilute phase and dispersed concentrated phase coexist. The formation and fragmentation of particle clusters in different local regions of the bed have time synchronization. There exists strong effect of solid back-mixing in MFB, leading to relatively low outlet solid flux. The gas back-mixing effect is, however, not so distinct. The axial porosity shows a monotonically increasing distribution with the bed height, but does not strictly follow the single exponential distribution. The solid volume fraction at the bottom of the bed is significantly lower than the correlated value in CFB. The axial heterogeneous distribution of the cross-sectional average porosity in the lower half of the bed is also weakened. The radial porosity shows a distribution pattern of higher in the central region and lower in the sidewall region. Compared with the correlation results in CFB, the porosity near the central region in MFB is relatively low, while the porosity near the sidewall region is relatively high. The geometric size of the container, especially the radial size, is much smaller than that of a CFB. The relative area of contact between the MFB wall and particles is larger, and the wall friction factor becomes more significant. This may be the main reason for the time synchronization of particle agglomeration, severe particle back-mixing, low outlet solid flux, monotonically increasing axial porosity and weak core-annular structure of radial porosity.

The purpose of this work was analysing the compaction of a cohesive material using different DEM simulators to determine the equivalent contact models and identify how some parameters of the simulations affect the compaction results (maximum force and compacts appearance) and computational costs. For that purpose, three cohesion contact models were tested (‘linear cohesion’ in EDEM; ‘SJKR’ and ‘SJKR2’ in LIGGGHTS). The influence of the particle size distribution (PSD) on the results was also investigated. Further assessments were performed on the effect of selecting different timesteps, using distinct conversion tolerances for exporting the 3D models to STL files and moving the punch with different speeds. Consequently, it was possible to determine that a timestep equal to a 10% Rayleigh timestep, a conversion tolerance of 0.01 mm and a punch speed of 0.2 m/s are adequate for simulating the compaction process using the contact models in this work. In addition, the results determined that the maximum force was influenced by the PSD because of the rearrangement of the particles. The PSD was also related to the computational cost because of the number of simulated particles and their sizes. Finally, an equivalence was found between the linear cohesion and SJKR2 contact models

Over the last decades, an important number of regional and global Digital Elevation Models (DEMs) have been released publicly. As a consequence, researchers need to choose between several of these models to perform their studies and to use these DEMs as third-party data to compute derived products (e.g., for orthorectification). However, the comparison of DEMs is not trivial. For most quantitative comparisons, DEMs require to be expressed in the same Coordinate Reference System (CRS), sampled over the same grid (i.e., be at the same ground sampling distance with the same pixel-is-area or pixel-is-point convention), and with heights relative to the same Vertical Reference System (VRS). Thankfully, many open tools allow us to perform these transformations precisely and easily. Despite these rigorous transformations, local or global planimetric displacements may still be observed from one DEM to another. These displacements, or disparities, may lead to significant biases in comparisons of DEM elevations or derived products such as slope, aspect or curvature. Therefore, the control of DEM planimetric accuracy is certainly the most important task to perform at first before any comparison. This paper presents the disparity analysis method enhanced to achieve a sub-pixel accuracy, by interpolating the linear regression coefficients computed within an exploration window. This new method is significantly faster than oversampling the input data because it uses the correlation coefficients that have been already computed in the disparity analysis. To demonstrate the robustness of this algorithm, artificial displacements have been performed in a 11x11 grid with a 0.1-pixel step in both directions. The Copernicus DEM GLO-30 has been interpolated by a “classical” bicubic to achieve these 11x11 displacements. This validation method has been applied in four 10 km x 10 km DEMIX tiles showing different height distributions. Globally, this new sub-pixel accuracy method is robust. Artificial displacements have been retrieved with with typical errors (e_b) ranging from 12 to 20 % of the pixel size (worst case in Croatia).These errors in displacement retrievals are not constant in the 11x11 grid and the overall error E_b depends on the roughness (altimetry variations) encountered in the different tiles. The second idea of this paper is to assess the impact of the bicubic parameter (slope of the weight function at a distance d=1 of the interpolated point) on the accuracy of the displacement retrieval. By considering E_b as a quality indicator, tests have been performed in the four DEMIX tiles making the bicubic parameter vary between -1.5 and 0.0 by step of 0.1. For each DEMIX tile, the best bicubic (BBC) parameter b* is interpolated from the four E_b minimal values. This BBC parameter b* is low for flat areas (around -0.95) and higher in mountainous areas (around -0.75). The roughness indicator is the standard deviation of the slope norms computed from all the pixels of a tile. A logarithmic regression analysis performed between the roughness indicator and the BBC parameter b* computed in 67 DEMIX tiles shows a high correlation (r=0.717). The logarithmic regression formula estimating the BBC parameter from the roughness indicator is generic and may be applied to estimate the displacements between two different DEMs. This formula may also be used to set up a future Adaptative Best BiCubic (ABBC) that will estimate the local roughness in a sliding window to compute a local BBC.

Seepage erosion is one of the main reasons for the local collapse or instability of embankments. To investigate the characteristics and mechanism of seepage erosion for cohesionless soils, the model tests by using an independently developed seepage erosion device and the numerical simulations based on the Discrete Element Method-Computational Fluid Dynamics (DEM-CFD) coupling model were carried out. The results show that the seepage erosion process of the cohesionless soil can be characterized by four stages: stable seepage stage, fine particles upward migration stage, sand samples boiling stage, and erosion damage stage. The skeleton structure of soil sample under seepage flow is continually changed due to the loss of fine soil particles, which results in a significant decrease in the sample strength and may ultimately lead to the failure of the sample. The results of this study can provide references and bases for the design, construction, and long-term service of embankments or earth dams under the complex seepage condition, reducing the risk of seepage erosion.

Archaeological topography identification from high-resolution DEMs is a current method that is used with high success in archaeological prospecting of wide areas. I present a methodology trough which burial mounds (tumuli) from LiDAR DEMS can be identified. This methodology uses geomorphometric and statistical methods to identify with high accuracy burial mound candidates. Peaks, defined as local elevation maxima are found as a first step. In the second step, local convexity watershed segments and their seeds are compared with positions of local peaks and the peaks that correspond or have in vicinity local convexity segments seeds are selected. The local convexity segments that correspond to these selected peaks are further feed to a Random Forest algorithm together with shape descriptors and descriptive statistics of geomorphometric variables in order to build a model for the classification. Multiple approaches to tune and selected the proper training dataset, settings and variables were tested. The validation of the model was performed on the full dataset where the training was performed and on an external dataset in order to test the usability of the method for other areas in a similar geomorphological and archaeological setting. The validation was performed against manually mapped and field checked burial mounds from two neighbor study areas of 100 km² each. The results show that by training the Random Forest on a dataset composed of between 75% to 100% of the segments corresponding to burial mounds and ten times more non-burial mounds segments selected using latin hypercube sampling, 93% of the burial mound segments from the external dataset are identified. There are 42 false positive cases that need to be checked, and there are two burial mound segments missed. The method shows great promise to be used for burial mound detection on wider areas by delineating a certain number of tumuli mounds for model training.

This paper addresses the study of the stress field in composites continua with the multiscale approach of the DECM (Discrete Element modeling with the Cell Method). The analysis focuses on composites consisting of a matrix with inclusions of various shapes, to investigate whether and how the shape of the inclusions changes the stress field. The purpose is to provide a numerical explanation for some of the main failure mechanisms of concrete, which is precisely a composite consisting of a cement-based matrix and aggregates of various shapes. Actually, while extensive experimental campaigns detailed the shape effect of concrete aggregates in the past, so far it has not been possible to model the stress field within the inclusions and on the interfaces accurately. The reason for this lies in the limits of the differential formulation, which is the basis of the most commonly used numerical methods. The Cell Method (CM), on the contrary, is an algebraic method that provides descriptions up to the micro-scale, independently of the presence of rheological discontinuities or concentrated sources. This makes the CM useful for describing the shape effect of the inclusions, on the micro-scale. When used together with a multiscale approach, it also models the macro-scale behavior of periodic composite continua, without losing accuracy on the micro-scale. The DECM uses discrete elements precisely to provide the CM with a multiscale approach.

Mostly populous city like Chennai is subjected to frequent flooding due to its complex nature of natural and man-made activities. From the analysis of the past records of flood events of 1943,1976,1985,2005 and 2008,it has been observed Adayar watershed is subjected to cataclysmic flooding in low-lying areas of the city and its suburbs because of inoperativeness of the local drainage system, rainfall associated with cyclonic activity, topography of the terrain, encroachments along the floodplain, hugh upstream flow discharge into the river and the highly impervious area which blocked the runoff to flow into the storm water drainage. After looking into these problems of flooding, a study have been conducted on Adayar watershed to develop a 2D hydrodynamic model for the two scenarios of existing condition of storm water drainage network and revised conditions of storm water drainage network using high resolution Lidar DEM to assess the volume of runoff with respect to time and duration on flood peaks for the two flood events of 2005 and 2015.Secondly to develop a 1D flood model to predict the river stages during peak floods using MIKE 11 for the Adayar watershed. Thirdly to integrate the coupled 1D and 2D model using MIKEFLOOD for assessing the extent of inundation in the floodplain area of Adayar river. Finally results from the integrated model have been validated and the results found satisfactory. As a part of mitigation measures, two flood mitigation measures have been adopted. One measure such as revised storm water drainage system which enhances the flood carrying capacity of the drains and results in less inundated area which solves the problem of urban flooding and second measure such as regrading the river bed which reduces the floodplain inundation around the adjacent area of the river. After adopting these measures, the river is free to flow into the sea without any blockades.

This paper presents a new numerical method for multiscale modeling of composite materials. The new numerical model, called DECM, consists in a DEM (Discrete Element Method) approach of the Cell Method (CM) and combines the main features of both the DEM and the CM. In particular, it offers the same degree of detail as the CM, on the microscale, and manages the discrete elements individually such as the DEM—allowing finite displacements and rotations—on the macroscale. Moreover, the DECM is able to activate crack propagation until complete detachment and automatically recognizes new contacts. Unlike other DEM approaches for modeling failure mechanisms in continuous media, the DECM does not require prior knowledge of the failure position. Furthermore, the DECM solves the problems in the space domain directly. Therefore, it does not require any dynamic relaxation techniques to obtain the static solution. For the sake of example, the paper shows the results offered by the DECM for axial and shear loading of a composite two-dimensional domain with periodic round inclusions. The paper also offers some insights into how the inclusions modify the stress field in composite continua.

Collisional mountain belts commonly develop intramontane basins from mechanical and isostatic subsidence during orogenic development. These frequently display a relict top surface, evidencing a change interval from basin infilling to erosion often via capture or overspill. Such surfaces provide markers that inform on orogenic growth patterns via climate and base level interplay. Here, we describe the top surface from the Sorbas Basin, a key intramontane basin within the Betic Cordillera (SE Spain). The surface is fragmentary comprising high elevation hilltops and discontinuous ridges developed onto the variably deformed final basin infill outcrop (Gochar Formation). We reconstruct surface configuration using DEM interpolation and apply ¹⁰Be/²⁶Al cosmonuclides to assess surface formation timing. The surface is a degraded Early Pleistocene erosional pediment developed via autogenic switching of alluvial fan streams under stable dryland climate and base level conditions. Base level lowering since the Middle Pleistocene focused headwards incision up interfan drainages, culminating in fan head capture and fan morphological preservation within the abandoned surface. Post abandonment erosion has lowered the basin surface by 31 m (average) and removed ~5.95 km³ of fill. Regional basin comparisons reveal a phase of Early Pleistocene surface formation, marking landscape stability following the most recent Pliocene-Early Pleistocene mountain building. Post-surface erosion rate quantification is low and in accordance with ¹⁰Be denudation rates typical of the low uplift Betic Cordillera.

The result of the multidisciplinary collaboration of researchers from different areas of knowledge to validate a solar radiation model is presented. The MAPsol is a 3D local-scale adaptive solar radiation model that allows to estimate direct, diffuse and reflected irradiance for clear sky conditions. The model includes the adaptation of the mesh to complex orography and albedo, and considers the shadows cast by the terrain and buildings. The surface mesh generation is based on surface refinement, smoothing and parameterisation techniques and allows the generation of high quality adapted meshes with a reasonable number of elements. Another key aspect of the paper is the generation of a high-resolution digital elevation model (DEM). This high-resolution DEM is constructed from LiDAR data and its resolution if two times more accurate than the publicly available DEMs. The validation process uses direct and global solar irradiance data obtained from pyranometers at the University of Salamanca located in an urban area affected by systematic shading from nearby buildings. This work provides an efficient protocol for studying solar resources, with particular emphasis on areas of complex orography and dense buildings where shadows can potentially make solar energy production facilities less efficient.

For Interferometry Synthetic Aperture Radar (InSAR), the normal baseline is one of the main factors that affect the accuracy of the ground elevation. For Gaofen-3 (GF-3) InSAR processing, the poor accuracy of the real-time orbit determination resulting in a large baseline error, leads to the modulation error in azimuth and the slope error in range for timely Digital Elevation Model (DEM) generation. In order to address this problem, a baseline estimation method based on external DEM is proposed in this paper. Firstly, according to the characteristic of the real-time orbit of GF-3 images, orbit fitting is executed to remove the non-linear error factor. Secondly, the height errors are obtained in slant-range plane between Shuttle Radar Topography Mission (SRTM) DEM and the GF-3 generated DEM after orbit fitting. At the same time, the height errors are used to estimate the baseline error which has a linear variation. In this way, the orbit error can be calibrated by the estimated baseline error. Finally, DEM generation is performed by using the modified baseline and orbit. This procedure is implemented iteratively to achieve a higher accuracy DEM. Based on the results of GF-3 interferometric SAR data for Hebei, the effectiveness of the proposed algorithm is verified and the accuracy of GF-3 real-time DEM products can be improved extensively.

Abstract：Ballast pick-up velocity is a crucial parameter in the design of a slurry system in excavating devices, it influences transportation stability, safety, and economy of the system. In this study, a vertical pipe hydraulic experiment with different ballast sizes under water condition was conducted to investigate the ballast pick-up velocity in the shield slurry system. Thereafter, the coupling method of computational fluid dynamics and the discrete element method (CFD-DEM) was employed to establish a numerical model of the ballast pick-up velocity in the shield slurry system, and its applicability was verified by comparing the experiment results. Additionally, the ballast pick-up process under actual working conditions was studied, and the influence of ballast size and concentration on pick-up velocity was investigated. Overall, the primary findings of the study are as follows: 1) The numerical method can be applied to predict both the single ballast and ballast group pick-up velocity. 2) The pick-up process of the ballast group can be divided into four stages. 3) The conveying velocity of the ballast should be higher than the pick-up velocity of the single ballast to facilitate smooth discharge. 4) For a certain increase in ballast concentration, the pick-up velocity tends to be stable.

Craters are the most prominent geomorphological features on the surface of celestial bodies, which is playing a crucial role in studying the formation and evolution of celestial bodies, as well as in landing and planning for surface exploration. Currently, the main automatic crater detection models and datasets focus on the detection of large and medium craters. In this paper, we created 23 small lunar crater datasets for model training based on the Chang’E-2 (CE-2) DOM, DEM, Slope, and integrated data with 7 kinds of visualization stretching methods. And then, we proposed the YOLO-Crater model for Lunar and Martian small crater detection by replacing EioU and VariFocal loss to solve the crater samples imbalance problem and introducing a CBAM attention mechanism to mitigate interference from the complex extraterrestrial environment. The results show that the accuracy (P = 87.86%, R = 66.04%, and F1 = 75.41%) of the Lunar YOLO-Crater model based on the DOM-MMS (Maximum-Minimum Stretching) dataset is the highest and better than that of YOLOX model. And the Martian YOLO-Crater, trained by the Martian dataset from the 2022 GeoAI Martian Challenge, achieves good performance with P = 88.37%, R = 69.25%, and F1 = 77.65%. It indicates that the YOLO-Crater model has a strong transferability and generalization capability, which can be applied to detect small craters on the Moon and other celestial bodies.

Structure from Motion with Multi-View Stereo photogrammetry (SfM) is increasingly utilised in geoscience investigations as a cost-effective method of acquiring high resolution (sub-meter) topographic data, but has not been thoroughly tested in gullied savanna systems. The aim of this study was to test the accuracy of topographic models derived from aerial (via an Unmanned Aerial Vehicle, ‘UAV’) and ground-based (via a handheld digital camera, ‘Ground’) SfM in modelling a hillslope gully system in dry-tropical savanna, and to assess the strengths and limitations of the approach at different scales. A UAV survey covered an entire hillslope gully system (0.715 km²), whereas a Ground survey covered a single gully within the broader system (650 m²). SfM topographic models, including Digital Surface Models (DSM) and dense point clouds, were compared against RTK-GPS point data and a pre-existing airborne LiDAR Digital Elevation Model (DEM). Results indicate UAV SfM can deliver topographic models with a resolution and accuracy suitable to define gully systems at a hillslope scale (e.g., 0.1 m resolution with ~ 0.5 – 1.3 m elevation error), while ground-based SfM is more capable of quantifying gully morphology (e.g., 0.01 m resolution with ~ 0.1 m elevation error). Key strengths of SfM for these applications include: the production of high resolution 3D topographic models and ortho-photo mosaics, low survey instrument costs (< $AUD 3,000); and rapid survey time (4 and 2 hours for UAV and Ground survey respectively). Current limitations of SfM include: difficulties in reconstructing vegetated surfaces; uncertainty as to optimal survey and processing designs; and high computational demands. Overall, this study has demonstrated great potential for SfM to be used as a cost-effective tool to aid in the mapping, modelling and management of hillslope gully systems at different scales, in tropical savanna landscapes and elsewhere.

Autonomous unmanned aerial vehicles (UAVs) require highly reliable navigation information. Generally, navigation systems with the inertial navigation system (INS) and global navigation satellite system (GNSS) have been widely used. However, the GNSS is vulnerable to jamming and spoofing. The terrain referenced navigation (TRN) technique can be used to solve this problem. In this study, to obtain reliable navigation information even if a GNSS is not available or the degree of terrain roughness is not determined, we propose a federated filter based INS/GNSS/TRN integrated navigation system. we also introduce a TRN system that combines batch processing and an auxiliary particle filter to ensure stable flight of UAVs even in a long-term GNSS-denied environment. As an altimeter sensor for the TRN system, we use an interferometric radar altimeter (IRA) to obtain reliable navigation accuracy in high altitude flight. In addition, a parallel computing technique with general-purpose computing on graphics processing units (GPGPU) is applied to process a high resolution terrain database and a nonlinear filter in real time on board. Finally, we verify the performance of the proposed system through software-in-the-loop (SIL) tests and captive flight tests in a GNSS unavailable environment.

Abstract：Objective: To compare the effects of fermentation on the chemical composition of Gastrodia Tuder Halimasch Powder, and to provide a basis for the processing, processing and clinical application of this medicinal material. Methods: UHPLC-Q-Orbitrap HRMS was used to identify compound structure by comparing excimer ions, characteristic fragment ions and reference substances, and referring to relevant literature and database search. The main differential components were screened by orthogonal partial least squares-discriminant analysis (OPLS-DA), and the differential components were quantitatively studied by high performance liquid chromatography, and the differences in the types and contents of chemical components before and after fermentation were compared. Results: A total of 77 compounds were identified or preliminarily derived after fermentation, including 21 terpenes, 18 organic acids, 9 flavonoids, 7 nucleosides, 6 amides, 4 amino acids, 2 pyrrolidone derivatives, 2 sterols, 2 alkaloids and 6 others.Through high performance liquid chromatography, it was found that the content of daidzein genistein and ergosterol was significantly increased after fermentation, and a new component ergosterol was produced. Compared with Gastrodia Tuder Halimasch Powder before fermentation, the content of genistin and other components was significantly reduced, and there may be a process of genistin conversion to genistein.Conclusion: The fermentation process will have a certain influence on the types and contents of the chemical components of Gastrodia Tuder Halimasch Powder, and ergosterol, genistein and .daidzein may be the main pharmacodynamic components, which can provide a useful reference for the clinical rational application, quality control and mechanism of action of Gastrodia Tuder Halimasch Powder.

This study examines the activities conducted at the archaeological site of Aptera in Crete, Greece. The research was part of the DIACHRONIC LANDSCAPES International Design Workshop, organized by the CAM (Center for Mediterranean Architecture), TUC (Technological University of Crete School of Architecture), and UNIFE (University of Ferrara Department of Architecture). This article outlines the methods used for data acquisition and processing on a territorial scale, which generated digital outputs necessary for the analysis and design phases of the workshop, as well as for further examination of the results. The collected data, obtained through low-cost aerial photogrammetric surveying and GNSS terrestrial coordinate detection, were integrated in a Structure from Motion workflow that led to the creation and exportation of various digital outputs, such as point clouds, DTM, DSM, orthophotos, and contour lines. An accuracy analysis was performed to evaluate the effectiveness and efficiency of the digital models compared to the implemented surveying strategies, including the Ground Control Point and Quality Check Point marker positioning strategy. The resulting digital models proved to be valuable assets for analysis and design within the workshop and provided insightful prospects for future research and territorial-scale projects.