
This research explored the life cycle analysis and environmental cost-benefit assessment of converting ash waste from hospital medical waste incineration into environmentally safe paving block raw materials. The growing concerns about medical waste disposal and its environmental impact necessitate innovative solutions for sustainable waste management. This research aimed to evaluate the feasibility and environmental implications of reusing hospital waste into raw materials for paving block mixtures. This research, a comprehensive life cycle analysis, examined the environmental impacts of medical waste collection for the production and use of paving blocks. Additionally, we conducted an environmental cost-benefit assessment to ascertain the economic feasibility and potential environmental impact forecasts of this recycling approach. The research results show that converting hospital medical waste ash into mixed raw materials for paving blocks not only immobilizes heavy metals but also provides a sustainable alternative for non-building materials. These findings highlight the potential for significant environmental and economic benefits, making this approach a promising strategy for waste management and sustainable construction practices. The cost of preventing environmental damage (eco-cost) in the process of converting ash from the incineration of medical waste into a mixture of raw materials for paving blocks is IDR 600,180.9 per cycle.
Citation: Siti Rachmawati, Syafrudin, Budiyono, Ellyna Chairani, Iwan Suryadi. Life cycle analysis and environmental cost-benefit assessment of utilizing hospital medical waste into heavy metal safe paving blocks[J]. AIMS Environmental Science, 2024, 11(5): 665-681. doi: 10.3934/environsci.2024033
[1] | Jing Liao, Jinling Wang, Jiarong Li, Xin Jiang . The dynamics and control of a multi-lingual rumor propagation model with non-smooth inhibition mechanism. Mathematical Biosciences and Engineering, 2024, 21(4): 5068-5091. doi: 10.3934/mbe.2024224 |
[2] | Yi Wang, Shicheng Zhong, Guo Wang . Dynamic selection of clarification channels in rumor propagation containment. Mathematical Biosciences and Engineering, 2023, 20(8): 14995-15017. doi: 10.3934/mbe.2023672 |
[3] | Zhenzhen Liao, Shujing Gao, Shuixian Yan, Genjiao Zhou . Transmission dynamics and optimal control of a Huanglongbing model with time delay. Mathematical Biosciences and Engineering, 2021, 18(4): 4162-4192. doi: 10.3934/mbe.2021209 |
[4] | Bevina D. Handari, Dipo Aldila, Bunga O. Dewi, Hanna Rosuliyana, Sarbaz H. A. Khosnaw . Analysis of yellow fever prevention strategy from the perspective of mathematical model and cost-effectiveness analysis. Mathematical Biosciences and Engineering, 2022, 19(2): 1786-1824. doi: 10.3934/mbe.2022084 |
[5] | Yi Wang, Shicheng Zhong, Guo Wang . Preventing online disinformation propagation: Cost-effective dynamic budget allocation of refutation, media censorship, and social bot detection. Mathematical Biosciences and Engineering, 2023, 20(7): 13113-13132. doi: 10.3934/mbe.2023584 |
[6] | Lili Liu, Xi Wang, Yazhi Li . Mathematical analysis and optimal control of an epidemic model with vaccination and different infectivity. Mathematical Biosciences and Engineering, 2023, 20(12): 20914-20938. doi: 10.3934/mbe.2023925 |
[7] | Mohammed H. Alharbi . Global investigation for an "SIS" model for COVID-19 epidemic with asymptomatic infection. Mathematical Biosciences and Engineering, 2023, 20(3): 5298-5315. doi: 10.3934/mbe.2023245 |
[8] | Jiajia Zhang, Yuanhua Qiao, Yan Zhang . Stability analysis and optimal control of COVID-19 with quarantine and media awareness. Mathematical Biosciences and Engineering, 2022, 19(5): 4911-4932. doi: 10.3934/mbe.2022230 |
[9] | Chenwei Song, Rui Xu, Ning Bai, Xiaohong Tian, Jiazhe Lin . Global dynamics and optimal control of a cholera transmission model with vaccination strategy and multiple pathways. Mathematical Biosciences and Engineering, 2020, 17(4): 4210-4224. doi: 10.3934/mbe.2020233 |
[10] | Liang Tian, Fengjun Shang, Chenquan Gan . Optimal control analysis of malware propagation in cloud environments. Mathematical Biosciences and Engineering, 2023, 20(8): 14502-14517. doi: 10.3934/mbe.2023649 |
This research explored the life cycle analysis and environmental cost-benefit assessment of converting ash waste from hospital medical waste incineration into environmentally safe paving block raw materials. The growing concerns about medical waste disposal and its environmental impact necessitate innovative solutions for sustainable waste management. This research aimed to evaluate the feasibility and environmental implications of reusing hospital waste into raw materials for paving block mixtures. This research, a comprehensive life cycle analysis, examined the environmental impacts of medical waste collection for the production and use of paving blocks. Additionally, we conducted an environmental cost-benefit assessment to ascertain the economic feasibility and potential environmental impact forecasts of this recycling approach. The research results show that converting hospital medical waste ash into mixed raw materials for paving blocks not only immobilizes heavy metals but also provides a sustainable alternative for non-building materials. These findings highlight the potential for significant environmental and economic benefits, making this approach a promising strategy for waste management and sustainable construction practices. The cost of preventing environmental damage (eco-cost) in the process of converting ash from the incineration of medical waste into a mixture of raw materials for paving blocks is IDR 600,180.9 per cycle.
With the rapid development of social networks, the spread of information becomes rapid and extensive, including the spread of rumors. The spread of rumors can cause significant economic losses and even bad political influence[1]. Research on spread of rumors has a long tradition. For decades, one of the most popular ideas has been the analysis of the dynamics of rumor propagation [2,3,4,5,6,7]. Moreover, the main purpose of most studies on rumor propagation has focused on the study of immunization strategies against rumors [8,9,10,11,12,13].
One of the main theoretical and conceptual frameworks used to model the spread of rumors is to model rumors as infectious diseases. This has been widely adopted in the field of spread of rumors research. The first serious discussion and analysis of rumor propagation appeared in the 1960s by Daley and Kendall, which is called DK model[14,15]. In this model, people in rumor spreading networks are divided into three states: Susceptible (S) are those who are not having been exposed to rumor; Infected (I) are those who are being exposed to rumor and believed to spread it. Removed(R) are those who have no interest to spread rumor. In 2001, Zanette et al.[16,17] studied the propagation rule of rumor in small-world network based on SIR model, and proved the existence of critical threshold for rumor propagation for the first time. Moreno et al.[18] in 2004 studied the propagation behavior of rumors based on SIR model on scale-free complex networks. This research approach was gradually extended to model the evolution of SIR. A number of authors [19,20] have considered the effects of influence of media on rumor spreading. Several previous studies [21,22,23] have also explored the relationship between rumor propagation and hesitation mechanisms. There is a large number of published studies [24,25] that describe the communication model of competitive information in multiplex social networks. Liu et al.[26,27] conducted a series of trials in which he mixed comments and rumor spreading in the network.
In terms of information dissemination modeling, it seems to be a common treatment to use information dissemination among friends as the main dissemination channel. As far as we know, few studies have focused on modeling information dissemination through non-friends. In fact, users in social networks not only rely on their friends to spread information, but also receive information by non-friend. In view of this situation, this paper discusses the rumor propagation model under two kinds of propagation modes (friend propagation and marketing account propagation). In some social networks, such as TikTok and other popular APPs, users can not only receive information posted by friends, but also receive a large number of non-friends' information, mainly marketing account information. In fact, the marketing account is a professional account for forwarding and commenting on hot issues in the society or a certain circle. These accounts are well-informed and have a large number of followers, so its information is more likely to be received by others, although it may not be his followers. Therefore, rumors once forwarded by marketing accounts will quickly ferment in the network. Based on the above, a main research problem of this paper is to establish a suitable rumor propagation model and control strategy.
The main contributions of this paper are as follows:
1. First, the rumor propagation model under the coexistence of two communication modes was established, and steady-state analysis was carried out.
2. Second, considering the negative losses caused by rumors and the cost of controlling rumor, the optimal control strategy to reduce rumor propagation density by suppressing individual accounts and marketing accounts is discussed.
3. Third, this paper compares the correlation between the spread rate of the rumor spread of friends accounts and marketing accounts when the network average degree is different.
4. Fourth, through the sensitivity analysis of parameters, the influences of different parameters including initial value on rumor propagation is analyzed and compared.
5. Finally, the superiority over the optimal control strategy is proved by simulation experiment.
The subsequent materials are organized in this fashion. In Section 2, we establish the Susceptible-Infected-Marketing-Removed (SIMR) model based on the coexistence of friend propagation and marketing account propagation, and conduct a steady state analysis. In Section 3, we perform theoretical analysis on the optimal control problem and a dynamic strategy for rumor immunity. In Section 4, the stability of SIMR model is verified by simulation experiments, and the influence of parameters and initial value on rumor propagation is discussed, and the influence of optimal control strategies for suppressing individual and marketing accounts on rumor propagation with different average degree is discussed. A brief conclusion is drawn in Section 5.
Inspired by the classic SIR model of infectious diseases, in order to describe the spread of rumors in a social network with a public marketing account (which can push information to a non-friend account), we propose the Susceptible-Infected-Marketing-Removed (SIMR) model. We divide the people in the network into four states: Susceptible (people who have not been exposed to rumors); Infected (people who believe in rumors and spread them through friends); Marketing (people who believe in rumors and push them through non-friends); Removed (people who have no interest to spread rumors). S(t),I(t),M(t),R(t) are the relative density of susceptible, infected, marketing, removed nodes at time t, respectively. (1−S(t)−I(t)−M(t)−R(t)) denotes the density of the empty nodes that can make the new users transplant into the online social networks with a certain constant rate b. For this purpose, we consider the following ordinary differential equations.
{dS(t)dt=b(1−S(t)−I(t)−M(t)−R(t))−λ1S(t)I(t)⟨k⟩−λ2M(t)S(t)−dS(t),dI(t)dt=λ1S(t)I(t)⟨k⟩+λ2M(t)S(t)−δI(t)−dI(t)−γI(t),dM(t)dt=γI(t)−βM(t)−dM(t),dR(t)dt=δI(t)+βM(t)−dR(t),S(0)≥0,I(0)≥0,M(0)≥0,R(0)≥0,0≤t≤T. | (2.1) |
Where ⟨k⟩=∑nk=1kP(k) implies the average degree in complex network and P(k) stands for a connective distribution function. As shown in Figure 1, the rumor spreading rules of the model can be summarized as follows:
1. We consider the non-closed nature of social networks. Assume that the rate of joining network and exiting network is b and d, respectively.
2. When a susceptible node contacts infected nodes, the susceptible node turns into an infected node with rate λ1. When the susceptible node contacts the marketing accounts, the susceptible node becomes the infected node with rate λ2.
3. Infected nodes are converted to marketing nodes at a certain rate γ to spread messages.
4. An infected node becomes removed node with the rate δ when contacts removed nodes or loses interest in spreading rumors.
5. A marketing account becomes removed node with the rate β when loses interest in spreading rumors.
Parameters | Description | Range(/days) |
b | Join networks rate | 0.04−0.4 |
d | Exit networks rate | 0.04−0.4 |
λ1 | Rate of a susceptible node contacts the infected node | 0.009−0.9 |
λ2 | Rate of a susceptible node contacts the Marketing node | 0.009−0.9 |
δ | Rate from infected node to removed node | 0.009−0.9 |
γ | Rate from infected node to Marketing node | 0.009−0.9 |
β | Rate from Marketing node to removed node | 0.009−0.9 |
The vector E(t)=(S(t),I(t),M(t),R(t)) represents the relative density of state in the network at time t. The system (2.1) can be written as
{dE(t)dt=f1(E(t)),E(0)≥0,0≤t≤T. | (2.2) |
Taking I(t)=M(t)=0 in system (2.1), the rumor-free equilibrium point is E0=(b/(b+d),0,0,0).
In general, the basic reproduction number R0 means the average number of infections in a purely susceptible population[28]. Let χ=(I,M,R,S)T, then the system (2.1) can be written as χ′=Ψ(χ)−Φ(χ), where
Ψ(χ)=(λ1SI⟨k⟩+λ2MS000),Φ(χ)=((δ+d+γ)I−γI+βM+dM−δI−βM+dR−b+b(S+I+M+R)+λ1⟨k⟩SI+λ2MS+dS). | (2.3) |
The Jacobian matrices of Ψ and Φ evaluated at the rumor-free equilibrium χ0=(I0,M0,R0,S0)=(0,0,0,b/(b+d)) are given by
J(Ψ|χ0)=(F000),J(Φ|χ0)=(V0V1V2), | (2.4) |
where
F=(λ1⟨k⟩S0λ2S000),V=(d+δ+γ0−γβ+d), | (2.5) |
V1=(−δ−βb+λ1⟨k⟩S0b+λ2S0),V2=(d0bb+d). | (2.6) |
Then
FV−1=(λ1⟨k⟩S0λ2S000)⋅(1d+δ+γ0γ(d+δ+γ)(β+d)1β+d)=(λ1⟨k⟩b(d+δ+γ)(b+d)+λ2bγ(d+δ+γ)(β+d)(b+d)λ2b(β+d)(b+d)00). | (2.7) |
Therefore, the basic reproduction number R0 defined by the spectral radius FV−1 [29] is
R0=ρ(FV−1)=λ1⟨k⟩b(d+δ+γ)(b+d)+λ2bγ(d+δ+γ)(β+d)(b+d). | (2.8) |
For the analytical and numerical computation of R0 in general structured population models see e.g., [30,31].
The rumor-free equilibrium point of SIMR model and the basic reproduction number of the model have been calculated. When the system reaches the rumor-free equilibrium point E0=(b/(b+d),0,0,0), the rumor disappears. And when rumor persistently spreads, it means the system has reached the rumor-prevailing equilibrium point E1=(S1,I1,M1,R1). The rumor-prevailing equilibrium should satisfy
{b(1−S(t)−I(t)−M(t)−R(t))−λ1S(t)I(t)⟨k⟩−λ2M(t)S(t)−dS(t)=0,λ1S(t)I(t)⟨k⟩+λ2M(t)S(t)−δI(t)−dI(t)−γI(t)=0,γI(t)−βM(t)−dM(t)=0,δI(t)+βM(t)−dR(t)=0. | (2.9) |
Then the rumor-prevailing equilibrium E1 has the following form
{S1=(β+d)(δ+d+γ)λ1(β+d)⟨k⟩+λ2γ,I1=bd(b+d)(γ+δ+d)−(β+d)dλ1(β+d)⟨k⟩+λ2γ,M1=γ⋅I1β+d,R1=(δd+βγd(β+d))I1. | (2.10) |
Since
I1=(β+d)dλ1(β+d)⟨k⟩+λ2γ(λ1(β+d)⟨k⟩+λ2γ(β+d)(b+d)(d+γ+δ)−1)=(β+d)dλ1(β+d)⟨k⟩+λ2γ(R0−1). |
Thus, the positive rumor-prevailing equilibrium point E1=(S1,I1,M1,R1) exists if R0>1, which satisfies (2.9). In another word, when the system (2.1) reaches stability, the rumor will disappear if R0<1. On the contrary, when R0>1, rumor continues to spread in the system.
Theorem 2.1. If R0<1, the rumor-free equilibrium point E0 of system (2.1) is locally asymptotically stable.
Proof. We can find the Jacobian matrix J(E0), which is a 4×4 matrix, is in the following form
J(E0)=(J11J12J13J14J21J22J23J24J31J32J33J34J41J42J43J44)=(−b−d−b−λ1b⟨k⟩b+d−b−λ2bb+d00λ1b⟨k⟩b+d−δ−d−γλ2bb+d00γ−β−d00δβ−d) | (2.11) |
→(−b−dJ′12−b−λ2bb+d00J′22λ2bb+d000−β−d0000−d). | (2.12) |
Where J′12,J′22 are obtained by the elementary transformation,
{J′12=−b(1+γβ+d)−bλ1⟨k⟩b+d−bγλ2(b+d)(β+d),J′22=b(λ1(β+d)⟨k⟩+λ2γ)(b+d)(β+d)−(δ+d+γ). | (2.13) |
We can easily conclude that J(E0) is an upper triangular matrix, and its four eigenvalues are −b−d,−β−d,−d,b(λ1(β+d)⟨k⟩+λ2γ)(b+d)(β+d)−(δ+d+γ) respectively. Obviously the first three eigenvalues are all negative. And then the last eigenvalue J′22<0 if and only if R0<1, since
b(λ1(β+d)⟨k⟩+λ2γ)(b+d)(β+d)−(δ+d+γ)=λ1⟨k⟩b(β+d)+λ2bγ−(d+δ+γ)(β+d)(b+d)(β+d)(b+d)<0⇔R0−1=λ1⟨k⟩b(d+δ+γ)(b+d)+λ2bγ(d+δ+γ)(β+d)(b+d)−1=λ1⟨k⟩b(β+d)+λ2bγ−(d+δ+γ)(β+d)(b+d)(d+δ+γ)(β+d)(b+d)<0. | (2.14) |
Hence, the rumor-free equilibrium E0 is locally asymptotically stable if R0<1 by Routh-Hurwitz criterion[32], completing the proof.
So E0 is locally asymptotically stable. If R0>1, the matrix J(E0) has a positive eigenvalue, so E0 is unstable. In the following, we consider the globally asymptotically stability of the rumor-free equilibrium E0, which is the important content of the stability analysis.
Lemma 2.1. Since there are S(0),I(0),M(0),R(0)≥0, E=(S,I,M,R) is a solution of the system (2.1) and N=S+I+M+R, then the solutions of the model are uniformly bounded.
Proof. From the equations of system (2.1), we can easily observed that
dN(t)dt=b−b(S(t)+I(t)+M(t)+R(t))−d(S(t)+I(t)+M(t)+R(t))=b−(b+d)N(t), |
It is easy to see that
0≤limt→∞supN(t)≤bb+d. |
It follows from the first equation of system (2.1) as
dSdt=b(1−N(t))−λ1S(t)I(t)⟨k⟩−λ2M(t)S(t)−dS(t)≥bdb+d−λ1S(t)I(t)⟨k⟩−λ2M(t)S(t)−dS(t), |
It is easy to see that
dS(t)dt+(λ1I(t)⟨k⟩+λ2M(t)+d)S(t)≥bdb+d,S(t)≥e−∫t0(λ1I(u)⟨k⟩+λ2M(u)+d)du[∫t0bdb+de∫t0(λ1I(u)⟨k⟩+λ2M(u)+d)dudt+S(0)]≥0. |
In the following, we proof that I(t),M(t)≥0 for all t>0. Assume there exists t∗1,t∗2 such that I(t∗1),M(t∗2) is negative, and t1=sup{t>0:I(t)>0,M(t)>0,∈[0,t]}. Thus, t1>0. There will be three cases as follows:
1) We have t1 such that I(t1)=M(t1)=0.
2) We have t1 such that I(t1)>0,M(t1)=0.
3) We have t1 such that I(t1)=0,M(t1)>0.
In 1), it easy to obtain I(t)=M(t)=0 for all t>t1.
In 2), for all t∈[0,t1), I(t)>0,M(t)>0 holds. By the third equation of system (2.1), we obtain
dM(t1)dt=γI(t1)>0. |
Thereby, there exists a sufficiently small positive constant ε such that for any t∈(t1−ε,t1), with M(t)<0 holds. This is in contradiction with case 2), then M(t)≥0 for all t>0. It is easy to obtain that the case 3) also contradicts the conditions known, then I(t)≥0 for all t>0. In similar fashion it can be shown that R(t)≥0 for all t>0.
So all solutions of the model are confined in the region {(S(t),I(t),M(t),R(t))∈R4+∪{0}:N(t)=bb+d}. Therefore, Lemma 2.1 has been proved. Therefore, the positive invariant set of the system (2.1) is Ω, where
Ω={(S,I,M,R)|(S(t),I(t),M(t),R(t))∈R4+∪{0}:0≤S(t),I(t),M(t),R(t)≤bb+d,t≥0}. |
Here, we use the method developed by Castillo-Chavez et al.[33]. we list two conditions that if met, also guarantee the global asymptotic stability of the rumor-free state. We rewrite the model system (2.1) as
{dxdt=F(x,H),dHdt=G(x,H),G(x,0)=0, | (2.15) |
where x∈Rm denotes (its components) the number of uninfected individuals including susceptible, removed, et al. and H∈Rn denotes (its components) the number of infected individuals including infected, et al. u0=(x∗,0) denotes the rumor-free equilibrium of this system.
Lemma 2.2. [33] If the equilibrium point u0=(x∗,0) of (2.1) is locally asymptotically stable when R0<1, the fixed point U0=(x∗,0) is a globally asymptotic stable equilibrium of (2.1) provided that R0<1 and that assumptions (H1) and(H2) are satisfied.
(H1) For dxdt=F(x,0), x∗ is globally asymptotically stable,
(H2) G(x,H)=AH−ˆG(x,H), ˆG(x,H)≥0 for (x,H)∈Ω, where A=DHG(x∗,0) is a M-matrix (the off diagonal elements of A are nonnegative) and Ω is the region where the model is biological sense.
Theorem 2.2. If R0<1, then the rumor-free equilibrium E0=(S0,I0,M0,R0) of system (2.1) is globally asymptotically stable.
Proof. Let x=(S,R),H=(I,M), u0=(S0,R0,I0,M0)=(x∗,0), where x∗=(bb+d,0), then
dxdt=(dSdt,dRdt)=(b−(b+d)S,−dR). |
When S(t) is equal to bb+d=S0, R(t) is equal to 0, we can obtain F(x,0)=0. As t→+∞, there are S(t)→bb+d, R(t)→0. Hence x∗=(S0,0) is globally asymptotically stable.Thus the conditions (H1) is satisfied. And G(x,H)=AH−ˆG(x,H), where
A=(λ1⟨k⟩S0−(δ+d+γ)λ2S0γ−β−d),ˆG(x,H)=(λ1⟨k⟩(S0−S)I+λ2(S0−S)M0). | (2.16) |
According to the Lemma 2.1, it can be obtained ˆG(x,H)≥0, and A is a M-matrix; the conditions (H2) is satisfied, and by Lemma 2.2, the rumor-free equilibrium E0 is globally asymptotically stable if R0<1.
This means that when R0<1, there are no more rumors in the network after the system is stabilized, i.e., the rumors disappear.
Theorem 2.3. Let R0>1, the rumor-prevailing equilibrium E1=(S1,I1,M1,R1) (2.10) of system (2.1) is asymptotically stable.
Proof. The Jacobian matrix of the model at E1(S1,I1,M1,R1) is given by
J(E1)=(−b−d−λ1⟨k⟩I1−λ2M1−b−λ1⟨k⟩S1−b−λ2S1−bλ1⟨k⟩I1+λ2M1λ1⟨k⟩S1−(δ+d+γ)λ2S100γ−β−d00δβ−d). | (2.17) |
The characteristics equation of J(E1) is
λ4+A1λ3+A2λ2+A3λ+A4=0, |
where
{M1=J11=−b−d−λ1⟨k⟩I1−λ2M1<0,M2=J22=λ1⟨k⟩S1−(δ+d+γ)=(δ+d+γ)(λ1⟨k⟩(β+d)λ1⟨k⟩(β+d)+λ2γ−1)<0,M3=J33=−β−d<0,M4=J44=−d<0, |
A1=−(J11+J22+J33+J44)>0. | (2.18) |
{M12=−(b+d+λ1⟨k⟩I1+λ2M1)(λ1S1⟨k⟩−δ−d−γ)+(b+λ1⟨k⟩S1)(λ1⟨k⟩I1+λ2M1)>0,M13=(b+d+λ1⟨k⟩+λ2M1)(β+d)>0,M14=(b+d+λ1⟨k⟩I1+λ2M1)d>0,M23=−[λ1⟨k⟩S1−(δ+d+γ)](β+d)−γλ2S1=−(δ+d+γ)(β+d)[λ1⟨k⟩(β+d)λ1⟨k⟩(β+d)+λ2γ−1]−γλ2S1=λ2γ(δ+d+γ)(β+d)λ1⟨k⟩(β+d)+λ2γ−γλ2S1=0,M24=[λ1⟨k⟩S1−(δ+d+γ)](−d)>0,M34=(β+d)d>0, |
A2=M12+M13+M14+M23+M24+M34>0. | (2.19) |
{M123=(b+d)[(β+d)(λ1⟨k⟩S1−(δ+d+γ))+λ2γS1]−(λ1⟨k⟩I1+λ2M1)[(b+δ+d+γ)(β+d)+bγ]=−(b+d)M23−(λ1⟨k⟩I1+λ2M1)[(b+δ+d+γ)(β+d)+bγ]<0,M124=(b+d+λ1⟨k⟩I1+λ2M1)d(λ1S1⟨k⟩−(δ+d+γ))−(λ1⟨k⟩I1+λ2M1)[d(b+λ1⟨k⟩S1)+bδ]<0,M234=−d[(β+d)(δ+d+γ−λ1⟨k⟩S1)−λ2γS1]=0, |
A3=M123+M124+M234>0. | (2.20) |
A4=det(J)=−(b+d+λ1⟨k⟩I1+λ2M1)M234−(λ1⟨k⟩I1+λ2M1)[−b(γβ+(β+d)δ)−d((b+λ1⟨k⟩S1)(β+d)+(b+λ2S1)γ)]>0. | (2.21) |
Therefore, according to The Routh-Hurwitz criterion, when R0>1, the system (2.1) is asymptotically stable.When R0>1, the system stabilizes at the rumor-prevailing equilibrium point, i.e., rumors in the network will continue to spread.
In this section, based on the SIMR model, there are two different ways of spreading rumors: infection between friends and marketing account push. Under the premise of considering the cost of rumor control, this paper proposes the optimal control strategy for controlling the spread of rumors by controlling individual and marketing accounts. Consider Θ(t)=(θ1(t),θ2(t)),0≤t≤T as the control variable of the SIMR problem, where θ1(t),θ2(t) represents the control variable for the individual account and marketing account, respectively. The upper limit of this strategy is determined as shown in the following remarks:
Remark 3.1. Suppose the feasible region of Θ(t),0≤Θ≤¯Θ,t∈(0,T], where ¯Θ are the upper bounds of Θ. The upper bound ¯Θ is determined the budgeted costs of immunization rumors.
So we assume the SIMR strategy is
{Θ=(θ1(t),θ2(t))∈L[0,T]2|0≤θ1(t)≤¯θ1,0≤θ2(t)≤¯θ2,0≤t≤T}, |
where L[0,T]2 represents the set of Lebesgue integrable functions defined on [0,T] [34]. In this paper, under the premise of controlling the cost of rumor, the expected cost effectiveness caused by rumor can be minimized. Assume that the loss caused by rumor in the network is J1, and the cost of rumor control in the network is J2. The following remark can be obtained.
Remark 3.2. Assuming that the unit loss caused by spreading rumors in individual accounts and marketing accounts per unit time is a constant c1 and c2, respectively. Then the total loss caused by the spreading of rumor in the time horizon [0,T] can be expressed as
J1(Θ)=∫T0c1I(t)+c2M(t)dt. | (3.1) |
Remark 3.3. Assuming that the unit cost of controlling individual accounts and marketing accounts per unit time is a constant c3 and c4, respectively. Then the total cost in the time horizon [0,T] can be expressed as
J2(Θ)=∫T0c3θ1(t)+c4θ2(t)dt. | (3.2) |
In summary, the expected cost effectiveness of rumor spread is
J(Θ)=J1(Θ)+J2(Θ)=∫T0c1I(t)+c2M(t)+c3θ1(t)+c4θ2(t)dt≡∫T0F(E(t),Θ(t))dt. | (3.3) |
Based on the above description, the rumor immunity problem of the SIMR model is established as the following optimal control problem
minθ∈Θ∫T0F(E(t),Θ(t))dtsubjectto | (3.4) |
{dS(t)dt=b(1−S(t)−I(t)−M(t)−R(t))−λ1S(t)I(t)⟨k⟩−λ2M(t)S(t)−dS(t),dI(t)dt=λ1S(t)I(t)(1−θ1(t))⟨k⟩+λ2M(t)(1−θ2(t))S(t)−δI(t)−dI(t)−γI(t),dM(t)dt=γI(t)−βM(t)−dM(t),dR(t)dt=δI(t)+βM(t)−dR(t)+λ1S(t)I(t)θ1⟨k⟩+λ2M(t)θ2(t)S(t),S(0)≥0,I(0)≥0,M(0)≥0,R(0)≥0,0≤t≤T. | (3.5) |
In order to solve the optimal control problem, we adopted the optimal principle of Pontryagin [35] and define the Lagrangian and Hamilton function of the optimal control problem as follows:
H=c1I(t)+c2M(t)+c3θ1(t)+c4θ2(t)+μ1dSdt+μ2dIdt+μ3dMdt+μ4dRdt=c1I(t)+c2M(t)+c3θ1(t)+c4θ2(t)+μ1[b(1−S(t)−I(t)−M(t)−R(t))−λ1S(t)I(t)⟨k⟩−λ2M(t)S(t)−dS(t)]+μ2[λ1S(t)I(t)(1−θ1(t))⟨k⟩+λ2M(t)(1−θ2(t))S(t)−δI(t)−dI(t)−γI(t)]+μ3[γI(t)−βM(t)−dM(t)]+μ4[δI(t)+βM(t)−dR(t)+λ1S(t)I(t)θ1⟨k⟩+λ2M(t)θ2(t)S(t)], | (3.6) |
where μ1,μ2,μ3,μ4 are the adjoint functions. We obtain the necessary conditions for optimal control of SIMR problems as follows.
Theorem 3.1. Suppose Θ(t)={θ1(t),θ2(t)} is an optimal control of the SIMR problem (3.4), E(t)=(S(t),I(t),M(t),R(t)) is the solution to the associated rumor spreading model (3.5). Then there exists an adjoint function μ(t)=(μ1(t),μ2(t),μ3(t),μ4(t)) such that the following equations hold.
{dμ1dt=μ1(b+λ1⟨k⟩I+λ2M+d)−μ2[λ1⟨k⟩I(1−θ1)+λ2M(1−θ2)]−μ4(λ1Iθ1⟨k⟩+λ2Mθ2),dμ2dt=−c1+μ1(b+λ1⟨k⟩S)−μ2(λ1⟨k⟩S(1−θ1)−δ−d−γ)−μ3γ−μ4(δ+λ1Sθ1⟨k⟩),dμ3dt=−c2+μ1(b+λ2S)−μ2(λ2S(1−θ2))+μ3(β+d)−μ4(β+λ2Sθ2),dμ4dt=μ1b+μ4d,0≤t≤T,μ1(T)=μ2(T)=μ3(T)=μ4(T)=0. | (3.7) |
We obtain the optimal control ˜Θ=(~θ1,~θ2) as follow:
~θ1={¯θ1,g1(θ1)<0,0,g1(θ1)>0, | (3.8) |
~θ2={¯θ2,g2(θ2)<0,0,g2(θ2)>0, | (3.9) |
where g1(θ1)=c3−λ1IS⟨k⟩μ2+μ4λ1IS⟨k⟩,g2(θ2)=c4−λ2MSμ2+μ4λ2MS.
Proof. According the Minimum principle of Pontryagin, there exists μ=(μ1,μ2,μ3,μ4) such that
{dμ1dt=−∂H∂S,0≤t≤T,dμ2dt=−∂H∂I,0≤t≤T,dμ3dt=−∂H∂M,0≤t≤T,dμ4dt=−∂H∂R,0≤t≤T. | (3.10) |
By the optimal conditions, we have
{∂H∂θ1=c3−λ1IS⟨k⟩μ2+μ4λ1IS⟨k⟩,∂H∂θ2=c4−λ2MSμ2+μ4λ2MS. | (3.11) |
According to the Pontryagin Minimum Principle, the optimal solution of the objective function in the time horizon [0,T] is Θ(t)=argmin˜θ∈ΘH(E(t),¯Θ,μ(t)).
In this section, we introduce the synthetic scale-free network. The degree distribution of scale-free network follows the power law property with P(k)∼ak−3, n=500, and the constant a satisfis ∑nk=1P(k)=1. By a simple calculation, we can conclude that the average degree of the scale-free network structure ⟨k⟩=∑nk=1kP(k)=1.367. The initial conditions are given by S(0)=0.7, I(0)=0.2, M(0)=0.05, R(0)=0. To verify the accuracy of the model, we have conducted numerical computation using the Runge-Kutta method and MATLAB by setting parameter value in the system[36,37].
Example 1: Stabilities of the equilibrium points.
In this part, we will verify the impact on R0 on the stability of system (2.1). When the parameter reference value is the first row of Table 2, we calculate the basic reproduction number R0=0.41<1. When the parameter reference value is the second row in Table 2, we calculate the basic reproduction number R0=1.43>1. Figures 2(a), (b) indicate the stability of individual ratios when R0<1 and R0>1, respectively. Figure 2(c) shows that the stability of nodes density when R0<1 (R0=0.41 and R0=0.89).
Parameters | b | λ1 | λ2 | d | δ | γ | β | R0 |
Data | 0.4 | 0.1 | 0.1 | 0.04 | 0.1 | 0.1 | 0.1 | 0.41 |
0.4 | 0.3 | 0.3 | 0.04 | 0.1 | 0.1 | 0.1 | 1.43 | |
0.4 | 0.12 | 0.1 | 0.04 | 0.1 | 0.1 | 0.1 | 0.89 |
According to Theorem 2.2, system (2.1) is globally asymptotically stable in E0. According to Theorem 2.3, system (2.1) is asymptotically stable in E1. Moreover, in the same case, we take R0<1 as an example, when R0 is larger, the transmission rate of infected individuals tending to 0 is slower. The density of the infected individuals and R0 has a positive linear relationship, and the greater the value of the R0, the greater the value of the density of the infected individuals. Finally, Figure 2 reveals the relative volatility of the early systems of rumor. As time goes by, the system tends to be stable, so controlling the rumors in the early stages often leads in better results.
Example 2: The effect of initial value on system stability.
Figures 3 and 4 indicate the influence of different initial S(0) and I(0) on rumor propagation results. The initial value is that the density of individual spreaders I(0) ranges from 0.05 to 0.45 and the density of susceptible S(0) ranges from 0.9 to 0.5, M(0)=0.05,R(0)=0. It can be seen that although the initial value of the system is different, the system (2.1) stabilized at the same value. In Figure 3, the infected individual is stable at 0, and in Figure 4 the system is stable at the rumor-prevailing equilibrium solution E1(S1,I1,M1,R1).
Figures 3 and 4 reveal that when the basic reproduction number R0 of the rumor system is fixed, the initial value (S(0),I(0),M(0),R(0)) of the individual in the network will not affect the stability of the rumor. Because the basic number of reproduction means the speed of rumor propagation, it is the decisive factor in the continued spread of rumors.
Example 3: The effect of parameters on rumor diffusion.
Figures 5–7 indicate the effect of the parameter λ1,λ2,γ,δ,β on rumor diffusion. Parameter λ1,λ2 respectively represents the conversion rate from susceptible individuals to infected individuals through the influence of friends and marketing accounts. Obviously, it can be seen from the Figures 5(a), (b) that the higher the value of the parameter, the higher the density of infection, and ultimately the higher the rate of infected individuals in the network.
Parameter γ represents the rate that the infected individual converted into a marketing account. It can be seen by Figures 6(a), (b) that both R0<1 or R0>1, the change of γ has little effect on rumors in the network.
Parameter δ and β respectively represent the rate that individual accounts and marketing accounts lose interest in rumors and thus become removed nodes. It can be seen from Figures 7(a), (b) that the rate of rumor spreading individuals decreases with the increase of parameters. In particular, when the parameter is small, the rate of rumor spreading individuals decreases significantly with the increase of parameter.
In conclusion, infected individuals increased with λ1 and λ2, and decreased with γ, δ, β. As can be seen from Figure 2, the relationship between the density of infected individuals and basic reproduction number R0 is a positive correlation. Then in the Example 4, the parameter sensitivity analysis on R0 is verified.
Example 4: Sensitivity analysis of parameter on R0.
In this part, we focus on the effect of parameters on R0. As can be seen by (2.8) and Figure 5, the density of infected individual increases when λ1 and λ2 increases. Figure 8 refers to the impact of λ1 and λ2 on the basic reproduction number R0 in the network of rumor spreading when the average degree ⟨k⟩ is different.
As can be seen from Figure 8, the larger the ⟨k⟩, the greater the number of friends in the network node, and the easier it is for rumors to spread among friends. This means that when the connection between nodes in the network is relatively small, the spread of rumors is more affected by marketing account push. On the contray, when the average degree of the networks is greater, rumors are easier to spread through friends. This is why marketing accounts encourage ordinary users to focus on their accounts.
Figure 9 shows that δ has the most significant effect on R0. Moreover, when δ is less than 0.4, the value of R0 decreases significantly as δ increase, and similarly, when β is less than 0.2, the value of R0 decreases significantly as β increases. Therefore, it is a very effective way to control rumors by enhancing the immune rate of friends accounts and marketing accounts.
Example 5: The effects of the optimal control strategies.
It can be concluded from Figure 9 that the spread of rumors can be effectively controlled by increasing the rate of rumor-infected individuals and marketing accounts converting to removed nodes. In this part, we discuss the control strategy of controlling rumors by suppressing individual accounts and marketing accounts. We simulated the effect of the rumor control strategy on three synthetic BA scale-free networks. The scale-free network with an average degree ⟨k⟩=2.04,4.06,7.92 were selected. The degree distribution of scale-free network is shown in Figure 10.
Figure 11 shows the comparison results of optimal control strategy and uniform control strategy under three networks with different average degree ⟨k⟩=2.04,4.06,7.92 when c1=100,c2=100, c3=10 and c4=10. The red line represents the value of J(Θpoc) and the black line represents the minimum value of J(Θp,q),p=0,0.1,...,1,q=0,0.1,...,1. It can be seen that the red function value is lower than the black function value. It is can conclude that Θpoc is superior to all the uniform control in three network in terms of the cost of rumor control.
Controlling cost changes in individual accounts and marketing accounts will have an impact on control. When ⟨k⟩=4.06, the changes of costs c3 and c4 are shown in Table 3, the corresponding value of J(minΘp,q) and J(Θpoc) are shown in Figure 12. In Figure 12, regardless of c3 and c4, the cost of optimal control is lower than the minimum cost of uniform control. And when c3=c4>40, the minimum value of uniform control remains unchanged, because when the control cost is too high, even using low-intensity control will greatly increase the total cost. Therefore, the minimum cost of unified control strategy is Θ0,0. On the contrary, the optimal control can dynamically adjust the control strategy according to the different cost of control.
Parameters | Value1 | Value2 | Value3 | Value4 | Value5 |
c3 | c3=10 | c3=20 | c3=30 | c3=40 | c3=50 |
c4 | c4=10 | c4=20 | c4=30 | c4=40 | c4=50 |
Figure 13 shows that under optimal control, the density of rumor propagation nodes, including individual accounts and marketing accounts, decreases by varying degrees. As can be seen from Figure 13, the optimal control strategy of rumors significantly reduces the density of infected individuals, so as to achieve a satisfactory control effect. Combined with Figures 11–13, the optimal control strategy is an ideal strategy to control the density of infected individuals. On the other hand, the optimal control strategy is also the best choice from the perspective of economic benefits. This also provides a reference value for controlling rumors.
In this paper, we have proposed a model of information propagation based on two kinds of propagation types. In the network, the state of the user was divided into S,I,M,R. Different from the traditional SIR Information communication model, it increased the marketing state, and could push information to non-friends. First, we calculated the basic reproduction number R0 by the method of the next generation matrix. In addition, we discussed the existence of the rumor-free equilibrium point and rumor-prevailing equilibrium point, and proved the globally asymptotically stability of the rumor-free equilibrium point when R0<1, and the asymptotically stability of rumor-prevailing equilibrium point when R0>1. More importantly, we proposed an optimal control strategy for rumors. Finally, the correctness of the above theory was verified by numerical simulation. Firstly, we verified the stability of the model and discussed the impact of initial value on the stability of the rumor. Secondly, we used sensitivity analysis to discuss the impact of parameters on R0 and draw two conclusions. On the one hand, the influence of λ1 on R0 increased with the increase of the average degree of the network. On the other hand, the rate δ of I to R was the most significant effect on R0. The rate β of I to R was the second effect of R0. Finally, based on the above conclusions, we proposed the optimal control strategy of the two kinds of immunity and verified the superiority of the optimal control strategy, and provided the reference for the control of rumors.
This work is supported by the National Natural Science Foundation of China under Grant Nos. 61807028, 61772449 and 61802332. The authors are grateful to valuable comments and suggestions of the reviewers.
The authors declare no potential conflict of interests.
[1] |
Ma Y, Lin X, Wu A, et al. (2020) Suggested guidelines for emergency treatment of medical waste during COVID-19: Chinese experience. Waste Dispos Sustain En 2: 81–84. https://doi.org/10.1007/s42768-020-00039-8 doi: 10.1007/s42768-020-00039-8
![]() |
[2] | Wang Y, Narayanan M, Shi X, et al. (2022) Plant growth-promoting bacteria in metal-contaminated soil: Current perspectives on remediation mechanisms. Front Microbiol 13: 966226. |
[3] |
Chen Y, Liu L, Feng Q, et al. (2012) Key issues study on the operation management of medical waste incineration disposal facilities. Procedia Environ Sci 16: 208–213. https://doi.org/10.1016/j.proenv.2012.10.029 doi: 10.1016/j.proenv.2012.10.029
![]() |
[4] |
Vieira DR, Calmon JL, Coelho FZ (2016) Life cycle assessment (LCA) applied to the manufacturing of common and ecological concrete: A review. Constr Build Mater 124: 656–666. https://doi.org/10.1016/j.conbuildmat.2016.07.125 doi: 10.1016/j.conbuildmat.2016.07.125
![]() |
[5] |
Unger S, Landis A (2016) Assessing the environmental, human health, and economic impacts of reprocessed medical devices in a Phoenix hospital's supply chain. J Clean Prod 112: 1995–2003. https://doi.org/10.1016/j.jclepro.2015.07.144 doi: 10.1016/j.jclepro.2015.07.144
![]() |
[6] |
Ingrao C, Lo Giudice A, Mbohwa C, et al. (2014) Life cycle inventory analysis of a precast reinforced concrete shed for goods storage. J Clean Prod 79: 152–167. https://doi.org/10.1016/j.jclepro.2014.05.030 doi: 10.1016/j.jclepro.2014.05.030
![]() |
[7] |
Sherwani AF, Usmani JA, Varun (2010) Life cycle assessment of solar PV based electricity generation systems: A review. Renew Sust Energy Rev 14: 540–544. https://doi.org/10.1016/j.rser.2009.08.003 doi: 10.1016/j.rser.2009.08.003
![]() |
[8] |
Kua HW, Kamath S (2014) An attributional and consequential life cycle assessment of substituting concrete with bricks. J Clean Prod 81: 190–200. https://doi.org/10.1016/j.jclepro.2014.06.006 doi: 10.1016/j.jclepro.2014.06.006
![]() |
[9] |
Dong YH, Ng ST, Kwan AHK, et al. (2015) Substituting local data for overseas life cycle inventories: A case study of concrete products in Hong Kong. J Clean Prod 87: 414–422. https://doi.org/10.1016/j.jclepro.2014.10.005 doi: 10.1016/j.jclepro.2014.10.005
![]() |
[10] |
Di Maria F, Beccaloni E, Bonadonna L, et al. (2020) Minimization of spreading of SARS-CoV-2 via household waste produced by subjects affected by COVID-19 or in quarantine. Sci Total Environ 15: 1–6. https://doi.org/10.1016/j.scitotenv.2020.140803 doi: 10.1016/j.scitotenv.2020.140803
![]() |
[11] |
Miao J, Li J, Wang F, et al. (2022) Characterization and Evaluation of the Leachability of Bottom Ash from a Mobile Emergency Incinerator of COVID-19 Medical Waste: A Case Study in Huoshenshan Hospital, Wuhan, China. J Environ Manage 303: 1–8. https://doi.org/10.1016/j.jenvman.2021.114161 doi: 10.1016/j.jenvman.2021.114161
![]() |
[12] |
Rachmawati S, Syafrudin, Budiyono (2023) Potential to be used for paving blocks. Jurnal Kesehatan Masyarakat 19: 312–318. https://doi.org/10.15294/kemas.v19i2.44392 doi: 10.15294/kemas.v19i2.44392
![]() |
[13] |
Makarichi L, Jutidamrongphan W, Techato KA (2018) The evolution of waste-to-energy incineration: A review. Renew Sustain Energy Rev 91: 812–821. https://doi.org/10.1016/j.rser.2018.04.088 doi: 10.1016/j.rser.2018.04.088
![]() |
[14] | Chen M, Yang H (2016) Current status of medical wastes disinfection and disposal technologies. Chin J Disinfect 33: 171–174. |
[15] |
Praveenkumar S, Sankarasubramanian G (2019) Mechanical and durability properties of bagasse ash-blended high-performance concrete. SN Appl Sci 19: 1–7. https://doi.org/10.1007/s42452-019-1711-x doi: 10.1007/s42452-019-1711-x
![]() |
[16] |
Memon BA, Khanzada GM, Oad M, et al. (2020) Tensile strength of concrete with biomedical waste. World J Eng 6: 81–90. https://doi.org/10.13140/RG.2.2.30128.79367 doi: 10.13140/RG.2.2.30128.79367
![]() |
[17] |
Deepak A, Sharma V, Kumar D (2022) Life cycle assessment of biomedical waste management for reduced environmental impacts. J Clean Prod 349: 131376. https://doi.org/10.1016/j.jclepro.2022.131376 doi: 10.1016/j.jclepro.2022.131376
![]() |
[18] |
Liu Y, Yao D, Xu Z, et al. (2023) Comparative analysis of life cycle water accounting of the Lurgi low-pressure methanol production process with biomass or coal as raw materials. Sci Total Environ 856: 159129. https://doi.org/10.1016/j.scitotenv.2022.159129 doi: 10.1016/j.scitotenv.2022.159129
![]() |
[19] |
Purwanto P, Citra ADP (2019) Recycling and processing of solid waste into products of the cosmetic packaging industry. J Phys 1295. https://doi.org/10.1088/1742-6596/1295/1/012042 doi: 10.1088/1742-6596/1295/1/012042
![]() |
[20] |
Rovira J, Nadal M, Schuhmacher M, et al. (2018) Concentrations of trace elements and PCDD/Fs around a municipal solid waste incinerator in Girona (Catalonia, Spain). Human health risks for the population living in the neighborhood. Sci Total Environ 630: 34–45. https://doi.org/10.1016/j.scitotenv.2018.02.175 doi: 10.1016/j.scitotenv.2018.02.175
![]() |
[21] |
Chen Y, Ding Q, Yang X, et al. (2013) Application countermeasures of non-incineration technologies for medical waste treatment in China. Waste Manag Res 31: 1237–1244. https://doi.org/10.1177/0734242X13507314 doi: 10.1177/0734242X13507314
![]() |
[22] |
Rozumová L, Motyka O, Čabanová K, et al. (2015) Stabilization of Waste Bottom Ash Generated from Hazardous Waste Incinerators. J Environ Chem Eng 3: 1–9. https://doi.org/10.1016/j.jece.2014.11.006 doi: 10.1016/j.jece.2014.11.006
![]() |
[23] |
Rachmawati S, Syarifuddin S, Budiyono B (2023) Quality of paving blocks soaking water made from medical waste incineration ash. IOP Conf Ser Earth Environ Sci 1–7. https://doi.org/10.1088/1755-1315/1268/1/012063 doi: 10.1088/1755-1315/1268/1/012063
![]() |
[24] |
Wu G, Kang H, Zhang X, et al. (2010) A critical review on the bio-removal of hazardous heavy metals from contaminated soils: Issues, progress, eco-environmental concerns and opportunities. J Hazard Mater 174: 1–8. https://doi.org/10.1016/j.jhazmat.2009.09.113 doi: 10.1016/j.jhazmat.2009.09.113
![]() |
[25] |
Pertiwi V, Joko T, Dangiran HL (2017) Evaluasi pengelolaan limbah bahan berbahaya dan beracun (B3) di Rumah Sakit Roemani Muhammadiyah Semarang. Jurnal Kesehatan Masyarakat 5: 420–430. https://doi.org/10.14710/jkm.v5i3.17260 doi: 10.14710/jkm.v5i3.17260
![]() |
[26] |
Demir AT, Moslem S (2024) A novel fuzzy multi-criteria decision-making for enhancing the management of medical waste generated during the coronavirus pandemic. Eng Appl Artif Intel 133: 108465. https://doi.org/10.1016/j.engappai.2024.108465 doi: 10.1016/j.engappai.2024.108465
![]() |
[27] | Intergovernmental Panel on Climate Change (IPCC) (2006), CHAPTER 1, 1–21. |
[28] | Directorate General of Pollution and Environmental Damage Control (2021) Guidelines for the Preparation of a Life Cycle Assessment (LCA) Report Ministry of Environment and Forestry, 1–82. |
[29] |
Zhao L, Zhang FS, Wang K, et al. (2008) Chemical properties of heavy metals in typical hospital waste incinerator ashes in China. Waste Manage 29: 1114–1121. https://doi.org/10.1016/j.wasman.2008.09.003 doi: 10.1016/j.wasman.2008.09.003
![]() |
[30] | LaGrega MD, Buckingham PL, dan Evans JC (2010) Hazardous Waste Management, New York: Waveland Press. |
[31] | Widayatno T, Yuliawati T, Susilo AA, et al. (2017) Adsorpsi logam berat (Pb) dari limbah cair dengan adsorben arang bambu aktif. Jurnal Teknologi Bahan Alam 1: 17–23. |
[32] |
Biswas B, Qi F, Biswas JK, et al. (2018). The fate of chemical pollutants with soil properties and processes in the climate change paradigm—a review. Soil Syst 2: 1–20. https://doi.org/10.3390/soilsystems2030051 doi: 10.3390/soilsystems2030051
![]() |
[33] |
Kochany EL (2018) Effect of climate change on humic substances and associated impacts on the quality of surface water and groundwater: A review. Sci Total Environ 640: 1548–1565. https://doi.org/10.1016/j.scitotenv.2018.05.376 doi: 10.1016/j.scitotenv.2018.05.376
![]() |
[34] |
Khan MAI (2018) Microbial diversity changes with rhizosphere and hydrocarbons in contrasting soils. Ecotox Environ Safe 156: 434–442. https://doi.org/10.1016/j.ecoenv.2018.03.006 doi: 10.1016/j.ecoenv.2018.03.006
![]() |
[35] |
Pechyen C, Ummartyotin S (2017) Development of isotactic polypropylene and stearic acid-modified calcium carbonate composite: A promising material for microwavable packaging. Polym Bull 74: 431–444. https://doi.org/10.1007/s00289-016-1722-3 doi: 10.1007/s00289-016-1722-3
![]() |
[36] |
Sydow M, Chrzanowski Ł, Hauschild MZ, et al. (2020) Influence of metal speciation on soil ecotoxicity impacts in life cycle assessment. J Environ Manage 266. https://doi.org/10.1016/j.jenvman.2020.110611 doi: 10.1016/j.jenvman.2020.110611
![]() |
[37] |
Plouffe G, Ceccile B, Louisse D (2015) Case study: Taking Zinc speciation into account in terrestrial ecotoxicity considerably impacts life cycle assessment results. J Clean Prod 108: 1002–1008. https://doi.org/10.1016/j.jclepro.2015.06.050 doi: 10.1016/j.jclepro.2015.06.050
![]() |
[38] |
Hong YJ, Liao W, Yan ZF, et al. (2020) Progress in the research of the toxicity effect mechanisms of heavy metals on freshwater organisms and their water quality criteria in China. J Chem 2020: 1–12. https://doi.org/10.1155/2020/9010348 doi: 10.1155/2020/9010348
![]() |
[39] |
Fu Q, Weng N, Fujii M, et al. (2018) Temporal variability in Cu speciation, phytotoxicity, and soil microbial activity of Cu-polluted soils as affected by elevated temperature. Chemosphere 194: 285–96. https://doi.org/10.1016/j.chemosphere.2017.11.183 doi: 10.1016/j.chemosphere.2017.11.183
![]() |
[40] |
Noyes PD, Seas CL (2015) Forecasting the impacts of chemical pollution and climate change interactions on the health of wildlife. Curr Zool 61: 669–89. https://doi.org/10.1093/czoolo/61.4.669 doi: 10.1093/czoolo/61.4.669
![]() |
[41] |
Pham N, Babcsányi I, Farsang A (2022) Ecological risk and enrichment of potentially toxic elements in the soil and eroded sediment in an organic vineyard (Tokaj Nagy Hill, Hungary). Environ Geochem Hlth 44: 1893–1909. https://doi.org/10.1007/s10653-021-01076-w doi: 10.1007/s10653-021-01076-w
![]() |
[42] |
Viveros IV, Levasseur A, Bulle C, et al. (2023) Modelling the influence of climate change on characterization factors for copper terrestrial ecotoxicity. J Clean Prod 414. https://doi.org/10.1016/j.jclepro.2023.137601 doi: 10.1016/j.jclepro.2023.137601
![]() |
[43] |
Pelesaraei AN, Mohammadkashi N, Naderloo L, et al. (2022) Principal of environmental life cycle assessment for medical waste during COVID-19 outbreak to support sustainable development goals. Sci Total Environ 25. https://doi.org/10.1016/j.scitotenv.2022.154416 doi: 10.1016/j.scitotenv.2022.154416
![]() |
[44] |
Kouassi HK, Murayama T, Ota M (2022) Life cycle analysis and cost-benefit assessment of the waste collection system in Anyama, Cote d'Ivoire. Sustainability 14. https://doi.org/10.3390/su142013062 doi: 10.3390/su142013062
![]() |
[45] |
Xing YF, Xu YH, Shi MH, et al. (2016) The impact of PM2.5 on the human respiratory system. J Thorac Dis 8. https://doi.org/10.3978/j.issn.2072-1439.2016.01.19 doi: 10.3978/j.issn.2072-1439.2016.01.19
![]() |
[46] |
Suryadi I, Nugraha AP, Fitriani, N, et al. (2022) The determinant of lung function disorders of the textile industry spinning section. Jurnal Kesehatan Masyarakat 17: 475–482. https://doi.org/10.15294/kemas.v17i4.25069 doi: 10.15294/kemas.v17i4.25069
![]() |
[47] |
Gnonsoro UP, Ake Assi YED, Sangare NS, et al. (2022) Health risk assessment of heavy metals (Pb, Cd, Hg) in hydroalcoholic gels of Abidjan, Côte d'Ivoire. Biol Trace Elem Res 200. https://doi.org/10.1007/s12011-021-02822-y doi: 10.1007/s12011-021-02822-y
![]() |
[48] |
Oase R, Nukpezah D, Darko DA, et al. (2023) Accumulation of heavy metals and human health risk assessment of vegetable consumption from a farm within the Korle lagoon catchment. Heliyon 9: 1–16. https://doi.org/10.1016/j.heliyon.2023.e16005 doi: 10.1016/j.heliyon.2023.e16005
![]() |
[49] |
Emmanuel UC, Chukwudi MI, Monday SS, et al. (2022) Human health risk assessment of heavy metals in drinking water sources in three senatorial districts of Anambra State, Nigeria. Toxicol Rep 9: 869–875. https://doi.org/10.1016/j.toxrep.2022.04.011 doi: 10.1016/j.toxrep.2022.04.011
![]() |
[50] |
Jaishankar M, Tseten T, Anbalagan N, et al. (2014) Toxicity, mechanism and health effects of some heavy metals. Interdiscip Toxicol 7: 60–72. https://doi.org/10.2478/intox-2014-0009 doi: 10.2478/intox-2014-0009
![]() |
[51] |
Gasia J, Fabiani C, Chàfer M, et al. (2021) Life cycle assessment and life cycle costing of an innovative component for refrigeration units. J Clean Prod 295. https://doi.org/10.1016/j.jclepro.2021.126442 doi: 10.1016/j.jclepro.2021.126442
![]() |
[52] |
Manda BMK, Bosch H, Karanam S, et al. (2016) Value creation with life cycle assessment: An approach to contextualize the application of life cycle assessment in chemical companies to create sustainable value. J Clean Prod 126: 337–351. https://doi.org/10.1016/j.jclepro.2016.03.020 doi: 10.1016/j.jclepro.2016.03.020
![]() |
[53] |
Indriyani D, Darundiati Y, Dewanti N (2017) Analisis risiko kesehatan lingkungan pajanan debu kayu pada pekerja Di industri mebel Cv. Citra Jepara Kabupaten Semarang. Jurnal Kesehatan Masyarakat 5: 571–580. https://doi.org/10.14710/jkm.v5i5.19179 doi: 10.14710/jkm.v5i5.19179
![]() |
[54] |
Peruzzini M, Michele G, Eugenia M (2013) Product-service sustainability assessment in virtual manufacturing enterprises. IFIP Adv Inform Commun Technol 408: 13–21. https://doi.org/10.1007/978-3-642-40543-3_2 doi: 10.1007/978-3-642-40543-3_2
![]() |
[55] |
Ji S, Lee B, Yi MY (2021) Building life-span prediction for life cycle assessment and life cycle cost using machine learning: A big data approach. Build Environ 205: 108267. https://doi.org/10.1016/j.buildenv.2021.108267 doi: 10.1016/j.buildenv.2021.108267
![]() |
[56] |
Adelfio L, Giallanza A, La Scalia G, et al. (2023) Life cycle assessment of a new industrial process for sustainable construction materials. Ecol Indic 148: 110042. https://doi.org/10.1016/j.ecolind.2023.110042 doi: 10.1016/j.ecolind.2023.110042
![]() |
[57] |
Patel P, Schwartz D, Wang X, et al. (2022) Technoeconomic and life-cycle assessment for electrocatalytic production of Furandicarboxylic Acid. ACS Sustain Chem Eng 10: 4206–4217. https://doi.org/10.1021/acssuschemeng.1c08602 doi: 10.1021/acssuschemeng.1c08602
![]() |
1. | Fulian Yin, Xinyi Tang, Tongyu Liang, Yanjing Huang, Jianhong Wu, External intervention model with direct and indirect propagation behaviors on social media platforms, 2022, 19, 1551-0018, 11380, 10.3934/mbe.2022530 | |
2. | Rui Shi, Xiaoyong Hou, Chang Liu, Praveen Kumar Reddy Maddikunta, Model of Negative Emotional Information Communication among Netizens under Corporate Negative Events, 2022, 2022, 1563-5147, 1, 10.1155/2022/3527980 | |
3. |
Rongrong Qiao, Yuhan Hu,
Dynamic analysis of a SI1I2ADSI1I2AD information dissemination model considering the word of mouth,
2023,
111,
0924-090X,
22763,
10.1007/s11071-023-09021-5
|
|
4. | Youjun Lu, Sen Wu, Jiayin Wei, Yanjun Liang, Li Deng, 2024, Dynamic Analysis and Optimal Control of SIWCR Rumor Propagation Model with Rumor Countering and Whistleblowing Mechanism, 979-8-3503-5359-4, 911, 10.1109/EEI63073.2024.10696185 | |
5. | 永缜 邹, SIMR Rumor Propagation Model of Considering Reporting Mechanism, 2024, 13, 2324-8696, 5131, 10.12677/mos.2024.135464 | |
6. | Sen Wu, Youjun Lu, Congli Chen, Shichen Liu, 2024, Dynamic analysis of SICDR rumor propagation model with time delay and psychological factors, 979-8-3503-7510-7, 746, 10.1109/CISAT62382.2024.10695275 | |
7. | Ru Wang, Zhengyou Xia, 2024, ISPR: A rumor propagation model based on dual transmission modes of friends and recommendation, 979-8-3503-6360-9, 551, 10.1109/ICIBA62489.2024.10868431 |
Parameters | Description | Range(/days) |
b | Join networks rate | 0.04−0.4 |
d | Exit networks rate | 0.04−0.4 |
λ1 | Rate of a susceptible node contacts the infected node | 0.009−0.9 |
λ2 | Rate of a susceptible node contacts the Marketing node | 0.009−0.9 |
δ | Rate from infected node to removed node | 0.009−0.9 |
γ | Rate from infected node to Marketing node | 0.009−0.9 |
β | Rate from Marketing node to removed node | 0.009−0.9 |
Parameters | b | λ1 | λ2 | d | δ | γ | β | R0 |
Data | 0.4 | 0.1 | 0.1 | 0.04 | 0.1 | 0.1 | 0.1 | 0.41 |
0.4 | 0.3 | 0.3 | 0.04 | 0.1 | 0.1 | 0.1 | 1.43 | |
0.4 | 0.12 | 0.1 | 0.04 | 0.1 | 0.1 | 0.1 | 0.89 |
Parameters | Value1 | Value2 | Value3 | Value4 | Value5 |
c3 | c3=10 | c3=20 | c3=30 | c3=40 | c3=50 |
c4 | c4=10 | c4=20 | c4=30 | c4=40 | c4=50 |
Parameters | Description | Range(/days) |
b | Join networks rate | 0.04−0.4 |
d | Exit networks rate | 0.04−0.4 |
λ1 | Rate of a susceptible node contacts the infected node | 0.009−0.9 |
λ2 | Rate of a susceptible node contacts the Marketing node | 0.009−0.9 |
δ | Rate from infected node to removed node | 0.009−0.9 |
γ | Rate from infected node to Marketing node | 0.009−0.9 |
β | Rate from Marketing node to removed node | 0.009−0.9 |
Parameters | b | λ1 | λ2 | d | δ | γ | β | R0 |
Data | 0.4 | 0.1 | 0.1 | 0.04 | 0.1 | 0.1 | 0.1 | 0.41 |
0.4 | 0.3 | 0.3 | 0.04 | 0.1 | 0.1 | 0.1 | 1.43 | |
0.4 | 0.12 | 0.1 | 0.04 | 0.1 | 0.1 | 0.1 | 0.89 |
Parameters | Value1 | Value2 | Value3 | Value4 | Value5 |
c3 | c3=10 | c3=20 | c3=30 | c3=40 | c3=50 |
c4 | c4=10 | c4=20 | c4=30 | c4=40 | c4=50 |