
The generation of large volumes of agricultural and agroindustrial waste in the state of Tabasco represents a significant waste management challenge. We aimed to determine the bioenergy potential of five types of biomasses: Banana rachis, coconut shell, cocoa pod husk, sugarcane bagasse, and palm kernel shell, generated in agricultural and agroindustrial processes. This research involved characterizing and evaluating the energy quality of these biomasses by determining their calorific values and assessing their viability as fuel alternative sources. Additionally, we explored these biomasses' calorific value potential to reduce the inadequate disposal of wastes, reduce environmental impact, and provide alternative uses for these materials, which are typically discarded or have limited added value in the southeast region. The yield of waste generation per amount of production was estimated, with cocoa pod husk biomass and sugarcane bagasse, banana rachis, coconut shell, and palm kernel shell generating 0.685, 0.283, 0.16, 0.135, and 0.0595 kg of biomass per kg of crop, respectively. The bioenergy potential was evaluated through direct measurements using a calorimeter bomb, and indirect measurements using stoichiometric calculations. Four stoichiometric methods based on predictive equations were employed to determine the energy content of the biomasses from their elemental composition (Dulong, Friedl, Channiwala, Boie). The biomasses with the highest calorific values were coconut shell and cocoa pod husk, with values of 16.47 ± 0.24 and 16.02 ± 1.54 MJ/kg, respectively. Moreover, banana rachis had the lowest calorific value at 13.68 ± 3.22 MJ/kg. The calorific values of the sugarcane bagasse and palm kernel shell were 13.91 ± 0.98 and 15.29 ± 1.02, respectively. The factorial experimental design and statistical analysis revealed trends and magnitudes in the evaluation of energy determination methods and types of waste. The predictive equation of Dulong showed the highest similarity to the experimental values, especially for coconut shell (16.02 ± 0.08 MJ/kg). The metal content in biomasses such as palm kernel shell and coconut shell were below the limits established in ISO 17225:2014. Finally, our results indicated that coconut shell has superior characteristics for potential use as an alternative fuel, whereas banana rachis requires exploring alternative utilization options.
Citation: Nathaly A. Díaz Molina, José A. Sosa Olivier, José R. Laines Canepa, Rudy Solis Silvan, Donato A. Figueiras Jaramillo. Evaluation of the bioenergy potential of agricultural and agroindustrial waste generated in southeastern Mexico[J]. AIMS Energy, 2024, 12(5): 984-1009. doi: 10.3934/energy.2024046
[1] | Akhil Kumar Srivastav, Pankaj Kumar Tiwari, Prashant K Srivastava, Mini Ghosh, Yun Kang . A mathematical model for the impacts of face mask, hospitalization and quarantine on the dynamics of COVID-19 in India: deterministic vs. stochastic. Mathematical Biosciences and Engineering, 2021, 18(1): 182-213. doi: 10.3934/mbe.2021010 |
[2] | Rahat Zarin, Usa Wannasingha Humphries, Amir Khan, Aeshah A. Raezah . Computational modeling of fractional COVID-19 model by Haar wavelet collocation Methods with real data. Mathematical Biosciences and Engineering, 2023, 20(6): 11281-11312. doi: 10.3934/mbe.2023500 |
[3] | Fang Wang, Lianying Cao, Xiaoji Song . Mathematical modeling of mutated COVID-19 transmission with quarantine, isolation and vaccination. Mathematical Biosciences and Engineering, 2022, 19(8): 8035-8056. doi: 10.3934/mbe.2022376 |
[4] | Yangyang Yu, Yuan Liu, Shi Zhao, Daihai He . A simple model to estimate the transmissibility of the Beta, Delta, and Omicron variants of SARS-COV-2 in South Africa. Mathematical Biosciences and Engineering, 2022, 19(10): 10361-10373. doi: 10.3934/mbe.2022485 |
[5] | Matthew Hayden, Bryce Morrow, Wesley Yang, Jin Wang . Quantifying the role of airborne transmission in the spread of COVID-19. Mathematical Biosciences and Engineering, 2023, 20(1): 587-612. doi: 10.3934/mbe.2023027 |
[6] | Qinghua Liu, Siyu Yuan, Xinsheng Wang . A SEIARQ model combine with Logistic to predict COVID-19 within small-world networks. Mathematical Biosciences and Engineering, 2023, 20(2): 4006-4017. doi: 10.3934/mbe.2023187 |
[7] | A. M. Elaiw, Raghad S. Alsulami, A. D. Hobiny . Global dynamics of IAV/SARS-CoV-2 coinfection model with eclipse phase and antibody immunity. Mathematical Biosciences and Engineering, 2023, 20(2): 3873-3917. doi: 10.3934/mbe.2023182 |
[8] | Zita Abreu, Guillaume Cantin, Cristiana J. Silva . Analysis of a COVID-19 compartmental model: a mathematical and computational approach. Mathematical Biosciences and Engineering, 2021, 18(6): 7979-7998. doi: 10.3934/mbe.2021396 |
[9] | Hamdy M. Youssef, Najat A. Alghamdi, Magdy A. Ezzat, Alaa A. El-Bary, Ahmed M. Shawky . A new dynamical modeling SEIR with global analysis applied to the real data of spreading COVID-19 in Saudi Arabia. Mathematical Biosciences and Engineering, 2020, 17(6): 7018-7044. doi: 10.3934/mbe.2020362 |
[10] | Sarafa A. Iyaniwura, Rabiu Musa, Jude D. Kong . A generalized distributed delay model of COVID-19: An endemic model with immunity waning. Mathematical Biosciences and Engineering, 2023, 20(3): 5379-5412. doi: 10.3934/mbe.2023249 |
The generation of large volumes of agricultural and agroindustrial waste in the state of Tabasco represents a significant waste management challenge. We aimed to determine the bioenergy potential of five types of biomasses: Banana rachis, coconut shell, cocoa pod husk, sugarcane bagasse, and palm kernel shell, generated in agricultural and agroindustrial processes. This research involved characterizing and evaluating the energy quality of these biomasses by determining their calorific values and assessing their viability as fuel alternative sources. Additionally, we explored these biomasses' calorific value potential to reduce the inadequate disposal of wastes, reduce environmental impact, and provide alternative uses for these materials, which are typically discarded or have limited added value in the southeast region. The yield of waste generation per amount of production was estimated, with cocoa pod husk biomass and sugarcane bagasse, banana rachis, coconut shell, and palm kernel shell generating 0.685, 0.283, 0.16, 0.135, and 0.0595 kg of biomass per kg of crop, respectively. The bioenergy potential was evaluated through direct measurements using a calorimeter bomb, and indirect measurements using stoichiometric calculations. Four stoichiometric methods based on predictive equations were employed to determine the energy content of the biomasses from their elemental composition (Dulong, Friedl, Channiwala, Boie). The biomasses with the highest calorific values were coconut shell and cocoa pod husk, with values of 16.47 ± 0.24 and 16.02 ± 1.54 MJ/kg, respectively. Moreover, banana rachis had the lowest calorific value at 13.68 ± 3.22 MJ/kg. The calorific values of the sugarcane bagasse and palm kernel shell were 13.91 ± 0.98 and 15.29 ± 1.02, respectively. The factorial experimental design and statistical analysis revealed trends and magnitudes in the evaluation of energy determination methods and types of waste. The predictive equation of Dulong showed the highest similarity to the experimental values, especially for coconut shell (16.02 ± 0.08 MJ/kg). The metal content in biomasses such as palm kernel shell and coconut shell were below the limits established in ISO 17225:2014. Finally, our results indicated that coconut shell has superior characteristics for potential use as an alternative fuel, whereas banana rachis requires exploring alternative utilization options.
In December 2019, a new coronavirus infection, called corona virus disease (COVID-19) by the World Health Organization, was found in Wuhan, Hubei, China [1] and soon spread outside of China. Although two rather effective antivirals have been approved by the FDA and monoclonal antibody treatment has been in use for some time and is extremely effective, this viral contamination continues to spread rapidly across globe [2,3]. In March 2020 [1,2,4,5,6,7] the emergency arrived at the pandemic stage as the overall effect was accelerated. The number of cases and deaths with COVID-19 have reached 506 million and 6.2 million worldwide, respectively, as of 21 April 2022. Currently, only vaccines can attenuate the COVID-19 pandemic. Several vaccines, such as Pffizer-BioNtech, Moderna, AstraZeneca, etc., had been developed and approved via Emergency Use Authorization (EUA) worldwide [8]. On December 2, 2020, UK became the first country to approve the Pffizer-BioNtech COVID-19 vaccine, and several vaccines were approved under EUA by different countries in early 2021 [8].
Mathematical modeling is a useful technique to understand the dynamics of infectious diseases. Mathematical models of evolution of COVID-19 pandemic have been proposed to predict the spread of Covid-19 pandemic such as the peaks, total numbers of infections, and the ending time of COVID-19 epidemic in different countries [9,10,11,12,13]. Balasubramaniam et al. [9] used a stochastic epidemiological model (SEM) along with a machine learning technique to compute the infection rates in populations. Yadav [10] employed the well known SIR model, consisting of deterministic ODEs, to predict the peak and ending time of COVID-19 pandemic in India. Hassen et al. [11] also employed the SIR model to predict the evolution of COVID-19 cases in the Maghreb Central countries. Bhardwaj [12] used the natural exponential function with a time dependent growth rate to estimate the peak dates and the number of total infections in several countries. Jamshidi et al. [13] proposed a computational algorithm to predict the number and locations of new epicenters.
Santosh [14] pointed out that complex models were needed to study the spread of COVID-19 which depends on many factors such as hospital capacity, testing rate, demographics, population density, poverty, etc. Furthermore, as COVID-19 continues to spread, mathematical models studying pharmaceutical or non-pharmaceutical factors, which can help with reducing effect of the spread and devising the control strategies for fighting COVID-19 pandemic, are as important as the models attempting to predict the future evolution of COVID-19 pandemic [15,16,17]. Dawoudi [18] proposed a computational algorithm to study coronavirus DNA sequences for discovery of new drugs. Based on the SIR model, Wanjau, Ndam, and Asempapa et al. [19,20,21] proposed mathematical models with more sub-populations such as susceptible, exposed, infectious, isolated, and recovered individuals. Asempapa et al. [21] further divided some of the above sub-populations into high-risk and low-risk sub-populations. These complex models contain many coupled ODEs with many parameters. However, daily information of these sub-populations, such as the number of isolated individuals, high-risk and low-risk sub-populations, or the number of exposed individuals, is not available in many countries/regions. Estimation of parameters and the plausibility of the model may be hindered by lack of information in many countries/regions.
Rojas-Vallejos [22] suggested that mathematical models should not be simpler than necessary, but also not more complex than needed to provide helpful insight. Recently, Alvarez et al. [23] have proposed a simple yet novel mathematical model of COVID-19 epidemic to study the effect of social distancing and detection of positive cases in urban areas. The mathematical model considered two populations, population of the cumulative confirmed cases and population of recovery/death. Daily reported data of the number of cumulative cases and the number of recovered or died individuals are available for most countries, making the model plausible for their study. The computational results are interesting. However, the mathematical model has not been correctly formulated. Motivated by the model in [23], the purpose of this paper is to build a mathematical model to study the impacts of pharmaceutical factors (vaccination, detection efficiency, and health care capacity) and non-pharmaceutical factors (lock-down and social distancing) on COVID-19 pandemic.
To provide details, the mathematical model is presented in Section 2. Parameter fitting and numerical simulation are performed in Section 3. Finally, a brief conclusion is made in Section 4.
In this paper, a two-stage model is to be constructed. The first-stage model contains two state variables, X and R, and the impacts of lock-down, social distancing, and testing intensity. The variable X represents the cumulative COVID-19 cases. The variable R represents the number of individuals who have been infected but removed from the population and thus do not contribute to the propagation of COVID-19. The construction of first-stage model is motivated by the model proposed by Alvarez et al. [23]. However, the mathematical model proposed by Alvarez et al. [23] was not correctly formulated nor consistent with the implementation of the solutions in Excel spreadsheets (Supplemental File F1 in [23]). The model given in [23] is as follows:
dX/dt=μ0(1−σ)(X−R)(1−X/P0), | (1) |
dR/dt=α∫t=t−delay_qt=0dX/dt+(1−α)∫t=t−delay_rt=0dX/dt, | (2) |
where P0 is the population of the region, which may be a country, state, or urban area, and μ0 is the intrinsic infection rate. Both σ and α were modeled by nondecreasing functions of t [23]. The function σ represents the effect of social distancing. The infection rate is reduced by σ due to social distancing effect. In [23], Eq (2) describes the rate at which infected patients are retrieved from the general population. It has been assumed in [23] that a fraction (α) of infected individuals is effectively removed from the general population with a delay of delay_q days from the onset of virus to positive diagnosis. These individuals do not contribute to the propagation of COVID-19. Another fraction (1−α) of infected individuals is not effectively removed from the general population until their recovery or death after delay_r days. It was stated in [23] that the second term was formulated by the integral of all infected individuals recovered or died from the onset of the epidemic episode in the region with a delay of delay_r days. Apparently, Eq (2) and the above statement do not describe the rate of change of R. Furthermore, Eqs (1) and (2) were converted into equations of differences in [23]:
ΔX=(μ0(1−σ)(X−R)(1−X/P0))Δt, | (3) |
ΔR={α∫t=t−delay_qt=0dX/dt+(1−α)∫t=t−delay_rt=0dX/dt}Δt. | (4) |
Using the fact that X(t) is nondecreasing, Eq (4) is equivalent to
ΔR=(αX(t−delay_q)+(1−α)X(t−delay_r))Δt≥(αX(t−delay_r)+(1−α)X(t−delay_r))Δt=X(t−delay_r)Δt | (5) |
if delay_r≥delay_q. The left hand side of Eq (5) represents the increment of the population R over a time span [t,t+Δt] while the right hand side uses the cumulative infected cases over the whole time span with a delay of a few days [0,t−delay_r], multiplied by Δt. This shows that Eq (2) was not correctly formulated in [23]. Nevertheless, the elements and terms considered in the mathematical model in [23] are novel and can be easily applied to construct a modified model to study the evolution of the desease.
Following the construction concept of Eqs (1) and (2) with modification and incorporation of the lock-down effect, the first-stage model proposed in this paper is given by a deterministic system of ODEs with time delays:
dX/dt=μ1(1−σ(t))(X(t)−R(t)+K(t))(1−X(t)/P0), | (6) |
dR/dt=α(t−τ1)dX(t−τ1)/dt+(1−α(t−τ2))dX(t−τ2)/dt, | (7) |
where μ1 is the intrinsic infection rate and K is the number of infected individuals who reside outside of the region and travel to the region. The functions σ(t) and α(t) represent the same effects as in Eqs (1) and (2). The notations delay_q and delay_r are replaced by shorter notations τ1 and τ2, respectively.
Equation (6) describes the rate of change in cumulative COVID-19 cases. Assume that it is of similar status of COVID-19 outbreak in neighboring regions. This assumption leads to K(t)=β1(X(t)−R(t)). Thus X(t)−R(t)+K(t)=(1+β1)(X(t)−R(t)) where β1 represents the unlock-down factor. At the stage of implementation of the lock-down orders, β1=0.
Social distancing involves the change of behavior [24] and it is difficult for individuals to stick to the social distancing guidelines [25]. Implementation of social distancing policy can become effective after people have learned the importance of social distancing [26]. This behavioral change can be related to the S-shaped learning curve, measuring an individual who is new to a task [27]. The bottom of the S curve indicates slow progressing initially as it requires more time to master new skills. The S curve also describes the strongest resistance to adoption of a change [28]. After the learner becomes proficient in the skills, it requires less time to complete the task. The curve eventually levels off, indicating a plateau or new challenges [27]. Alvarez et al. [23] reported that using a linear ramp between two values over a time period for σ could reproduce the actual progression of COVID-19. The behavioral change of social distancing described above and the formulation of σ given in [23] suggest that a sigmoid function, which is regularly used in mathematical modeling, is a good fit for σ.
The effect of social distancing σ(t), shown in Figure 1(a), is modeled as follows:
σ(t)=L11+e−k1(t−t1),0≤L1≤1, | (8) |
where k1 represents the steepness of the curve, t1 is the time point at which the learner reaches the maximum growth rate, and L1 represents the maximum performance measure. When a learner develops new skills required for a new task through previous experience, the bottom of the curve might not start with zero performance. In this case, Eq (8) is modified as follows:
σ(t)=M1+L11+e−k1(t−t1),0≤M1+L1≤1, | (9) |
where M1 is the previous-experience factor. Note that Eq (9) becomes Eq (8) when M1=0. Thus Eq (9) is used in this paper as it represents a more general case than Eq (8). Figures 1(a), (b) give examples of S-shaped learning curves without and with prior experience, respectively.
Equation (7) describes the rate of change of the removed population R. At the initial stage of COVID-19 outbreak, most countries faced a shortage of testing kits [29]. As a result, infected individuals were tested only after the onset of symptoms. The function α has the same meaning as in Eq (2). Recall that τ1 is the delay from the onset of symptoms to positive diagnosis. The delay takes about 2–5 days [23]. We use τ1=4 in this paper. It has been reported that the median time to recovery is 13.5 days [30]. We assume that these individuals are removed from the general population after their recovery or death with a delay of τ2=14 days. From Eq (7), 0≤α(t)≤1, X(t)=0 for t<0, and dX(t)/dt=0 for t<0, we have
R(t)=∫t0(dR(ξ)/dξ)dξ=∫t0[α(ξ−τ1)dX(ξ−τ1)/dξ+(1−α(ξ−τ2))dX(ξ−τ2)/dξ]dξ=∫t−τ10α(ξ)(dX(ξ)/dξ)dξ+∫t−τ20(1−α(ξ))(dX(ξ)/dξ)dξ=∫t−τ20α(ξ)(dX(ξ)/dξ)dξ+∫t−τ1t−τ2α(ξ)(dX(ξ)/dξ)dξ+∫t−τ20(1−α(ξ))(dX(ξ)/dξ)dξ=X(t−τ2)+∫t−τ1t−τ2α(ξ)(dX(ξ)/dξ)dξ≤X(t−τ1)≤X(t). | (10) |
Equation (10) shows that R(t) includes the cumulative infected individuals over the time span [0,t−τ2] and the removed individuals over [t−τ2,t−τ1]. The expression X(t)−R(t) represents the infected individuals who are capable of infecting susceptible individuals in the region. Equation (10) also shows that X(t)−R(t) includes the individuals contracted the disease over [t−τ1,t] and the individuals contracted the disease over [t−τ2,t−τ1] who have not been removed from the general population. Since X(t) is a differentiable, nondecreasing function and bounded above by P0, the limit limt→∞X(t) exists and limt→∞(dX/dt)=0. Thus, Eq (10) implies limt→∞R(t)=limt→∞X(t); all the infected individuals will be included in the population R eventually.
The outbreak of COVID-19 is an unprecedented public health crisis affecting all countries worldwide. Testing and treatment strategies are new tasks for health-care workers. Thus, it is reasonable to expect an improvement in effective diagnosis with time. Here, the removal factor α(t) is also modeled by an S-curve,
α(t)=M2+L21+e−k2(t−t2),0≤M2+L2≤1. | (11) |
The first-stage model becomes
dX/dt=μ1(1+β1)(1−σ(t))(X(t)−R(t))(1−X(t)/P0), | (12) |
dR/dt=α(t−τ1)dX(t−τ1)/dt+(1−α(t−τ2))dX(t−τ2)/dt, | (13) |
where σ(t) and α(t) are given as in Eqs (9) and (11), respectively. Figure 2 gives a schematic of Eqs (12) and (13).
COVID-19 vaccines had been approved by many countries in early 2021 [8]. In the second-stage model, four state variables X1, X2, Y1, and R are considered. The variables X1 and X2 represent the cumulative COVID-19 cases in unvaccinated and vaccinated sub-populations, respectively. The variable Y1 represents the unvaccinated uninfected sub-population. The practical meaning of R is the same as in Eqs (12) and (13). Many countries have experienced the second wave of COVID-19, and the second wave has caused a dramatic increase in COVID-19 cases which exceed health care capacity [31,32,33,34]. The second-stage model further considers the factor of health care capacity and is given as follows:
dX1/dt=μ1(1+β1)(1−σ(t))(X1(t)+X2(t)−R(t))(Y1(t)/P0)−γ(t)X1(t)h1+X1(t), | (14) |
dX2/dt=μ2(1+β1)(1−σ(t))(X1(t)+X2(t)−R(t))(Y2(t)/P0)+γ(t)X1(t)h1+X1(t), | (15) |
dY1/dt=−μ1(1+β1)(1−σ(t))(X1(t)+X2(t)−R(t))(Y1(t)/P0)−γ(t)Y1(t)h1+Y1(t) | (16) |
dR/dt=(1−H(t))(α(t−τ1))dX(t−τ1)/dt+(1−α(t−τ2))dX(t−τ2)/dt+H(t−(τ2−τ1))α(t−τ2)dX(t−τ2)/dt | (17) |
where X(t)=X1(t)+X2(t) represents cumulative infected cases and Y2(t)=P0−X1(t)−X2(t)−Y1(t) represents the vaccinated uninfected sub-population. The function
γ(t)=β2(M3+L31+e−k3(t−t3)) | (18) |
is a factor related to vaccination rate with the maximum vaccination rate β2. Note that γ(t) is modeled by an S-curve. Vaccination plays an important role in controlling the spread of the disease. Due to rapid development of COVID-19 vaccines, doubts and worries of vaccination safety lead to hesitancy or refusal initially [35]. Many people preferred to wait until others have taken the vaccine [36]. However, vaccination willingness changed over time, where vaccination refusal was reported from >70% in March 2020 to <50% in October 2020 [36]. Another research study reported an increase in vaccine acceptance during a study period from February to September 2020 [37]. The behavioral response to vaccination acceptance can be related to adoption and diffusion of technology which is often described by the S-curve [28]. The vaccination rate depends on the size of unvaccinated sub-populations X1 and Y1. A factor X1/(h1+X1) is incorporated into the term of vaccination rate. This factor describes the situations that the rate of change of vaccination is determined by vaccination willingness, γ(t), if a considerable population is unvaccinated (X1>>h1 or X1/(h1+X1)≈1). The rate of change of vaccination decreases and becomes small as majority of the population are vaccinated (X1<<h1 or X1/(h1+X1)≈0). Furthermore, the factor X1/(h1+X1) can also prevent the solutions from becoming negative.
Let τ3 be the average number of days that a COVID-19 patient is hospitalized for. The function H(t)∈[0,1] represents the utilization factor of health care resources. Due to limited health care resources available for COVID-19 patients, we assume that the portion of COVID-19 patients who can be admitted to hospital depends on H(t). When H(t)<<1 (1−H(t)≈1), ample health care resources are available for COVID-19 patients; the patients with positive test results can be admitted to hospital. In contrast, H(t)≈1 represents the situation when available health care resources have been fully employed. In this case, few or none COVID-19 patients can be admitted to hospital. Let Hτ3(t) represent the number of patients being hospitalized. For simplicity, we shall assume that these patients do not contribute to the propagation of the disease. Since (1−H(t))(α(t−τ1))(dX(t−τ1)/dt) represents the rate of change of these individuals,
Hτ3(t)=∫tt−τ3(1−H(ξ)α(ξ−τ1)(dX(ξ−τ1)/dξ)dξ | (19) |
represents the number of patients being hospitalized at t. We model the utilization factor of health care resources using H(t)=Hτ3(t)h2+Hτ3(t) which depends on the number of patients being hospitalized at t. Thus
Hτ3(t)=∫tt−τ3(1−Hτ3(ξ)h2+Hτ3(ξ))α(ξ−τ1)(dX(ξ−τ1)/dξ)dξ. | (20) |
Since τ2 represents the average time period it requires to recover from the disease, we use τ3=τ2 in this paper.
The parameter values in the mathematical model are fitted using the daily reported data of Tamil Nadu of India. The dataset is taken from [38], a repository supported by Center for Systems Science and Engineering (CSSE) at Coronavirus Resource Center of John Hopkins University (JHU) through GitHub. The features, in daily reported data, related to Eqs (12) and (13) are "Confirmed", "Deaths", "Recovered" and "Active", where Confirmed equals the sum of Deaths, Recovered, and Active. The data for the feature Confirmed is used to fit the cumulative cases X(t). Figure 4 shows that the sum of Deaths and Recovered in daily reported data follows the same progression as Confirmed with a delay of roughly 10 days. Note that the mathematical analysis (Eq (10)) agrees with this property. The sum of Deaths and Recovered is used to fit R(t). In the second-stage model, X2(t)+Y2(t) represents the number of individuals who received at least one dose of vaccine in the data. The vaccination data for Tamil Nadu is not available. To work around this, the daily vaccination data for India is downloaded from [39]. Consider populations of 7.8×107 persons and 1.4×109 persons in Tamil Nadu and India, respectively [40]. The India Express reported that over 76% of eligible population have received at least one dose of COVID vaccine on 25 November 2021. Daily vaccination data for India shows that 776, 515, 911 individuals (about 55.5% of the population) in India had received at least one dose of COVID vaccine by 25 November 2021. Let f=(0.76⋅7.8×107)/(0.555⋅1.4×109)=0.076, and let N(t) be the number of individuals in India that had been vaccinated. Then, fN(t) is used to estimate the number of individuals in Tamil Nadu who had been vaccinated. The lock-down dates are taken from lock-down orders posted in https://www.tn.gov.in/go_view/dept/26 by Tamil Nadu government. Two periods of lock-down orders took place in 2020, 8 April 2020 through 15 September 2020 and 28 October 2020 through 30 November 2020.
The daily reported data of COVID-19 cases from 6 March 2020 to 28 February 2021 (359 days) along with the lock-down information is used to estimate the parameter values for the first-stage model, Eqs (12) and (13). Let T=359 be the number of days from 6 March 2020 to 28 February 2021. Let ˜Xi and ˜Ri, i=0,1,⋯,T be Confirmed and the sum of Deaths and Recovered from the daily reported data, respectively. Equations (12) and (13) are integrated over [0,T] using Euler's method with Δt=0.01 day. Ten parameters listed in Table 1 are to be fitted using 360 data points. The parameter values are determined by solving the minimization problem min√∑Ti=0(X(i)−˜Xi)2+√∑Ti=0(R(i)−˜Ri)2. The Nelder Mead simplex method [41] is employed to solve the least square error of the fitting. The Nelder Mead simplex method is an iterative method. The Matlab build-in functions fminsearch can also fulfill the tasks. In the fitting process, we set 200 for the maximum number of iteration. The fitting process is executed a few times until the evolution of X(t) and R(t) agrees with the daily reported data. Figure 5 shows the fitting results with √∑Ti=0(X(i)−˜Xi)2+√∑Ti=0(R(i)−˜Ri)2=1.3×105 (persons) and with the parameter values determined and summarized in Table 1. The parameter fitting result shows that the intrinsic infection rate in Tamil Nadu is 0.2053 day−1 which is between the growth rates initially and after measures for the regions studied by Alvarez et al. [23].
Parameter | Brief Description | Value | Units | Reference |
P0 | 7.8×107 | Population of Tamil Nadu State, India | individual | [40] |
μ1 | 0.2053 | Intrinsic infection rate | Day−1 | fitted with data [38] |
β1 | 0.0484 | Unlock-down factor | Dimensionless | fitted with data [38] |
M1 | 0 | Pre-experience in σ | Dimensionless | fitted with data [38] |
t1 | 29 | Time parameter in σ | Day | fitted with data [38] |
k1 | 0.1935 | Steepness of σ | Day−1 | fitted with data [38] |
L1 | 0.4533 | Maximum increase in σ | Dimensionless | fitted with data [38] |
M2 | 0 | Pre-experience in α | Dimensionless | fitted with data [38] |
t2 | 105 | Time parameter in α | Day | fitted with data [38] |
k2 | 0.0248 | Steepness of α | Day−1 | fitted with data [38] |
L2 | 0.6405 | Maximum increase in α | Dimensionless | fitted with data |
τ1 | 4 | Delay in detection | Day | [23] |
τ2 | 14 | Recovery time | Day | [30] |
Parameter estimation of dynamical systems is often associated with parameter uncertainty due to data errors and computational errors in parameter fitting processes. The analysis of parameter uncertainty is important as it determines how reliable the model outputs are. Bootstrap method is a popular technique for assessment of parameter uncertainty due to its simplicity. A bootstrapping sampling was proposed for assessment of parameter uncertainty of a mathematical model of population dynamics [42]. This method can be directly applied to the analysis of parameter uncertainty in this paper. A random sample with replacement of size T+1 is drawn over {0,1,⋯,T}. After discarding the repeated numbers, we have 0≤t1<t2<⋯<tI≤T. The parameter values are estimated using the datasets ˜Xti and ˜Rti, i=1,⋯,I with the parameter values listed in Table 1 as the initial guess. The Nelder Mead simplex method is performed for 500 iterations to obtain the bootsrtap estimates of parameter values. The above parameter estimation process with resampling is repeated 1000 times. The bootstrap statistics are shown in Table 2. As about 95% of the data values of a distribution usually fall within 2 standard deviations from the mean. The intervals [Mean-2Std, Mean + 2Std] are given in Table 2. These narrow intervals and small standard deviations show that the bootstrap estimates of parameter values closely distributed around the mean values. The parameter values used for the first-stage model (Table 1) are rather close to these mean values of bootstrap estimates and lie within the intervals of Mean ± 2 × Std. Further investigation of correlations between parameters has shown that the correlation coefficients of μ1 and L1 and of t2 and L2 are 0.77 and –0.96, respectively. Figure 6 shows the plots of their distributions. These plots also show that these bootstrap estimates are densely distributed around the mean values. All other parameters show no correlation or weak correlation with each other.
Parameter | μ1 | M1 | t1 | k1 | L1 |
Mean | 0.2064 | 0 | 29.3 | 0.1786 | 0.4565 |
Std | 0.004 | 0 | 0.9 | 0.0233 | 0.0081 |
Mean ± 2Std | [0.1984, 0.2144] | [0, 0] | [27.5, 31.1] | [0.1320, 0.2252] | [0.4403, 0.4727] |
Parameter | β1 | M2 | t2 | k2 | L2 |
Mean | 0.0485 | 0 | 105.7 | 0.0249 | 0.6395 |
Std | 0.0017 | 0 | 2 | 0.0004 | 0.0076 |
Mean ± 2Std | [0.0451, 0.0519] | [0, 0] | [101.7, 109.7] | [0.0241, 0.0257] | [0.6243, 0.6547] |
The parameter β1 is related to the lock-down effect. The infection rate is reduced by 1−1/(1+β1)=4.6% during lock-down, where β1=0.0484. The effect for lock-down is considered insignificant in this paper. The functions σ(t) and α(t) are plotted in Figure 7 to show the effect of social distancing and efficiency of detection of COVID-19 positive cases, respectively. Figure 7(a) shows that the social distancing factor starts with 0 initially, implying social distancing is considered a new concept in Tamil Nadu during the first wave. The social distancing factor reaches the plateau of 0.4533 around Day 50. Alvarez et al. [23] used social distancing factors increasing form 0 on March 10 to 0.75 on March 23, 2000, for New York City (NYC) in their simulation. They also showed that a social distancing factor of about 70% in Mexico City. Compared with the social distancing effect in NYC and Mexico City, our simulation result shows that the social distancing effect in Tamil Nadu increases relatively slowly and levels off at a relatively low level. Chakraborty et al. [43] have pointed out that the first wave of COVID-19 infection in India took off in March 2020, and the lack of the use of face covering and compliance with social distancing measures are among significant factors that cause the fast spread of COVID-19 in India.
Figure 7(b) shows that detection of positive cases does not start with zero performance meaning that although COVID-19 is a novel disease, health care professionals adopt new skills through prior experience. The function α(t) starts with a low level of 0.05 and slowly increases to a plateau of 0.64 over several months. In the simulation example for NYC, Alvarez et al. [23] used α values from 0 to 0.76 in two weeks, from March 20 to April 7, 2020. They also reported that South Korea implemented widespread testing with α≈0.98 in controlling the spread of the disease whereas Mexico City suffered from a low level of testing with α=0.1. Compared with NYC and South Korea, the detection of positive cases in Tamil Nadu is relatively ineffective.
To study the effect on the evolution of COVID-19 cases caused by the change in σ, we enhance L1 by 10%. Let L1=0.498 and all parameter values are kept the same as in Table 1. Figure 8(a) shows the evolution of cumulative cases X(t); the spread of COVID-19 can be effectively controlled by improving social distancing behaviors. Next, we enhance L2 by 10%. Let L2=0.704 which increases the detection rate. Figure 8(b) shows that the number of cumulative cases X(t) decreases dramatically. Comparing Figures 8(a) and (b), the effectiveness in controlling the spread of COVID-19 by enhancing detection rate is not as significant as that by improving social distancing measures. The simulation suggests that implementing non-pharmaceutical interventions may be an effective method to slow the spread of COVID-19.
It has been reported that the efficacy of all COVID-19 vaccines exceeds 70% [44], and the RNA based vaccines are highly efficacious, with an efficacy over 92% [44,45]. The vaccine efficacy varies over time [45]. The vaccination data used in this study contains individuals who have received one or two doses. Based on the above reported efficacy of COVID-19 vaccines, we assume a collective vaccine efficacy of 80% and let μ2=0.2μ1 in the parameter fitting process. We also assume that the vaccination willingness is the main factor determining the vaccination rate; a relatively small value h1=104 (persons) is used in this paper, compared with the population P0 in Tamil Nadu. For the utilization factor of health care resources, h2=5×105 (persons) is assumed, which is about one percent of P0. The daily reported data of COVID-19 cases and vaccination from 1 March 2021 to 4 August 2021 (156 days) are used to estimate the parameter values for the second-stage model.
Let ˜Xi, i=0,1,2,⋯,156, be the cumulative COVID-19 cases from daily reported data. We fit X1(t)+X2(t) to ˜Xi. Note that X1(t)+X2(t) represents the cumulative COVID-19 cases in the mathematical model. Let ˜Vi, i=0,1,2,⋯,156, be the number of vaccinated individuals from the vaccination dataset. We fit P0−(X1(t)+Y1(t)), which represents the number of vaccinated individuals in the mathematical model, to ˜Vi. Let ˜Ri, i=0,1,2,⋯,156, be the sum of Recovery and Deaths from daily reported data. The removed population R(t) is fitted to ˜Ri. The fitting process is the same as that used for the first-stage model. Equation (20) is approximated by the left-hand sum with dξ=0.01 day. Fourteen parameters listed in Table 3, including μ1, β2, Li, Mi, ki, ti, i=1,2,3, are to be fitted. The parameter values are determined by minimizing √∑Ti=0(X1(i)+X2(i)−˜Xi)2+√∑Ti=0(R(i)−˜Ri)2+0.3×√∑Ti=0(V(i)−˜Vi)2, where T=156. Figure 9 shows the fitting results with √∑Ti=0(X1(i)+X2(i)−˜Xi)2+√∑Ti=0(R(i)−˜Ri)2+0.3×√∑Ti=0(V(i)−˜Vi)2=3.6×106 (persons) and with the parameter values determined and summarized in Table 3. The infection rate μ1=0.227 day−1 is higher in the second wave than that in the first wave. India scientists have reported mutant variants which are more infectious in the second wave, and India experienced the exponential growth of COVID-19 infection beginning in March 2021 [43,46]. Figure 4 shows that in the first and second waves, the number of COVID-19 cases was growing exponentially in March 2020 and March 2021, respectively.
Parameter | Value | Brief Description | Units | Reference |
P0 | 7.8×107 | Population of Tamil Nadu State, India | individual | [40] |
μ1 | 0.227 | Intrinsic infection rate | Day−1 | fitted with data [38,39] |
μ2 | 0.0454 | Infection rate after vaccinated | Day−1 | assumed, based on [44,45] |
β1 | 0.0484 | Unlock-down factor | Dimensionless | fitted with data [38] |
h1 | 104 | HSC of vaccination | Individual | assumed |
M1 | 0.43 | Pre-experience in σ | Dimensionless | fitted with data [38,39] |
t1 | 83.9 | Time parameter in σ | Day | fitted with data [38,39] |
k1 | 0.2421 | Steepness of σ | Day−1 | fitted with data [38,39] |
L1 | 0.265 | Maximum increase in σ | Dimensionless | fitted with data [38,39] |
h2 | 5×105 | HSC of HCR | Individual | assumed |
M2 | 0.3491 | Pre-experience in α | Dimensionless | fitted with data [38,39] |
t2 | -5.8882 | Time parameter in α | Day | fitted with data [38,39] |
k2 | -0.6178 | Steepness of α | Day−1 | fitted with data [38,39] |
L2 | 0.5406 | Maximum increase in α | Dimensionless | fitted with data [38,39] |
M3 | 0.0216 | Pre-experience in γ | Dimensionless | fitted with data [38,39] |
t3 | 282.2 | Time parameter in γ | Day | fitted with data [38,39] |
k3 | 0.0163 | Steepness of γ | Day−1 | fitted with data [38,39] |
L3 | 0.7209 | Maximum increase in γ | Dimensionless | fitted with data [38,39] |
β2 | 3155833 | Maximum vaccination rate | Individual Day−1 | fitted with data [38,39] |
τ1 | 4 | Delay in detection | Day | [23] |
τ2 | 14 | Recovery time | Day | [30] |
τ3 | 14 | Hospitalization time | Day | [30] |
The bootstrap statistics for the parameters of the second-stage model are shown in Table 4. Similarly, small standard deviations for bootstrap estimates are found for the second-stage model. The parameter values used for the second-stage model (Table 3) are rather close to these mean values. The bootstrap statistics show no correlation or weak correlation between all parameters except for μ1 and t1 with a correlation coefficient of –0.9 (Figure 10).
Parameter | μ1 | M1 | t1 | k1 | L1 |
Mean | 0.2339 | 0.4304 | 82.2 | 0.2363 | 0.2658 |
Std | 0.01 | 0.0011 | 2.6 | 0.0391 | 0.0082 |
Mean ± 2Std | [0.2139, 0.2539] | [0.4282, 0.4326] | [77, 87.4] | [0.1581, 0.3145] | [0.2494, 0.2822] |
Parameter | M2 | t2 | k2 | L2 | |
Mean | 0.3574 | -5.53 | -0.6341 | 0.5625 | |
Std | 0.0219 | 0.52 | 0.0803 | 0.0318 | |
Mean ± 2Std | [0.3136, 0.4012] | [-6.57, -4.49] | [-0.7947, -0.4735] | [0.4989, 0.6261] | |
Parameter | M3 | t3 | k3 | L3 | β2 |
Mean | 0.0221 | 282 | 0.0165 | 0.7264 | 3.1757×106 |
Std | 0.0015 | 7 | 0.0007 | 0.0528 | 1.872×105 |
Mean ± 2Std | [0.0191, 0.0251] | [268,296] | [0.0151, 0.0179] | [0.6208, 0.832] | [2.8013,3.5501]×106 |
During the second wave, major hospitals faced a shortage of hospital beds, medicines, oxygen and ventilator supplies, leading to high mortality [43,46]. The detection rate α(t), shown in Figure 11(a), is lower compared with that in the first wave (Figure 7(b)). The parameters k2 and t2 are negative indicating that the detection rate α(t) is decreasing and α(t), for t≥0, contains only the lower part of the curve. The detection rate decreases quickly initially, levels off soon, and remains at a low level. Combining this effect with shortage of health care resources, Figure 11(b) shows that the overall rate of hospitalization of the positive cases continue to decrease over time. Shortage of health care resources is one of the factors responsible for fast increased number of cases in Tamil Nadu [46]. Jain et al. [46] suggested that prevention was better than cure for this disease and the implementation of social distancing and the use of face masks was able to help with controlling the disease. Figure 11(c) shows that during the first two months of the second wave, implementation of social distancing is not effective and remains at the same level as that in the first wave. Then, the social distancing effect increases during May and Jun 2021 and levels off at about 0.7 which is much higher than that during the first wave. Chakraborty et al. [43] pointed out that social distancing was not strictly implemented during the first wave in India due to crowded style of living arrangement. Our study shows that people in Tamil Nadu may have learned their lessons and follow social distancing measures after some time in the second wave. Nevertheless, a factor of 0.7 for social distancing measures is not effective compared with other cities studied in [23].
In this paper, mathematical models using deterministic ODEs with time delays are proposed for the study of impacts of pharmaceutical and non-pharmaceutical factors on the COVID-19 pandemic for the state of Tamil Nadu, India. The first-stage model, designed for the first wave in the region, considers two populations, the population of infected individuals and the population of removed individuals. The second-stage model, designed for the second wave, considers four populations, unvaccinated infected population, vaccinated infected population, unvaccinated uninfected population, and removed population. Some factors such as social distancing, detection efficiency, and vaccination willingness are modeled by learning curves. Daily reported data supported by CSSE at Coronavirus Resource Center of JHU are used to fit the model [38,39] and to determine the parameter values. The model produces a close match to the evolution of COVID-19 in Tamil Nadu (Figures 5 and 9).
The preliminary results obtained from the simulation indicate that the lock down contributed limited effects on controlling the spread of the disease in Tamil Nadu. The social distancing implementation and detection of positive cases are relative ineffective compared with some big cities studied by Alvarez et al. [23] in the first wave. These observations are in agreement with the research report that social distancing was not strictly implemented during the first wave in India [43]. The implementation of social distancing remains at a low level during the first half period of the second wave but is improved after some time in the second wave (Figure 11(c)). The detection efficiency and hospitalization rate decrease rapidly and remains at low levels over the second wave (Figure 11(a), (b)). Jain et al. [46] pointed out that Tamil Nadu faced shortage of health care resources during the second wave of COVID-19. The impacts of other factors will also need to be considered and is currently under future research directions of the authors.
A Shankaranarayanan was supported by the Ministry of Science and Technology (MOST) of Taiwan under the grant no. MOST 110-2222-E-035-002; H.-C. Wei was supported by the MOST of Taiwan under the grant no. MOST 110-2115-M-035-004.
The authors declare there is no conflict of interest.
[1] | FAO (2022) FAOSTAT, Land use, 2022. Available from: https://www.fao.org/faostat/en/#data/RL. |
[2] | FAO (2022) FAOSTAT, Crops and livestock products, 2022. Available from: https://www.fao.org/faostat/en/#data/QCL. |
[3] | FAO (2022) FAOSTAT, Value of agricultural production, 2022. Available from: https://www.fao.org/faostat/en/#data/QV. |
[4] | FAO (2022) FAOSTAT, Macro indicators, 2022. Available from: https://www.fao.org/faostat/en/#data/MK. |
[5] | SIAP (2022) Statistical yearbook of agricultural production, 2022. Available from: https://nube.siap.gob.mx/cierreagricola/. |
[6] | SIAP (2023) Performance of the agri-food GDP in the third quarter of 2023 (2022: Ⅲ–2023 Ⅲ). |
[7] | SEDEC (2023) State profile, 2023. Available from: https://tabasco.gob.mx/sites/default/files/users/sdettabasco/Perfil%20del%20Estado.pdf. |
[8] | INEGI Quarterly indicator of State economic activity (ITAEE), Tabasco. |
[9] |
Chauhan A, Upadhyay S, Saini G, et al. (2022) Agricultural crop residue-based biomass in India: Potential assessment, methodology and key issues. Sustainable Energy Technol Assess 53: 102552.https://doi.org/10.1016/j.seta.2022.102552 doi: 10.1016/j.seta.2022.102552
![]() |
[10] |
Tauro R, García C, Skutsch M, et al. (2018) The potential for sustainable biomass pellets in Mexico: An analysis of energy potential, logistic costs and market demand. Renewable Sustainable Energy Rev 82: 380–389.https://doi.org/10.1016/j.rser.2017.09.036 doi: 10.1016/j.rser.2017.09.036
![]() |
[11] |
Peñaloza DF, Laiton LJ, Caballero DF, et al. (2012) Sciencimetric study of trends in the utilization of cocoa by-products (Theobroma cacao L.). Rev Espacio I+D, Innov Desarro 10: 83–94.https://doi.org/10.31644/IMASD.27.2021.a05 doi: 10.31644/IMASD.27.2021.a05
![]() |
[12] | Balladares CA (2016) Physicochemical characterization of cocoa and coffee leachates from the Ecuadorian coast, as potential sources of bioethanol production. Available from: http://hdl.handle.net/10553/22931. |
[13] |
Lu F, Rodriguez J, Van Damme I, et al. (2018) Valorisation strategies for cocoa pod husk and its fractions. Curr Opin Green Sustainable Chem 14: 80–88.https://doi.org/10.1016/j.cogsc.2018.07.007 doi: 10.1016/j.cogsc.2018.07.007
![]() |
[14] |
Vásquez ZS, de Carvalho Neto DP, Pereira GVM, et al. (2019) Biotechnological approaches for cocoa waste management: A review. Waste Manage 90: 72–83.https://doi.org/10.1016/j.wasman.2019.04.030 doi: 10.1016/j.wasman.2019.04.030
![]() |
[15] |
Campos R, Nieto KH, Oomah BD (2018) Cocoa (Theobroma cacao L.) pod husk: Renewable source of bioactive compounds. Trends Food Sci Technol 81: 172–184.https://doi.org/10.1016/j.tifs.2018.09.022 doi: 10.1016/j.tifs.2018.09.022
![]() |
[16] |
Fernandes ERK, Marangoni C, Souza O, et al. (2013) Thermochemical characterization of banana leaves as a potential energy source. Energy Convers Manage 75: 603–608.https://doi.org/10.1016/j.enconman.2013.08.008 doi: 10.1016/j.enconman.2013.08.008
![]() |
[17] |
Espina R, Barroca R, Abundo MLS (2022) The optimal high heating value of the torrefied coconut shells. Eng Technol Appl Sci Res 12: 8605–8610.https://doi.org/10.48084/etasr.4931 doi: 10.48084/etasr.4931
![]() |
[18] | Trujillo AF, Arias LS (2013) Coconut, a renewable resource for the design of green materials. Entre Cienc Ing 7: 93–100. Available from:https://revistas.ucp.edu.co/index.php/entrecienciaeingenieria/article/view/637. |
[19] |
Irawan A, US Latifah, DIP Meity (2017) Effect of torrefaction process on the coconut shell energy content for solid fuel, AIP Conf Proc 1826: 020010.https://doi.org/10.1063/1.4979226 doi: 10.1063/1.4979226
![]() |
[20] | Triana O, León TS, Céspedes MI, et al. (2014) Characterization of sugarcane harvest residues stored in bulk. ICIDCA Sobre Deriv Caña Azúc 48: 65–70. Available from:https://www.redalyc.org/articulo.oa?id = 223131337010. |
[21] |
Kalifa MA, Habtu NG, Jembere AL, et al. (2024) Characterization and evaluation of torrefied sugarcane bagasse to improve the fuel properties. Curr Res Green Sustain Chem 8: 100395.https://doi.org/10.1016/j.crgsc.2023.100395 doi: 10.1016/j.crgsc.2023.100395
![]() |
[22] |
Schmitt CC, Moreira R, Neves RC, et al. (2020) From agriculture residue to upgraded product: The thermochemical conversion of sugarcane bagasse for fuel and chemical products. Fuel Proc Technol 197: 106199.https://doi.org/10.1016/j.fuproc.2019.106199 doi: 10.1016/j.fuproc.2019.106199
![]() |
[23] |
Liew RK, Nam WL, Chong MY, et al. (2018) Oil palm waste: An abundant and promising feedstock for microwave pyrolysis conversion into good quality biochar with potential multi-applications. Proc Saf Environ Prot 115: 57–69.https://doi.org/10.1016/j.psep.2017.10.005 doi: 10.1016/j.psep.2017.10.005
![]() |
[24] |
Ohimain EI, Izah SC (2014) Energy self-sufficiency of smallholder oil palm processing in Nigeria. Renewable Energy 63: 426–431.https://doi.org/10.1016/j.renene.2013.10.007 doi: 10.1016/j.renene.2013.10.007
![]() |
[25] |
Syamsiro M, Saptoadi H, Tambunan BH, et al. (2012) A preliminary study on use of cocoa pod husk as a renewable source of energy in Indonesia. Energy Sustainable Dev 16: 74–77.https://doi.org/10.1016/j.esd.2011.10.005 doi: 10.1016/j.esd.2011.10.005
![]() |
[26] |
Vásquez ZS, de Carvalho Neto DP, Pereira GVM, et al. (2019) Biotechnological approaches for cocoa waste management: A review. Waste Manage 90: 72–83.https://doi.org/10.1016/j.wasman.2019.04.030 doi: 10.1016/j.wasman.2019.04.030
![]() |
[27] |
Redondo C, Rodríguez M, Vallejo S, et al. (2020) Biorefinery of biomass of agro-industrial banana waste to obtain high-value biopolymers. Molecules 25: 3829.https://doi.org/10.3390/molecules25173829 doi: 10.3390/molecules25173829
![]() |
[28] |
da Silva JCG, Alves JLF, de Araujo WV, et al. (2019) Pyrolysis kinetics and physicochemical characteristics of skin, husk, and shell from green coconut wastes. Energy Ecol Environ 4: 125–132.https://doi.org/10.1007/s40974-019-00120-x doi: 10.1007/s40974-019-00120-x
![]() |
[29] |
Sarkar JK, Wang Q (2020) Different pyrolysis process conditions of South Asian waste coconut shell and characterization of gas, bio-char, and bio-oil. Energies 13: 1970.https://doi.org/10.3390/en13081970 doi: 10.3390/en13081970
![]() |
[30] |
Nadzri SNIHA, Sultan MTH, Shah AUM, et al. (2022) A comprehensive review of coconut shell powder composites: Preparation, processing, and characterization. J Thermoplast Compos Mater 35: 2641–2664.https://doi.org/10.1177/089270572093080 doi: 10.1177/089270572093080
![]() |
[31] |
Dewajani H, Zamrudy W, Irfin Z, et al. (2023) Utilization of Indonesian sugarcane bagasse into bio asphalt through pyrolysis process using zeolite-based catalyst. Mater Today Proc 87: 383–389.https://doi.org/10.1016/j.matpr.2023.04.171 doi: 10.1016/j.matpr.2023.04.171
![]() |
[32] |
Sales A, Lima SA (2010) Use of Brazilian sugarcane bagasse ash in concrete as sand replacement. Waste Manage 30: 1114–1122.https://doi.org/10.1016/j.wasman.2010.01.026 doi: 10.1016/j.wasman.2010.01.026
![]() |
[33] |
Abnisa F, Daud WMAW, Husin WNW, et al. (2011) Utilization possibilities of palm shell as a source of biomass energy in Malaysia by producing bio-oil in pyrolysis process. Biomass Bioenergy 35: 1863–1872.https://doi.org/10.1016/j.biombioe.2011.01.033 doi: 10.1016/j.biombioe.2011.01.033
![]() |
[34] |
Sosa JA, Laines JR, García DS, et al. (2022) Activated carbon: A review of residual precursors, synthesis processes, characterization techniques, and applications in the improvement of biogas. Environ Eng Res 28: 220100.https://doi.org/10.4491/eer.2022.100 doi: 10.4491/eer.2022.100
![]() |
[35] |
Pérez M de L, Hernández JC, Bideshi DK, et al. (2020) Agave: a natural renewable resource with multiple applications. J Sci Food Agric 100: 5324–5333.https://doi.org/10.1002/jsfa.10586 doi: 10.1002/jsfa.10586
![]() |
[36] |
Borrega M, Hinkka V, Hörhammer H, et al. (2022) Utilizing and valorizing oat and barley straw as an alternative source of lignocellulosic fibers. Materials 15: 7826.https://doi.org/10.3390/ma15217826 doi: 10.3390/ma15217826
![]() |
[37] |
Ali AH, Wanderlind EH, Almerindo GI (2024) Activated carbon obtained from malt bagasse as a support in heterogeneous catalysis for biodiesel production. Renewable Energy 220: 119656.https://doi.org/10.1016/j.renene.2023.119656 doi: 10.1016/j.renene.2023.119656
![]() |
[38] |
Chung WJ, Shim J, Ravindran B (2022) Application of wheat bran-based biomaterials and nano-catalyst in textile wastewater. J King Saud Univ Sci 34: 101775.https://doi.org/10.1016/j.jksus.2021.101775 doi: 10.1016/j.jksus.2021.101775
![]() |
[39] |
Armynah B, Tahir D, Tandilayuk M, et al. (2019) Potentials of biochars derived from bamboo leaf biomass as energy sources: Effect of temperature and time of heating. Int J Biomater 2019: 1–9.https://doi.org/10.1155/2019/3526145 doi: 10.1155/2019/3526145
![]() |
[40] |
Sagastume A, Cabello JJ, Hens L, et al. (2020) The energy potential of agriculture, agroindustrial, livestock, and slaughterhouse biomass wastes through direct combustion and anaerobic digestion. The case of Colombia. J Clean Prod 269: 122317.https://doi.org/10.1016/j.jclepro.2020.122317 doi: 10.1016/j.jclepro.2020.122317
![]() |
[41] |
Mdhluli FT, Harding KG (2021) Comparative life-cycle assessment of maize cobs, maize stover and wheat stalks for the production of electricity through gasification vs traditional coal power electricity in South Africa. Clean Environ Syst 3: 100046.https://doi.org/10.1016/j.cesys.2021.100046 doi: 10.1016/j.cesys.2021.100046
![]() |
[42] |
Appiah NB, Li J, Rooney W, et al. (2019) A review of sweet sorghum as a viable renewable bioenergy crop and its techno-economic analysis. Renewable Energy 143: 1121–1132.https://doi.org/10.1016/j.renene.2019.05.066 doi: 10.1016/j.renene.2019.05.066
![]() |
[43] |
Mazurkiewicz J, Marczuk A, Pochwatka P, et al. (2019) Maize straw as a valuable energetic material for giogas plant feeding. Materials 12: 3848.https://doi.org/10.3390/ma12233848 doi: 10.3390/ma12233848
![]() |
[44] |
Niju S, Swathika M, Balajii M (2020) Pretreatment of lignocellulosic sugarcane leaves and tops for bioethanol production. Lignocellulosic Biomass Liquid Biofuels, 301–324.https://doi.org/10.1016/B978-0-12-815936-1.00010-1 doi: 10.1016/B978-0-12-815936-1.00010-1
![]() |
[45] |
Chala B, Oechsner H, Latif S, et al. (2018) Biogas potential of coffee processing waste in Ethiopia. Sustainability 10: 2678.https://doi.org/10.3390/su10082678 doi: 10.3390/su10082678
![]() |
[46] |
Mekuria D, Diro A, Melak F, et al. (2022) Adsorptive removal of methylene blue dye using biowaste materials: Barley bran and enset midrib leaf. J Chem 2022: 1–13.https://doi.org/10.1155/2022/4849758 doi: 10.1155/2022/4849758
![]() |
[47] |
Kosheleva RI, Mitropoulos AC, Kyzas GZ (2019) Synthesis of activated carbon from food waste. Environ Chem Lett 17: 429–438.https://doi.org/10.1007/s10311-018-0817-5 doi: 10.1007/s10311-018-0817-5
![]() |
[48] |
Prieto García F, Canales-Flores RA, Prieo-Méndez J, et al. (2022) Evaluation of three lignocellulose biomass materials (barley husk, corn cobs, agave leaves) as precursors of activated carbon. Rev Fac Cienc 11: 17–39.https://doi.org/10.15446/rev.fac.cienc.v11n1.97719 doi: 10.15446/rev.fac.cienc.v11n1.97719
![]() |
[49] |
Araújo L, Machado AR, Pintado M, et al. (2023) Toward a circular bioeconomy: Extracting cellulose from grape stalks. Eng Proc 86.https://doi.org/10.3390/ECP2023-14746 doi: 10.3390/ECP2023-14746
![]() |
[50] |
Martinelli FRB, Ribeiro FRC, Marvila MT, et al. (2023) A Review of the use of coconut fiber in cement composites. Polymers 15: 1309.https://doi.org/10.3390/polym15051309 doi: 10.3390/polym15051309
![]() |
[51] |
Elnagdy NA, Ragab TIM, Fadel MA, et al. (2024) Bioethanol production from characterized pre-treated sugarcane trash and Jatropha agrowastes. J Biotechnol 386: 28–41.https://doi.org/10.1016/j.jbiotec.2024.02.015 doi: 10.1016/j.jbiotec.2024.02.015
![]() |
[52] | Darmayanti R, Wika Amini H, Fitri Rizkiana M, et al. (2019) Lignocellulosic material from main indonesian plantation commodity as the feedstock for fermentable sugar in biofuel production. ARPN J Eng Appl Sci 14: 3524–3534. |
[53] |
Perea MA, Manzano F, Hernandez Q, et al. (2018) Peanut shell for energy: Properties and its potential to respect the environment. Sustainability 10: 3254.https://doi.org/10.3390/su10093254 doi: 10.3390/su10093254
![]() |
[54] |
Daud WMAW, Ali WSW (2004) Comparison on pore development of activated carbon produced from palm shell and coconut shell. Bioresour Technol 93: 63–69.https://doi.org/10.1016/j.biortech.2003.09.015 doi: 10.1016/j.biortech.2003.09.015
![]() |
[55] |
Abnisa F, Arami A, Daud WMAW, et al. (2013) Characterization of bio-oil and bio-char from pyrolysis of palm oil wastes. Bioenergy Res 6: 830–840.https://doi.org/10.1007/s12155-013-9313-8 doi: 10.1007/s12155-013-9313-8
![]() |
[56] |
Sable H, Kumar V, Mishra R, et al. (2024) Bamboo stem derived biochar for biosorption of Cadmium (Ⅱ) ions from contaminated wastewater. Environ Nanotechnol Monit Manag 21: 100936.https://doi.org/10.1016/j.enmm.2024.100936 doi: 10.1016/j.enmm.2024.100936
![]() |
[57] |
Manríquez A, Sierra J, Muñoz P, et al. (2020) Analysis of urban agriculture solid waste in the frame of circular economy: Case study of tomato crop in integrated rooftop greenhouse. Sci Total Environ 734.https://doi.org/10.1016/j.scitotenv.2020.139375 doi: 10.1016/j.scitotenv.2020.139375
![]() |
[58] |
Soares IS, Perrechil F, Grandis A, et al. (2024) Cassava waste (stem and leaf) analysis for reuse. Food Chem Adv 4: 100675.https://doi.org/10.1016/j.focha.2024.100675 doi: 10.1016/j.focha.2024.100675
![]() |
[59] |
Gómez E, Roriz CL, Heleno SA, et al. (2021) Valorisation of black mulberry and grape seeds: Chemical characterization and bioactive potential. Food Chem 337: 127998.https://doi.org/10.1016/j.foodchem.2020.127998 doi: 10.1016/j.foodchem.2020.127998
![]() |
[60] |
Al Afif R, Pfeifer C, Pröll T (2020) Bioenergy recovery from cotton stalk. Advances Cotton Res. https://doi.org/10.5772/intechopen.88005 doi: 10.5772/intechopen.88005
![]() |
[61] |
Stavropoulos GG, Zabaniotou AA (2005) Production and characterization of activated carbons from olive-seed waste residue. Microporous Mesoporous Mater 82: 79–85.https://doi.org/10.1016/j.micromeso.2005.03.009 doi: 10.1016/j.micromeso.2005.03.009
![]() |
[62] |
Dominguez EL, Uttran A, Loh SK, et al. (2020) Characterisation of industrially produced oil palm kernel shell biochar and its potential as slow release nitrogen-phosphate fertilizer and carbon sink. Mater Today Proc 31: 221–227.https://doi.org/10.1016/j.matpr.2020.05.143 doi: 10.1016/j.matpr.2020.05.143
![]() |
[63] |
Toro JL, Carrillo ES, Bustos D, et al. (2019) Thermogravimetric characterization and pyrolysis of soybean hulls. Bioresour Technol Rep 6: 183–189.https://doi.org/10.1016/j.biteb.2019.02.009 doi: 10.1016/j.biteb.2019.02.009
![]() |
[64] |
Marafon AC, Salomon KR, Amorim ELC, et al. (2020) Use of sugarcane vinasse to biogas, bioenergy, and biofertilizer production. Sugarcane Biorefinery, Technol Perspectives, 179–194.https://doi.org/10.1016/B978-0-12-814236-3.00010-X doi: 10.1016/B978-0-12-814236-3.00010-X
![]() |
[65] |
Sindhu R, Binod P, Pandey A, et al. (2019) Biofuel production from biomass. Current Dev Biotechnol Bioeng, 79–92.https://doi.org/10.1016/B978-0-444-64083-3.00005-1 doi: 10.1016/B978-0-444-64083-3.00005-1
![]() |
[66] | Solís JA, Morales M, Ayala RC, et al. (2012) Obtaining activated carbon from agro-industrial waste and its evaluation in the removal of color from sugarcane juice. Tecnol, Cienc, Educ 27: 36–48. Available from:https://www.redalyc.org/articulo.oa?id = 48224413006. |
[67] |
Medhat A, El-Maghrabi HH, Abdelghany A, et al. (2021) Efficiently activated carbons from corn cob for methylene blue adsorption. Appl Surface Sci Advances 3: 100037.https://doi.org/10.1016/j.apsadv.2020.100037 doi: 10.1016/j.apsadv.2020.100037
![]() |
[68] |
Abbey CYB, Duwiejuah AB, Quianoo AK (2023) Removal of toxic metals from aqueous phase using cacao pod husk biochar in the era of green chemistry. Appl Water Sci 13: 57.https://doi.org/10.1007/s13201-022-01863-5 doi: 10.1007/s13201-022-01863-5
![]() |
[69] | Marín FJ, García RM, Barrezueta SA (2020) Results of the application of biochar obtained from banana and cocoa residues in corn cultivation. Rev Cient Agroecosist 8: 83–88. Available from:https://aes.ucf.edu.cu/index.php/aes/article/view/404. |
[70] |
Pinzon DA, Adarme CA, Vargas LY, et al. (2022) Biochar as a waste management strategy for cadmium contaminated cocoa pod husk residues. Int JRecycl Org Waste Agric 11: 101–115.https://doi.org/10.30486/ijrowa.2021.1920124.1192 doi: 10.30486/ijrowa.2021.1920124.1192
![]() |
[71] |
Nosratpour MJ, Karimi K, Sadeghi M (2018) Improvement of ethanol and biogas production from sugarcane bagasse using sodium alkaline pretreatments. J Environ Manage 226: 329–339.https://doi.org/10.1016/j.jenvman.2018.08.058 doi: 10.1016/j.jenvman.2018.08.058
![]() |
[72] | Santos CV, Lourenzani AEBS, Mollo M, et al. (2021) Study of the biogas potential generated from residue: peanut shells. Rev Bras Ciênc Ambient (Online) 56: 318–326.https://doi.org/10.5327/Z21769478765 |
[73] |
Tsai WT, Lee MK, Chang YM (2006) Fast pyrolysis of rice straw, sugarcane bagasse and coconut shell in an induction-heating reactor. J Anal Appl Pyrolysis 76: 230–237.https://doi.org/10.1016/j.jaap.2005.11.007 doi: 10.1016/j.jaap.2005.11.007
![]() |
[74] |
Zheng J, Yi W, Wang N (2008) Bio-oil production from cotton stalk. Energy Convers Manage 49: 1724–1730.https://doi.org/10.1016/j.enconman.2007.11.005 doi: 10.1016/j.enconman.2007.11.005
![]() |
[75] |
Kim SW, Koo BS, Lee DH (2014) Catalytic pyrolysis of palm kernel shell waste in a fluidized bed. Bioresour Technol 167: 425–432.https://doi.org/10.1016/j.biortech.2014.06.050 doi: 10.1016/j.biortech.2014.06.050
![]() |
[76] |
Robak K, Balcerek M (2018) Review of second-generation bioethanol production from residual biomass. Food Technol Biotechnol 56: 174–187.https://doi.org/10.17113/ftb.56.02.18.5428 doi: 10.17113/ftb.56.02.18.5428
![]() |
[77] |
Guerrero AB, Ballesteros I, Ballesteros M (2018) The potential of agricultural banana waste for bioethanol production. Fuel 213: 176–185.https://doi.org/10.1016/j.fuel.2017.10.105 doi: 10.1016/j.fuel.2017.10.105
![]() |
[78] | Kaur M, Kaur M (2012) A review on utilization of coconut shell as coarse aggregates in mass concrete. Int J Appl Eng Res 7: 2063–2065. |
[79] |
Tomar R, Kishore K, Singh Parihar H, et al. (2021) A comprehensive study of waste coconut shell aggregate as raw material in concrete. Mater Today Proc 44: 437–443.https://doi.org/10.1016/j.matpr.2020.09.754 doi: 10.1016/j.matpr.2020.09.754
![]() |
[80] |
Chávez V, Valencia A, Córdova C, et al. (2017) Banana stem leachates: Obtaining and potential uses. Cuad Biodivers 53: 1–8.https://doi.org/10.14198/cdbio.2017.53.01 doi: 10.14198/cdbio.2017.53.01
![]() |
[81] |
Abdul Wahid FRA, Saleh S, Abdul Samad NAF (2017) Estimation of higher heating value of torrefied palm oil wastes from proximate analysis. Energy Proc 138: 307–312.https://doi.org/10.1016/j.egypro.2017.10.102 doi: 10.1016/j.egypro.2017.10.102
![]() |
[82] |
Florian TDM, Villani N, Aguedo M, et al. (2019) Chemical composition analysis and structural features of banana rachis lignin extracted by two organosolv methods. Ind Crops Prod 132: 269–274.https://doi.org/10.1016/j.indcrop.2019.02.022 doi: 10.1016/j.indcrop.2019.02.022
![]() |
[83] |
Granados DA, Velásquez HI, Chejne F (2014) Energetic and exergetic evaluation of residual biomass in a torrefaction process. Energy 74: 181–189.https://doi.org/10.1016/j.energy.2014.05.046 doi: 10.1016/j.energy.2014.05.046
![]() |
[84] | Jirón EG, Rodríguez K, Bernal C (2020) Cellulose nanofiber production from banana rachis. IJESC 10: 24683–24689. |
[85] |
Balogun AO, Lasode OA, McDonald AG (2018) Thermochemical and pyrolytic analyses of Musa spp. residues from the rainforest belt of Nigeria. Environ Prog Sustainable Energy 37: 1932–1941.https://doi.org/10.1002/ep.12869 doi: 10.1002/ep.12869
![]() |
[86] |
Meramo SI, Ojeda KA, Sanchez E (2019) Environmental and safety assessments of industrial production of levulinic acid via acid-catalyzed dehydration. ACS Omega 4: 22302–22312.https://doi.org/10.1021/acsomega.9b02231 doi: 10.1021/acsomega.9b02231
![]() |
[87] |
Guerrero AB, Aguado PL, Sánchez J, et al. (2016) GIS-Based assessment of banana residual biomass potential for ethanol production and power generation: A case study. Waste Biomass Valor 7: 405–415.https://doi.org/10.1007/s12649-015-9455-3 doi: 10.1007/s12649-015-9455-3
![]() |
[88] |
Ozyuguran A, Akturk A, Yaman S (2018) Optimal use of condensed parameters of ultimate analysis to predict the calorific value of biomass. Fuel 214: 640–646.https://doi.org/10.1016/j.fuel.2017.10.082 doi: 10.1016/j.fuel.2017.10.082
![]() |
[89] |
Titiloye JO, Abu Bakar MS, Odetoye TE (2013) Thermochemical characterisation of agricultural wastes from West Africa. Ind Crops Prod 47: 199–203.https://doi.org/10.1016/j.indcrop.2013.03.011 doi: 10.1016/j.indcrop.2013.03.011
![]() |
[90] |
Ghysels S, Acosta N, Estrada A, et al. (2020) Integrating anaerobic digestion and slow pyrolysis improves the product portfolio of a cocoa waste biorefinery. Sustainable Energy Fuels 4: 3712–3725.https://doi.org/10.1039/D0SE00689K doi: 10.1039/D0SE00689K
![]() |
[91] |
Akinola AO, Eiche JF, Owolabi PO, et al. (2018) Pyrolytic analysis of cocoa pod for biofuel production. Niger J Technol 37: 1026.https://doi.org/10.4314/njt.374.1866 doi: 10.4314/njt.374.1866
![]() |
[92] |
Londoño-Larrea P, Villamarin-Barriga E, García AN, et al. (2022) Study of cocoa pod husks thermal decomposition. Appl Sci 12: 9318.https://doi.org/10.3390/app12189318 doi: 10.3390/app12189318
![]() |
[93] |
Tsai C, Tsai W, Liu S, et al. (2018) Thermochemical characterization of biochar from cocoa pod husk prepared at low pyrolysis temperature. Biomass Convers Biorefin 8: 237–243.https://doi.org/10.1007/s13399-017-0259-5 doi: 10.1007/s13399-017-0259-5
![]() |
[94] |
Adjin M, Asiedu N, Dodoo D, et al. (2018) Thermochemical conversion and characterization of cocoa pod husks a potential agricultural waste from Ghana. Ind Crops Prod 119: 304–312.https://doi.org/10.1016/j.indcrop.2018.02.060 doi: 10.1016/j.indcrop.2018.02.060
![]() |
[95] |
Channiwala SA, Parikh PP (2002) A unified correlation for estimating HHV of solid, liquid and gaseous fuels. Fuel 81: 1051–1063.https://doi.org/10.1016/S0016-2361(01)00131-4 doi: 10.1016/S0016-2361(01)00131-4
![]() |
[96] |
Kabir R, Anwar S, Yusup S, et al. (2022) Exploring the potential of coconut shell biomass for charcoal production. Ain Shams Eng J 13: 101499.https://doi.org/10.1016/j.asej.2021.05.013 doi: 10.1016/j.asej.2021.05.013
![]() |
[97] |
Borel LDMS, de Lira TS, Ataíde CH, et al. (2021) Thermochemical conversion of coconut waste: Material characterization and identification of pyrolysis products. J Therm Anal Calorim 143: 637–646.https://doi.org/10.1007/s10973-020-09281-y doi: 10.1007/s10973-020-09281-y
![]() |
[98] |
Said M, John G, Mhilu C, et al. (2015) The study of kinetic properties and analytical pyrolysis of coconut shells. J Renewable Energy 2015: 1–8.https://doi.org/10.1155/2015/307329 doi: 10.1155/2015/307329
![]() |
[99] |
Rout T, Pradhan D, Singh RK, et al. (2016) Exhaustive study of products obtained from coconut shell pyrolysis. J Environ Chem Eng 4: 3696–3705.https://doi.org/10.1016/j.jece.2016.02.024 doi: 10.1016/j.jece.2016.02.024
![]() |
[100] |
Gani A, Erdiwansyah, Desvita H, et al. (2024) Comparative analysis of HHV and LHV values of biocoke fuel from palm oil mill solid waste. Case Stud Chem EnvironEng 9: 100581.https://doi.org/10.1016/j.cscee.2023.100581 doi: 10.1016/j.cscee.2023.100581
![]() |
[101] |
Nizamuddin S, Jayakumar NS, Sahu JN, et al. (2015) Hydrothermal carbonization of oil palm shell. Korean J Chem Eng 32: 1789–1797.https://doi.org/10.1007/s11814-014-0376-9 doi: 10.1007/s11814-014-0376-9
![]() |
[102] |
Kim SJ, Jung SH, Kim JS (2010) Fast pyrolysis of palm kernel shells: Influence of operation parameters on the bio-oil yield and the yield of phenol and phenolic compounds. Bioresour Technol 101: 9294–9300.https://doi.org/10.1016/j.biortech.2010.06.110 doi: 10.1016/j.biortech.2010.06.110
![]() |
[103] |
Uemura Y, Omar WN, Tsutsui T, et al. (2011) Torrefaction of oil palm wastes. Fuel 90: 2585–2591.https://doi.org/10.1016/j.fuel.2011.03.021 doi: 10.1016/j.fuel.2011.03.021
![]() |
[104] |
Chang G, Huang Y, Xie J, et al. (2016) The lignin pyrolysis composition and pyrolysis products of palm kernel shell, wheat straw, and pine sawdust. Energy Convers Manage 124: 587–597.https://doi.org/10.1016/j.enconman.2016.07.038 doi: 10.1016/j.enconman.2016.07.038
![]() |
[105] |
Marrugo G, Valdés CF, Chejne F (2016) Characterization of colombian agroindustrial biomass residues as energy resources. Energy Fuels 30: 8386–8398.https://doi.org/10.1021/acs.energyfuels.6b01596 doi: 10.1021/acs.energyfuels.6b01596
![]() |
[106] |
Ma Z, Chen D, Gu J, et al. (2015) Determination of pyrolysis characteristics and kinetics of palm kernel shell using TGA–FTIR and model-free integral methods. Energy Convers Manage 89: 251–259.https://doi.org/10.1016/j.enconman.2014.09.074 doi: 10.1016/j.enconman.2014.09.074
![]() |
[107] |
Liew RK, Chong MY, Osazuwa OU, et al. (2018) Production of activated carbon as catalyst support by microwave pyrolysis of palm kernel shell: a comparative study of chemical versus physical activation. Res Chem Intermed 44: 3849–3865.https://doi.org/10.1007/s11164-018-3388-y doi: 10.1007/s11164-018-3388-y
![]() |
[108] |
Athira G, Bahurudeen A, Appari S (2021) Thermochemical conversion of sugarcane Bagasse: composition, reaction kinetics, and characterisation of by-products. Sugar Tech 23: 433–452.https://doi.org/10.1007/s12355-020-00865-4 doi: 10.1007/s12355-020-00865-4
![]() |
[109] |
Kanwal S, Chaudhry N, Munir S, et al. (2019) Effect of torrefaction conditions on the physicochemical characterization of agricultural waste (sugarcane bagasse). Waste Manage 88: 280–290.https://doi.org/10.1016/j.wasman.2019.03.053 doi: 10.1016/j.wasman.2019.03.053
![]() |
[110] |
Chen WH, Ye SC, Sheen HK (2012) Hydrothermal carbonization of sugarcane bagasse via wet torrefaction in association with microwave heating. Bioresour Technol 118: 195–203.https://doi.org/10.1016/j.biortech.2012.04.101 doi: 10.1016/j.biortech.2012.04.101
![]() |
[111] |
Iryani DA, Kumagai S, Nonaka M, et al. (2017) Characterization and production of solid biofuel from sugarcane bagasse by hydrothermal carbonization. Waste Biomass Valor 8: 1941–1951.https://doi.org/10.1007/s12649-017-9898-9 doi: 10.1007/s12649-017-9898-9
![]() |
[112] | Beta San Miguel Beta San Miguel, Responsabilidad BSM. Available from: https://www.bsm.com.mx/resp_ambiental.html. |
[113] | Oleopalma (2022) 2022 Sustainable ability report. |
[114] | American Society of Testing Methods (2002) ASTM D-2974, standard test methods for moisture, ash, and organic matter of peat and other organic soils, 2002. Available from: https://www.astm.org/d2974-14.html. |
[115] |
Sheng C, Azevedo JLT (2005) Estimating the higher heating value of biomass fuels from basic analysis data. Biomass Bioenergy 28: 499–507.https://doi.org/10.1016/j.biombioe.2004.11.008 doi: 10.1016/j.biombioe.2004.11.008
![]() |
[116] |
Friedl A, Padouvas E, Rotter H, et al. (2005) Prediction of heating values of biomass fuel from elemental composition. Anal Chim Acta 544: 191–198.https://doi.org/10.1016/j.aca.2005.01.041 doi: 10.1016/j.aca.2005.01.041
![]() |
[117] |
Wang C, Deng X, Xiang W, et al. (2020) Calorific value variations in each component and biomass-based energy accumulation of red-heart Chinese fir plantations at different ages. Biomass Bioenergy 134: 105467.https://doi.org/10.1016/j.biombioe.2020.105467 doi: 10.1016/j.biombioe.2020.105467
![]() |
[118] |
Yin CY (2011) Prediction of higher heating values of biomass from proximate and ultimate analyses. Fuel 90: 1128–1132.https://doi.org/10.1016/j.fuel.2010.11.031 doi: 10.1016/j.fuel.2010.11.031
![]() |
[119] | Environmental Protection Agency USA (1996) EPA, METHOD 6010B, Inductively coupled plasma-atomic emission spectrometry, 1996. |
[120] |
Husain Z, Zainac Z, Abdullah Z (2002) Briquetting of palm fibre and shell from the processing of palm nuts to palm oil. Biomass Bioenergy 22: 505–509.https://doi.org/10.1016/S0961-9534(02)00022-3 doi: 10.1016/S0961-9534(02)00022-3
![]() |
[121] |
Sosa JA, Laines JR, Guerrero D, et al. (2022) Bioenergetic valorization of Sargassum fluitans in the Mexican Caribbean: The determination of the calorific value and washing mechanism. AIMS Energy 10: 45–63.https://doi.org/10.3934/energy.2022003 doi: 10.3934/energy.2022003
![]() |
[122] |
Demirbas A (2002) Relationships between heating value and lignin, moisture, ash and extractive contents of biomass fuels. Energy Explor Exploit 20: 105–111.https://doi.org/10.1260/014459802760170420 doi: 10.1260/014459802760170420
![]() |
[123] |
López MI, Sosa JA, Laines JR, et al. (2023) Aerobic biotransformation of Sargassum fluitans in combination with sheep manure: optimization of control variables. Chem Ecol 39: 823–842.https://doi.org/10.1080/02757540.2023.2263427 doi: 10.1080/02757540.2023.2263427
![]() |
[124] | Vergara GR (2022) Physical and energetic characterization of the fibrous residue from the processing of African palm by varying the percentage of humidity for bioenergy use in the Quevepalma company. Available from: http://repositorio.espe.edu.ec/handle/21000/28641. |
[125] | Huaraca JN (2022) Determination of the calorific value of agroindustrial waste from wheat husk and barley straw enriched with cellulose nanoparticles as an energy alternative. Available from: http://dspace.espoch.edu.ec/handle/123456789/20236. |
1. | Karan Thakkar, Julia Regazzini Spinardi, Jingyan Yang, Moe H. Kyaw, Egemen Ozbilgili, Carlos Fernando Mendoza, Helen May Lin Oh, Impact of vaccination and non-pharmacological interventions on COVID-19: a review of simulation modeling studies in Asia, 2023, 11, 2296-2565, 10.3389/fpubh.2023.1252719 |
Parameter | Brief Description | Value | Units | Reference |
P0 | 7.8×107 | Population of Tamil Nadu State, India | individual | [40] |
μ1 | 0.2053 | Intrinsic infection rate | Day−1 | fitted with data [38] |
β1 | 0.0484 | Unlock-down factor | Dimensionless | fitted with data [38] |
M1 | 0 | Pre-experience in σ | Dimensionless | fitted with data [38] |
t1 | 29 | Time parameter in σ | Day | fitted with data [38] |
k1 | 0.1935 | Steepness of σ | Day−1 | fitted with data [38] |
L1 | 0.4533 | Maximum increase in σ | Dimensionless | fitted with data [38] |
M2 | 0 | Pre-experience in α | Dimensionless | fitted with data [38] |
t2 | 105 | Time parameter in α | Day | fitted with data [38] |
k2 | 0.0248 | Steepness of α | Day−1 | fitted with data [38] |
L2 | 0.6405 | Maximum increase in α | Dimensionless | fitted with data |
τ1 | 4 | Delay in detection | Day | [23] |
τ2 | 14 | Recovery time | Day | [30] |
Parameter | μ1 | M1 | t1 | k1 | L1 |
Mean | 0.2064 | 0 | 29.3 | 0.1786 | 0.4565 |
Std | 0.004 | 0 | 0.9 | 0.0233 | 0.0081 |
Mean ± 2Std | [0.1984, 0.2144] | [0, 0] | [27.5, 31.1] | [0.1320, 0.2252] | [0.4403, 0.4727] |
Parameter | β1 | M2 | t2 | k2 | L2 |
Mean | 0.0485 | 0 | 105.7 | 0.0249 | 0.6395 |
Std | 0.0017 | 0 | 2 | 0.0004 | 0.0076 |
Mean ± 2Std | [0.0451, 0.0519] | [0, 0] | [101.7, 109.7] | [0.0241, 0.0257] | [0.6243, 0.6547] |
Parameter | Value | Brief Description | Units | Reference |
P0 | 7.8×107 | Population of Tamil Nadu State, India | individual | [40] |
μ1 | 0.227 | Intrinsic infection rate | Day−1 | fitted with data [38,39] |
μ2 | 0.0454 | Infection rate after vaccinated | Day−1 | assumed, based on [44,45] |
β1 | 0.0484 | Unlock-down factor | Dimensionless | fitted with data [38] |
h1 | 104 | HSC of vaccination | Individual | assumed |
M1 | 0.43 | Pre-experience in σ | Dimensionless | fitted with data [38,39] |
t1 | 83.9 | Time parameter in σ | Day | fitted with data [38,39] |
k1 | 0.2421 | Steepness of σ | Day−1 | fitted with data [38,39] |
L1 | 0.265 | Maximum increase in σ | Dimensionless | fitted with data [38,39] |
h2 | 5×105 | HSC of HCR | Individual | assumed |
M2 | 0.3491 | Pre-experience in α | Dimensionless | fitted with data [38,39] |
t2 | -5.8882 | Time parameter in α | Day | fitted with data [38,39] |
k2 | -0.6178 | Steepness of α | Day−1 | fitted with data [38,39] |
L2 | 0.5406 | Maximum increase in α | Dimensionless | fitted with data [38,39] |
M3 | 0.0216 | Pre-experience in γ | Dimensionless | fitted with data [38,39] |
t3 | 282.2 | Time parameter in γ | Day | fitted with data [38,39] |
k3 | 0.0163 | Steepness of γ | Day−1 | fitted with data [38,39] |
L3 | 0.7209 | Maximum increase in γ | Dimensionless | fitted with data [38,39] |
β2 | 3155833 | Maximum vaccination rate | Individual Day−1 | fitted with data [38,39] |
τ1 | 4 | Delay in detection | Day | [23] |
τ2 | 14 | Recovery time | Day | [30] |
τ3 | 14 | Hospitalization time | Day | [30] |
Parameter | μ1 | M1 | t1 | k1 | L1 |
Mean | 0.2339 | 0.4304 | 82.2 | 0.2363 | 0.2658 |
Std | 0.01 | 0.0011 | 2.6 | 0.0391 | 0.0082 |
Mean ± 2Std | [0.2139, 0.2539] | [0.4282, 0.4326] | [77, 87.4] | [0.1581, 0.3145] | [0.2494, 0.2822] |
Parameter | M2 | t2 | k2 | L2 | |
Mean | 0.3574 | -5.53 | -0.6341 | 0.5625 | |
Std | 0.0219 | 0.52 | 0.0803 | 0.0318 | |
Mean ± 2Std | [0.3136, 0.4012] | [-6.57, -4.49] | [-0.7947, -0.4735] | [0.4989, 0.6261] | |
Parameter | M3 | t3 | k3 | L3 | β2 |
Mean | 0.0221 | 282 | 0.0165 | 0.7264 | 3.1757×106 |
Std | 0.0015 | 7 | 0.0007 | 0.0528 | 1.872×105 |
Mean ± 2Std | [0.0191, 0.0251] | [268,296] | [0.0151, 0.0179] | [0.6208, 0.832] | [2.8013,3.5501]×106 |
Parameter | Brief Description | Value | Units | Reference |
P0 | 7.8×107 | Population of Tamil Nadu State, India | individual | [40] |
μ1 | 0.2053 | Intrinsic infection rate | Day−1 | fitted with data [38] |
β1 | 0.0484 | Unlock-down factor | Dimensionless | fitted with data [38] |
M1 | 0 | Pre-experience in σ | Dimensionless | fitted with data [38] |
t1 | 29 | Time parameter in σ | Day | fitted with data [38] |
k1 | 0.1935 | Steepness of σ | Day−1 | fitted with data [38] |
L1 | 0.4533 | Maximum increase in σ | Dimensionless | fitted with data [38] |
M2 | 0 | Pre-experience in α | Dimensionless | fitted with data [38] |
t2 | 105 | Time parameter in α | Day | fitted with data [38] |
k2 | 0.0248 | Steepness of α | Day−1 | fitted with data [38] |
L2 | 0.6405 | Maximum increase in α | Dimensionless | fitted with data |
τ1 | 4 | Delay in detection | Day | [23] |
τ2 | 14 | Recovery time | Day | [30] |
Parameter | μ1 | M1 | t1 | k1 | L1 |
Mean | 0.2064 | 0 | 29.3 | 0.1786 | 0.4565 |
Std | 0.004 | 0 | 0.9 | 0.0233 | 0.0081 |
Mean ± 2Std | [0.1984, 0.2144] | [0, 0] | [27.5, 31.1] | [0.1320, 0.2252] | [0.4403, 0.4727] |
Parameter | β1 | M2 | t2 | k2 | L2 |
Mean | 0.0485 | 0 | 105.7 | 0.0249 | 0.6395 |
Std | 0.0017 | 0 | 2 | 0.0004 | 0.0076 |
Mean ± 2Std | [0.0451, 0.0519] | [0, 0] | [101.7, 109.7] | [0.0241, 0.0257] | [0.6243, 0.6547] |
Parameter | Value | Brief Description | Units | Reference |
P0 | 7.8×107 | Population of Tamil Nadu State, India | individual | [40] |
μ1 | 0.227 | Intrinsic infection rate | Day−1 | fitted with data [38,39] |
μ2 | 0.0454 | Infection rate after vaccinated | Day−1 | assumed, based on [44,45] |
β1 | 0.0484 | Unlock-down factor | Dimensionless | fitted with data [38] |
h1 | 104 | HSC of vaccination | Individual | assumed |
M1 | 0.43 | Pre-experience in σ | Dimensionless | fitted with data [38,39] |
t1 | 83.9 | Time parameter in σ | Day | fitted with data [38,39] |
k1 | 0.2421 | Steepness of σ | Day−1 | fitted with data [38,39] |
L1 | 0.265 | Maximum increase in σ | Dimensionless | fitted with data [38,39] |
h2 | 5×105 | HSC of HCR | Individual | assumed |
M2 | 0.3491 | Pre-experience in α | Dimensionless | fitted with data [38,39] |
t2 | -5.8882 | Time parameter in α | Day | fitted with data [38,39] |
k2 | -0.6178 | Steepness of α | Day−1 | fitted with data [38,39] |
L2 | 0.5406 | Maximum increase in α | Dimensionless | fitted with data [38,39] |
M3 | 0.0216 | Pre-experience in γ | Dimensionless | fitted with data [38,39] |
t3 | 282.2 | Time parameter in γ | Day | fitted with data [38,39] |
k3 | 0.0163 | Steepness of γ | Day−1 | fitted with data [38,39] |
L3 | 0.7209 | Maximum increase in γ | Dimensionless | fitted with data [38,39] |
β2 | 3155833 | Maximum vaccination rate | Individual Day−1 | fitted with data [38,39] |
τ1 | 4 | Delay in detection | Day | [23] |
τ2 | 14 | Recovery time | Day | [30] |
τ3 | 14 | Hospitalization time | Day | [30] |
Parameter | μ1 | M1 | t1 | k1 | L1 |
Mean | 0.2339 | 0.4304 | 82.2 | 0.2363 | 0.2658 |
Std | 0.01 | 0.0011 | 2.6 | 0.0391 | 0.0082 |
Mean ± 2Std | [0.2139, 0.2539] | [0.4282, 0.4326] | [77, 87.4] | [0.1581, 0.3145] | [0.2494, 0.2822] |
Parameter | M2 | t2 | k2 | L2 | |
Mean | 0.3574 | -5.53 | -0.6341 | 0.5625 | |
Std | 0.0219 | 0.52 | 0.0803 | 0.0318 | |
Mean ± 2Std | [0.3136, 0.4012] | [-6.57, -4.49] | [-0.7947, -0.4735] | [0.4989, 0.6261] | |
Parameter | M3 | t3 | k3 | L3 | β2 |
Mean | 0.0221 | 282 | 0.0165 | 0.7264 | 3.1757×106 |
Std | 0.0015 | 7 | 0.0007 | 0.0528 | 1.872×105 |
Mean ± 2Std | [0.0191, 0.0251] | [268,296] | [0.0151, 0.0179] | [0.6208, 0.832] | [2.8013,3.5501]×106 |