Co-Combustion Characteristics of Typical Biomass and Coal Blends by Thermogravimetric Analysis

In this study, the co-combustion characteristics of coal and biomass blends (20, 40, 60, 80, and 100 wt%) were investigated by thermogravimetric analysis. All the samples were operated under an oxidative atmosphere, with a heating rate of 20 C/min. The reaction stages, ignition and burnout temperature, maximum weight loss rate, and different combustion indices were determined. When the percentage of biomass in the blends was increased, the maximum mass loss rate was enhanced in the second region, and the ignition and burnout temperature was lowered, indicating the higher reactivity and better combustion performance of the samples. The comprehensive performance index presented an N shape with the increasing biomass blending ratio. Based on various combustion indices, 20% was an optimum percentage for the co-utilization of coal-biomass blends. A significant promoting interaction was observed between corn straw and rice straw blends, while inhibiting effects occurred between rice husk and coal. The kinetic parameters of the blends were evaluated by the Coats and Redfern method using the nth-order reaction model. The value of activation energy and the pre-exponential factor increased with the decreasing biomass percentage in the blends.


INTRODUCTION
In recent years, the rapid economic development has caused increasing energy consumption, but the limited amount of nonrenewable energy resources and the depletion in the future pose a profound challenge to energy demand. In addition, the combustion of fossil fuels aggravates environmental pollution because of the emissions of nitrogen and sulfur pollutants as well as CO 2 . Thus, the utilization of renewable energy has attracted considerable critical attention. Among these renewable energies, biomass may play an essential role due to its advantages of fuel flexibility, high combustion efficiency, low pollutant emission, and carbon neutrality (Jayaraman et al., 2017).
Numerous technologies have been extensively investigated for potential utilization of biomass, including combustion, pyrolysis, gasification, and liquefaction (Gil et al., 2010). Co-combustion is considered the state-of-the-art technology to utilize biomass for replacing fossil fuels, such as coal, to meet the stationary energy demand. Several reasons account for mixing biomass with coal before burning. The co-combustion of coal and biomass blends would directly help to cut down the consumption of fossil fuels. At the same time, a biofuel product is added to coal to obtain better burning performance during the combustion reaction (Wang et al., 2009). During the co-combustion process, a higher volatile matter content (normally more than 35%) is considered to provide a stable flame (Sahu et al., 2014), attained by the addition of biomass. The slagging and fouling problems of heating surfaces commonly faced in biomass combustion could be mitigated or eliminated through co-combusting of coal and biomass (Haykiri-Acma and Yaman, 2008). Additionally, there is no need for existing power plants to undergo many modifications for continuous use. In conclusion, the co-utilization of biomass and coal blends in present coal-fired power plants possibly leads to a mass of benefits in the environment, technology, and economy (Kastanaki and Vamvuka, 2006).
Knowledge of combustion characteristics and kinetic mechanisms is crucial for comprehension and modeling in industrial-scale furnaces, co-firing of biomass and coal, and burning of coal alone (Hu et al., 2019). Thermogravimetric (TG) analysis is commonly used to obtain a rapid investigation and comparison in thermal characteristics and kinetic parameters during the combustion of various fuels, including coal and biomass. Many researchers (Gil et al., 2010;Yanfen and Xiaoqian, 2010;Su et al., 2013;Lin et al., 2015;Lu and Chen, 2015;Jayaraman et al., 2017;Chen et al., 2018;Hu et al., 2019;Wang et al., 2019;Wang et al., 2020) have studied the cocombustion behaviors of different biomass and coal/biomass blends. A significant decrease in ignition and burnout temperature was observed by Wang et al. (2019) and Konwar et al. (2019). According to Chen et al. (2018) and Liu et al. (2015), the combustion indices, such as maximum mass loss, ignition performance, and the comprehensive combustion index, were evaluated. Guo et al. (2020) and Shi et al. (2019) found the interaction between coal and biomass, including promoting synergy and inhibiting effects. The combustion kinetic parameters, including activation energy and the preexponential factor, using the Coats and Redfern method (Gil et al., 2010;Wang et al., 2012;Lin et al., 2015;Jayaraman et al., 2017;Chen et al., 2018) and the Flynn-Wall-Ozawa method (Xie and Ma, 2013;Liu et al., 2015;Konwar et al., 2019;Wang et al., 2019;Guo et al., 2020) were determined.
Crop straws and rice husk are the most fundamental agricultural residues in quantitative terms (Sathitruangsak and Madhiyanon, 2017). As the world's largest grain manufacturer, China owns plentiful biomass resources, and the amount of collectable straw was approximately 674 million tons in 2017 . Also, China is the leading country to cultivate rice, producing the most rice paddy (Hossain et al., 2018). Compared with other biomass, straw and rice husk char contain more high alkali and alkaline earth metals, which could catalyze the co-combustion process with coal and enhance combustion performance. Therefore, it was necessary to study the combustion and co-combustion characteristics and kinetic mechanisms to provide a theoretical basis and an optimal blending ratio of coal and biomass fuel.
In this research, co-combustion characteristics and kinetic parameters of mixed coal and three biomass samples (corn straw, rice straw, and rice husk) and their blends are examined under an oxidative atmosphere with TG equipment. The main objective of this study is to analyze the thermal behaviors and kinetic mechanisms of coal and biomass blends. This study also provides data on blending fuels in their combustion applicability. Finally, these results can thoroughly control the coal/biomass characteristics according to the demand.

Materials
The corn straw (CS), rice straw (RS), and rice husk (RH) used in this study were collected from Jilin Province, China. Coal was obtained from Inner Mongolia, which is a mixture of brown and bituminous coal. The samples were milled into powder with a particle size of 75-150 μm. Then, different blend ratios were prepared with biomass of 20, 40, 60, and 80 wt% and named, for example, 20% CS, 40% RS, and 60% RH, depending on the proportion of the biomass in the blends. The prepared mixtures of coal and biomass were mixed manually into a homogeneous distribution. The photographs of raw, ground, and mixed samples are shown in Figure 1, and the ultimate and proximate analyses of biomass and coal samples are presented in Table 1.

Experimental Setup and Methods
Non-isothermal experiments were carried out employing a NETZSCH analyzer (model STA 449 F3) to investigate cocombustion behaviors and kinetic parameters of the samples. The sample mass for each test was approximately 5 mg. The experiment temperature was set from ambient temperature (20 C) to 1,000 C with a 20 C/min heating rate and a flow rate of 60 ml/min under air conditions. All experiments were reduplicated to ensure the reproducibility and consistency of the results. The weight loss and derivative thermogravimetry (DTG) of the samples were analyzed as a function of the temperature.

Combustion Parameters
The characteristic parameters of the combustion process were obtained from TG and DTG profiles. The ignition temperature (T i ) was determined using the intersection method (IM) (Lu and Chen, 2015), and the burnout temperature (T b ) was defined when the conversion reached 98% (Wang et al., 2019). In order to have a comprehensive analysis of the combustion characteristics of the samples, the flammability index, C, and the comprehensive combustion index, S, were introduced to compare the cocombustion performance and are shown as follows (Hu et al., 2019): where DTG max and DTG mean are the maximum and average mass loss rates, respectively. The higher index S represents the more vigorous burning of the samples and the faster burnout of char.

Kinetic Analysis
The reaction of coal, biomass, and their blends is generally considered as a heterogeneous solid-phase reaction (Gil et al., 2010;Sahu et al., 2014), which follows the Arrhenius law, and can be described as where α m 0 − m t m 0 − m ∞ is the conversion degree or mass conversion ratio obtained from TG curves, of which m 0 and m ∞ are the initial mass and final masses of the samples, respectively, and m t represents the sample mass at time t. From the Arrhenius equation, the reaction rate constant can be expressed as where A, R, and E are the pre-exponential factor, universal gas constant, and activation energy, respectively. For nonisothermal reactions, the heating rate remains constant (β dT dt ). Therefore, the following could be inferred combining the equation: In this study, the Coats and Redfern (CR) method (Coats and Redfern, 1964) is used to determine the activation energy and reaction mechanism, which is widely used to estimate kinetic parameters for biomass combustion. The hypothetical model of the reaction is based on the nth reaction model Chen et al., 2018); therefore, f(α) (1 − α) n . The detailed form is shown as In For the basis of F RT ≫ 1, 1 − 2RT E ≈ 1, and the final form takes the following equation: In By plotting In − ln(1−α) T 2 and vs. 1/T, E could be obtained from the slope of a straight line. Different reaction order values were chosen to determine the best fitting ones.

Fuel Properties of a Single Sample
The thermochemical behaviors demonstrated by TG analysis could offer vital insights into the combustibility and transforming patterns of fuels (Luo et al., 2014). As shown in Figure 2, the ignition temperature for coal, CS, RS, and RH was 401°C, 241°C, 256°C, and 270°C, respectively, and the burnout temperature for the four samples was 557°C, 486°C, 481°C, and 479°C, respectively. These results indicated that the biomass fuels had a much lower ignition temperature and were easier to achieve burnout due to the higher volatile content in biomass, resulting in better thermal reactivity. Besides, the higher carbon content and the lower oxygen content in coal also contributed to the higher decomposition temperature as the chemical bond energy for breaking the C-C bond was higher than that for the C-O and C-H bonds (Konwar et al., 2019). Therefore, the destruction of the macromolecular structure and chemical bonds delayed the burnout of coal, and the reaction occurred at a higher temperature. The proximate analysis of coal and biomass in Table 1 also confirmed all these results shown above. The temperature ranges of the main stages, together with the weight loss, are shown in Table 2. It could be inferred that the combustion process of three biomass samples was partitioned into three stages (Jayaraman et al., 2017;Wang et al., 2019): the dehydration process (stage A), the release and combustion of the volatiles (stage B), and the burning and burnout of residues and fixed carbon (stage C). The demarcation of each stage was clarified by the inflection point in the DTG curves (Lin et al., 2015). In the case of the coal sample, there was only one stage during its combustion process after the moisture evaporation, which was the primary carbonizing stage accompanying the released carbon dioxide and hydrogen (Jayaraman et al., 2017). A single weight-loss peak with a wider temperature range was observed, corresponding to a longer overlapping reaction process with the combustion of the volatiles and char oxidation (Vamvuka and Sfakiotakis, 2011).
In all samples, the initial stage was observed from indoor temperature to 120°C. In this stage, the mass loss depended on the moisture content of the fuel properties, and a small DTG peak occurred in the DTG curve. The second stage for CS, RS, and RH was extended from 116 to 408°C, from 119 to 401°C, and from 156 to 372°C, and the mass loss was 69.6, 61.35, and 53.28%, respectively. It was seen that the weight loss during this stage was almost close to the volatile content in Table 1, accounting for  the majority of the total weight mass loss. This result confirmed that the predominant form of combustion in the second stage is via volatile releasing and burning, and a higher DTG peak was observed at around 300 C for three samples. It was attributed to the decomposition and combustion of celluloses (hemicellulose and cellulose), and a small account of lignin also combusted in this stage. The hemicellulose consists of a low degree of the amorphous structure, resulting in bad thermal stability, and the decomposition temperature extended from 220 to 315°C. However, the macromolecular structure of the cellulose is composed of glucose without branching, leading to a higher decomposition temperature (315-390 C). However, the thermal behavior of lignin occupied the whole decomposition process with a wider temperature range. This highly cross-linked polymer contained various chemical bonds, such as C-C bonds and ether bonds, causing a wide decomposition temperature (160-627 C). The third stage occupied the remaining small part of the total loss with a higher temperature range, as shown in Table 2. During this stage, the DTG peak for three biomass samples appeared at around 420°C, and CS showed a sharp and narrow peak, while RH and RS presented a smaller and smoother peak. A deduction was that the decomposition process Frontiers in Energy Research | www.frontiersin.org October 2021 | Volume 9 | Article 753622 5 changed the physicochemical properties of the char produced after the burnout of volatiles, and the combustion performance in the char combustion process significantly varied. Liu et al. (2015) pointed out that the formation or the structure of char differed greatly during the carbonization process, directly affecting the DTG peaks.

Co-Combustion Properties of the Biomass and Coal Blends
The mass loss and DTG profiles of coal and biomass blends are presented in Figure 3. The combustion process of the blends also contained three stages: moisture evaporation, the decomposition and combustion of volatile matters, and the oxidation of char. It can be seen that for RS and RH, the DTG curves of their blends were located between individual fuels, but for CS, there existed obvious overlapping in its DTG curve, corresponding to the evident decline in the TG curve at approximately 430°C for CS and 80% CS samples. It was easy to interpret that the intensity of the third DTG peak for CS varied distinctly with RS and RH owing to the higher fixed carbon content and the significantly lower ash content within CS in Table 1, generating a higher amount of combustible char, that is, more flammable substances would be produced and combusted in Stage C. Besides, the higher carbon and lower oxygen content required more oxygen during combustion, so the combustion of certain matters in Stage B was delayed, contributing to a higher DTG value in Stage C. When 20% of coal was added to the biomass sample, the combustion performance of the mixed sample improved. This was owing to the fact that the large amount of heat generated during the combustion of biomass volatiles would facilitate the combustion of coal and biomass char. While more coal was mixed with CS, the combustion of char needed more heat, but the combustion of biomass char could not provide sufficient heat; thus, the third DTG peak moved forward to a higher-temperature range. In a similar case in corn and wheat straw, Li et al. (2014) and Wang et al. (2019) identified the same phenomenon.
To compare combustion characteristics and quantitatively analyze, various combustion parameters are shown in Table 3. The main observation was that when the percentage of biomass in the blends increased, the ignition temperature immediately decreased, but a relatively slow reduction in the burnout temperature was found. It was noteworthy that the ignition temperature of the blends was extremely close to that of individual biomass samples, revealing that the ignition of two fuels occurred independently, and biomass played a predominant role at low temperatures. A possible explanation for this might be that many volatile matters were obtained and rapidly ignited. The combustion of the volatiles was an exothermic reaction, giving off a large amount of heat and increasing the surface temperature (Liu et al., 2015) and thereby the reaction rate of blends. In contrast, when the percentage of biomass was within 60%, the burnout temperature for blends was closer to that of the coal sample. This result suggested that the coal dominated the burnout process for the blends, which was probably related to the higher carbon content in coal, requiring a high temperature to burn out.
Another important index obtained for the curves was the maximum DTG value; DTG max was deemed to be proportionate to sample reactivity. DTG max1 and DTG max2 represented the maximum loss rates for Stages B and C, respectively. In the blends, DTG max1 increased with the percentage of biomass contained in the blends, which suggested that the higher the amount of biomass in the mixtures, the faster the rate of mass loss for the second stage or in other words, the higher the reactivity of the sample in this stage. This result may be explained by the fact that the biomass had a lower ignition temperature and almost all volatiles generated in Stage B would be combusted, contributing to the higher DTG value for the first mass loss peak. In stage C, this was not the case due to the addition of coal, and DTG max2 increased with the rising percentage of coal in the blends, suggesting that coal could improve the reactivity of the samples in the later stage of combustion. It was easy to understand that the higher the ratio of coal, the higher the fixed carbon in the blends, and the combustion of flammable substances including char enhanced the sample's reactivity. Interestingly, there was a turning point for blends where DTG max2 was bigger than DTG max1 when 40% of coal was mixed with biomass. This observation may support the hypothesis that biomass controlled the second combustion stage (Stage B), while coal dominated the third stage (Stage C).
Finally, to achieve a comprehensive evaluation of ignition and combustion performance of coal and biomass blends, two combustion indices as a function of biomass blending ratio are shown in Figure 4. The trend of S and C increased first, then declined, and rose again when the percentage of biomass exceeded 80%, overall exhibiting an "N" curve shape. What stands out in the figure is the marked value of C and S for 20% biomass mixed in blends, implying the improvement of combustion performance by the addition of biomass. Considering the stability of boiler operation and the limitation of the feeding system, the amount of biomass in the blends is suggested to be kept at around 20%, according to the results in this study. Previous research (Zhou et al., 2016) had confirmed that the blending proportion of 20% was regarded as the optimum blend for the co-combustion of corn stalk and bituminous coal. This conclusion was consistent with the review by Sahu et al. (2014), which reported that a biomass content of 20% favored the best additive effects.

Synergistic Effect Analysis
To identify whether interactions between coal and biomass happened during the co-combustion process, the calculated TG curves of the blends based on the experimental curves of individual fuels were compared with experimental curves. Assuming that there were no interactions between the two components, the calculated TG curves could be expressed as follows (Wang et al., 2019;Guo et al., 2020): TG cal x bio TG bio + x coaI TG coal , where TG bio and TG coal are the mass losses of individual fuels and x bio and x coaI are the percentages of biomass and coal, respectively.
The calculated TG curves of three blends are shown in Figure 5. To further elucidate the differences between the two curves of the mixed samples, the deviation was introduced as follows (Wang et al., 2019): where Δw was defined to explicitly demonstrate the interaction between biomass and coal during the co-fire process. Figure 6 exhibits the deviation value between TG cal and TG exp , and Δw > 0 and <0 represented the synergetic and suppressive interactions, respectively. It could be found that calculated and experimental curves were similar in the early stage, but deviation occurred with the proceeding co-combustion process. This implied that synergetic interactions existed between coal and biomass. As shown in Figure 5 and Figure 6, CS and RS exhibited similar patterns for calculated and experimental curves. From the initial stage to around 250°C (the ignition temperature of biomass), good agreement was observed between TG cal and TG exp for the biomass ratios of 20, 40, and 60% in the blends, which suggested that no interactions happened. With the increase of conversion, when the temperature rose to 300°C, the Δw curve slightly dropped for all blends.
This may be attributed to the thermal resistance effect, in which the biomass component of the blends softened, flowed, and then adhered to the burning surface. Besides, a higher temperature was needed for the devolatilization of the high coalification degree of coal, and the heat released by biomass combustion after ignition was used for the endothermic reaction of the coal's ignition process. However, when the temperature was over 350°C, it can be seen from Figure 6 that TG exp markedly excessed TG cal , which indicated that the inhibition turned into a promoting effect. The Δw reached a peak value at around 450-500°C, which was the burnout temperature of individual biomass. It could be explained by the fact that the ignition temperature of biomass char was low, and vast heat generated after combustion heated the fixed carbon in coal. At the same time, the alkaline matters in biomass ash, including oxides of potassium, sodium, and magnesium, catalyzed the co-combustion of the blends and then accelerated the reaction rate. This assumption was also confirmed by other researchers Wang et al., 2019;Guo et al., 2020), who found the promoting reactions between alkaline and alkaline earth metals in biomass ash and solid carbon in coal. It was noteworthy that the interaction between 60% biomass and 40% coal lagged behind other group blends; for CS, the interaction functioned as inhibition. This phenomenon may be driven by the fact that the devolatilization of chemical constituents from biomass blocked the pores of coal (Konwar et al., 2019), which could weaken heat-transfer efficiency and oxygen diffusion, therefore hindering the combustion of fixed carbon, inhibiting the devolatilization process and exerting resistance on gas diffusion  October 2021 | Volume 9 | Article 753622 8 . Besides, the generated fly ash would hinder the diffusion of the flame during 40% coal char combustion, thereby exhibiting a distinct flame retardant effect. Under such a blending ratio, the catalytic effects of biomass char were dramatically restricted by the blocking pores, while for 80% biomass samples, the reaction of biomass in different stages dominated the whole process due to the higher biomass blending ratio. Cong et al. (2020) confirmed this result in evaluating the performance of semi-coke and corn stover, and such a trend was also observed.
As for RH, which is shown in Figure 6C, the calculated and experimental TG profiles almost coincided for 20 and 60% RH blends. It was noted that for 40% RH, the calculated values were higher than experimental values and Δw > 0, which represented the synergetic interaction between two components in the sample, but inhibiting interactions were detected for 80% RH since an opposite trend was found between the two values. These results showed that the interactions in the blends switched from promotion to suppression, with the biomass ratio increased from 20 to 80%. A probable reason for the change was the highest content of SiO 2 in RH ash compared with straw biomass. Different from the catalytic effects of alkali matters in biomass ash, SiO 2 and Al 2 O 3 could hinder the reaction process. According to earlier studies, nearly 80% of RH ash was SiO 2 and this significantly inhibited the combustion process. This assumption corresponded to Wang et al. (2016) and Zhou et al. (2014), who demonstrated that the existence of SiO 2 in rice husk ash had obvious inhibiting effects during the co-combustion process; the experimental data were lower than calculated values in the temperature range of 575-800 K.

Kinetic Analysis
The kinetic parameters were calculated on the assumption that the reaction was separated into a single stage during the thermal conversion (Gil et al., 2010). Thus, two independent reactions were necessary for coal/biomass blends to describe the co-combustion process, and the conversion in each stage was recalculated for each reaction. In this study, several reaction order values (n 0.3, 0.6, 0.9, 1, 1.2, 1.5, 2, and 3) were applied to the TG data; the n value that gave the highest correlation coefficient represents the kinetic mechanisms for each reaction. Table 4 shows that in this research, reaction orders of 0.9-1.2 and 1-1.5 were the best fitting models for Stage B and Stage C, respectively. Similarly, Chen et al. (2018) also reported that 1.2 and 1.5 were the most suitable model functions for the two stages during the co-combustion of coal/biomass blends. Figure 7 shows the curve fitting of CR plots for all samples, which had the best linear results. Most of the resulting correlation coefficients (R2) were higher than 0.99 ( Table 4), indicating that the calculation profiles were reliable. A decreasing order of activation energy was found in Stage B, followed by RH > RS > CS, and the E values for the three biomasses were 83.31 kJ/mol, 79.42 kJ/mol, and 51.89 kJ/mol, respectively. As the reaction progressed, what stood out in Stage C was the sharp increase for the E values, which were 292.88 kJ/mol, 176.11 kJ/mol, and 125.72 kJ/mol, with the opposite order as RH < RS < CS. The E values were considered to be the minimum energies required for the combustion reactions, with the high value representing that the reaction started slowly. These results meant that CS was easier to ignite in the early stage, but more energy was needed to break down the barriers in the later combustion process. This was due to the difference in the constitution of the three biomasses; the decomposition and burning of volatile matter and fixed carbon greatly influenced the combustion process. A single or overlapping process of manifold reaction mechanisms, including diffusion and interface, also accounted for the different values (Huang et al., 2018). The difference in the E values between the two stages implied that a higher temperature was required for the char-burning stage, revealing that this reaction with higher activation energies was more temperature-dependent (Masnadi et al., 2014;Liu et al., 2015).
It was noted that in Stage B for all blends, the activation energy and pre-exponential factor values increased with increasing biomass ratio in the blends, while a significant opposite trend was observed in Stage C. Besides, when the biomass ratio was under 60%, there was no striking difference in E values between the coal sample and the blends in Stage C; for RH, the blending ratio even reached 80%. This finding also accorded with our earlier observations, which showed that coal dominated the oxidation of char in the later stage. The activation energy value of coal and biomass blends acquired in this study was comparable to those of previous studies. Zhou et al. (2016) reported that the activation energy of corn stalk/ bituminous coal blends was in the range of 64.5-94.8 kJ/mol. Xie and Ma (2013) demonstrated 97 kJ/mol of E value for rice straw using the Friedman method. Wang et al. (2020) specified an activation energy of 83.35 kJ/mol for rice husk with a second-order reaction mechanism.

CONCLUSIONS
In the air atmosphere, the combustion process of biomass and the blends could be portioned into two stages after the initial moisture evaporation compared to a single stage for the coal sample. Considering the pattern of characteristic combustion parameters (T i , T b , DTG max1 , DTG max2 , T max1 , T max2 , C, and S), as the biomass ratio increased, the combustion performance of the blends was enhanced. Based on the combustion characteristics, 20% biomass was the optimum percentage for the co-combustion process. The synergetic interaction of the co-combustion process for CS and RS blends was observed due to the catalytic effect of the biomass ash. However, an inhibiting effect for RH and coal in a higher biomass ratio was found, attained by a large amount of SiO 2 content in ash, hindering the combustion of fixed carbon. The reaction orders of 1.2 and 1.5 were the best fitting results for Stages B and C and sought to be the reaction mechanisms for biomass and the blends. In Stage B, the activation energy values for three biomasses showed a decreasing order as RH > RS > CS, and the E values were 83.31 kJ/ mol, 79.42 kJ/mol, and 51.89 kJ/mol, respectively. However, as the reaction was processed, an opposite trend was reported. Also, the E and A values increased with the increasing biomass percentage in Stage B, while the oxidation of char dominated Stage C and presented a decreasing pattern.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary material; further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
YY drafted the work and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. YH revised it critically for important intellectual content. JT and YW helped the acquisition, analysis, or interpretation of data for the work. ZW and SK provided approval for publication of the content.