
The milling of wheat flour is a process that has existed since ancient times. In the course of history, the techniques have improved, the equipment modernized. The interest of the miller in charge of the mill is still to ensure that a mill is functional and profitable, as well as to provide a consistent quality of flour. The production of organic flour means that methods of adding chemicals and unnatural agents are not possible. In organic flour production, it is necessary to work with the raw material. A grain of wheat is a living material, and its quality varies according to a multitude of factors. Challenges are therefore present at each stage of the value chain. The use of artificial intelligence techniques offers solutions and new perspectives to meet the different objectives of the miller. A literature review of artificial intelligence techniques developed at each stage of the value chain surrounding the issues of quality and yield is conducted. An analysis of a large number of variables, including process factors, process parameters and wheat grain quality from data collected on the value chain enables the development and training of artificial intelligence models. From these models, it is possible to develop decision support tools and optimize the wheat flour milling process. Several major research directions, other than constant quality, are to be studied to optimize the process and move towards a smart mill. This includes energy savings, resource optimization and mill performance.
Citation: Loïc Parrenin, Christophe Danjou, Bruno Agard, Robert Beauchemin. Future trends in organic flour milling: the role of AI[J]. AIMS Agriculture and Food, 2023, 8(1): 48-77. doi: 10.3934/agrfood.2023003
[1] | Sungwon Kim, Meysam Alizamir, Youngmin Seo, Salim Heddam, Il-Moon Chung, Young-Oh Kim, Ozgur Kisi, Vijay P. Singh . Estimating the incubated river water quality indicator based on machine learning and deep learning paradigms: BOD5 Prediction. Mathematical Biosciences and Engineering, 2022, 19(12): 12744-12773. doi: 10.3934/mbe.2022595 |
[2] | Cong Zhao, Na Deng . An actor-critic framework based on deep reinforcement learning for addressing flexible job shop scheduling problems. Mathematical Biosciences and Engineering, 2024, 21(1): 1445-1471. doi: 10.3934/mbe.2024062 |
[3] | Ruiping Yuan, Jiangtao Dou, Juntao Li, Wei Wang, Yingfan Jiang . Multi-robot task allocation in e-commerce RMFS based on deep reinforcement learning. Mathematical Biosciences and Engineering, 2023, 20(2): 1903-1918. doi: 10.3934/mbe.2023087 |
[4] | Shixuan Yao, Xiaochen Liu, Yinghui Zhang, Ze Cui . An approach to solving optimal control problems of nonlinear systems by introducing detail-reward mechanism in deep reinforcement learning. Mathematical Biosciences and Engineering, 2022, 19(9): 9258-9290. doi: 10.3934/mbe.2022430 |
[5] | Zilong Zhuang, Zhiyao Lu, Zizhao Huang, Chengliang Liu, Wei Qin . A novel complex network based dynamic rule selection approach for open shop scheduling problem with release dates. Mathematical Biosciences and Engineering, 2019, 16(5): 4491-4505. doi: 10.3934/mbe.2019224 |
[6] | Yangjie Sun, Xiaoxi Che, Nan Zhang . 3D human pose detection using nano sensor and multi-agent deep reinforcement learning. Mathematical Biosciences and Engineering, 2023, 20(3): 4970-4987. doi: 10.3934/mbe.2023230 |
[7] | Zichen Wang, Xin Wang . Fault-tolerant control for nonlinear systems with a dead zone: Reinforcement learning approach. Mathematical Biosciences and Engineering, 2023, 20(4): 6334-6357. doi: 10.3934/mbe.2023274 |
[8] | Hong Zhang, Zhenchao Xu, Yunxiang Wang, Yupeng Shen . An innovative parameter optimization of Spark Streaming based on D3QN with Gaussian process regression. Mathematical Biosciences and Engineering, 2023, 20(8): 14464-14486. doi: 10.3934/mbe.2023647 |
[9] | Jianguo Duan, Mengting Wang, Qinglei Zhang, Jiyun Qin . Distributed shop scheduling: A comprehensive review on classifications, models and algorithms. Mathematical Biosciences and Engineering, 2023, 20(8): 15265-15308. doi: 10.3934/mbe.2023683 |
[10] | Van Dong Nguyen, Dinh Quoc Vo, Van Tu Duong, Huy Hung Nguyen, Tan Tien Nguyen . Reinforcement learning-based optimization of locomotion controller using multiple coupled CPG oscillators for elongated undulating fin propulsion. Mathematical Biosciences and Engineering, 2022, 19(1): 738-758. doi: 10.3934/mbe.2022033 |
The milling of wheat flour is a process that has existed since ancient times. In the course of history, the techniques have improved, the equipment modernized. The interest of the miller in charge of the mill is still to ensure that a mill is functional and profitable, as well as to provide a consistent quality of flour. The production of organic flour means that methods of adding chemicals and unnatural agents are not possible. In organic flour production, it is necessary to work with the raw material. A grain of wheat is a living material, and its quality varies according to a multitude of factors. Challenges are therefore present at each stage of the value chain. The use of artificial intelligence techniques offers solutions and new perspectives to meet the different objectives of the miller. A literature review of artificial intelligence techniques developed at each stage of the value chain surrounding the issues of quality and yield is conducted. An analysis of a large number of variables, including process factors, process parameters and wheat grain quality from data collected on the value chain enables the development and training of artificial intelligence models. From these models, it is possible to develop decision support tools and optimize the wheat flour milling process. Several major research directions, other than constant quality, are to be studied to optimize the process and move towards a smart mill. This includes energy savings, resource optimization and mill performance.
In the case of the species are under the threat of extinction, scientific researches have been focused on evaluating the possibility on the population of the species will grow back to its previous size. Allee effect for social mammals was first introduced by [1], but since then the concept has been extended to a wide range of organisms. The Allee effect explains a situation in which low numbers of species are influenced by a positive interaction between population growth rate and size, which raises their probability of extinction. On this subject, a lot of attention has been paid to the corresponding spatial dispersal versions of predator-prey systems with strong Allee effect, and the spatiotemporal pattern formation has been studied in such systems. Under the influence of the Allee effect, the dynamics of the system can significantly change as chaotic patches or traveling population pulses [2,3]. Recently, pattern formation of a diffusive predator-prey model under the effect of strong Allee effect [4], pattern dynamics in prey-predator system is also focused by [5,6]. Most of the current trend of literature leaves the issue of the process of prevention when the population experiences extinction. Petrovskii et al. [7] also studied the Allee effect as a possible lagging mechanism between introduction and establishment. They found that, for a particular set of parameters, there are three forms of regimes when the invasive prey affected by the Allee effect. One of them is ‘regimes of anomalous extinction’. As a comparison to this extinction, it is shown in this study that extinction is delayed by a gradient that contributes to the linear term of predator mortality.
It can be mentioned in some studies about population dynamics that highlight the possibility of extinction or invasion of a particular species [8,9,10] and/or possible changes in the population structure [11,12] were established. With the development of the gradient analysis [13,14], spatial patterns have been reported in some literature [15,16,17,18,19,20]. In another study, this issue is considered by taking the diffusion coefficient as the spatial gradient for the predator-prey system [19]. In particular, herbivore dynamics are examined under the effect of gradient [20]. Taking the spatial gradient as one of the system parameters, i.e., the rate of growth of the prey, is described by [21]. However, the literature covers considerable progress on the issue of the spatial gradient. The effect of the spatial gradient on the system parameters to avoid an approaching extinction is ignored. Consequently, this issue remains somewhat unclear when the environmental conditions are not acceptable for sustainable species. Hence, the species move through space from unsuitable conditions to suitable one to survive. From this point of view, the increase in temperature caused by climate change and therefore the habitat becomes inappropriate for the organisms' survival forces them to migrate and adapt to new environments. This ecological issue constitutes the main purpose of this study.
Species challenge against to get rid of extinction. For this reason, species prefer to adapt instinctively or move to a new habitat to improve resistance against extinction due to their survival instincts. We can safely say that habitat change and beachhead are considered as a response to climate change. The new discovery by Jonkers et al. [22] has revealed the results of climate change on the aquamarine life in the oceans. This research has shown that climate change leads to migration of marine zooplankton communities in consequence of 12 years of research by [23]. However, we can also give the paper studied by Pecl et al. [24] as an evidence. They suggested that the first response of species is to shift their location to climate change to the cooler or warmer locations where the species last their life in a more suitable conditions for preventing extinction.
A large number of theoretical models have been proposed to study the spatial pattern of population dynamics and interactions [4,5,6,25,26,27,28]. Although considerable progress has been made in understanding the spatial distribution of many species, some of the important issues remain obscure. One of these important issues is the gradient effect on systems' spatial distributions. There are several works in the literature on the gradient effect of spatial patterns. The continuous distribution of species through environmental gradient is described with the development of the gradient analysis [13,14]. It is stated that spatially heterogeneous patterning can be described by the spatial gradient response of individuals [16,17,18]. Pascual found the spatial gradient to be taken into account both in terms of prey and predator diffusion in [19]. At the same rate of diffusion, the distribution of species is concentrated along the spatial gradient [29]. In particular, Post [20] focuses on the spatial gradient effect on herbivorous population dynamics. Geographical patterns and ecological speciation processes across environmental gradients of spatially structured populations are described in [15]. The spatial gradient of the prey growth rate is studied to demonstrate the stability of the system in [21]. Apart from the several apparent works, the spatial gradient continues to be literally an issue of vital importance.
The current literature on the gradient impact on system dynamics does not pay enough attention to the event that species face extinction, while environmental conditions for sustainable species survival are not sufficient. For this reason, this work revisits the prey-predator model where prey growth is dampened by the strong Allee effect introduced in [7] with the intention of further establishing the species persistence/extinction issue. The present paper extends that study to the case in which the mortality rate of the predator is spatially varying; that is, the mortality rate is given by a spatial function as a response to climate change.
In particular, the biological relevance here is that if individuals threatened by unsuitable environmental conditions are foreseen to result in extinction, whether species continue to diffuse along a spatial gradient in the direction of unsuitable conditions, e.g., more resources, suitable climate conditions, and less predation, etc. Therefore, the focus of this work about concerning this issue, the following scenarios are concentrated: whether the extinction of the species can be prevented or at least postponed to a distant future. In the case of extinction, if it happens, is there any warning before extinction? Obtained numerical results show that when the system is on the verge of extinction the pattern of spatio-temporal oscillations becomes completely regular with fixed period and amplitude. This can serve as an early warning [30,31,32].
The issues mentioned above are addressed by mathematical modelling. This begins with the question, even if the decreasing level of predator mortality hampers the growth of prey whether the spatial gradient has the power to prevent extinction. Therefore, the starting point of the predator mortality rate is taken where the system has an extinction without gradient case. Importantly, as an early warning signal, the system produces a regular distribution when the system is near to its extinction point which is why the obtained results are compatible with our previous climate change work performed on the oxygen-plankton model system [31,33]. On the other hand, one of the most important stages of species movement is the ‘beachhead’, i.e., the establishment of an invader as a self-sustaining community [34]. In this sense, the prey predator system [7] is redefined by expanding the system with gradient and beachhead, taking this change in the system as a result of rising temperatures, i.e., climate change.
In brief, in this work a diffusion-reaction model describing predator-prey interactions are studied. Gradient and beachhead are used as precautions when the species faces extinction by adding a spatial gradient to the mortality rate of the predator. Based on extensive numerical simulations and analytical observations, it is revealed that the spatial structure gives some early warning signals, i.e., regular distribution when the predator mortality rate is closer to the extinction point.
In light of the above, the paper is structured as follows. In Section 2, the predator-prey model system is revisited, suggesting that the predator mortality rate decreases linearly. In Section 3, the spatial system is analyzed in one-dimensional and two-dimensional case through comprehensive numerical simulations focused on the system parameter space. The obtained numerical results show that species' extinction is followed by regular pattern formation in the spatial system; see a similar pattern as a response to climate change in [31,32,35]. System persistence/extinction map is given in parameter plane (x1,ω) for different domains correspond to species' extinction and persistence. In particular, it is established that for the different choices of initial distributions for prey and predator, the spatial distribution of the system does not change and retains its extinction point where the entire system is extinct. The most important observation of this work is that the movement of species as a precaution to change in climate, i.e., addition of a spatial gradient and beachhead, can not stop the extinction. But this movement in space slow down the approaching extinction so that the species can continue to exist for a while. Finally, in the last section, the potential significance of this work is discussed.
We consider one dimension predator prey model, studied earlier in [7,36,37,38]:
∂h(X,T)∂T=d1∂2h∂X2+F(h)−f(h,p), | (2.1) |
∂p(X,T)∂T=d2∂2p∂X2+κf(h,p)−μp. | (2.2) |
where h is the densities of prey and p the densities of predator at moment T and position X [39,40,41,42]. d1 and d2 are the diffusion coefficients and the value of κ is the food utilisation. The functions F(h) and f(h,p) describe prey multiplications and predations, respectively. The last term μp describes predator mortality. Holling type II functional response [7,26,43,44] is used for predator response:
f(h,p)=Ahph+b | (2.3) |
where A is intensity for predation and b is half saturation value. The prey population is assumed to be damped by the Allee effect and parameterized as follows [45]:
F(h)=(4r(k−h0)2)h(h−h0)(k−h) | (2.4) |
where k is the carrying capacity, r is the per capita growth rate for the intensity of Allee effect defined by h0. Systems' dynamical structure is studied in detail with the stress of Allee effect in [7]. For more detailed explanations on the system parameter, see [7]. The system (2.1-2.2) is in a dimensional form. The dimensionless form of parameters is as
u=hk,v=pκk,t=aT,x=X√ad1 |
Eqs. (2.5 and 2.6) includes five dimensionless parameters instead of nine with the choice
We consider one dimension predator prey model, studied earlier in [7,36,38]:
α=kb,β=h0k,γ=4rbkAκ(k−h0)2,δ=mua,ϵ=d2d1. |
For more details on system parametrizations and its dimensionless version see [7].
The nonspatial counterpart of the system becomes:
dudt=γu(u−β)(1−u)−uv1+αu, | (2.5) |
dvdt=uv1+αu−δv. | (2.6) |
The steady state of the system (2.5-2.6) is a solution of the following equations:
0=γu(u−β)(1−u)−uv1+αu, | (2.7) |
0=uv1+αu−δv. | (2.8) |
One of the states is extinction E1=(0,0) is always stable for any values of parameters. Other is semi-trivial equilibrium E2=(β,0), E3=(1,0), extinction of predator case. Stability of the predator-free steady state is saddle. For more information on the stability of these boundary states and the phase plane of the system see [7]. The other steady state of the system is the coexistence one E4=(u,v). The steady state of the coexistence system is:
v=γ(u−β)(1−u)(1+αu),u=δ1−δα. |
Jacobian matrix of the system (2.5-2.6) is as follows:
J=(2uγ−3u2γ−γβ+2uγβ−v(1+uα)2−u(1+uα)v(1+uα)2u(1+uα)−δ). | (2.9) |
The eigenvalues are the solutions for the characteristic equation for each of the steady states:
det(Ji−λI)=0, |
where I is the unit matrix and Ji is the matrix (2.9) with the components calculated at the steady state Ei, i=1,2,3,4. For further details on the stability analysis can be found in [7]. The system is addressed by numerical simulations with different parameter values and initial conditions, species extinction thresholds with and without spatial gradient is studied.
The spatially extended version of the model is as follows:
∂u∂t=d1∂2u∂x2+γu(u−β)(1−u)−uv1+αu, | (2.10) |
∂v∂t=d2∂2v∂x2+uv1+αu−δv. | (2.11) |
Here u and v are the densities of prey and predator, respectively, at time t and position x. It should be emphasized here that the diffusion coefficients are chosen equal, i.e., d1/d2=1 in order to reduce the number of system parameters. However, it should be clarified that this was not a primary restriction and the specific value d1=d2 was chosen. Preliminary numerical simulations for the range of 0.5−2 for ϵ have shown that all the obtained regimes still exist [7]. Arbitrary choice of the initial conditions is not consistent with the inherent properties of the system [46]. So the choice of the system initial conditions should, therefore, correspond to its biological equivalents. In this case, the initial distribution of the species can be represented as a finite support function based on the assumption that foreign species are attacked by an exotic species. The initial conditions, therefore, can be taken into account [7]:
u(x,0)=u0,for−Δu<x<Δu,otherwiseu(x,0)=0, | (2.12) |
v(x,0)=v0,for−Δv<x<Δv,otherwisev(x,0)=0 | (2.13) |
where u0, v0 are the initial population densities and Δu, Δv correspond the initially invaded domain radius. This type of initial conditions in Eqs. (2.12 and 2.13) also refers to the biological control issue [47,48,49]. On the other hand to see the response of spatial distribution to different initial conditions, another type of condition used in our earlier work for the oxygen-plankton system [31,33,35] is taken into account. For this case, it is assumed that at the beginning of the spread, the domain already contains both populations at the density level corresponding to stationary state (˜u,˜v) with the introduction of a linear disruption to the predator distribution. We therefore have the following constant gradient form, i.e.,
u(x,0)=˜u,v(x,0)=˜v+ϵ(x−x0), | (2.14) |
where ˜u, ˜v are the steady states of the system dynamics, given by Eq. (2.15), with ϵ=0.001 and x0=200, i.e., auxiliary parameters [50].
˜u=δ1−αδ,˜v=γ(1+α˜u)(˜u−β)(1−˜u),. | (2.15) |
To identify the behavior of the system, extensive numerical simulations performed in [7] and the systems' parametric dependencies are sketched perfectly. In the light of this work, Eqs. (2.10 and 2.11) is set numerically on the 0<x<L domain by finite difference method using zero–flux boundary conditions. The steps of the numerical mesh are chosen as Δx=0.5 and Δt=0.01. Since we are mainly interested in spatial gradient effect on system dynamics, we keep the system parameters as in [7] (see Figure 5 in [7]) to use the bifurcation diagram of the systems shown for α=0.5, β=0.28, γ=7, Δu=5, Δv=3, u0=1, v0=1 and L=400.
As it is detailed in the introduction, when the surrounding environmental conditions are not favor of species and even also these conditions bring the species to the edge of extinction, species propagate in space to sustain their existence to avoid this environment. In [7], the prey-predator system is addressed on the basis that prey development is damped by a strong Allee effect, but with no attention to species persistence/extinction scenario under the effect of a spatial gradient. The rate of zooplankton mortality, i.e., δ, is quantified differently in this study. The mortality rate is specified as a function of space to determine the species precautions for potential unsuitable environmental conditions. Suppose that the species move away to get rid of unsuitable environmental conditions. We then consider δ=δ(x), but keep the other parameters unchanged.
Animal movements are affected by some external forces, one of which is the spatial forces, such as the fragmented habitat [51], slope [52], beachhead [34], etc. One of the important stages of the invasion is the ‘beachhead’, i.e., the establishment of an invader as a self-sufficient population. At least these variables should be included in the ‘realistic’ δ(x) function, but it becomes quite complicated. Since the purpose of this study is to reveal species' precautions, in this case, it is spatial movement, we consider the mortality rate to be a linearly decreasing function, assuming that spatial movement facilitates predator abundance, i.e., mortality rate of predator decrease in space. black The reason behind choosing predator mortality rate to decrease beyond a critical point threaten the system with extinction is that this choice increases the predator size to the point that the prey is pushed below the Allee effect threshold which results in extinction for both prey and predator. Hence the main purpose here is to show the dynamics become persistence under the influence of the spatial gradient in predator mortality rate despite the Allee effect. The possible simplest choice of δ is given as follows:
δ=δLforω=0,δ=δL−ω(x−x1)forω>0. | (2.16) |
Here, x1 is beachhead, δL is predator mortality rate before the spatial movement starts, and the ω quantifies the slope of the spatial gradient. Here this study aims to understand species' persistence/extinction rather than species' migration in detail, the simplest scenario in this context is to consider spatial gradients considering a decrease in predator mortality in space.
Therefore, we compare two possible approaches. The first is without a spatial gradient scenario, i.e., (ω=0). Second is with the case of the spatial gradient, i.e., with the effect of the spatial gradient, the mortality rate decreases (ω>0), and taking into account the beachhead as well. In particular, for the initial death rate of the predator, we find ω from appropriate range corresponding to the phase plane structure in (Figure 14 in [7]) and keep the system parameters as in [7] to see the species distributions under the effect of this new approach. It should be emphasized here that, contrary to [7], δ is taken as a function of space for the system (2.10-2.11).
To analyze the spatial stability of equilibrium points based on the knowledge of temporal system linearization method is used [53].
Theorem 2.1. The trivial equilibrium E1 is stable if w(x−x1)<l2d2+δL. The semitrivial equilibrium E2 is stable if d1>~d1, where ~d1=γβ(1−β)l2 and with given condition 2.21. The boundary equilibrium E3 is stable if d1>d⋆1, where d⋆1=γ(β−1)l2 and given stability condition 2.22 is satisfied.
Proof. Let u=ˆu+M, v=ˆv+N, where M and N are perturbations around an arbitrary equilibrium point (ˆu,ˆv). The linearized form of the system (2.10-2.11) can be given as
∂M∂t=a11M+a12N+d1∂2M∂x2, | (2.17) |
∂N∂t=a21M+a22N+d2∂2N∂x2, | (2.18) |
where a11=2γˆu−3γˆu2−γβ+2γˆuβ−ˆv(1+αˆu)2, a12=−ˆu1+αˆu, a21=ˆv(1+αˆu)2, a22=ˆu1+αˆu−δ. The eigen function is ∑∞r=0(arbr)exp(λt+ilx) for the wave number l>0 and the growth rate of perturbation λ>0, x is the spatial coordinate for reals ar,br. General uniform steady state of the spatial system (2.10-2.11) by (ˆM,ˆN), we have the perturbation in the following form,
(M(t,x)N(t,x))=(ˆMˆN)+(ϵ1ϵ2)eλt+ilx
ϵi are assumed to be small enough and the corresponding amplitudes for i=1,2. Then the problem is reduced to eigenvalue problem in the parameter λ for the following matrix at (ˆM,ˆN).
J(ˆu,ˆv)=[2uγ−3u2γ−γβ+2uγβ−v(1+uα)2−l2d1−u(1+uα)v(1+uα)2u(1+uα)−δ−l2d2] | (2.19) |
The eigenvalues evaluated at E1 are −γβ−l2d1<0 and −δ−l2d2. Since the main point of this work is to reveal the species extinction and persistence issue by taking predator mortality rate as a linear dependence in space, the term −(δL−ω(x−x1))−l2d2 should be less than zero for stable E1. Hence, −(δL−ω(x−x1))−l2d2<0
ω(x−x1)<l2d2+δL | (2.20) |
It means that if the condition 2.20 is satisfied the system will then collapse and the all species extinct. The eigenvalues evaluated at E2 are γβ(1−β)−l2d1 and β1+αβ−δ−l2d2. This boundary state is spatially stable if d1>~d1, where ~d1=γβ(1−β)l2 and
ω(x−x1)<l2d2+δL−β1+αβ | (2.21) |
is satisfied. The eigenvalues evaluated at E3 are γβ−γ−l2d1 and 11+α−δ−l2d2. This boundary state is spatially stable if d1>d⋆1, where d⋆1=γ(β−1)l2 and with the condition
ω(x−x1)<l2d2+δL−11+α | (2.22) |
is satisfied. The coexistence steady state characteristic equation is of the form
λ2+Q1λ+Q2=0 | (2.23) |
where Q1=−a11+l2(d1+d2), Q2=a21+l4(d1+d2)−l2d2a11. The stability of coexistence state and the sustainability of the species is detailed in numerical simulations. Note that in numerical simulations taking d1=d2=1 is a technical consequence of being dimensionless form.
Under the light of above analytical calculations, we start with the response of the systems to different initial conditions. Figure 1a, b show the spatial distributions of predator and prey for the obtained initial conditions given by Eqs. (2.12 and 2.13) and Eq. (2.14), respectively. There is no important qualitative difference between distributions provided by Eqs. (2.12 and 2.13) (shown in Figure 1a) and given by Eq. (2.14) (shown in Figure 1b). It indicates that taking Eqs. (2.12 and 2.13) as shown in Figure 1a and using it as the initial conditions does not allow any qualitative difference to the inherent spatial structures. The following simulations will, therefore, operate with the initial condition provided by Eqs. (2.12 and 2.13). Red solid lines in all numerical simulations show δ value to see the changes of δ in space. In general, the red line is shown in the first figures for the same parameter choice.
The purpose of this study is to describe the effect of spatial gradient and beachhead on species when the species are faced with unsuitable environmental conditions and these conditions lead to extinction. It is assumed that δ can change in space with the help of the gradient. Let us begin with the case where δ is taken from a safe region (i.e., Figure 14 in [7]). Since δ is predator mortality a small quantity of it places a great deal of pressure on its prey. As a consequence, low δ leads to the extinction of species. That's the point this research has been activated. Therefore, however, the issue is whether the spatial gradient and the beachhead prevent or postpone the eventual extinction even if its population is driven by its predator which is inevitable without a gradient case. For this purpose, the numerical simulations follow the order; first, observing the spatial distribution without a gradient case with the aim of finding the critical value of extinction. Second, beginning from this extinction point and introducing gradient and beachhead to the system to detail species persistence/extinction issue despite the Allee effect. Then some numerical observations on two dimensional space are performed.
Figure 2 shows the distribution of prey and predator at different times. We note the species invades the whole domain with irregular spatiotemporal oscillations of fixed values of δ. System dynamics is sustainable and extinction is not observed even for longer periods (i.e., t=20000 equivalent simulations are not given for brevity reasons). Note that there is no smooth pattern due to lack of gradient, i.e., ω. Figures 2 and 3 are provided to display the correspondence of the results obtained by a slight decrease in δ, which is detailed and observed in [7].
Figure 3 shows snapshots of prey and predator spatial distributions obtained at different times for ω=0.0001, δL=0.49 and beachhead x1=133. For this value of δ, there may be two somewhat different scenarios with invasive species with irregular spatiotemporal oscillations and moving patches (cf. [7] how invasive species can be established and general invasive behavior for the corresponding domain). The first scenario is shown in Figure 3a. With a slight slope, i.e., ω, the distributions of the species in δ do not change at different times. Such irregular patterns are connected to the condition of coexistence (smooth pattern) behind it. The second scenario is shown in Figure 3b, with the same value as δ but for larger time limits the irregular distribution appears like a series of separate patches. At larger times (t>5000) the system tends to give way to distinctly divided patches but this scenario is not seen here for the sake of brevity.
We now proceed to simulations to obtain an extinction threshold of species without a gradient case. The extinction value is meager in the phase plane map in [7], to observe that extensive numerical simulations performed. Note that the obtained numerical simulations indicate that the pattern of distribution when the system approaches its extinction value is regular, i.e., the spatial system gives certain early warning signals for potential extinction. Such findings support the regular patterning of critical transition approaches such as in the oxygen-plankton model [31] and in other systems [28,30,54,55,56,57,58].
Figure 4a, b show the prey and predator distributions obtained for different values of δ and Figure 4c is for the power spectrum analysis for the prey distribution of Figure 4b. It is observed that species' extinction happens when the system reaches a critical value of δcr=0.44. Therefore the δ values that are used for Figure 4a, b bring the system close to its extinction value so the system is invaded by regular patterning. Contrary to [46], regular oscillations continue for longer periods of time. Figure 4c shows the analysis of the power spectrum as evidence that the system (2.10-2.11) is capable of developing regular spatial patterns. The question here is what is called ‘regular spatial distribution’ is actually regular or irregular. In this case, the analysis of the power spectrum should be taken into account, and the system is accepted as a regular if the analysis of the power spectrum has some leading frequencies. For the regular structure see Figure 4b, some of the leading frequencies are readily distinguishable, see Figure 4c. This regular structure emerges as an early warning signal and says the system is the edge of extinction, see the earlier work on the early warning signal and the regime shift [31,56,58]. Note that taking δ values from a given range of (0.44,0.45] always results in regular species spatial distributions.
Since the issue of species extinction/persistence is interested, the gradient and the beachhead have been taken into account here. This implies that δcr (where the prey and predator species is extinct) is used as a starting point to see whether the gradient and the beachhead can prevent extinction.
The above results are demonstrated to show the spatial distribution of the system with and without ω=0 and x1. The following simulations are considered when the system is affected by the gradient of space and the beachhead.
Figure 5 shows the population density snapshots. The vertical blue dashed line indicates the assumed beachhead value. It is important to note here again that, without a spatial gradient, the species is extinct, but for the given parameter ω=0.001and provided beachhead x1=100, the species propagate through traveling pulses with irregular wavemakers (see system Figure 5 in [7]). Traveling pulses die when the system is closer to the critical value of δ. When the traveling pulse reaches this critical point in space (importantly this critical value of δ is totally different from δcr, when we call a new settlement of δ as δcr⋆, it is obvious that δcr⋆<δcr, the system results in the first extinction of prey followed by the extinction of predators. Note that, δcr⋆ is to be the place along with the space-axis beyond which neither prey nor predator exists. The system (2.10-2.11) is maintained for a sufficiently long time (t=20000) and the population retains its spatial distributions and does not get extinct. With a decrease in ω (i.e., ω=0.0006), in Figure 6 the beachhead stays the same as before (Figure 5, i.e., (x1=100)), the region with a smooth spatial distribution disappears and the area is occupied by non-stationary pulses.
There seems to be a common behavior in all species' spatial distribution that traveling pulses disappear where the system reaches δcr⋆. It may be worth looking at details. It is, therefore, determined to detect the positions of the traveling pulse when it comes to near extinction, i.e., when the prey density is sufficiently small (u>ϵ where ϵ=0.02). The general tendency of predators to extinction is qualitatively close to prey, therefore, we display only distribution of prey to avoid repetition. Figure 7a, b show the extinction position of the traveling pulses in space (i.e., xbar) for Figure 5 and Figure 6 for t=6000, respectively. As can be easily seen in Figure 7a, the extinction range can be limited to xbar≈180 and xbar≈270. The restriction is clear for ω=0.001, but there are some lower peaks for ω=0.0006. It implies that, in some situations, the traveling pulses extinct earlier in space than the general tendency of the system. Interestingly extinction exceeds the value of the beachhead for both slope values so it happens beyond the vertical blue line, i.e., the beachhead.
Figure 8 shows the population density snapshots for the given parameter ω=0.001 and different beachhead from above, i.e., x1=150. In Figure 8, traveling pulses are extinct beyond the beachhead and also the scenario is the same more or less. Figure 9 shows snapshots for ω=0.0006 in the same manner. Figure 10a, b show the extinction position of traveling pulses in space of different values of ω. Once, extinction exceeds the value of the beachhead for both slope values so extinction happens beyond the vertical blue line, i.e., the beachhead. This issue will be discussed in the last section.
Applying a different spatial gradient with the aid of the slope, i.e., ω, and the beachhead, i.e., x1, an important question emerges as to whether there is any relation between the beachhead and the slope. Figure 11 shows the data point set, where the diamonds represent the extinction position for the corresponding ω collected as a result of numerous numerical simulations. There is therefore a relationship between the critical point (beachhead) and the spatial gradient defined by simple rational function, i.e., 1ωα where α=1. The dashed line in Figure 11 plays a role as a separator of extinction and persistence regions. Thus, the area below the dashed line refers to the species extinction, while above it is the case where the population survives. Therefore, the persistence of the species can be determined by the simultaneous effect of the gradient slope and beachhead. If the slope rises, the extinction region of the species will decline. In this case, even if the beachhead is taken big enough, the species is not extinct.
Figure 12 is given to provide more details on temporal variations of population. Figure 12a–c indicate the local dynamics at fixed spatial positions typical for different parameter domains, i.e., extinction and persistence. Figure 12d–f show the phase plane of the local population densities obtained at a fixed position at space x1=50, i.e., the system below the dashed line in the extinction region, x1=94, i.e., the system on the slightly above extinction threshold. As it is observed in spatial distribution, the temporal distribution at closer to critical point results in regular distribution. Note that for ω=6×10−4, the threshold value is x1=93, see Figure 11 for more detail. x1=150, i.e., the system in the persistence region above the dashed line.
Now attempts are being made to create the predator-prey system's spatial dynamics, applied to two spatial dimensions, which are more realistic than its one spatial dimension counterpart. The corresponding two-dimensional simulations coincide with its one-dimensional counterparts, the boundary and the initial conditions are selected accordingly. Therefore, all two dimensional simulations are performed in a rectangular domain: (x,y)∈[0,Lx]×[0,Ly] at the domain boundaries Neumann boundary conditions are used,
∂(u,v)∂x|x=0,Lx=∂(u,v)∂y|x=0,Ly=0. | (3.1) |
The spatial distribution of prey density is given in a two-dimensional system. Predator density demonstrates qualitatively similar behavior to prey, which is why the distribution of prey is provided only for the sake of brevity. To observe the distribution of prey, different types of initial conditions are used. The first one (see Figure 13) given by Eqs. (3.2 and 3.3) defines spatially homogeneous prey distribution at steady state with constant gradient predator distribution for ω=0.0006, δL=0.44 and x1=100 where u0=v0=1, ϵ2=ϵ3=3.10−5.
u(x,y,0)=u0, | (3.2) |
v(x,y,0)=v0−ϵ2(x−200)−ϵ3(y−200), | (3.3) |
In the latter case, the system's initial distribution is as follows:
u(x,y,0)=u0−ϵ1(x−240)(x−294)−ϵ2(y−120)(y−280), | (3.4) |
v(x,y,0)=v0−ϵ2(x−200)−ϵ3(y−200). | (3.5) |
The two-dimensional system dynamics follow more or less the same scenario as the corresponding one-dimensional system. At an early stage (see Figure 13a), a smooth distribution of prey appears. This smooth distribution is followed by stripes corresponding to the traveling pulses in one dimension case, see Figure 13b. Eventually, the stripes are broken and give way to irregular structure in Figure 13c. The irregular structures swell and produce some circular shape and some broken irregular stripes, see Figure 13d. Observed that even for a larger time limit the upper part of the domain is not invaded regardless of the different initial condition, which ensures that the prey distribution stays close to zero.
The second one (see Figure 14) given by Eqs. (3.4 and 3.5) describes the spatially homogeneous distribution of prey at its steady state with a constant gradient of predators. ω=0.0006, δL=0.44 and x1=100 with u0=v0=1, ϵ1=2.10−7, ϵ2=3.10−5, ϵ3=ϵ4=6.10−5. Note that the initial predator mortality rate, i.e., δL=0.44, is based on its one–dimensional case. The only difference between Figure 13 and Figure 14 is the choice of the initial distributions. It can be inferred that there is no qualitative difference between the different initial conditions and the general the tendency of the spread of species. In both cases, the typical sequence of events is the smooth distribution that produces stripes, the stripes are followed by their broken versions and then the domain is invaded by irregular structures, like circular shapes, parts of broken stripes, etc.
It is noticed that the spatial patterns shown in Figures 13 and 14 exist in the parameter range where the species persistence would not be possible in the corresponding one dimensional case (i.e., for δL=0.44). It is, therefore, consistently observed that the dynamics of the two-dimensional system is persistent/sustainable in a broader range [35].
In this work a predator-prey reaction diffusion model is studied with an Allee effect on the prey. The model is based on earlier work Petrovskii et al. [7], in which the conditions under which the predator-prey system could invade an environment were studied. The present manuscript extends that study to the case in which the mortality rate of the predator is spatially varying; that is, the mortality rate is given in Eq. (2.16) by the function δ=δL−ω(x–x1) where δL is a constant, ω is the slope, x is distance along the x-axis, and x1 is what is termed the ‘beachhead' that the invading predator-prey system can establish. Therefore, in relation to changing climate, the dynamics of the prey-predator system is assumed to reflect the spatial gradient of predator mortality rate. The objective of this study is to determine the joint effect of the slope of the gradient and the ‘beachhead’, that is, the points along the spatial axis at which the system, despite the Allee effect, can establish a sustainable population as an adaptation mechanism to new environment.
Some interesting results are observed. One is that there is a simple inverse algebraic relationship between the slope of the gradient and critical point. This is found by allowing δL to take the value 0.44, which is the threshold point for extinction if there is no slope, i.e., ω=0. Another is that, under circumstances when the system is on the verge of extinction, for example, δL slightly above 0.44 and ω=0, the pattern of spatiotemporal oscillations becomes completely regular with a fixed period and amplitude. This can serve as an ‘early warning signals' of changing climate [31] that the system is on the verge of extinction if δL decreases further.
Gradient and beachhead are used as a precaution as species face extinction by applying a spatial gradient to the predator mortality rate. This predator prey system has a well-structured parameter plane is labeled in [7]. The structure of the system's spatial dynamics with linearly decreasing values of δ in space is described by comprehensive numerical simulations. The difference of this study from the literature is that the rate of predator mortality changes with space to account for the response to changing climate. It is observed that despite the Allee effect, the species can establish a sustainable population. Such spatial structure provides some early warning signals (i.e., regular distribution) when the rate of predator mortality is closer to the extinction point. As observed in the oxygen-plankton system [31,33] when the system is close to extinction, transitions from irregular distribution to regular one which is used for the early warning signal of climate change are also observed for this system.
Note that, the emergence of pulses due to the spatial decrease in δ corresponds to its non-gradient case phase plane map; see Domain 2 and 2∗ in Figures 13 and 14 [7]. Here, it can be seen that the traveling pulses die when the system becomes closer to the critical value of predator mortality rate, i.e., δ. When the traveling pulse reaches the critical point in time, the species are extinct.
It is observed that the persistence of the species can be interpreted by the slope (steepness) of the spatial gradient and the length of the beachhead. If the slope is raised, the persistence of the system is expanded. In this case, if the beachhead becomes sufficiently large, the species becomes extinct (see the systems' extinction tendency in Figure 11, i.e., a further increase in ω results in extinction). On the other hand, the gradient addition does not prevent extinction but at least postpones it and allows the species a chance to exist for a while. One of the most interesting results of this work is that the spatial gradient can not prevent the extinction of species, at least for the spatial gradient on prey mortality but it has the capacity to slow down the approaching extinction. Interestingly, when the system reaches its critical value resulting in extinction beyond this value, its spatial distribution becomes regular and compatible with the structure which was obtained in the oxygen-plankton system [31].
The one dimension to two dimensions spatially explicit extension is considered. Here, it is shown that the dying traveling pulse in 1-D is represented as broken propagating stripes in 2-D. It should also be emphasized that different choice of initial distributions for both 1-D and 2-D case, spatial system patterns remains qualitatively similar.
In conclusion, the most interesting question here is whether the spatial gradient of different slopes can protect the species from extinction as climate change adversely affects species. This question is verified in this work with extensive numerical simulations, with the choice of large enough ω, and for specific beachhead values, while leaving certain parameters set for convenience. It is observed that, at least for this system (2.10-2.11), a gradient can not stop an inevitable ecological disaster, but has capability to slow it down, which is why the species has a chance to survive for a while.
The author is thankful to Sergei Petrovskii (University of Leicester) for his stimulating discussions at the all stages of this study. Useful comments from editor and reviewers are appreciated.
The author has no conflict of interest.
[1] | Déragon F (2016) La ferme biologique, un espace d'éducation relative à l'éco-alimentation et de construction du rapport à la terre[Thèse de Maitrise]: Université du Québec à Montréal. |
[2] |
Hughner RS, McDonagh P, Prothero A, et al. (2007) Who are organic food consumers? A compilation and review of why people purchase organic food. J Consum Behav 6: 94–110. https://doi.org/10.1002/cb.210 doi: 10.1002/cb.210
![]() |
[3] | Willer H, Trávníček J, Meier C, et al. (2021) The World of Organic Agriculture Statistics and Emerging Trends 2021 Research Institute of Organic Agriculture FiBL, Frick, and IFOAM—Organics International, Bonn. |
[4] |
Cappelli A, Oliva N, Cini E (2020) Stone milling versus roller milling: A systematic review of the effects on wheat flour quality, dough rheology, and bread characteristics. Trends Food Sci Technol 97: 147–155. https://doi.org/10.1016/j.tifs.2020.01.008 doi: 10.1016/j.tifs.2020.01.008
![]() |
[5] |
Doblado-Maldonado AF, Pike OA, Sweley JC, et al. (2012) Key issues and challenges in whole wheat flour milling and storage. J Cereal Sci 56: 119–126. https://doi.org/10.1016/j.jcs.2012.02.015 doi: 10.1016/j.jcs.2012.02.015
![]() |
[6] | Campbell GM (2007) Chapter 7—Roller milling of wheat. In: Salman AD, Ghadiri M, Hounslow MJ, Handbook of Powder Technology, Elsevier Science B.V., 383–419. https://doi.org/10.1016/S0167-3785(07)12010-8 |
[7] | Boudreau A, Ménard G (1992) Le Blé: éléments fondamentaux et transformation. |
[8] |
Mateos-Salvador F, Sadhukhan J, Campbell GM (2011) The normalised Kumaraswamy breakage function: A simple model for wheat roller milling. Powder Technol 208: 144–157. https://doi.org/10.1016/j.powtec.2010.12.013 doi: 10.1016/j.powtec.2010.12.013
![]() |
[9] | FAO (1985) Standard for wheat flour. |
[10] | IAOM (2018) Fundamentals of Flour Milling. Lenexa: International Association of Operative Millers, 422. |
[11] | Steffan P (2012) An optimization model: Minimizing flour millers' costs of production by blending wheat and additives[Thèse de Maitrise]: Kansas State University. |
[12] |
Vrček Ⅳ, Čepo DV, Rašić D, et al. (2014) A comparison of the nutritional value and food safety of organically and conventionally produced wheat flours. Food Chem 143: 522–529. https://doi.org/10.1016/j.foodchem.2013.08.022 doi: 10.1016/j.foodchem.2013.08.022
![]() |
[13] |
Borghi B, Giordani G, Corbellini M, et al. (1995) Influence of crop rotation, manure and fertilizers on bread making quality of wheat (Triticum aestivum L.). Eur J Agron 4: 37–45. https://doi.org/10.1016/S1161-0301(14)80015-4 doi: 10.1016/S1161-0301(14)80015-4
![]() |
[14] |
Johansson E, Svensson G (1998) Variation in bread-making quality: Effects of weather parameters on protein concentration and quality in some Swedish wheat cultivars grown during the period 1975–1996. J Sci Food Agric 78: 109–118. https://doi.org/10.1002/(SICI)1097-0010(199809)78:1%3C109::AID-JSFA92%3E3.0.CO; 2-0 doi: 10.1002/(SICI)1097-0010(199809)78:1<109::AID-JSFA92>3.0.CO;2-0
![]() |
[15] |
Triboi E, Abad A, Michelena A, et al. (2000) Environmental effects on the quality of two wheat genotypes: 1. quantitative and qualitative variation of storage proteins. Eur J Agron 13: 47–64. https://doi.org/10.1016/S1161-0301(00)00059-9 doi: 10.1016/S1161-0301(00)00059-9
![]() |
[16] |
Köksal G, Batmaz İ, Testik MC (2011) A review of data mining applications for quality improvement in manufacturing industry. Expert Syst Appl 38: 13448–13467. https://doi.org/10.1016/j.eswa.2011.04.063 doi: 10.1016/j.eswa.2011.04.063
![]() |
[17] |
Misra NN, Dixit Y, Al-Mallahi A, et al. (2020) IoT, big data and artificial intelligence in agriculture and food industry. IEEE Int Things J 9: 6305–6324. https://doi.org/10.1109/JIOT.2020.2998584 doi: 10.1109/JIOT.2020.2998584
![]() |
[18] | Kelleher JD (2019) Deep Learning (MIT Press Essential Knowledge series). The MIT Press, Illustrated edition, 296. |
[19] |
Fisher OJ, Watson NJ, Escrig JE, et al. (2020) Considerations, challenges and opportunities when developing data-driven models for process manufacturing systems. Comput Chem Eng 140: 106881. https://doi.org/10.1016/j.compchemeng.2020.106881 doi: 10.1016/j.compchemeng.2020.106881
![]() |
[20] |
Branlard G, Dardevet M, Saccomano R, et al. (2001) Genetic diversity of wheat storage proteins and bread wheat quality. Euphytica 119: 59–67. https://doi.org/10.1023/A:1017586220359 doi: 10.1023/A:1017586220359
![]() |
[21] |
Maphosa L, Langridge P, Taylor H, et al. (2014) Genetic control of grain yield and grain physical characteristics in a bread wheat population grown under a range of environmental conditions. Theor Appl Genet 127: 1607–1624. https://doi.org/10.1007/s00122-014-2322-y doi: 10.1007/s00122-014-2322-y
![]() |
[22] |
Wrigley CW, Blumenthal C, Gras PW, et al. (1994) Temperature variation during grain filling and changes in wheat-grain quality. Funct Plant Biol 21: 875–885. https://doi.org/10.1071/PP9940875 doi: 10.1071/PP9940875
![]() |
[23] |
Mefleh M, Conte P, Fadda C, et al. (2020) From seed to bread: Variation in quality in a set of old durum wheat cultivars. J Sci Food Agric 100: 4066–4074. https://doi.org/10.1002/jsfa.9745 doi: 10.1002/jsfa.9745
![]() |
[24] |
Aydin N, Sayaslan A, Sönmez M, et al. (2020) Wheat flour milling yield estimation based on wheat kernel physical properties using artificial neural networks. Int J Intell Syst Appl Eng 8: 78–83. https://doi.org/10.18201/ijisae.2020261588 doi: 10.18201/ijisae.2020261588
![]() |
[25] |
Dowell FE, Maghirang EB, Xie F, et al. (2006) Predicting wheat quality characteristics and functionality using near-infrared spectroscopy. Cereal Chem 83: 529–536. https://doi.org/10.1094/CC-83-0529 doi: 10.1094/CC-83-0529
![]() |
[26] |
Caporaso N, Whitworth MB, Fisk ID (2018) Near-Infrared spectroscopy and hyperspectral imaging for non-destructive quality assessment of cereal grains. Appl Spectrosc Rev 53: 667–687. https://doi.org/10.1080/05704928.2018.1425214 doi: 10.1080/05704928.2018.1425214
![]() |
[27] |
Assadzadeh S, Walker CK, McDonald LS, et al. (2022) Prediction of milling yield in wheat with the use of spectral, colour, shape, and morphological features. Biosyst Eng 214: 28–41. https://doi.org/10.1016/j.biosystemseng.2021.12.005 doi: 10.1016/j.biosystemseng.2021.12.005
![]() |
[28] |
Delwiche SR, Souza EJ, Kim MS (2013) Limitations of single kernel near-infrared hyperspectral imaging of soft wheat for milling quality. Biosyst Eng 115: 260–273. https://doi.org/10.1016/j.biosystemseng.2013.03.015 doi: 10.1016/j.biosystemseng.2013.03.015
![]() |
[29] |
Unlersen MF, Sonmez ME, Aslan MF, et al. (2022) CNN–SVM hybrid model for varietal classification of wheat based on bulk samples. Eur Food Res Technol 248: 2043–2052. https://doi.org/10.1007/s00217-022-04029-4 doi: 10.1007/s00217-022-04029-4
![]() |
[30] |
Zhu J, Li H, Rao Z, et al. (2023) Identification of slightly sprouted wheat kernels using hyperspectral imaging technology and different deep convolutional neural networks. Food Control 143: 109291. https://doi.org/10.1016/j.foodcont.2022.109291 doi: 10.1016/j.foodcont.2022.109291
![]() |
[31] |
Gucbilmez CM, Sahin M, Akcacik AG, et al. (2019) Evaluation of GlutoPeak test for prediction of bread wheat flour quality, rheological properties and baking performance. J Cereal Sci 90: 102827. https://doi.org/10.1016/j.jcs.2019.102827 doi: 10.1016/j.jcs.2019.102827
![]() |
[32] |
Sisman CB, Ergin AS (2011) The effects of different storage buildings on wheat quality. J Appl Sci 11: 2613–2619. https://doi.org/10.3923/jas.2011.2613.2619 doi: 10.3923/jas.2011.2613.2619
![]() |
[33] |
Jia C, Sun D-W, Cao C (2001) Computer simulation of temperature changes in a wheat storage bin. J Stored Prod Res 37: 165–177. https://doi.org/10.1016/S0022-474X(00)00017-5 doi: 10.1016/S0022-474X(00)00017-5
![]() |
[34] |
Kibar H (2015) Influence of storage conditions on the quality properties of wheat varieties. J Stored Prod Res 62: 8–15. https://doi.org/10.1016/j.jspr.2015.03.001 doi: 10.1016/j.jspr.2015.03.001
![]() |
[35] | Posner E, Deyoe C (1986) Changes in milling properties of newly harvested hard wheat during storage. Cereal Chem 63: 451–456. |
[36] |
Campbell GM, Sharp C, Wall K, et al. (2012) Modelling wheat breakage during roller milling using the Double Normalised Kumaraswamy Breakage Function: Effects of kernel shape and hardness. J Cereal Sci 55: 415–425. https://doi.org/10.1016/j.jcs.2012.02.002 doi: 10.1016/j.jcs.2012.02.002
![]() |
[37] | Posner ES, Hibbs AN (2005) Wheat flour milling: American Association of Cereal Chemists, Inc. https://doi.org/10.1094/1891127403 |
[38] | Yoon BS, Brorsen BW, Lyford CP (2002) Value of increasing kernel uniformity. J Agric Resour Econ 27: 481–494. |
[39] |
González-Torralba J, Arazuri S, Jarén C, et al. (2013) Influence of temperature and r.h. during storage on wheat bread making quality. J Stored Prod Res 55: 134–144. https://doi.org/10.1016/j.jspr.2013.10.002 doi: 10.1016/j.jspr.2013.10.002
![]() |
[40] | Catterall P (1998) Flour milling. In: Cauvain SP, Young LS, Technology of Breadmaking, Boston, MA: Springer US., 296–329. https://doi.org/10.1007/978-1-4615-2199-0_12 |
[41] |
Tibola CS, Fernandes JMC, Guarienti EM (2016) Effect of cleaning, sorting and milling processes in wheat mycotoxin content. Food Control 60: 174–179. https://doi.org/10.1016/j.foodcont.2015.07.031 doi: 10.1016/j.foodcont.2015.07.031
![]() |
[42] |
Magyar Z, Véha A, Pepó P, et al. (2019) Wheat cleaning and milling technologies to reduce DON toxin contamination. Acta Agraria Debreceniensis 0: 89–95. https://doi.org/10.34101/actaagrar/2/3684 doi: 10.34101/actaagrar/2/3684
![]() |
[43] | Bettge A, Rubenthaler G, Pomeranz Y (1989) Air‐aspirated wheat cleaning in grading and in separation by functional properties. Cereal Chem 66: 15–18. |
[44] |
Hayta M, Çakmakli Ü (2001) Optimization of wheat blending to produce breadmaking flour. J Food Process Eng 24: 179–192. https://doi.org/10.1111/j.1745-4530.2001.tb00539.x doi: 10.1111/j.1745-4530.2001.tb00539.x
![]() |
[45] | Elevi B, Öztürk H, Kacir Z (2017) Optimization of wheat and flour blending for cost minimization by using mathematical modelling. Kastamonu Üniversitesi İktisadi ve İdari Bilimler Fakültesi Dergisi, 264–272. |
[46] |
Kweon M, Martin R, Souza E (2009) Effect of tempering conditions on milling performance and flour functionality. Cereal Chem 86: 12–17. https://doi.org/10.1094/CCHEM-86-1-0012 doi: 10.1094/CCHEM-86-1-0012
![]() |
[47] |
Warechowska M, Markowska A, Warechowski J, et al. (2016) Effect of tempering moisture of wheat on grinding energy, middlings and flour size distribution, and gluten and dough mixing properties. J Cereal Sci 69: 306–312. https://doi.org/10.1016/j.jcs.2016.04.007 doi: 10.1016/j.jcs.2016.04.007
![]() |
[48] |
Hook SCW, Bone GT, Fearn T (1982) The conditioning of wheat. The effect of increasing wheat moisture content on the milling performance of uk wheats with reference to wheat texture. J the Sci Food Agric 33: 655–662. https://doi.org/10.1002/jsfa.2740330711 doi: 10.1002/jsfa.2740330711
![]() |
[49] |
Doblado-Maldonado AF, Flores RA, Rose DJ (2013) Low moisture milling of wheat for quality testing of wholegrain flour. J Cereal Sci 58: 420–423. https://doi.org/10.1016/j.jcs.2013.08.006 doi: 10.1016/j.jcs.2013.08.006
![]() |
[50] |
Noort MWJ, van Haaster D, Hemery Y, et al. (2010) The effect of particle size of wheat bran fractions on bread quality—Evidence for fibre-protein interactions. J Cereal Sci 52: 59–64. https://doi.org/10.1016/j.jcs.2010.03.003 doi: 10.1016/j.jcs.2010.03.003
![]() |
[51] |
Cappelli A, Guerrini L, Parenti A, et al. (2020) Effects of wheat tempering and stone rotational speed on particle size, dough rheology and bread characteristics for a stone-milled weak flour. J Cereal Sci 91: 102879. https://doi.org/10.1016/j.jcs.2019.102879 doi: 10.1016/j.jcs.2019.102879
![]() |
[52] |
Parrenin L, Danjou C, Agard B, et al. (2022) Predicting the moisture content of organic wheat in the first stage of tempering. IFAC-PapersOnLine 55: 678–683. https://doi.org/10.1016/j.ifacol.2022.09.484 doi: 10.1016/j.ifacol.2022.09.484
![]() |
[53] |
Fang C, Campbell GM (2003) On predicting roller milling performance Ⅴ: Effect of moisture content on the particle size distribution from first break milling of wheat. J Cereal Sci 37: 31–41. https://doi.org/10.1006/jcrs.2002.0476 doi: 10.1006/jcrs.2002.0476
![]() |
[54] |
Lin S, Gao J, Jin X, et al. (2020) Whole-wheat flour particle size influences dough properties, bread structure and in vitro starch digestibility. Food Funct 11: 3610–3620. https://doi.org/10.1039/C9FO02587A doi: 10.1039/C9FO02587A
![]() |
[55] |
Tóth Á, Prokisch J, Sipos P, et al. (2006) Effects of particle size on the quality of winter wheat flour, with a special focus on macro‐ and microelement concentration. Commun Soil Sci Plant Anal 37: 2659–2672. https://doi.org/10.1080/00103620600823117 doi: 10.1080/00103620600823117
![]() |
[56] | Pagani MA, Marti A, Bottega G (2014) Chapter 2 : Wheat milling and flour quality evaluation. In: Zhou W, Hui YH, De Leyn I, et al., Bakery Products Science and Technology, John Wiley & Sons, Ltd. 20–49. https://doi.org/10.1002/9781118792001.ch2 |
[57] |
Campbell GM, Webb C (2001) On predicting roller milling performance: Part Ⅰ. The breakage equation. Powder Technol 115: 234–242. https://doi.org/10.1016/S0032-5910(00)00348-X doi: 10.1016/S0032-5910(00)00348-X
![]() |
[58] |
Campbell GM, Bunn PJ, Webb C, et al. (2001) On predicting roller milling performance: Part Ⅱ. The breakage function. Powder Technol 115: 243–255. https://doi.org/10.1016/S0032-5910(00)00349-1 doi: 10.1016/S0032-5910(00)00349-1
![]() |
[59] |
Cappelli A, Mugnaini M, Cini E (2020) Improving roller milling technology using the break, sizing, and reduction systems for flour differentiation. Lwt-Food Sci Technol 133: 110067. https://doi.org/10.1016/j.lwt.2020.110067 doi: 10.1016/j.lwt.2020.110067
![]() |
[60] |
Dal-Pastro F, Facco P, Bezzo F, et al. (2015) Data-based multivariate modeling of a grain comminution process. Comput Aided Chem Eng 37: 2219–2224. https://doi.org/10.1016/B978-0-444-63576-1.50064-9 doi: 10.1016/B978-0-444-63576-1.50064-9
![]() |
[61] |
Dal-Pastro F, Facco P, Bezzo F, et al. (2016) Data-driven modeling of milling and sieving operations in a wheat milling process. Food Bioprod Process 99: 99–108. https://doi.org/10.1016/j.fbp.2016.04.007 doi: 10.1016/j.fbp.2016.04.007
![]() |
[62] |
Fang Q, Biby G, Haque E, et al. (1998) Neural network modeling of physical properties of ground wheat. Cereal Chem 75: 251–253. https://doi.org/10.1094/CCHEM.1998.75.2.251 doi: 10.1094/CCHEM.1998.75.2.251
![]() |
[63] | Stefan E-M, Voicu G, Constantin G, et al. (2018) Effects of wheat seeds characteristics on roller milling process—a review. |
[64] |
Fang C, Campbell GM (2003) On predicting roller milling performance Ⅳ: Effect of roll disposition on the particle size distribution from first break milling of wheat. J Cereal Sci 37: 21–29. https://doi.org/10.1006/jcrs.2002.0475 doi: 10.1006/jcrs.2002.0475
![]() |
[65] |
Campbell GM, Fang C, Muhamad Ⅱ (2007) On predicting roller milling performance Ⅵ: Effect of kernel hardness and shape on the particle size distribution from first break milling of wheat. Food Bioprod Process 85: 7–23. https://doi.org/10.1205/fbp06005 doi: 10.1205/fbp06005
![]() |
[66] |
Hook SCW, Bone GT, Fearn T (1982) The conditioning of wheat. The influence of roll temperature in the bühler laboratory mill on milling parameters. J Sci Food Agric 33: 639–644. https://doi.org/10.1002/jsfa.2740330709 doi: 10.1002/jsfa.2740330709
![]() |
[67] |
Fang Q, Haque E, Spillman CK, et al. (1998) Energy requirements for size reduction of wheat using a roller mill. Trans ASAE 41: 1713–1720. https://doi.org/10.13031/2013.17314 doi: 10.13031/2013.17314
![]() |
[68] |
Kalitsis J, Minasny B, Quail K, et al. (2021) Application of response surface methodology for optimization of wheat flour milling process. Cereal Chem 98: 1215–1226. https://doi.org/10.1002/cche.10474 doi: 10.1002/cche.10474
![]() |
[69] |
Oliver J, Blakeney A, Allen H (1993) The colour of flour streams as related to ash and pigment contents. J Cereal Sci 17: 169–182. https://doi.org/10.1006/jcrs.1993.1017 doi: 10.1006/jcrs.1993.1017
![]() |
[70] |
Banu I, Stoenescu G, Ionescu V, et al. (2011) Estimation of the baking quality of wheat flours based on rheological parameters of the Mixolab Curve. Czech J Food Sci 29: 35–44. https://doi.org/10.17221/40/2009-CJFS doi: 10.17221/40/2009-CJFS
![]() |
[71] |
Dowell FE, Maghirang EB, Pierce RO, et al. (2008) Relationship of bread quality to kernel, flour, and dough properties. Cereal Chem 85: 82–91. https://doi.org/10.1094/CCHEM-85-1-0082 doi: 10.1094/CCHEM-85-1-0082
![]() |
[72] |
Parenti O, Guerrini L, Zanoni B (2020) Techniques and technologies for the breadmaking process with unrefined wheat flours. Trends Food Sci Technol 99: 152–166. https://doi.org/10.1016/j.tifs.2020.02.034 doi: 10.1016/j.tifs.2020.02.034
![]() |
[73] |
Parenti O, Guerrini L, Cavallini B, et al. (2020) Breadmaking with an old wholewheat flour: Optimization of ingredients to improve bread quality. Lwt-Food Sci Technol 121: 108980. https://doi.org/10.1016/j.lwt.2019.108980 doi: 10.1016/j.lwt.2019.108980
![]() |
[74] |
Torbica A, Blazek KM, Belovic M, et al. (2019) Quality prediction of bread made from composite flours using different parameters of empirical rheology. J Cereal Sci 89: 102812. https://doi.org/10.1016/j.jcs.2019.102812 doi: 10.1016/j.jcs.2019.102812
![]() |
[75] |
Cappelli A, Guerrini L, Cini E, et al. (2019) Improving whole wheat dough tenacity and extensibility: A new kneading process. J Cereal Sci 90: 102852. https://doi.org/10.1016/j.jcs.2019.102852 doi: 10.1016/j.jcs.2019.102852
![]() |
[76] |
Ktenioudaki A, Butler F, Gallagher E (2010) Rheological properties and baking quality of wheat varieties from various geographical regions. J Cereal Sci 51: 402–408. https://doi.org/10.1016/j.jcs.2010.02.009 doi: 10.1016/j.jcs.2010.02.009
![]() |
[77] |
Oliver JR, Allen HM (1992) The prediction of bread baking performance using the farinograph and extensograph. J Cereal Sci 15: 79–89. https://doi.org/10.1016/S0733-5210(09)80058-1 doi: 10.1016/S0733-5210(09)80058-1
![]() |
[78] | Abbasi H, Emam-Djomeh Z, Seyedin S (2011) Application of artificial neural network and genetic algorithm for predicting three important parameters in bakery industries. Int J Agric Sci Res 2: 51–64. |
[79] |
Różyło R, Laskowski J (2011) Predicting bread quality (bread loaf volume and crumb texture). Pol J Food Nutr Sci 61: 61–67. https://doi.org/10.2478/v10222-011-0006-8 doi: 10.2478/v10222-011-0006-8
![]() |
[80] | Gómez Sarduy J, Viego P, Diaz Torres Y, et al. (2018) A new energy performance indicator for energy management system of a wheat mill plant. Int J Energy Econ Policy 8: 324–330. |
1. | Wei Li, Wenyin Gong, Differential evolution with quasi-reflection-based mutation, 2021, 18, 1551-0018, 2425, 10.3934/mbe.2021123 | |
2. | Luka Grbčić, Lado Kranjčević, Siniša Družeta, Ivana Lučin, Efficient Double-Tee Junction Mixing Assessment by Machine Learning, 2020, 12, 2073-4441, 238, 10.3390/w12010238 | |
3. | Chengyu Hu, Rui Qiao, Zhe Zhang, Xuesong Yan, Ming Li, Dynamic Scheduling Algorithm Based on Evolutionary Reinforcement Learning for Sudden Contaminant Events Under Uncertain Environment, 2022, 2, 2096-9929, 213, 10.23919/CSMS.2022.0014 | |
4. | Chengyu Hu, Qiuming Wang, Wenyin Gong, Xuesong Yan, Multi-objective deep reinforcement learning for emergency scheduling in a water distribution network, 2022, 14, 1865-9284, 211, 10.1007/s12293-022-00366-9 | |
5. | Li Liang, Ahmed Farouk, Water Pollution Prediction Based on Deep Belief Network in Big Data of Water Environment Monitoring, 2021, 2021, 1875-919X, 1, 10.1155/2021/8271950 | |
6. | Zhen Yan, Shuijia Li, Wenyin Gong, An adaptive differential evolution with decomposition for photovoltaic parameter extraction, 2021, 18, 1551-0018, 7363, 10.3934/mbe.2021364 | |
7. | Demetrios G. Eliades, Stelios G. Vrachimis, Alireza Moghaddam, Ioannis Tzortzis, Marios M. Polycarpou, Contamination event diagnosis in drinking water networks: A review, 2023, 13675788, 10.1016/j.arcontrol.2023.03.011 | |
8. | Albira Acharya, Jia Liu, Sangmin Shin, Evaluating the multi-dimensional resilience of water distribution networks to contamination events, 2023, 23, 1606-9749, 1416, 10.2166/ws.2023.058 | |
9. | Xingwang Shen, Shimin Liu, Bin Zhou, Tao Wu, Qi Zhang, Jinsong Bao, Digital Twin-Driven Reinforcement Learning Method for Marine Equipment Vehicles Scheduling Problem, 2024, 21, 1545-5955, 2173, 10.1109/TASE.2023.3289915 | |
10. | Xingwang Shen, Shimin Liu, Bin Zhou, Yu Zheng, Jinsong Bao, 2022, Digital twin-based scheduling method for marine equipment material transportation vehicles, 978-1-6654-9042-9, 100, 10.1109/CASE49997.2022.9926428 | |
11. | Roberto Ortega, Dana Carciumaru, Alexandra D. Cazares-Moreno, Reinforcement learning for watershed and aquifer management: a nationwide view in the country of Mexico with emphasis in Baja California Sur, 2024, 6, 2624-9375, 10.3389/frwa.2024.1384595 | |
12. | Tymoteusz Miller, Grzegorz Michoński, Irmina Durlik, Polina Kozlovska, Paweł Biczak, Artificial Intelligence in Aquatic Biodiversity Research: A PRISMA-Based Systematic Review, 2025, 14, 2079-7737, 520, 10.3390/biology14050520 | |
13. | Yu Zheng, Qianyue Hao, Jingwei Wang, Changzheng Gao, Jinwei Chen, Depeng Jin, Yong Li, A Survey of Machine Learning for Urban Decision Making: Applications in Planning, Transportation, and Healthcare, 2025, 57, 0360-0300, 1, 10.1145/3695986 |