Discrete distribution estimation is a fundamental statistical tool, which is widely used to perform data analysis tasks in various applications involving sensitive personal information. Due to privacy concerns, individuals may not always provide their raw information, which leads to unpredictable biases in the final results of estimated distribution. Local Differential Privacy (LDP) is an advanced technique for privacy protection of discrete distribution estimation. Currently, typical LDP mechanisms provide same protection for all items in the domain, which imposes unnecessary perturbation on less sensitive items and thus degrades the utility of final results. Although, several recent works try to alleviate this problem, the utility can be further improved. In this paper, we propose a novel notion called Item-Oriented Personalized LDP (IPLDP), which independently perturbs different items with different privacy budgets to achieve personalized privacy protection. Furthermore, to satisfy IPLDP, we propose the Item-Oriented Personalized Randomized Response (IPRR) based on the observation that the sensitivity of data shows an inverse relationship with the population size of respective individuals. Theoretical analysis and experimental results demonstrate that our method can provide fine-grained privacy protection and improve data utility simultaneously.

With the increasing concern of environment, the energy-efficiency scheduling of manufacturing industry is becoming urgent and popular. In turning processes, both spindle speed and processing time can affect the final energy consumption and thus the spindle speed and scheduling scheme need to be optimized simultaneously. Since the turning workshop can be regarded as the flexible flow shop, this paper formulates a mixed integer linear programming model for energy-efficient scheduling of flexible flow shop. Accordingly, a new decoding method is developed by considering of the optimization of spindle speed and scheduling scheme simultaneously, and an estimation of distribution algorithm adopting the new decoding method is proposed to solve large-size problems. The parameters of this algorithm are determined by statistical technique with a simplified practical case. In order to validate the proposed method, a case from practical factories is studied, in which the makespan can be shortened by 25.22%, and the consumed energy can be saved by 5.48%. These results demonstrate the effectiveness of the proposed mathematical model and algorithm.

The power function distribution is a flexible waiting time model that may provide better fit for some failure data. This paper presents the comparison of the maximum likelihood estimates and the Bayes estimates of two-parameter power function distribution. The Bayes estimates are obtained, using conjugate priors, under five loss functions consist of square error, precautionary, weighted, LINEX and DeGroot loss function. The Gibbs sampling algorithm is proposed to generate samples from posterior distributions and in result the Bayes estimates are computed. The comparison of the maximum likelihood estimates and the Bayes estimates are done through the root mean squared errors. One real-life data set is analyzed to illustrate the evaluation of proposed methods of estimation. Finally, results from the simulation are discussed to assess the performance behavior of the maximum likelihood estimates and the Bayes estimates.

Establishing global high-rate secure communications is a potential application of continuous-variable quantum key distribution (CVQKD) but also challenging for long-distance transmissions in metropolitan areas. The discrete modulation(DM) can make up for the shortage of transmission distance that has a unique advantage against all side-channel attacks, however its further performance improvement requires source preparation in the presence of noise and loss. Here, we consider the effects of photon catalysis (PC) on the DM-involved source preparation for lengthening the maximal transmission distance of the CVQKD system. We address a zero-photon catalysis (ZPC)-based source preparation for enhancing the DM-CVQKD system. The statistical fluctuation due to the finite length of data is taken into account for the practical security analysis. Numerical simulations show that the ZPC-based DM-CVQKD system can not only achieve the extended maximal transmission distance, but also contributes to the reasonable increase of the secret key rate. This approach enables the DM-CVQKD to tolerate lower reconciliation efficiency, which may promote the practical implementation solutions compatible with classical optical communications using state-of-the-art technology.

Forest ecosystems make a greater contribution to carbon (C) stocks than any other terrestrial ecosystem. To understand the role of regional forest ecosystems in global climate change and carbon exchange, forest C stock and its spatial distribution within the small (2,300 km2) Liuxihe River basin were analyzed to determine the different contributors to the C stock. Forest C stocks were quantified by measuring the biomass of trees, understory vegetation, litter and roots, as well as soil organic C, using data from field samples and laboratory experiments. The results showed that forests stored 38.04 Tg C in the entire basin, with secondary and planted forests accounting for 89.82% and 10.18%, respectively, of the stored C. Five types of forests, a subtropical evergreen broad-leaved forest, a subtropical coniferous and broad-leaved mixed forest, a subtropical coniferous forest, a timber forest, and a non-wood forest, stored 257.55 ± 15.01, 218.92 ± 9.59, 195.24 ± 18.29, 177.42 ± 17.55, and 117.86 ± 6.04 Mg C ha−1, respectively. In the forest ecosystem C stocks of the basin, soils averagely contribute about 73.78%, not including root underground biomass. It provides a comprehensive method for forest ecosystem carbon investigation and forest management in small basin scale.

In recent decades, there have been numerous endeavors to develop a novel category of survival distributions possessing enhanced flexibility through the extension of existing distributions. This article constructs and validates the statistical properties of a novel survival distribution in order to obtain an alternative distribution that is suitable for analyzing survival data by presenting the novel mixture of the Fréchet distribution along with statistical properties such as the probability density function (PDF), cumulative distribution function (CDF), rth ordinary moment, skewness, kurtosis, moment-generating function, mean, variance, mode, survival function, hazard function, and asymptotic behavior, as well as constructing the estimators of the unknown parameter by employing the expectation-maximization (EM) algorithm, and simulated annealing. Additionally, the performance of the proposed estimators was compared with bias, mean squared errors (MSE), and simulated variances, and given an illustrative example of the proposed distribution to the survival data set in order to show that the proposed distribution is appropriate for the right-skewed data. This will be extremely advantageous in survival analysis.

To assure the opacity, an optimal problem is to retain the occurring events sequences as many as possible. Contrary to such problem, the other optimal problem is to preserve the minimal occurring events sequences. In the paper, based on the choosing cost, an optimal opacity-enforcing model is presented to minimize the sum of discount choosing cost, which is not only to preserve the opacity, but also to retain the secret to the maximum. To solve the model, two scenarios on opacity are considered. For the two scenarios, some algorithms are presented to achieve the optimal solution for the model by using the method of dynamic programming. And then, the solutions produced by the algorithms are proved correctly by theoretical proof. Finally, some illustrations for the algorithms are given.

The paper presents an error model of a measurement chain containing link performing discrete wavelet transform algorithm which is most often the last stage of measurement signal processing. The goal is to determine the uncertainty budget of the input quantities of the wavelet transform algorithm. The error model takes into account parts of analog, analog-to-digital and digital processing, describing the properties of subsequent fragments of the chain using their transmittance and processing functions. The proposed model enables the description of both deterministic and non-deterministic properties of signal errors. The proposed model was validated using an example measurement chain created for this purpose.

The accompanying mechanical behavior during the flow of ore-rock dispersion is an important factor leading to the instability and failure of the orepass in the ore storage section. It is of great significance to accurately understand the instability characteristics of the well wall in the ore section, during the ore drawing process by understanding the flow characteristics and internal mechanical transfer mechanism of ore-rock dispersion. Using the discrete element method to analyze the flow characteristics, contact compactness, stress distribution characteristics and contact force probability distribution of the ore-rock dispersion, it realizes the quantitative characterization of the damage degree of ore-rock flow, and reveals the damage mechanism of mine storage wall. The results show that: (1) The friction between ore-rock particles, between ore-rock and well wall is an important reason for the change of the macroscopic flow form of ore-rock from “one” to “V” type. (2) The compactness and stress concentration of ores are increasing with the increase of storage depth, which is shown that with the increase of storage depth, the greater the probability of well wall instability. (3) The probability distribution of contact force strength between particles decreases exponentially in the whole process of ore drawing, in which the strong chain plays a major role in the stability of the dispersion system. (4) The number of overpressure and the overpressure coefficient can be used to quantitatively characterize the wall damage degree under the action of ore-rock flow. The dynamic side pressure decreases periodically in the form of exponential fluctuations, and the dynamic load formed by the ore-rock flow mainly played a role in the lower area of the ore storage section.

This paper proposed an optimal tracking control scheme through adaptive dynamic programming(ADP) for a class of partially unknown discrete-time nonlinear systems based on radial basis function neural network(RBF-NN). In order to acquire the unknown system dynamics, we use two RBF-NNs, the one is used to construct the identifier, and the another is used to directly approximate the steady-state control input, where a novel adaptive law is proposed to update neural network weights. While the optimal feedback control and the cost function are derived via feedforward neural networks approximating, it is proposed to regulate the tracking error, the critic network and the actor network are then trained online to obtain the solution of the associated Hamilton–Jacobi–Bellman (HJB) equation being built under the ADP framework. Simulations verify the effectiveness of the optimal tracking control technique using the neural networks.

This paper presents a method to detect and classify the high impedance fault that occur in the medium voltage distribution network using discrete wavelet transform (DWT) and adaptive neuro-fuzzy inference system (ANFIS). The network is designed using Matlab software and various faults such as high impedance, symmetrical and unsymmetrical fault have been applied to study the effectiveness of the proposed ANFIS classifier method. This is achieved by training the ANFIS classifier using the features (standard deviation values) extracted from the three phase fault current signal by DWT technique for various cases of fault with different values of fault resistance in the system. The success and discrimination rate obtained for identifying and classifying the high impedance fault from the proffered method is 100% whereas the values are 66.7% and 85% respectively for conventional fuzzy based approach. The results indicate that the proposed method is more efficient to identify and discriminate the high impedance fault accurately from other power system faults in the system.

Modelling diameter distribution is an integral part of forest management. It requires the use of a suitable probability density function or cumulative distribution function with the appropriate fitting method. In this study, we compared the suitability of eight probability density functions fitted by derivative and optimization methods used for diameter distribution estimation of forest stands. The derivation and optimization were used on A Charlier, beta, generalized beta, gamma, Gumbel, Johnson’s SB, and Weibull (2 and 3-parameter). We used data from 167 permanent sample plots from Atlantic forest (Quercus robur) and 59 temporary sample plots from tropical forests (Tectona grandis). The quality of the fits was evaluated with different indices such as Kolmogorov–Smirnov, Cramer-von Mises, mean absolute error, bias and mean squared error. The results showed Johnson’s SB function was more suitable for describing the diameter distribution of the stands. Johnson’s SB, 3-parameter Weibull and generalized beta had consistent performance regardless of the fitting methods. The quality of fits produced by gamma, Gumbel and 2-parameter Weibull was poor.

The causes of algal blooms in reservoirs are often complexly intertwined with chemical, physical, and biological factors such as the supply of nutrients. Observation of phytoplankton distribution with high spatiotemporal resolution is necessary to track the nutrient sources that cause algal blooms and to understand their behavior in response to wind and water temperature stratification. Observation from a UAV, which has excellent temporal and spatial resolution, is considered to be an effective method to obtain water quality information comprehensively. On the other hand, it is not only the growth of plankton that affects the color of the water surface but also turbidity. Furthermore, since the brightness value of passive sensors such as optical cameras changes depending on the amount of insolation, it is necessary to perform analysis after making corrections for this. In this study, we attempted to develop a method for estimating chlorophyll concentration from aerial images taken from UAVs using machine learning that takes into account brightness correction based on insolation and the spatial distribution of turbidity evaluated by satellite image analysis.

The paper deals with the Compound Inverse Rayleigh distribution, shown to constitute a proper model for the characterization of the probability distribution of extreme values of wind-speed, a topic which is gaining growing interest in the field of renewable generation assessment, both in view of wind power production evaluation and the wind-tower mechanical reliability and safety. The first part of the paper illustrates such model starting from its origin as a generalization of the Inverse Rayleigh model - already proven to be a valid model for extreme wind-speeds - by means of a continuous mixture generated by a Gamma distribution on the scale parameter, which gives rise to its name. Moreover, its validity to interpret different field data is illustrated, also by means of numerous numerical examples based upon real wind speed measurements. Then, a novel Bayes approach for the estimation of such extreme wind-speed model is proposed. The method relies upon the assessment of prior information in a practical way, that should be easily available to system engineers. In practice, the method allows to express one’s prior beliefs both in terms of parameters, as customary, and/or in terms of probabilities. The results of a large set of numerical simulations – using typical values of wind-speed parameters - are reported to illustrate the efficiency and the accuracy of the proposed method. The validity of the approach is also verified in terms of its robustness with respect to significant differences compared to the assumed prior information.

Linear state estimation (SE) formulation under a rectangular coordinate system has been proved to be applicable for real-time distribution network management. Micro phasor measurements’ model can be accommodated into this kind of SE easily. However, voltage magnitude, active power and reactive power measurements are transformed to linear measurements with large node voltage error. To cope with this issue, a linear state estimation under a polar coordinate system is adopted at the first stage to obtain accurate enough complex node voltage, and then nonlinear measurements are transformed to be linear with complex node voltage. At the second stage, linear SE under a rectangular coordinate system can be adopted to satisfy more strict network constraints. The proposed two-stage linear SE is validated on balanced 14, 33, 70,84, 119, 135 nodes network and IEEE 13, 34, 37,123 unbalanced test feeders.

In wireless sensor networks (WSNs), each sensor node can estimate the global parameter from the local data in distributed manner. This paper proposed a robust diffusion estimation algorithm based on minimum error entropy criterion with self-adjusting step-size, which are referred to as diffusion MEE-SAS (DMEE-SAS) algorithm. The DMEE-SAS algorithm has fast speed of convergence and is robust against non-Gaussian noise in the measurements. The detailed performance analysis of the DMEE-SAS algorithm is performed. By combining the DMEE-SAS with diffusion minimum error entropy (DMEE) algorithms, an Improving DMEE-SAS algorithm is proposed, in non-stationary environment where tracking is very important. The Improving DMEE-SAS algorithm can avoid insensitivity of the DMEE-SAS algorithm due to the small effective step-size near the optimal estimator, and obtain a fast convergence speed. Numerical simulations are given to verify the effectiveness and advantages of these proposed algorithms.

This article presents a concise survey of basic discrete and semi-discrete nonlinear models which produce two- and three-dimensional (2D and 3D) solitons, and a summary of main theoretical and experimental results obtained for such solitons. The models are based on the discrete nonlinear Schrödinger (DNLS) equations and their generalizations, such as a system of discrete Gross-Pitaevskii (GP) equations with the Lee-Huang-Yang corrections, the 2D Salerno model (SM), DNLS equations with long-range dipole-dipole and quadrupole-quadrupole interactions, a system of coupled discrete equations for the second-harmonic generation with the quadratic (χ⁽²⁾) nonlinearity, a 2D DNLS equation with a superlattice modulation opening mini-gaps, a discretized NLS equation with rotation, a DNLS coupler and its PT-symmetric version, a system of DNLS equations for the spin-orbit-coupled (SOC) binary Bose-Einstein condensate, and others. The article presents a review of basic species of multidimensional discrete modes, including fundamental and vortex solitons, their bound states, gap solitons populating mini-gaps, symmetric and asymmetric solitons in the conservative and PT-symmetric couplers, cuspons in the 2D SM, discrete SOC solitons of the semi-vortex and mixed-mode types, 3D discrete skyrmions, and some others.

Panel data of our interest consist of a moderate number of panels, while the panels contain a small number of observations. An estimator of common breaks in panel means without a boundary issue for this kind of scenario is proposed. In particular, the novel estimator is able to detect a common break point even when the change happens immediately after the first time point or just before the last observation period. Another advantage of the elaborated change point estimator is that it results in the last observation in situations with no structural breaks. The consistency of the change point estimator in panel data is established. The results are illustrated through a simulation study. As a by-product of the developed estimation technique, a theoretical utilization for correlation structure estimation, hypothesis testing, and bootstrapping in panel data is demonstrated. A practical application to non-life insurance is presented as well.

In this study, we present different methods of estimating fuzzy reliability of a two-parameter Rayleigh distribution via the maximum likelihood estimator, median first-order statistics estimator, quartile estimator, L-moment estimator, and mixed Thompson-type estimator. The mean-square error MSE as a measurement for comparing the considered methods using simulation through different values for the parameters and unalike sample sizes is used. The results of simulation show that the fuzziness values are better than the real values for all sample sizes, as well as the fuzzy reliability at the estimation of the Maximum

Ambulatory heart rate (HR) monitors that acquire electrocardiogram (ECG) or/and photoplethysmographm (PPG) signals from the torso, wrists, or ears are notably less accurate in tasks associated with high levels of movement compared to clinical measurements. However, reliable estimation of HR can be obtained through data fusion from different sensors. These methods are especially suitable for multimodal hearable devices, where heart rate can be tracked from different modalities, including electrical ECG, optical PPG, and sounds (heart tones). Combined information from different modalities can compensate for single source limitations. In this paper, we evaluate the possible application of data fusion methods in hearables. We assess data fusion for heart rate estimation from simultaneous in-ear ECG and in-ear PPG, recorded on eight subjects while performing 5-minute sitting and walking tasks. Our findings show that data fusion methods provide a similar level of mean absolute error as the best single-source heart rate estimation but with much lower intra-subject variability, especially during walking activities. We conclude that data fusion methods provide more robust HR estimation than a single cardiovascular signal. These methods can enhance the performance of wearable devices, especially multimodal hearables, in heart rate tracking during physical activity.

The variation of aerosols, especially dust aerosol, in time and space plays an important role in climate forcing studies. Aerosols can effectively reduce land surface longwave emission and re-emit energy at a colder temperature, making estimation of downwelling surface longwave radiation (DSLR) with satellite data difficult. Using the latest atmospheric radiative transfer code (MODTRAN 5.0), we simulate the outgoing longwave radiation (OLR) and DSLR under different land surface and atmospheric profile conditions. The results show that dust aerosol has an obvious “warming” effect to longwave radiation compared with other aerosols, that aerosol longwave radiative forcing (ALRF) increased with increasing aerosol optical depth (AOD), and that the atmospheric water vapor content (WVC) is critical to the understanding of ALRF. A method is proposed to improve the accuracy of DSLR estimation from satellite data for the skies under heavy dust aerosols. The AOD and atmospheric WVC under cloud-free conditions with a relatively simple satellite-based radiation model that yields the high accurate DSLR under heavy dust aerosol are used explicitly as model input to reduce the effects of dust aerosol on the estimation of DSLR. Validations of the proposed model with satellites data and field measurements show that it estimates the DSLR accurately under heavy dust aerosol skies. The root mean square errors (RMSEs) are 20.4 W/M² and 24.2 W/M² for Terra and Aqua satellites, respectively, at the Yingke site, and the biases are 2.7 W/M² and 9.6 W/M², respectively. For the Arvaikheer site, the RMSEs are 23.2 W/M² and 19.8 W/M² for Terra and Aqua, respectively, and the biases are 7.8 W/M² and 10.5 W/M², respectively. The proposed method is especially applicable to acquire relatively high accurate DSLR under heavy dust aerosol using MODIS data with available WVC and AOD data.

The Active Distribution Networks (ADNs) are Multi-Input Multi-Output (MIMO) systems with coupled dynamics which cause interactions among the control loops of the Distributed Energy Resources (DERs). This undesired effect lead to performance degradation in the voltage control. To mitigate the effects of these unavoidable coupling, the present paper proposes a systematic design procedure based on the analysis of the interaction's sources. In details, each DER is equipped with a double loop PI to control the active and reactive powers output by the Voltage Source Converter (VSC) which connects the DER to the network's node. Furthermore, to guarantee Ancillary Services (ASs), the two loops are coupled by a simple mechanism of cooperation of the active power to voltage regulation realized by a filtered droop law. To achieve voltage regulation with reduced loop interactions, the PI parameters and the filter's pulse are designed according to a two sequential steps procedure based on the (Internal Model Control) IMC technique. Simulation studies are finally presented to demonstrate that the proposed design method achieves both reduction of the loop interaction and robust voltage control in presence of model parameter uncertainty in the MIMO model.

An effective empirical statistical method is developed to improve the process of mineral resource estimation of seabed polymetallic nodules and is applied to analyse the abundance of seabed polymetallic nodules in the Clarion Clipperton Zone (CCZ). The newly proposed method is based on three hypotheses as the foundation for a model of “Idealized Nodules”, which was validated by analysing nodule samples collected from the seabed within the Tonga Offshore Mining Limited (TOML) exploration contract. Once validated, the “Idealized Nodule” model was used to deduce a set of empirical formulae for predicting the nodule resources, in terms of Percentage Coverage and Abundance. The formulae were then applied to analysing a total of 188 sets of nodule samples collected across the TOML areas, comprising box-core samples and towed camera images collected by one of the authors and detailed in [4]. The analysis also relies upon detailed box-core sample measurements from other areas reported by [7]. Numerical results for resource prediction were compared with field measurements, and reasonable agreement has been achieved. The new method has the potential to achieve more accurate mineral resource estimation with reduced sample numbers and sizes. They may also have application in improving the efficiency of design and configuration of mining equipment.

Daily maximum surface air temperature (Tamax) is a crucial factor for understanding complex land surface processes under rapid climate change. Remote detection of Tamax has widely relied on the empirical relationship between air temperature and land surface temperature (LST), a product derived from remote sensing. However, little is known about how such a relationship is affected by the high heterogeneity in landscapes and dynamics in seasonality. This study aims to advance our understanding of the roles of land cover and seasonal variation in the estimation of Tamax using the MODIS (Moderate Resolution Imaging Spectroradiometer) LST product. We developed statistical models to link Tamax and LST in the middle and lower reaches of the Yangtze River in China for five major land-cover types (i.e., forest, shrub, water, impervious surface, cropland, and grassland) and two seasons (i.e., growing season and non-growing season). Results show that the performance of modeling the Tamax-LST relationship was highly dependent on land cover and seasonal variation. Estimating Tamax over grasslands and water bodies achieved superior performance; while uncertainties were high over forested lands that contained extensive heterogeneity in species types, plant structure, and topography. We further found that all the land-cover specific models developed for the plant non-growing season outperformed the corresponding models developed for the growing season. Discrepancies in model performance mainly occurred in the vegetated areas (forest, cropland, and shrub), suggesting an important role of plant phenology in defining the statistical relationship between Tamax and LST. For impervious surfaces, the challenge of capturing the high spatial heterogeneity in urban settings using the low-resolution MODIS data made Tamax estimation a difficult task, which was especially true in the growing season.

Moso bamboo is famous for fast growing and biomass accumulation, as well as high annual output for timber and bamboo shoots. These high outputs require high nutrient inputs to maintain and improve stand productivity. Soil nitrogen (N), phosphorus (P) and potassium (K) are important micronutrients for plant growth and productivity. Due to high variability of soils, analysing spatial patterns of soil N, P and K stocks is necessary for scientific nutrient management in Moso bamboo forests. In this study, soils were sampled from 138 locations across Yong’an City and ordinary kriging was applied for spatial interpolation of soil N, P and K stocks. Soil N stock showed a strong spatial dependence while soil N and P stocks presented a moderate spatial dependence, indicating soil N was mainly controlled by intrinsic factors while soil N and P stocks were controlled by both intrinsic and extrinsic factors. Different spatial patterns were observed for soil N, P and K stocks across the whole study area, indicating that fertilizations with different ratios of N:P:K should be applied for different sites to maintain and improve stand productivity. The total soil N, P and K stocks within 0-60 cm were 0.624, 0.020 and 0.583 Tg, respectively.

The interest in virtual fencing, as opposed to traditional physical barriers, has increased since it makes it easier to move the boundaries of animal enclosures for agricultural and rewilding purposes. Despite this, the implementation of virtual fence systems, such as the Nofence© system, requires large financial investments. The necessity of individual collars also makes it more difficult to implement this system on a larger scale. If a herd of cattle tend to maintain cohesive social groups around dominant individuals, it might be possible to use fewer collars thereby making the system more cost efficient. This study investigates the pattern in spatial distribution in a herd of 17 Angus cows on Fanø in Denmark with GPS locations, using a Nofence© system. The aim of this paper is to determine how individuals position themselves in a herd, spatially, and identify a pattern in ranks. The method used in this study examines the distances between an individual to the rest of the herd mates using nearest neighbour and distribution plots. Contrary to the expectations, this study did not show a distinct pattern in herd distribution. While some tendencies in spatial distribution patterns were observed, only a low concordance could be found (W=0.15,p=3.02e−16). Therefore, interactions among individuals should be considered in future studies, rather than only investigating the herd as one unit. A cumulative curves of the ranks estimated over the entire periods allowed however a rough estimation of the hierarchy and allowed the identification of the highest ranked cows. Therefore it might be possible to use fewer collars if the high ranked cows are discovered, but the estimation of the hierarchy will require a period of observation overall several days as the hierarchical positions are fluctuating too much from days to days. Therefore the cumulative rank curves seems to be a valid solution.

Revealing the spatial distribution mechanism of the coupling effect between soil water and salt is essential for managing saline-alkali farmland. This study utilized the multifractal method to analyze soil water and salt under varying water content levels in arid saline-alkali farmland. Soil samples were collected on the second (S1), fifth (S2), eighth (S3), eleventh (S4), and fourteenth (S5) days after a rainfall event. The findings showed that a significant decrease in soil water content, contrasting with an increase in soil salt content throughout the soil layers post-rainfall. As water content decreased, the spatial variability of soil water initially increased from S1-3 and then decreased, while the spatial variability of soil salt decreased. The spatial distribution between soil water and salt exhibited a high correlation at S3-4 (with relative water content of soil ranging from 0.52 to 0.75) due to their coupling effect. However, soil salt was unevenly leached by rainfall at high water content levels (S1-2) and precipitated at low water content levels (S5), resulting in low spatial variability correlations between soil water and salt. This study elucidated the coupling process of soil water and salt, identifying their spatiotemporal distribution mechanism in dryland agricultural areas.

This article presents a review of the distribution system state estimation (DSSE) algorithms from the model-based approach to the recent data-driven methodologies. The insufficient measurements which result in low observability have motivated the need to shift from the conventional algorithm to data-driven approaches that can successfully estimate states despite the DSSE challenges. This article discusses the nonlinearity in the DSSE problem formulation and how different model-based methods have been proposed to mitigate the problems of robustness, ill-conditioning, and complexity. Moreover, approximate DSSE to obviate nonlinearity were discussed: complex linearization, small angle approximation, convexification and compressed sensing. Furthermore, probabilistic DSSE methods were also discussed in the need to quantify the uncertainty associated with the state estimation results. Also, data-driven methods applicable to DSSE, pseudo measurement generations, and topology identification were also discussed using machine and deep learning methods. Lastly, a recent approach that employs a hybrid of model-based and data-driven methods using matrix and tensor completion is surfacing because they work in a low observable condition of the network and can estimate the states satisfactorily. With this review, researchers can look further into developing and improving on the model-based data-driven methods less susceptible to the barriers in conventional DSSE algorithms.

In this paper we describe a nonparametric maximum likelihood (NPML) method for estimating multivariate mixing distributions. Given $N$ independent observations, convexity theory shows that the NPML estimator is discrete with at most $N$ support points. The original infinite NPML problem then becomes the finite dimensional problem of finding the location and probability of the support points. The probability of the support points is found by a Primal-Dual Interior-Point method; the location of the support points is found by an Adaptive Grid method. Our method is able to handle high-dimensional and complex multivariate mixture models.An important application is discussed for the problem of population pharmacokinetics and a non-trivial example is treated.Our algorithm has been successfully applied in hundreds of published pharmacometric studies. In addition to population pharmacokinetics, this research also applies to empirical Bayes estimation and many other areas of applied mathematics.

Infrared thermography has already been proven to be a significant method in non-destructive evaluation since it gives information with immediacy, rapidity and low cost. However, the thorniest issue for wider application of IRT is the quantification. In this work, we proposed a specific depth quantifying technique by employing the Gated Recurrent Units (GRU) in composite material samples via pulsed thermography (PT). Carbon Fiber Reinforced Polymer(CFRP) embedded with flat bottom holes were designed via Finite Element Method (FEM) modeling in order to precisely control the depth and geometrics of the defects. The GRU model automatically quantified the depth of defects presented in the CFRP material. The proposed method evaluated the accuracy and performance of synthetic CFRP data from FEM for defect depth predictions.

Accurate power curve modelling is essential to continuously evaluate the performance of a wind turbine (WT). In this work, we characterize the wind power curves using SCADA data acquired at a frequency of 5 minutes in a wind farm (WF) consisting of 5 WTs. Regarding the non-parametric methods, we select artificial neural networks (ANNs) to make curve estimations. Given that, we have the curves provided by the manufacturer of the WTs given by some very precisely measured pair of wind speed and power points. We can evaluate the difference between the manufacturer characterization and the ones estimated with the data provided by the SCADA system. Before the estimation, we propose a method of filtering the anomalies based on the characteristics provided by the manufacturer. We use three-quarters of the available data for curve estimation and one-quarter for the test. One WT suffered a break in the test part, so we can check how the test estimates reflect this problem in its wind-power curve compared to the estimations obtained in the WTs that worked adequately.

As the adoption of renewable energy sources grows, ensuring a stable power balance across various time frames has become a central challenge for modern power systems. In line with the "dual carbon" objectives and the seamless integration of renewable energy sources, harnessing the advantages of various energy storage resources and coordinating the operation of long-term and short-term storage have become pivotal directions for future energy storage deployment. To address the complexities arising from the coupling of different time scales in optimizing energy storage capacity, this paper proposes a method for energy storage planning that accounts for power imbalance risks across multiple time scales. Initially, the Seasonal and Trend decomposition using the Loess (STL) decomposition method is utilized to temporally decouple actual operational data. Subsequently, power balance computations are performed based on the obtained data at various time scales to optimize the allocation of different types of energy storage capacities and assess the associated imbalance risks. Finally, the effectiveness of the proposed approach is validated through hourly applications using real-world data from a province in China over recent years.

Urban rail transit hub platform is the most important area for passenger flow distribution. In order to calculate passenger flow volume in platform and evaluate platform service level during rush hours, this paper presents a method for modeling and simulation of passenger flow distribution in platform. Passenger flow distribution model (PFDM) is proposed based on the basic analysis and the superposition principle of passenger flow. Simulation design for PFDM is proposed by Anylogic, which contains simulation process and simulation model. Experiment results show that PFDM and simulation design are effective and accordant with the reality scenario, and the simulation precision is comparatively ideal. This research could provide a beneficial reference for train scheduling and operation management under the viewpoint of traffic safety and service level.

Accelerated degradation tests (ADT) are widely used in the manufacturing industry to obtain information on the reliability of components and materials, through degrading the lifespan of the product by applying an acceleration factor which causes damage to the material. The main objective is to obtain fast information which is modeled to estimate the characteristics of the material life under normal conditions of use and to save time and expenses. The purpose of this work is to estimate the lifespan distribution of gold nanoparticles stabilized with lipoic acid (GNPs@LA) through accelerated degradation tests applying sodium chloride (NaCl) as an acceleration factor. For this, the synthesis of GNPs@LA was carried out, a constant stress ADT (CSADT) was applied, and the non-linear Wiener process was proposed with random effects, error measures and different covariability for the adjustment of the degradation signals. The information obtained with the test and analysis allows us to obtain the life distribution in GNPs@LA, the results make possible to determine the guaranteed time for a possible commercialization and successful application based on the stability of the material. In addition, for the evaluation and selection of the model, the Akaike and Bootstraping criteria were used.

This paper studies the duration patterns of unemployment spells for immigrants and the determinants of unemployment’s completion into one of a number of studied labour market states in Finland. We estimate a duration model for unemployment with competing risks of its terminating into employment, labour market training or economic inactivity. Taking into account the wide period of observation and opportunities to analyse processes of labour market integration during various periods of economic development in Finland, in combination with the individualistic character of the labour careers of immigrants, this research is beneficial owing to the many various findings concerning labour market integration of immigrants. The approach undertaken in this research has a dualistic “descriptive-dynamic” character under which integration is understood as a never-ending process, which is conditioned by a time period of long-term existence and a context of solitary action. We find that transitions out of unemployment spells have a cyclical character; after every new “cycle” in unemployment, the probability of terminating unemployment decreases further. We also find that ascriptive factors make sense in the process of job-placement of immigrants from unemployment. Therefore, the gender, education and age of immigrants, as well as the effect of the period in which first unemployment occurred, potentially predict transitions out of unemployment and further labour market integration of immigrants.

Fire detection from satellite sensors relies on an accurate estimation of the unperturbed state of a target pixel, from which an anomaly can be isolated. Methods for estimating the radiation budget of a pixel without fire depend upon training data derived from the location's recent history of brightness temperature variation over the diurnal cycle, which can be vulnerable to cloud contamination and the effects of weather. This study proposes a new method that utilises the common solar budget found at a given latitude in conjunction with an area's local solar time to aggregate a broad-area training dataset, which can be used to model the expected diurnal temperature cycle of a location. This training data is then used in a temperature fitting process with the measured brightness temperatures in a pixel, and compared to pixel-derived training data and contextual methods of background temperature determination. Results of this study show similar accuracy between clear-sky medium wave infrared upwelling radiation and the diurnal temperature cycle estimation compared to previous methods, with demonstrable improvements in processing time and training data availability. This method can be used in conjunction with brightness temperature thresholds to provide a baseline for upwelling radiation, from which positive thermal anomalies such as fire can be isolated.

Motivated by empirical observations, we propose a possible extension of Gibrat's law. By applying it into the random growth theory of income distribution (Gabaix, 2009), we find that the income distribution is described by a generalized Pareto-type distribution (GPD) with three parameters. We observe that there is a parameter /elta in the GPD that plays a key role in determining the shape of income distribution. By using the Kolmogorov-Smirnov test, we empirically show that, for typical market-economy countries, /elta is close to 0 significantly, such that the income distribution is characterized by a two-class pattern in which the bottom 90% of the population is approximated by an exponential distribution and the richest 1%~3% is approximated by an asymptotic power law. However, we empirically find that, for China both in planned economy period and in early stages of market reformation (from 1978 to 1990), /elta is significantly deviated from 0, such that the bottom of the population no longer conforms to an exponential distribution.

The main objective of this work is the assessment of the annual number of hospital admissions for respiratory disease (HARD) due to the exposure to in-healable particulate matter (PM10), within the greater Athens area (GAA), Greece. Towards this aim, the time series of the particulate matter with aerodynamic diameter less than 10μm (PM10) recorded in six monitoring stations located in the GAA, for a 13-year period 2001-2013, is used. Initially, a descriptive statistical treatment of PM10 concentrations took place. Furthermore, the AirQ2.2.3 software developed by the WHO was used to evaluate adverse health effects by PM10 in the GAA during the examined period. The results show that, during the examined period PM10 concentrations present a significant decreasing trend. Also, the mean annual HARD cases per 100,000 inhabitants ranged between 20 (suburban location) and 40 (city centre location). Approximately 70% of the annual HARD cases are due to city centre residents. In all examined locations, a declining trend in the annual number of HARD cases is appeared. Moreover, a strong relation between the annual number of HARD cases and the annual number of days exceeding the European Union daily PM10 threshold value was found.

Employing remote sensing techniques to evaluate crop residue cover (CRC) enables rapid mapping of tillage practices over large cropland areas, which has significant practical implications for monitoring the implementation of conservation tillage practices. This research utilized Sentinel-2 MSI and Landsat-8 OLI remotely sensed data, various spectral indices were calculated and the datasets were correlated with the field-measured CRC data. Selected spectral indices, highly significantly correlated, were chosen to build a correlation model between the indices and the CRC. Considering Sentinel-2 MSI and Landsat 8 OLI remotely sensed data variations in spectral and spatial scales, the accuracy of CRC models was evaluated using R2 and RMSE values. From the results, it's evident that all six spectral indices, including STI, NDTI, SRNDR, NDI5, NDI7 and NDSVI, registered correlation coefficients with the CRC exceeded 0.3 in the study area. STI and NDTI showed a stronger correlation to CRC, with coefficients of 0.810 and 0.803, respectively. For the Landsat 8 data, the lowest correlation obtained was for NDSVI with the Landsat 8 data with a correlation coefficient of 0.358. Applying the linear regression technique, STI and NDTI obtained higher accuracy compared to other indices in the Sentinel-2 MSI data. For NDTI and STI, The R2 values were 0.650 and 0.665, respectively, while their corresponding RMSE values were 3.96% and 4.03%. The highest model accuracy was obtained for the STI and NDTI constructed from the Landsat 8 OLI data (R2 values 0.41 and 0.40; RMSE values 5.25% and 5.27%, respectively). The CRC model built using partial least squares regression (PLSR) exhibited greater precision when applied to the Sentinel-2 MSI data (R2 value 0.894; RMSE value 2.18%) in comparison to the Landsat 8 OLI data (R2 value 0.695; RMSE value 3.70%). PLSR has the potential to enhance the precision of the model. In addition, the study area shows that the estimation of the CRC is better in the one-dimensional linear regression model and the PLSR model created with the Sentinel-2 MSI optical data compared to the Landsat 8 OLI optical remotely sensed data. The validation suggests that the CRC estimation in the study area is more precise when using the model relying on Sentinel-2 satellite imagery rather than the Landsat 8 OLI imagery.

Remote sensing datasets are increasingly being used to provide spatially explicit large scale evapotranspiration (ET) estimates. The focus of this study was to estimate and thematically map on a pixel-by-pixel basis, the actual evapotranspiration (ET_a) of the Wonji Shoa Sugarcane Estate using the Surface Energy Balance Algorithm for Land (SEBAL), Simplified Surface Energy Balance (SSEB) and Operational Simplified Surface Energy Balance (SSEBop) algorithms. The results obtained revealed that the ranges of the daily ETa estimated on January 25, February 26, September 06 and October 08, 2002 using SEBAL were 0.0 - 6.85, 0.0 – 9.36, 0.0 – 3.61, 0.0 – 6.83 mm/day; using SSEB 0.0 - 6.78, 0.0 – 7.81, 0.0 – 3.65, 0.0 – 6.46 mm/day, and SSEBop were 0.05 - 8.25, 0.0 – 8.82, 0.2 – 4.0, 0.0 – 7.40 mm/day, respectively. The Root Mean Square Error (RMSE) values between SSEB and SEBAL, SSEBop and SEBAL, and SSEB and SSEBop were 0.548, 0.548, and 0.99 for January 25, 2002; 0.739, 0.753, and 0.994 for February 26, 2002;0.847, 0.846, and 0.999 for September 06, 2002; 0.573, 0.573, and 1.00 for October 08, 2002, respectively. The standard deviation of ET_a over the sugarcane estate showed high spatio-temporal variability perhaps due to soil moisture variability and surface cover. The three algorithm results showed that well watered sugarcane fields in the mid-season growing stage of the crop had higher ET_a values compared with the other dry agricultural fields confirming that they consumptively use more water. Generally during the dry season, ET_a is limited to water surplus areas only and in wet season, ET_a was high throughout the entire sugarcane estate. The evaporation fraction (ETrF) results also followed the same pattern as the daily ET_a over the sugarcane estate. The total crop and irrigation water requirement and effective rainfall estimated using the Cropwat model were 2468.8, 2061.6 and 423.8 mm/yr for January 2001 planted and 2281.9, 1851.0 and 437.8 mm/yr for March 2001 planted sugarcanes, respectively. The mean annual ET_a estimated for the whole estate were 107 Mm³, 140 Mm³, and 178 Mm³ using SEBAL, SSEB, and SSEBop, respectively. Even though the algorithms should be validated through field observation, they have potential to be used for effective estimation of ET in the sugarcane estate.

Different variants of thinning for discrete random variables are studied. The thinning procedure allows to introduce an analog of scale parameter for positive integer-valued random variables. Sufficient and necessary conditions for the existence of such a scale are given.

Introduction: Hypertension is a serious public health problem in Morocco. The objective of this pilot study was to estimate the salt intake of Moroccan adults by measuring 24-hour urinary sodium excretion. Methods: 132 participants were recruited from the central region of Morocco. Sodium, potassium and creatinine levels were measured using the 24-hour urine samples. Additional data included anthropometric measurements and socio-demographic questionnaire. Results: The daily sodium, potassium and creatinine excretions were 2779.1±1334.9 mg, 1350.0±642.8 mg and 820.3 ± 381.2 mg, respectively. Data analysis revealed that 71.2% of the participants had a daily sodium intake more than 2000 mg (5g of salt) recommended by the WHO. The mean of potassium excretion was 1350.0±642.8 mg /day, this average was lower than 3.51g per day recommended by the WHO. Conclusion: Public health measures to reduce sodium and increase potassium consumption in order to decrease the population’s risk of high blood pressure and heart disease are recommended.

The aim of this study is to obtain the Bayes estimators and the maximum likelihood estimators (MLEs) for the unknown parameters of the Rayleigh Weibull (RW) distribution based on progres-sive type-II censored samples. The approximate Bayes estimators are calculated using the idea of Lindley and Tierney-Kadane's approximation method under the squared-error loss function when the Bayes estimators are not handed in explicit forms. In this study, the approximate Bayes esti-mates are compared with the maximum likelihood estimates in the aspect of the estimated risks (ERs) using Monte Carlo simulation. In addition, the coverage probabilities of the parametric bootstrap estimates are calculated. Real lifetime data sets belonging to the cancer types as bladder cancer, head and neck cancer, and leukemia are used to illustrate the emprical results belonging to the approximate Bayes estimates, the maximum likelihood estimates, and the parametric bootstrap intervals.

This article presents a new methodology for modeling and analyzing discrete dynamic systems through the construction of inverse algebraic models, giving rise to the Theory of Inverse Discrete Dynamical Systems. Key concepts such as inverse modeling, structural analysis in inverse algebraic trees, and the establishment of topological equivalences for the analytical transfer of properties between the canonical system and its inverted counterpart are introduced. Central theorems on homeomorphic invariance and topological transport are demonstrated, consolidating the validity of transferring cardinal attributes between both equivalent dynamic representations. As a concrete application, an alternative proof of the historic Collatz Conjecture is presented through the rigorous construction of the associated inverse model and the analytical transfer of properties exhibited in the inverted tree structure. Furthermore, evidence has been presented on Conway's Conjecture in the Game of Life through meticulous inverse modeling in bounded cases, overcoming initial limitations of the theory in combinatorial explosions. Although a complete proof of the conjecture is not achieved, this positively expands the modeling and analysis capacity of the methodology. The impact of properly applying this novel theory to expand understanding and solve open problems in discrete dynamical systems seems vast and profound.

Hasselmann and coauthors proposed the discrete interaction approximation (DIA) as the best tool replacing the nonlinear evolution term in a numerical wind-wave model. Much later, Polnikov and Farina radically improved the original DIA by means of location all the interacting four wave vectors, used in the DIA configuration, exactly at the nodes of the numerical frequency-angular grid. This provides nearly two-times enhancing the speed of numerical calculation for the nonlinear evolution term in a wind-wave model. For this reason, the proposed version of the DIA was called as the fast DIA (FDIA). In this paper we demonstrate all details of the FDIA concept for several frequency-angular numerical grids of high resolution, with the aim of active implementation the FDIA in modern versions of world-wide used wind-wave models.

We formulate the statistics of the discrete multicomponent fragmentation event using a methodology borrowed from statistical mechanics. We generate the ensemble of all feasible distributions that can be formed when a single integer multicomponent mass is broken into fixed number of fragments and calculate the combinatorial multiplicity of all distributions in the set. We define random fragmentation by the condition that the probability of distribution be proportional to its multiplicity and obtain the partition function and the mean distribution in closed form. We then introduce a functional that biases the probability of distribution to produce in a systematic manner fragment distributions that deviate to any arbitrary degree from the random case. We corroborate the results of the theory by Monte Carlo simulation and demonstrate examples in which components in sieve cuts of the fragment distribution undergo preferential mixing or segregation relative to the parent particle.

This algorithm is designed to perform Discrete Fourier Transforms (DFT) to convert temporal data into spectral data. This algorithm obtains the Fourier Transforms by studying the Coefficient of Determination of a series of artificial sinusoidal functions with the temporal data, and normalizing the variance data into a high-resolution spectral representation of the time-domain data with a finite sampling rate. What is especially beneficial about this DFT algorithm is that it can produce spectral data at any user-defined resolution.

In order to investigate the size distributions and seasonal variations of carbonaceous aerosols (OC and EC), the carbonaceous species were collected and then analyzed by using a 9-stage Anderson-type aerosol sampler and DRI Model 2001A Thermal/Optical Carbon Analyzer on the typical industrial city Nanjing in Yangtze River Delta, China in the summer, autumn and winter of 2013 and spring of 2014. OC, EC, SOC and POC exhibited obvious seasonal variations, with the highest level in winter (39.1±14.0, 5.7±2.1, 23.6±11.7 and 14.1±5.7 μg•m^-3) and the lowest level in summer (20.6±6.7, 3.3±2.0, 12.2±3.8 and 8.4±4.1 μg•m^-3), and were mainly centralized in PM_1.1 in four seasons. The concentrations of OC in PM_1.1 varied in the order of winter > autumn > spring > summer, while EC ranked in the order of autumn > winter > summer > spring. In the PM_1.1-2.1 and PM_2.1-10, the concentrations of OC and EC decreased in the sequence of winter > spring > autumn > summer. The size spectra of OC, EC and SOC had bimodal distributions in four seasons, except for EC with four peaks in summer. The size spectra of POC varied greatly with seasons, exhibiting bimodal distribution in winter, trimodal distribution in spring and summer, and four peaks in autumn. The OC/EC ratios were 7.0, 6.3, 7.6 and 6.9 in spring, summer, autumn and winter, respectively, which demonstrated the abundance of secondary organic aerosols in Nanjing. The sources of carbonaceous aerosol varied significantly with seasons, and were dominated by vehicle exhaust, coal and biomass burning in PM_2.1, and dominant by dust, coal and biomass burning in PM_2.1-10.

Abstract: Global ecological concern regarding the transfer of fouling organisms to ship hulls is increasing. This study investigated the species composition, dominant species, distribution patterns, community structure, and life-cycle differences of hull-fouling macroinvertebrates on five research vessels (R/Vs: Isabu, Onnuri, Eardo, Jangmok 1, and Jangmok 2) operated by the Korea Institute of Ocean Science and Technology (KIOST). Hull-fouling macroinvertebrates were collected three to five times on quadrats from the upper and middle sectors of the hull sides, bottom, and niche areas (the propellers and shafts and the thrusters). A total of 47 macroinvertebrate species were identified, represented by 8,519 individuals(ind.)/m2 and a biomass of 1,967 gWWt/m2 on the five vessels. The number of species, density, and biomass were greater on the coastal vessels the Eardo, Jangmok 1, and Jangmok 2 than on the ocean-going vessels the Isabu and Onnuri. Among the coastal vessels, barnacles were the most abundant and had the greatest density, while mollusks had the highest biomass. Differences between hull sectors showed that the highest species abundance and density appeared on all hulls in ports and bays where the Jangmok 1 operated, while the highest species abundance, density, and biomass were identified in the niche areas of the Eardo, which operated farther from the coast. The hull-fouling macroinvertebrates that exceeded 1% of all organisms were the barnacles Amphibalanus amphitrite, Balanus trigonus, and Amphibalanus improvisus; the polychaete Hydroides ezoensis; the bivalves Magallana gigas and Mytilus galloprovincialis; and the amphipod Jassa slatteryi. The dominant species were cosmopolitan and globally distributed, and many of them were cryptogenic. Six native species were identified: M. gigas, H. ezoensis, the amphipod Melita koreana, the isopod Cirolana ko-reana, and the barnacle B. trigonus and F. kondakovi. Eight non-indigenous species (NIS) were de-tected: barnacles A. amphitrite and A. improvisus, the bivalve M. galloprovincialis, the polychaete Perinereis nuntia, amphipods J. slatteryi and Caprella californica, and bryozoans Bugulina californica and Bugula neritina. Of the fouling macroinvertebrates found on the vessel hulls, 13% were native and 17% were NIS. The nMDS ordination analysis of the five vessels showed four clusters of fouling macroinvertebrate communities: two ocean clusters (the Isabu and Onnuri), one coastal cluster (the Eardo and Jangmok 2), and one local cluster (the Jangmok 1). More diverse communities developed on the hulls of vessels that operated locally rather than globally or in deep oceans. Spearman rank correlation showed that biological indices such as number of species, density, and biomass correlated highly and positively with the numbers of days on coastal operations and negatively with the number days on ocean operations. The species diversity index correlated positively with the total number of anchoring days and coastal operation days and negatively with the total number of operation days and ocean operation days. The macroinvertebrates differed by the area of operation, the port of anchorage, the number of days in operation and at anchor, and the hull sectors. There is no previous research data on hull-fouling macroinvertebrates in the Republic of Korea, and this study provides a basis for future studies to identify introduced species and their differences based on operation area.

Recycling waste products is an environmental-friendly activity that can bring benefits to accompany, saving manufacturing costs and improving economic efficiency. For the beer industry, recycling bottles can reduce manufacturing costs and reduce the industry's carbon footprint. This paper presents a model for a multi-objective collection-distribution center location and allocation problem in a closed loop supply chain for the beer industry, in which the objective is to minimize total costs and transportation pollution. Uncertainties in the form of randomness and fuzziness are jointly handled in this paper to ensure a more practical problem solution, for which returned bottle sand unusable bottles are considered fuzzy random variables. A heuristic algorithm based on priority-based global-local-neighbor particle swarm optimization (pb-glnPSO) is applied to ensure reliable solutions for this NP-hard problem. A case study on a beer operation company is conducted to illustrate the application of the proposed model and demonstrate the priority-based global-local-neighbor particle swarm optimization.

With the progression of AI, embedding advanced AI technologies into small robotics and mobile devices has become essential, driving research towards lightweighting AI models. Our study enhances the EPINET depth estimation model for light field images, aiming for compactness and faster inference while preserving accuracy. We conducted two-step experiments aimed at enhancing inference efficiency: Initially, by adjusting input streams and convolution layers, we simplified the CNN model, achieving faster inference times at the cost of reduced accuracy. To address this reduction in accuracy, we then applied knowledge distillation, allowing the simplified model to learn from the original model’s more complex patterns. In our quantitative experiments using two error metrics, MSE (Mean Squared Error) and BadPix, we identified optimal knowledge positions and evaluated the required complexity for the student model. As a result, our method improved MSE by 21% and BadPix by 14% compared to training without it. Furthermore, the student model achieved an inference speed 13% faster than the teacher model and surpassed its accuracy by 10% in MSE. Additionally, we demonstrated that repeatedly applying our approach could further enhance both model compactness and accuracy.

Xinjiang, as a climate-sensitive region in northwest China, holds significant importance in studying extreme weather events for agricultural production and socioeconomic development. Based on data from 52 meteorological stations and 23 extreme climate indicators in Xinjiang from 1960 to 2021, the 5a moving average, linear trend estimation, climate tendency rate, and GIS spatial analysis methods are used to analyze the distribution and change characteristics of extreme climate. The results revealed that from 1960 to 2021, TXx, TNx, TXn, TNn, FD, ID, TN10p, TX10p, and WSDI showed a significant upward trend, while SU, TR, TN90p, TX90p, and CSDI showed a significant downward trend. What is more, the range of changes in the extreme temperature indexes from 1990 to 2021 is higher than that of 1960 to 1989. The high-value areas of extreme indexes are primarily concentrated in the northern part, while the high-value areas of cold event indexes are mainly distributed in the southern part. The upward/downward trends all accounting for over 80.00% of the entire region. The Rx1day, Rx5day, PRCPTOT, R95p, R99p, R10, R20, and R25 showed a significant growth trend from 1960 to 2021, and the range of change in the extreme precipitation indexes from 1990 to 2021 was lower than that from 1960 to 1989. Moreover, their high-value areas and high climate tendency areas are mainly concentrated in the northern and western parts of Xinjiang, with upward trends all accounting for over 71.00% of the entire region. The research results provide theoretical foundations for formulating climate risk strategies in the northwest region of China.

Contemporary financial systems, despite their inherent complexity and resilience, face significant challenges related to transparency, efficiency, and security. Notable deficiencies in transparency and emerging security vulnerabilities underscore the pressing need for innovative solutions in current financial practices. Hence, this paper introduces a financial distribution system based on a consortium blockchain. Given organizations’ preference for keeping sensitive data private, the partially decentralized nature of consortium blockchains plays a pivotal role in the proposed system’s architecture. Moreover, within this system, the proposed smart contract serves a dual role. It functions not only as a transactional tool but also as a specialized solution designed to improve transactional procedures. Its multifaceted capabilities include ensuring data accuracy, facilitating swift execution, and ensuring adherence to regulatory requirements within the dynamic financial distribution landscape. Through rigorous testing and empirical analysis utilizing the Hyperledger Besu platform, it has been observed that the performance of the proposed system surpasses traditional blockchain frameworks. The results unequivocally demonstrate the proposed system’s proficiency in handling high-frequency financial transactions with minimal computational burden.

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal “core” population plus a tenuous energetic “halo” population. At first, such a feature was empirically fitted with the Kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the Kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously-emitted thermal fluctuations and the background turbulence leads to the formation of non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the Kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.

Differences are introduced as outputs of linear systems called differencers, being considered two classes: shift and scale-invariant. Several types are presented, namely: nabla and delta, bilateral, tempered, bilinear, stretching, and shrinking. Both continuous and discrete-time are described. ARMA type systems based on differencers are introduced and exemplified. In passing, the incorrectness of the usual delta difference is shown.

The basis of the discrete-impulse energy supply (DIES) concept is the efficient use of supplied energy. The references describe in detail the general principles of DIES, examine the energy and thermodynamic aspects and the main mechanisms of intensification that can be initiated on the basis of this principle. DIES mechanisms conveniently can be divided into hard and soft ones. The former should be used to stimulate hydromechanical processes, and the latter to accelerate the processes of phase heat and mass transfer, or for the purpose of intensive mixing of multicomponent media. The authors have studied the possibility of using DIES to intensify the hydromechanical processes, in particular emulsification of milk fat (homogenization of milk, preparation of spreads), processing of cream cheese masses. Objects of research were whole non-homogenized milk, fat emulsions, cream cheese mass. In order to evaluate the efficiency of milk homogenization the homogenization coefficient change was studied, which was determined by centrifugation method as the most affordable and accurate one. Emulsions were evaluated according to the degree of destabilization, resistance and dispersion of the fat phase. The rheological characteristics of cheese masses were evaluated by the effective viscosity change.

In this paper we show how to efficiently implement parallel discrete simulations on Multi-Core and GPU architectures through a real example of application: a cellular automata model of laser dynamics. We describe the techniques employed to build and optimize the implementations using OpenMP and CUDA frameworks. We have evaluated the performance on two different hardware platforms that represent different target market segments: high-end platforms for scientific computing, using an Intel Xeon Platinum 8259CL server with48cores and also an NVIDIA Tesla6V100 GPU, both running on Amazon Web Server (AWS) Cloud, and on a consumer-oriented platform, using an Intel Core i9 9900k CPU and an NVIDIA GeForce GTX 1050 TI GPU. Performance results are compared and analysed in detail. We show that excellent performance and scalability can be obtained in both platforms, and we extract some important issues that imply a performance degradation for them. We also found that current Multi-Core CPUs with large core numbers can bring a performance very near to that of GPUs, even similar in some cases.

Discrete Double Hilbert Exponential sums along polynomials. This is preprint version (not peer reviewed) And we now check the proof of (ii) of Lemma 3.2 whether it is really true.

In this paper, we introduce two discrete degenerate random variables, namely the degenerate binomial and degenerate Poisson random variables. We deduce the expectations of the degenerate binomial random variables. We compute the generating function of the moments of the degenerate Poisson random variables, which leads us to define the new type degenerate Bell polynomials, and hence obtain explicit expressions for the moments of those random variables in terms of such polynomials. We also get the variances of the degenerate Poisson random variables. Finally, we illustrate two examples of the degenerate Poisson random variables.

Let be a power function $f_{r,M}(s)$ defined for every $s$ within the finite set $M$ as follows $$f_{r,M}(s)= \begin{cases} s^r, \ &s\in M,\\ 0, \ &\mathrm{otherwise}. \end{cases} $$ Let a discrete convolution of $f_{r,M}(s)$ be denoted as follows $\mathrm{Conv}_{r,M}[n]=(f_{r,M}*f_{r,M})[n]$. Let a real coefficients $A_{m,j}$ be given by the following recurrence $$ A_{m,j} = \begin{cases} 0, & \mathrm{if } \ j<0 \ \mathrm{or } \ j>m, \\ (2j+1)\binom{2j}{j} \sum_{d=2j+1}^{m} A_{m,d} \binom{d}{2j+1} \frac{(-1)^{d-1}}{d-j} B_{2d-2j}, & \mathrm{if } \ 0 \leq j < m, \\ (2j+1)\binom{2j}{j}, & \mathrm{if } \ j=m. \end{cases} $$ In this paper we show that for every $n>0$ the following odd-power identities involving coefficients $A_{m,j}$ and convolution transform $\mathrm{Conv}_{r,M}[n]$ hold $$ \begin{split} n^{2m+1}+1&=\sum_{r=0}^{m}A_{m,r}\mathrm{Conv}_{r,\mathbb{N}}[n],\\ n^{2m+1}-1&=\sum_{r=0}^{m}A_{m,r}\mathrm{Conv}_{r,\mathbb{Z}_{>0}}[n],\\ n^{2m+1}&=\sum_{r=0}^{m}A_{m,r}\sum_{k=1}^{n} k^r(n-k)^r\\ &=\sum_{r=0}^{m}A_{m,r}\sum_{k=0}^{n-1} k^r(n-k)^r. \end{split} $$

We study certain function algebras and their operator algebra completions on r-discrete abelian groupoids, the corresponding conditional expectations, maximal abelian subalgebras (masa) and eigen-functionals.

Understanding the dynamics of Organic Carbon mineralization is fundamental in forecasting biosphere to atmosphere Net Carbon Ecosystem Exchange (NEE). With this perspective, we developed 3D-CMCC-PSM, a new version of the hybrid Process Based Model 3D‐CMCC FEM where also heterotrophic respiration (R_h) is explicitly simulated. The aim was to quantify NEE as a forward problem, by subtracting Ecosystem Respiration (R_eco) to Gross Primary Productivity (GPP). To do so, we developed a simplification of the Soil Carbon dynamics routine proposed in DNDC [1]. The method calculates decomposition as a function of soil moisture, temperature, state of the organic compartments, and relative abundance of microbial pools. Given the pulse dynamics of soil respiration, we introduced modifications in some of the principal constitutive relations involved in phenology and littering sub-routines. We quantified the model structure related uncertainty in NEE, by running our training simulations over 1000 random parameter-sets extracted from parameters distributions expected from literature. 3D-CMCC-PSM predictability was tested on independent time series for 6 Fluxnet sites. The model resulted in daily and monthly estimations highly consistent with the observed time series. It showed lower predictability in Mediterranean ecosystems, suggesting that it may need further improvements in addressing evapotranspiration and water dynamics.

Although Shannon mutual information has been widely used, its effective calculation is often difficult for many practical problems, including those in neural population coding. Asymptotic formulas based on Fisher information sometimes provide accurate approximations to the mutual information but this approach is restricted to continuous variables because the calculation of Fisher information requires derivatives with respect to the encoded variables. In this paper, we consider information-theoretic bounds and approximations of the mutual information based on Kullback--Leibler divergence and Rényi divergence. We propose several information metrics to approximate Shannon mutual information in the context of neural population coding. While our asymptotic formulas all work for discrete variables, one of them has consistent performance and high accuracy regardless of whether the encoded variables are discrete or continuous. We performed numerical simulations and confirmed that our approximation formulas were highly accurate for approximating the mutual information between the stimuli and the responses of a large neural population. These approximation formulas may potentially bring convenience to the applications of information theory to many practical and theoretical problems.

Inverse Sumudu transform multiple shifting properties are used to design methodology for solving ordinary differential equations. Then algorithm applied to solve Whittaker and Zettl equations to get their new exact solutions and profiles which shown through Maple complex graphicals. Table of inverse Sumudu transforms for elementary functions given for supporting the differential equations solving using inverse Sumudu transform.

In recent years, unmanned aerial vehicles (UAVs) have become a valuable platform for many applications, including communication networks. UAV-enabled wireless communication faces challenges in complex urban and dynamic environments. UAVs can suffer from power limitations and path losses caused by non-line-of-sight connections, which may hamper better communication performance. To address these issues, reconfigurable intelligent surfaces (RIS) have been proposed as a helpful tool to enhance UAV communication networks. However, due to the high mobility of UAV, complex channel environments, and dynamic RIS configurations, it is challenging to estimate the link quality of ground users. In this paper, we propose a link quality estimation model using a gated recurrent unit (GRU) to assess the link quality of ground users for a multi-user RIS-assisted UAV-enabled wireless communication system. Our proposed framework uses a time series of user channel data and RIS phase shift information to estimate the quality of the link for each ground user. Simulation results showed that the proposed GRU model effectively and accurately estimates the link quality of the ground users in the RIS-assisted UAV-enabled wireless communication network.

Open source software such as OpenDSS has given a lot of help to distribution network researchers and educators. With high penetration of distributed renewable energy resources into distribution network, tradition distribution steady state analysis software such as OpenDSS is faced with difficulty in handling distributed generators. Three-phase distributed generators are often modeled in sequence frame while unbalanced distribution network are usually modeled in phase frame. So a load flow in sequence-phase coupled frame is proposed to handle models described in both frames. Voltage controlled DGs which are difficult to cope with in OpenDSS are handled in proposed program. The steady state analysis platform is programmed with open source Modelica language and the main aim of this paper is to introduce an open source platform for active distribution network steady analysis include load flow and short circuit analysis which can be easily adopted and improved by other educators and researchers.

In this research article we solve the problem of synchronization and anti-synchronization of chaotic systems described by discrete and time-delayed variable fractional order differential equations. To guarantee the synchronization and anti-synchronization of these systems, we use the well-known PID control theory and the Lyapunov-Krasovskii stability theory for discrete systems of variable fractional order.We illustrate the results obtained through simulation with examples, in which it can be seen that our results are satisfactory, thus achieving synchronization and anti-synchronization of chaotic systems of variable fractional order with discrete time delay.

Abstract Using the Pauli algebra version of the Dirac equation six symmetric but non-equivalent copies of the QED Lagrangian are constructed. The group of symmetries is the group of outer automorphisms of the finite quaternion group, isomorphic to the group of permutations on three letters. Solutions to corresponding Dirac equations form six sectors of fermion fields. The intersections between different sectors can only contain massless fields, implying that the sectors are distinct. If the symmetry could be extended to the full Standard Model Lagrangian it would predict the ratio of dark matter to ordinary matter of 5 : 1 which is close to the ratio of 5.2 : 1 according to current observations.

A stochastic dominance (SD) relation can be defined by two different perspectives: One from the view of distributions, and the other one from the view of expected utilities. In early days, Fishburn investigated SD from the view of distributions and we refer this perspective as Fishburn’s SD. One of his many results was the development of fractional-order SD for continuous distributions. However, discrete fractional-order SD may not be generalized directly since some properties of fractional calculus do not have a discrete counterpart. In this paper, we develop a discrete analogue of fractional-order SD from the view of distributions. We generalize the order of SD by Lizama’s fractional delta operator, show the preservation of SD hierarchy, and formulate the utility classes that are congruent with our SD relations. This work brings a message that some results of discrete SD cannot be generalized directly from continuous SD. We characterize the difference between discrete and continuous fractional-order SD, as well as the way to handle them.

The author demonstrates a stable Lagrangian solid modeling method, tracking the interactions of solid mass particles rather than using a meshed grid. This numerical method avoids the problem of tensile instability often seen with smooth particle applied mechanics by having the solid particles apply stresses expected with Hooke's law, as opposed to using a smoothing function for neighboring solid particles. This method has been tested successfully with a bar in tension, compression, and shear, as well as a disk compressed into a flat plate, and the numerical model consistently matched the analytical Hooke's law as well as Hertz contact theory for all examples. The solid modeling numerical method was then built into a 2-D model of a pressure vessel, which was tested with liquid water particles under pressure and simulated with smoothed particle hydrodynamics. This simulation was stable, and demonstrated the feasibility of Lagrangian specification modeling for fluid–solid interactions.

3D printing facilitates the straightforward construction of microchannels with complex three-dimensional architectures. Here, we demonstrate 3D-printed modular mixing components that operate on the basis of splitting and recombining fluid streams to decrease interstream diffusion length. These are compared to helical mixers that operate on the principle of chaotic advection.

In this paper we introduce a computational model for the simulation of hypervelocity impact (HVI) phenomena which is based on the Discrete Element Method (DEM). Our paper constitutes the first application of DEM to the modeling and simulating of impact events for velocities beyond 5 kms⁻¹. We present here the results of a systematic numerical study on HVI of solids. For modeling the solids, we use discrete spherical particles that interact with each other via potentials. In our numerical investigations we are particularly interested in the dynamics of material fragmentation upon impact. We model a typical HVI experiment configuration where a sphere strikes a thin plate and investigate the properties of the resulting debris cloud. We provide a quantitative computational analysis of the resulting debris cloud caused by impact and a comprehensive parameter study by varying key parameters of our model. We compare our findings from the simulations with recent HVI experiments performed at our institute. Our findings are that the DEM method leads to very stable, energy–conserving simulations of HVI scenarios that map the experimental setup where a sphere strikes a thin plate at hypervelocity speed. Our chosen interaction model works particularly well in the velocity range where the local stresses caused by impact shock waves markedly exceed the ultimate material strength.

An investigation of pair partial radial distribution functions and molecular pair potentials within a system has established that the existing potential functions are rooted in the assumption of a static arrangement of atoms, overlooking their distribution and vibration. In this study, Hill’s proposed radial distribution function polynomials are applied for the pure gaseous state to a binary liquid alloy system to derive the pair potential energy. The partial radial distribution functions of 36 binary liquid alloy systems reported in literature were examined and then validated using four thermodynamic models. Results show that the molecular interaction volume model (MIVM) and regular solution model (RSM) outperform other models when an asymmetric method is used to calculate the partial radial distribution function. The MIVM exhibits a mean standard deviation (SD) of 0.095 and a mean absolute relative difference (ARD) of 32.2%. Similarly, the RSM demonstrates a mean SD of 0.078 and an ARD of 32.2%. The Wilson model yields a mean SD of 0.124 and a mean ARD of 226%. The nonrandom two-liquid (NRTL) model exhibits a mean SD of 0.225 and a mean ARD of 911%. On applying the partial radial distribution function symmetry method, the MIVM and RSM outperform the other approaches, with a mean SD of 0.143 and a mean ARD of 165.9% for the MIVM. The RSM yields a mean SD of 0.117 and an ARD of 208.3%. The Wilson model exhibits the mean values of 0.133 and 305.6% for SD and ARD, respectively. The NRTL model shows an average SD of 0.200 and an average ARD of 771.8%. Based on this result, the influence of the symmetry degree on the thermodynamic model is explored by examining the symmetry degree as defined by the experimental activity curves of the two components.

The Fruit Harvest Helper, a mobile application developed by Northwest Nazarene University’s (NNU) Robotics Vision Lab, aims to assist farmers in estimating fruit yield for apple orchards. Currently, farmers manually estimate the fruit yield for an orchard, which is a laborious task. Fruit Harvest Helper seeks to simplify their process. While prior research efforts at NNU concentrated on developing an iOS app for blossom detection, this current research aims to adapt that smart farming application for apple detection across multiple platforms, iOS and Android. The old and new applications were designed with an intuitive user interface that is easy for farmers to use, allowing for quick image selection and processing. Unlike before, the adapted app utilizes a color ratio-based image segmentation algorithm implemented in OpenCV C++ to detect apples in apple tree images selected for processing. The results of testing the algorithm with a dataset of images indicate an 8.52% Mean Absolute Percentage Error (MAPE) and a Pearson correlation coefficient of 0.6 between detected and actual apples on the trees. These findings were obtained by evaluating the images from both the east and west sides of the trees, which was the best method to reduce the error of this algorithm. Although the Fruit Harvest Helper shows promise, there are many opportunities for improvement. These opportunities include exploring alternative machine-learning approaches for apple detection, conducting real-world testing without any human assistance, and expanding the app to detect various types of fruit. The Fruit Harvest Helper mobile application is among the many mobile applications contributing to precision agriculture, nearing readiness for farmers to use in yield monitoring and farm management.

Packaged water consumption has grown rapidly in urban areas of many low and middle income countries, but particularly in Ghana. However, the sources of water used by this growing packaged water industry and its implications for water resource management and transport-related environmental impacts have not been described. This study aimed to assess the spatial distribution of regulated packaged water production in Ghana, both in relation to demand and for natural mineral water, to hydrogeological characteristics. 764 addresses for premises licenced to produce packaged water from 2009-2015 were mapped and compared to regional sachet water consumption and beverage import/export data examined. We find evidence to suggest packaged water is transported shorter distances in Ghana than in developed countries. For natural mineral waters, producers should be able to address the most widespread water quality hazards (including high salinity, iron and nitrates) in aquifers used for production through reverse osmosis treatment. The study suggests there is scope to integrate beverage product and groundwater regulatory databases to support groundwater management.

Eutrophication, a global concern, impacts water quality, ecosystems, and human health. It’s crucial to monitor algal blooms in freshwater reservoirs, as they indicate the trophic condition of a waterbody through Chlorophyll-a (Chla) concentration. Traditional monitoring methods, however, are expen-sive and time-consuming. Addressing this hindrance, we developed models using remotely sensed data from the Sentinel-2 satellite for large-scale coverage, including its bands and spectral indexes, to estimate the Chla concentration on 149 freshwater reservoirs in Ceará, Brazil. Several machine learning models were trained and tested, including k-nearest neighbours, random forests, extreme gradient boosting, the least absolute shrinkage, group method of data handling (GMDH), and sup-port vector machine models. A stepwise approach determined the best subset of input parameters. Using a 70/30 split for the training and testing datasets, the best-performing model was the GMDH, achieving an R2 of 0.91, MAPE of 102.34%, and RMSE of 20.38 g/L, which are values consistent with the ones found in the literature. Nevertheless, the predicted Chla concentration values were most sensitive to the red, green, and near infra-red bands.

This study advances the field of solar irradiance nowcasting by introducing a discrete-level classification approach, diverging from traditional continuous measurement methods. Grounded in a need for cost-effective and modular solutions, the research employs high-resolution computer vision coupled with a deep learning framework to predict Direct Normal Irradiance (DNI). By harnessing the capabilities of a Logitech C900 1080p HD camera and an NVIDIA Jetson module, the study achieves real-time data processing, pivotal for CSP systems' operational efficiency. The core of the methodology is the ResNet-50 convolutional neural network, refined via transfer learning on a bespoke dataset, culminating in a predictive accuracy of 85.78%. This discrete classification model contrasts with conventional, costly instruments like the MS-57, offering a novel and accessible alternative for DNI estimation. Such innovation not only demonstrates high predictive accuracy but also signifies a shift towards less resource-intensive and more adaptable solar energy forecasting tools, contributing a significant leap towards optimizing renewable energy systems.

By applying a running average (with a window-size= d), we could transform Discrete data to broad-range, Continuous values. When we have more than 2 columns and one of them is containing data about the tags of classification (Class Column), we could compare and sort the features (Non-class Columns) based on the R2 coefficient of the regression for running averages. The parameters tuning could help us to select the best features (the non-class columns which have the best correlation with the Class Column). “Window size” and “Ordering” could be tuned to achieve the goal. this optimization problem is hard and we need an Algorithm (or Heuristics) for simplifying this tuning. We demonstrate a novel heuristics, Called Simulated Distillation (SimulaD), which could help us to gain a somehow good results with this optimization problem.

Pit craters, and pit crater chains, are now recognised as being an important part of the surface morphology and structure of many planetary bodies, and are particularly remarkable on Mars. Pit craters do not possess the elevated rims, ejecta deposits, or other features that are typically associated with impact craters. They are thought to arise from the drainage/collapse of a relatively weak surficial material into an open (or widening) void in a much stronger material below. The creation of such voids has been suggested to be due to extensional fracturing/dilational faulting, shallow dike intrusion, lava tube collapse amongst other hypotheses. These craters have a very distinctive expression, often presenting funnel, cone, or bowl-shaped geometries. Analogue models of pit crater formation provide a map-view picture of their initiation and evolution but give little insight into their internal structure or geometry, but produce pits that typically have steep, nearly conical cross sections. Numerical modelling studies of their formation have been limited and have produced some quite interesting, but nonetheless puzzling, results whereby the simulated pit craters had generally convex (steepening downward) slope profiles with no distinct rim; quite unlike many of those observed on Earth or on Mars. To address these issues, I present here a high-resolution, 2D discrete element model of weak cover (regolith) collapse into either a static or a widening underlying void. I use frictional and frictional-cohesive discrete elements to represent a range of probable cover rheologies. Under Martian gravitational conditions, frictional-cohesive and frictional materials produce cone, bowl and scoop-shaped pit craters. For a given cover thickness, the specific crater shape depends on the amount of underlying void space created for drainage. When void space is small relative to cover thickness, craters have bowl or scoop-shaped geometries. In contrast, when void space is large relative to cover thickness, craters have cone-shaped geometries with essentially planar (nearing angle of repose) slope profiles. Frictional-cohesive materials exhibit more distinct rims than simple frictional materials and thus may reveal some stratigraphic layering on the pit crater walls. In the limit, when drainage from the overlying cover is insufficient to fill the underlying void, ´skylights´ into the deeper structure are created. Implications of these results for the interpretation of pit craters on Earth, Mars, other planets and moons are discussed.

In this study, we consider delta fractional Sturm--Liouville (DFSL) initial value problems in the sense of delta Caputo and delta Riemann-Liouville (R--L) operators. One of the properties of delta fractional difference operators which makes it different from nabla counterpart is to shift its domain. This feature makes it more complex than the nabla fractional operator. We obtain sum representation of solutions for DFSL initial value problems with the help of $\mathcal{Z}-$ transformation. Moreover, we get analytical solutions of homogeneous DFSL problem within Riemann-Liouville (R--L) and Caputo sense, discrete Sturm--Liouville (DSL) problem, continuous fractional Sturm--Liouville (FSL) problem in the sense of R--L and Caputo operators, and continuous Sturm--Liouville (SL) differential problem. From this point of view, we compare all the solutions with each other. Consequently, we show that all results for these four eigenvalue problems are compatible with each other and approach to each other while the orders tends to one, i.e. $\Delta^{\mu }\left( \Delta x\left( t-\mu \right) \right)\cong D_{0^{+}}^{\mu }\left( x^{\prime }\left( t\right) \right)\cong \Delta^2x(n-1) \cong x^{\prime \prime }\left( t\right) =\lambda x\left( t\right),\ \mu\rightarrow1 $ . We support our results comparatively by tables and simulations in detail.

Optimal dividend payment under a ruin constraint is a two objective control problem which—in simple models—can be solved numerically by three essentially different methods. One is based on a modified Bellman equation and the policy improvement method (see Hipp C., 2003). In this paper we use explicit formulas for running allowed ruin probabilities which avoid a complete search and speed up and simplify the computation. The second is also a policy improvement method, but without the use of a dynamic equation (see Hipp C., 2003). It is based on closed formulas for first entry probabilities and discount factors for the time until first entry (see Hipp C., 2016). Third a new, faster and more intuitive method which uses appropriately chosen barrier levels and a closed formula for the corresponding dividend value. Using the running allowed ruin probabilities, a simple test for admissibility—concerning the ruin constraint—is given. All these methods work for the discrete De Finetti model and are applied in a numerical example. The non stationary Lagrange multiplier method suggested in (see Hipp C., 2016), Section 2.2.2 does also yield optimal dividend strategies which differ from those in all other methods, and Lagrange gaps are present here. These gaps always exist in De Finetti models, see (see Hipp C., 2017).

High Efficiency Video Coding (HEVC) uses an 8-point filter and a 7-point filter, which are based on the discrete cosine transform (DCT), for the 1/2-pixel and 1/4-pixel interpolations, respectively. In this paper, discrete sine transform (DST)-based interpolation filters (IF) are proposed. The first proposed DST-based IFs (DST-IFs) use 8-point and 7-point filters for the 1/2-pixel and 1/4-pixel interpolations, respectively. The final proposed DST-IFs use 12-point and 11-point filters for the 1/2-pixel and 1/4-pixel interpolations, respectively. These DST-IF methods are proposed to improve the motion-compensated prediction in HEVC. The 8-point and 7-point DST-IF methods showed average BD-rate reductions of 0.7% and 0.3% in the random access (RA) and low delay B (LDB) configurations, respectively. The 12-point and 11-point DST-IF methods showed average BD-rate reductions of 1.4% and 1.2% in the RA and LDB configurations for the Luma component, respectively.

Time estimation was investigated in 24 healthy adults, 12 women and 12 men, before and after an exhaustive exercise. We compared the capability of estimating time intervals in the range from 1 to 5 s using tasks requiring mental count and tasks that did not allow the mental count. Time estimation and blood lactate levels were evaluated before and at the end of the exercise. We found that the perception of time intervals between 1 and 5 s was affected at the end of the exercise. The observed effects, associated with a significant increase in blood lactate levels, were different in the two types of time estimation used in the present study. In experimental conditions in which participants had to evaluate duration of the time interval by mental counting, at the end of exercise a significant reduction in the overestimation of time made at rest was observed. On the other hand, when participants had to assess the difference in duration between two events without the possibility to mentally count, at the end of the exercise a significant deterioration in performance was observed. In both cases, no gender differences were seen.

This paper is concerned with the optimal linear quadratic Gaussian (LQG) control problem for discrete time-varying system with multiple input delays and multiplicative noises. The main contributions are two-fold. Firstly, when the state variables can be observed exactly, we obtain a necessary and suffcient condition for the multiple-delays system in terms of the non-homogeneous relationship between the state and costate, which is the solution to the coupled forward and backward stochastic difference equations. Secondly, when the state variables are partially observed, we derive a suboptimal linear output feedback controller for the discrete-time system based on the obtained results of the optimal LQG control. Numerical examples are shown to illustrate the proposed algorithm.

The celebrated geodesic congruence equation of Raychaudhuri, together with the resulting singularity theorems of Penrose and Hawking that it enabled, yield a highly general set of conditions under which a spacetime (or, more generically, a pseudo-Riemannian manifold) is expected to become geodesically incomplete. However, the proofs of these theorems traditionally depend upon a collection of assumptions about the continuum spacetime (and, in the physical case, the stress-energy distribution defined over it), including its global structure, its energy conditions, the existence of trapped null surfaces and the various volume/intersection properties of geodesic congruences, that are inherently difficult to translate to the case of discrete spacetime formalisms, such as causal set theory or the Wolfram model. Some, such as the discrete analog of Raychaudhuri's equation for the volumes of geodesic congruences, are subtle to formulate due to intrinsic differences in the behavior of discrete vs. continuous geodesics; for others, such as the definition of a trapped null surface, no appropriate translation is known for the general discrete case due to a lack of a priori coordinate information. It is therefore a non-trivial question to ask whether (and to what extent) there exist equivalently general conditions under which one expects discrete spacetimes to become geodesically incomplete, and how these conditions might differ from those in the continuum. This article builds upon previous work, in which the conformal and covariant Z4 (CCZ4) formulation of the Cauchy problem for the Einstein field equations, with constraint-violation damping, was defined in terms of Wolfram model evolution over discrete (spatial) hypergraphs for the case of vacuum spacetimes, and proceeds to consider a minimal extension to the non-vacuum case by introducing a massive scalar field distribution, defined in either spherical or axial symmetry. Under appropriate assumptions, this scalar field distribution admits a physical interpretation as a collapsing (and, in the axially-symmetric case, uniformly rotating) dust, and we are able to show, through a combination of rigorous mathematical analysis and explicit numerical simulation, that the resulting discrete spacetimes converge asymptotically to either non-rotating Schwarzschild black hole solutions or maximally-rotating (extremal) Kerr black hole solutions, respectively. Although the assumptions used in obtaining these preliminary results are very strong, they nevertheless offer hope that a more general, perhaps ultimately ``Penrose-like'', singularity theorem may be provable in the discrete spacetime case too.

This addendum contains clarifications and a sharpening of some of the results on the VIT transform framework developed in [1]. The focus is on the right-coefficient and left- coefficient forms of the transform, the extraction of a first-order term from a left polynomial fraction, and the application to linear time-varying systems.

High levels of meat consumption are increasingly being criticised for ethical, environmental, and social reasons. Plant-based meat substitutes have been identified as healthy sources of protein that, in comparison to meat, offer a number of social, environmental and health benefits and may play a role in reducing meat consumption. However, there has been a lack of research on the role they can play in the policy agenda and how specific meat substitute attributes can influence consumers to replace partially replace meat in their diets.

1) Landsat operational land imager (OLI) data and consequent laboratory measurements were used to predict Chlorophyll-a (Chl-a) concentration and the trophic states for an inland lake within the East Kolkata Wetland, India; 2) The most suitable band ratio was identified by performing Pearson correlation analysis between Chl-a concentrations and possible OLI band and band ratios from the study points; 3) The results showed highest correlation coefficient from the band ratio OLI5/OLI4 with an R value of 0.85. The prediction model was then developed by applying regression analysis between the band ratio OLI5/OLI4 and Chl-a concentration of the study points. The reflectance ratios of the validation points were given as input on the prediction model and the model output was considered as predicted Chl-a values of the validation points to check the efficiency of the prediction model. The regression model between laboratory-derived Chl-a value and model-fitted Chl-a value of the validation points revealed a high correlation with an R² value of 0.78. Trophic State Index (TSI) of the lake was also calculated from laboratory-derived Chl-a value and model-fitted Chl-a value of the validation points. The study presented a high correlation of TSI determined from predicted data with TSI from laboratory reference data (R = 0.88). The TSI values of the lake ranged from 65 to 75 which indicate that the lake is appeared to be eutrophic to hypereutrophic conditions. 4) This empirical study showed that Landsat 8 OLI imagery can be effectively applied to estimate Chl-a levels and trophic states for inland lakes.

Medical image segmentation is a critical step in various healthcare applications, aiding in diagnosis, treatment planning, and disease monitoring. In this study, we investigate the efficacy of a segmentation approach based on the 2D wavelet transform. Leveraging a dataset comprising 10 diverse medical images, we evaluate the performance of our segmentation method using three key metrics: accuracy, precision, and recall. Our findings demonstrate that the proposed approach enhances segmentation accuracy, offering promising results compared to existing methods. By harnessing the multi-resolution feature extraction capabilities of the 2D wavelet transform, our method achieves improved delineation of medical image structures, paving the way for more accurate and efficient healthcare interventions.

Covariant integral quantization with rotational symmetry SO(3) is established for the quantum motion on this group manifold, or, alternatively, for Gabor signal analysis on this group. We revisit the action of the related (non-unitary) Weyl-Gabor operator on the Hilbert space of square integrable functions on SO(3) and disclose a set of various properties. By selecting weight functions on the corresponding discrete-continuous phase space, one derives related coherent states and Wigner transform.

The deployment of a manufacturing execution system (MES) holds promising potential in facilitating the accumulation of a substantial amount of inspection data. Low quality levels in discrete manufacturing environments are the result of multi-factor coupling and failure to find quality issues in a timely manner within manufacturing settings may trigger the propagation of defects downstream. Currently, most of the inspection quality methods are direct measurement followed by manual judgment. The integration of deep learning methods provides a feasible way to identify defects in a timely manner, thus improving the acceptance rate of factories. This paper focuses on the design of a data-driven quality prediction and control model around discrete manufacturing characteristics, and use fuzzy theory to evaluate the quality level of production stages. Building Multivariate Long and short-term memory learning hidden quality representations to extract predictions from multi-level inspection data in manufacturing systems. Finally, by validating the data of actual produced water dispensers according to three evaluation indexes, RMSE, MAE, MAPE, the results show that Multivariate Long and short-term memory has better prediction performance.

Stenting is a common method for treating atherosclerosis. A metal or polymer stent is deployed to open the stenosed artery or vein. After the stent is deployed, the blood flow dynamics influence the mechanics by compressing and expanding the structure. If the stent does not respond properly to the resulting stress, vascular wall injury or re-stenosis can occur. In this work, Discrete Multiphysics is used to study the mechanical deformation of the coronary stent and its relationship with the blood flow dynamics. The major parameters responsible for deforming the stent are sort in terms of dimensionless numbers and a relationship between the elastic forces in the stent and pressure forces in the fluid is established. The blood flow and the stiffness of the stent material contribute significantly to the stent deformation and affect the rate of deformation. The stress distribution in the stent is not uniform with the higher stresses occurring at the nodes of the structure.

We present progress in fast, high-resolution imaging, material classification, and fault detection using hyperspectral X-ray measurements. Classical X-ray CT approaches rely on data from many projection angles, resulting in long acquisition and reconstruction times. Additionally, conventional CT cannot distinguish between materials with similar densities. However, in additive manufacturing, the majority of materials used are known a priori. This knowledge allows to vastly reduce the data collected and increase the accuracy of fault detection. In this context, we propose an imaging method for non-destructive testing of materials based on the combination of spectral X-ray CT and discrete tomography. We explore the use of spectral X-ray attenuation models and measurements to recover the characteristic functions of materials in heterogeneous media with piece-wise uniform composition. We show by means of numerical simulation that using spectral measurements from a small number of angles, our approach can alleviate the typical deterioration of spatial resolution and the appearance of streaking artifacts.

With the widespread use of new equipment such as distributed photovoltaics, distributed energy storage, electric vehicles, and distributed wind power, the control of low-voltage distribution networks (LVDNs) has become increasingly complex. Acquiring the most recent topological structure is essential for conducting accurate analysis and real-time control of LVDNs. The signal-injection-based topology identification algorithm is favored for its speed and efficiency. This study proposes a new signal-injection-based topology identification algorithm aimed at solving the issues of reduced completeness and correctness of identification caused by missing feature signal records (FSRs). By analyzing the correlations between FSRs of different positions, the algorithm employs vertical and horizontal completion methods to effectively address missing data, combined with inclusion detection algorithms, to identify the topology of LVDNs. Based on the study of actual LVDN data, the results indicate that the algorithm not only significantly improves the completeness and correctness of topology identification for LVDNs but can also increase the number of devices in operation by evaluating the status of topological equipment, further enhancing the topology identification performance.

Wireless communication has proliferated due to the colossal augmentation of smart phones, mobile phones, tablets and laptops etc. Multimedia applications such as live videos, audios, still images, animated graphics and live communications all requires accompaniments of quality of service (QoS) as well as reliable end-to-end transfer of data. Information about the availability of resources in channel is must in order to improve the QoS. QoS can be determined in terms of capacity of channel, available bandwidth (ABW), and bulk transfer capacity (BTC). Bandwidth availability is a key check for improving the QoS in a network. The performance of a multimedia application is directly affected by the bandwidth availability. One of the most important QoS characteristic is ABW at a wireless route and it can be demarcated as least unused capacity of links instituting a network route. Since now there have been many bandwidth estimation techniques are available in the literature to increase the network performance. Bandwidth estimation techniques have been arranged into three leading classifications: i) Active probing bandwidth estimation techniques, ii) Passive bandwidth estimation technique and iii) Model based bandwidth estimation techniques. Each of these techniques are briefly discussed in this paper.

The theory of the continuous two-dimensional (2D) Fourier Transform in polar coordinates has been recently developed but no discrete counterpart exists to date. In this paper, we propose and evaluate the theory of the 2D discrete Fourier Transform (DFT) in polar coordinates. This discrete theory is shown to arise from discretization schemes that have been previously employed with the 1D DFT and the discrete Hankel Transform (DHT). The proposed transform possesses orthogonality properties, which leads to invertibility of the transform. In the first part of this two-part paper, the theory of the actual manipulated quantities is shown, including the standard set of shift, modulation, multiplication, and convolution rules. Parseval and modified Parseval relationships are shown, depending on which choice of kernel is used. Similar to its continuous counterpart, the 2D DFT in polar coordinates is shown to consist of a 1D DFT, DHT and 1D inverse DFT.

The stock beta coefficient literature extensively discusses the proper methods for the estimation of beta as well as its use in asset valuation. However, there are relatively few references with respect to the appropriate time horizon that investors should utilize when evaluating the risk-return relationship of a stock. We examine the appropriate time horizon for beta estimation differentiating our results by sector according to the Industry Classification Benchmark. We employ data from the NYSE and we estimate varying lengths of beta employing data from 30 to 250 trading days. The constructed beta series is then examined for the presence of breaks using the endogenous structural break literature. Results show evidence against the use of betas that employ more than 90 trading days of data provisional to the sector under study.

A three-year survey was conducted to estimate the incidence of grapevine trunk diseases (GTDs) in Greece and identify fungi associated with the disease complex. In total, 310 vineyards in different geographical regions in northern, central, and southern Greece were surveyed, and 533 fungal strains were isolated from diseased vines. Morphological, physiological and molecular (5.8S rRNA gene-ITS sequencing) analyses revealed that isolates belonged to 35 distinct fungal genera, including well-known (e.g. Botryosphaeria sp, Diaporthe spp., Eutypa sp., Diplodia sp., Fomitiporia sp., Phaeoacremonium spp., Phaeomoniella sp.) and lesser-known (e.g. Neosetophoma sp., Seimatosporium sp., Didymosphaeria sp., Kalmusia sp.) grapevine wood inhabitants. GTDs-inducing population structure differed significantly among the discrete geographical zones. Phaeomoniella chlamydospora (26.62%, n=70), Diaporthe spp. (18.25%, n=48) and F. mediterranea (10.27%, n=27) were the most prevalent in Heraklion, whereas D. seriata, Alternaria spp., P. chlamydospora and Fusarium spp. were predominant in Nemea (central Greece). In Amyntaio and Kavala (northern Greece), D. seriata was the most frequently isolated species (&gt;50% frequency). Multi-genes (rDNA-ITS, LSU, tef1-α, tub2, act) sequencing of selected isolates, followed by pathogenicity tests revealed that Neosetophoma italica, Seimatosporium vitis, Didymosphaeria variabile and Kalmusia variispora caused wood infection with the former being the most virulent. To the best of our knowledge, this is the first report of N. italica associated with GTDs worldwide. This is also the first record of K. variispora, S. vitis and D. variabile associated with wood infection of grapevine in Greece. The potential associations of disease indices with vine age, cultivar, GTD-associated population structure and the prevailed meteorological conditions in different viticultural zones in Greece are presented and discussed.

Paratemporal methods based on tree expansion have proven to be effective in efficiently generating trajectories of stochastic systems. However, combinatorial explosion of branching arising from multiple choice points presents a major hurdle that must be overcome to implement such techniques. In this paper we tackle this scalability problem by developing a systems theory-based framework covering both conventional and proposed tree expansion algorithms for speeding up discrete event system stochastic simulations while preserving desired accuracy. An example is discussed to illustrate the tree expansion framework in which a discrete event system specification (DEVS) Markov stochastic model takes the form of a tree isomorphic to a free monoid over the branching alphabet. We derive the computation times for baseline, non-merging, and merging tree expansion algorithms to compute the distribution of output values at any given depth. The results show the remarkable reduction from exponential to polynomial dependence on depth effectuated by node merging. We relate these results to the similarly reduced computation time of binomial coefficients underlying Pascal’s triangle. Finally, we discuss application of tree expansion to estimating temporal distributions in stochastic simulations involving serial and parallel compositions with potential real world use cases.

To investigate the variation patterns of spectral characteristics of main branches in fruit trees due to uneven distribution of growth, pruning, and fruit quality, as well as to determine the optimal fruit vi-bration harvesting, a vibration differential model of tree branches was constructed using the concen-trated mass method. A three-dimensional model of the fruit tree was built using finite element software. The spectral characteristics of each branch were obtained using the linear sweep frequency method. The method of changing the mass of each branch was employed to simulate the variation in branch mass. The frequency of occurrence, mean acceleration, and coefficient of variation of the inherent frequencies were used as evaluation indicators. The results revealed that when the vibration frequency of the fruit tree was at the peak frequency in the tree trunk spectrum, the frequency of occurrence and mean accel-eration were optimal. And the coefficient of variation of branch acceleration was less than 1, indicating a moderate dispersion. Additionally, reducing the branch mass resulted in an increase in the peak fre-quency of the tree trunk spectrum. Therefore, the growth, pruning, and uneven distribution of the branches in fruit trees can alter the spectral characteristics. When harvesting fruit by vibration, selecting the peak frequency in the trunk vibration spectrum as the excitation force frequency yields the best fruit detachment.