Abstract
Honey is a natural sweetener valued for its nutritional and potential health benefits. However, honey’s essential minerals and toxic metals can significantly impact its quality and safety. A total of 12 (n = 12) honey samples, including both local and imported, were assessed for concentrations of essential minerals (iron, zinc, and manganese) and toxic heavy metals (cadmium, chromium, and lead) using Atomic Absorption Spectrophotometry (AAS). The findings revealed significant variability in mineral content, with imported honey generally containing higher levels of iron and manganese. Zinc concentrations were similar between local and imported honey. Concerning toxic metals, cadmium levels in 4 local samples out of 6 and all imported (6) honey samples exceeded the maximum allowable limit. Health risk assessment using Target Hazard Quotient (THQ) and Hazard Index (HI) indicated no immediate non-carcinogenic risk (HI <1), though cadmium contributed the highest proportion to the HI, particularly in local honey. Chromium and lead levels were within safe limits for all samples. The study emphasizes the importance of regular monitoring and stringent quality control to ensure the safety of honey. Further research should investigate the sources of cadmium in local honey and develop mitigation strategies.
Introduction
Honey, a natural sweetener produced by bees from flower nectar, is globally valued for its unique flavor, nutritional properties, and potential health benefits. It is a complex mixture of sugars, enzymes, organic acids, vitamins, minerals, and other bioactive compounds.1,2 The mineral and metal content of honey is an essential aspect of its quality and safety. These elements are naturally present in honey and derive from the soil, water, and plants in the area where the bees forage. Consequently, beekeeping products, including honey, can serve as bioindicators of environmental contamination, reflecting the presence of heavy metals in the ecosystem.3 -5 Honey produced in different regions can exhibit distinct mineral profiles, reflecting the unique environmental characteristics of their production areas. 6
Conversely, the presence of toxic metals, such as lead (Pb), cadmium (Cd), and Chromium (Cr), can pose severe health hazards, including neurotoxicity, carcinogenicity, kidney damage, and other chronic conditions.7 -9 These toxic elements can enter honey through environmental contamination, primarily from industrial activities, vehicular emissions, and agricultural practices.8,10 A study conducted in Poland found variations in the levels of essential and toxic elements in honey from different regions, highlighting the impact of environmental factors on honey composition. 11 Similarly, research in Turkey revealed that honey samples contained trace amounts of heavy metals, with some exceeding permissible limits, indicating potential health risks. 12 These findings underscore the necessity of regular monitoring and stringent quality control measures to prevent contamination.
In Bangladesh, honey is a popular natural product, valued for its perceived medicinal and nutritional benefits. Dhaka, the capital of Bangladesh, is a rapidly growing megacity with a dense population and a significant demand for honey.13,14 Additionally, the Dhaka market is flooded with both locally produced and imported honey, raising concerns about the comparative quality and safety of these products.15,16 Previous studies on food contamination in Bangladesh have primarily focused on vegetables, fruits, grains, and animal products, with limited research on honey.17 -19 Moreover, some studies have investigated the mineral content of honey from different regions of Bangladesh. However, there is a lack of comprehensive data on the levels of essential minerals and heavy metals in local and imported honey from Dhaka, particularly with quantitative health risk assessment.20,21
Several analytical methods are used to determine the chemical composition of honey, primarily spectroscopic and spectrometric techniques. These include atomic absorption spectrometry (AAS), inductively coupled plasma mass spectrometry (ICP-MS), and inductively coupled plasma optical emission spectrometry (ICP-OES). 22 This low-temperature, wet-acid digestion-based AAS method ensures accurate detection and quantification of both beneficial minerals and potentially harmful heavy metals in honey. This study aims to evaluate the quality and safety of honey available in Dhaka by analyzing the mineral and metal content of both local and imported honey samples. By determining the concentrations of essential minerals (iron, zinc, and manganese) and potentially toxic metals (cadmium, chromium, and lead), and by conducting a quantitative health risk assessment, this research provides valuable insights into the safety of consuming honey in Dhaka. Additionally, the study highlights the need for regular monitoring and stringent quality control measures to ensure the safety of honey on the market.
Materials and Methods
Reagents
High-concentration (1000 mg/l) standards of Fe, Zn, Mn, Cd, Cr, and Pb were obtained from Scharlau Chemie, S.A. (Barcelona, Spain) through Labtex Bangladesh. Trace analytical-grade nitric acid (65% pure) and hydrogen peroxide (30% pure) were purchased from Merck Darmstadt, Germany, through Kuri & Company (Pvt.) Limited, Dhaka, Bangladesh. Ultrapure-grade distilled water was produced using the Biobase Ultrapure Water Purification System (BK-UP-20, China).
Sample Collection
A total of 12 (n = 12) honey samples were systematically collected between January and February 2024 for comprehensive analysis. The sampling strategy employed a stratified convenience approach to ensure representation of both local and imported honey categories available to Dhaka consumers.
Local honey samples: Six local honey samples were collected in duplicate (two 500 g jars per source) to ensure analytical reproducibility. Three unprocessed samples (LH1-LH3) were obtained directly from beekeepers affiliated with the Department of Entomology, Sher-e-Bangla Agricultural University. Three commercially processed local samples (LH4-LH6) were purchased from 3 different local brands available in the neighborhood markets of Mohammadpur, with 1 replicate jar collected per brand.
Imported honey samples: Six imported honey samples (IH1-IH6) were collected from 6 different international brands to capture product diversity in the markets of Mohammadpur, Dhaka. For each brand, 2 replicate jars (500 g each) were purchased from 2 separate retail outlets to account for batch variability.
This diverse collection facilitated comparisons between imported and locally produced honey, as well as between processed and unprocessed varieties, within Bangladesh.
Sample Preparation
Honey samples were processed using the existing method, with modifications. 23 Honey samples were placed in petri dishes and dried in a laboratory oven at 60°C to 70°C for 72 hours to achieve a constant dry weight. Approximately 10 g of each oven-dried honey sample was weighed and transferred into digestion tubes. A 15 ml diacid mixture (HNO3: H2O2 in a 2:1 ratio) was added to each tube. The tubes were placed in a hot water bath at around 90°C until the production of red NO2 fumes ceased. The remaining liquid evaporated to approximately 3 to 5 ml. After cooling, the tube contents were diluted with deionized water, filtered through Whatman No. 1 filter paper into a 100 ml volumetric flask, and then adjusted to volume with additional deionized water. Finally, the digested samples were transferred to clean plastic bottles for storage until further analysis by atomic absorption spectrophotometry (AAS). The chemical analysis of these samples was conducted at the food safety laboratory of the Department of Agricultural Chemistry at Sher-e-Bangla Agricultural University, Dhaka-1207, Bangladesh. This low-temperature, wet-acid digestion method was chosen for its simplicity, effectiveness in decomposing the high sugar content of honey, and suitability for the available laboratory infrastructure.
Instrumental Analysis
An Analytik Jena NovAA 400P Atomic Absorption Spectrophotometer was used to determine the total content of chromium, lead, cadmium, manganese, iron, and zinc in the honey samples (Analytik Jena NovAA 400P, 2012, country of origin: Germany). The hollow cathode lamps (HCL) in AAS were used for estimations in various situations, depending on the element being analyzed. The Aspect LS software was used to manage the AAS instrument and generate the absorbance readout for each element in each honey digest sample. The instrumental operating conditions of AAS are summarized in Table 1.
Instrumental Conditions of AAS.
Quality Control (QC) of the Analytical Method
Specific wavelengths were selected for each metal to minimize interference from other elements under study. Glassware and plastic containers were cleaned by soaking them in 10% extra-pure nitric acid (HNO3) for at least 24 hours, followed by a thorough rinse with ultra-pure water. The purity of the reagents and laboratory equipment was further assessed by analyzing reagent blanks containing ultra-pure water. As with honey samples, we handled standards and blanks. The calibration curves were obtained using at least 5 diluted solutions at concentrations ranging from 0 to 1.00 mg/l, prepared from single-elemental stock solutions at 1000.00 mg/l. This sample was appropriately diluted with deionized water prior to AAS analysis to bring the concentration within the calibration range. The linearity and determination coefficients (R2) were derived using the calibration curves of each element. The method’s accuracy and precision were determined by calculating the mean recovery (%) and relative standard deviation (RSD) of the quality control (QC) sample. The limit of detection (LOD) and limit of quantification (LOQ) were calculated using the following equations:
Statistical Data Analysis
All analytical analyses were carried out in triplicate, and results were expressed as mean ± standard deviation (SD). Statistical analyses were performed using RStudio IDE and Microsoft Office Excel 2019. Independent samples t-tests were conducted to compare the mean concentrations of minerals and metals between local and imported honey samples. Statistical significance was set at P < .05. Pearson’s correlation analysis was performed to assess relationships between mineral and metal concentrations within each honey category. The reliability of correlation results was interpreted with caution given the limited sample size.
Human Health Risk Assessment
To quantitatively assess potential health risks associated with heavy metal exposure through honey consumption, the Estimated Daily Intake (EDI), Target Hazard Quotient (THQ), and Hazard Index (HI) were calculated following USEPA methodology. 24
where C is the metal concentration in honey (mg/kg, fresh weight), IR is the daily honey ingestion rate (assumed as 20 g/person/day for the Bangladeshi population based on national consumption patterns), and BW is the average body weight (70 kg for adults).
where RfD is the oral reference dose (mg/kg/day). The RfD values used were: Cd = 0.001, Cr = 1.5, and Pb = 0.0035. 24 THQ values <1 indicate no significant health risk, while THQ >1 suggests potential non-carcinogenic effects.
HI <1 indicates no significant risk, while HI >1 suggests potential cumulative non-carcinogenic effects.
Results and Discussion
Quality Control of the Instrumental Method
The linearity of the analytical method was evaluated using standard solutions of known concentrations, yielding excellent linearity, with R2 values ranging from .995 to .997 for all elements analyzed (Table 2). This demonstrates a strong correlation between analyte concentration and instrument response, ensuring accurate quantification of elements in honey samples. The LODs and LOQs are sufficiently low to detect trace levels of these elements in honey, ensuring the instrumental method’s sensitivity. The average recovery percentages for the analyzed minerals and metals are 95.80% for Fe, 104.30% for Zn, 101.70% for Mn, 95.60% for Cd, 104.90% for Cr, and 96.40% for Pb, confirming the accuracy of instrumental quantification. Precision, measured as the relative standard deviation (RSD), ranged from 1.20% for Zn to 6.80% for Cr and Pb. The low RSD values for most elements indicate that the method yields consistent and reproducible results, which are essential for reliable analytical performance (Table 2).
Quality Control Parameters of the Analytical Method.
LOQ values were converted from mg/l to mg/kg using a dilution factor of 10 (10 g sample digested to 100 ml final volume).
Essential Mineral Contents of Honey Samples
The concentrations of essential minerals, specifically iron (Fe), zinc (Zn), and manganese (Mn), were measured in various local and imported honey samples collected from Dhaka City. The results are detailed in Table 3, presenting the mean concentrations and standard deviations (SD) for each mineral in both local and imported honey samples. The Fe content in local honey samples ranged from 5.789 (LH6) to 16.654 mg/kg (LH1), with notable variability. LH1 had the highest Fe concentration, indicating potential differences in environmental or botanical sources. Imported honey samples showed a broader and higher range of Fe concentrations, from 3.995 (IH1) to 30.650 mg/kg (IH6). IH6 had a significantly higher Fe content than other samples, suggesting considerable variation in the source regions of imported honey. Statistical analysis using independent sample t-tests revealed that the difference in mean Fe concentrations between local (9.08 ± 4.33 mg/kg) and imported (15.59 ± 9.28 mg/kg) honey was not statistically significant (P = .148). Zn concentrations in local honey samples varied from 1.350 (LH4) to 4.346 mg/kg (LH3). LH3 exhibited the highest Zn content, reflecting differences in soil composition and botanical sources. Imported honey samples had Zn concentrations ranging from 1.276 (IH4) to 3.532 mg/kg (IH6). IH6 showed the highest Zn content, indicating variability among different imported sources. No significant difference was observed between local and imported honey for Zn (P = .887). Mn levels in local honey samples ranged from 0.241 (LH5) to 1.638 mg/kg (LH2). The imported honey samples had Mn concentrations ranging from 0.871 (IH4) to 1.794 mg/kg (IH6). Like Fe and Zn, IH6 had the highest Mn content among imported samples. Like Fe and Zn, the difference in Mn concentrations between local and imported honey was not statistically significant (P = .369). The significant variability in mineral concentrations, particularly the high Fe content in IH6 and high Zn in LH3, may reflect differences in botanical origin as well as environmental factors. Previous studies have demonstrated that unifloral honeys (eg, heather, manuka) can have distinctly different mineral profiles compared to multifloral varieties.6,25
Essential Mineral Contents in the Studied Honey Samples.
Toxic Metal Contents of the Honey Sample
The concentrations of toxic metals, specifically cadmium (Cd), chromium (Cr), and lead (Pb), were measured in 6 local honey (LH) samples and 6 imported honey (IH) samples collected from Dhaka City. The results are detailed in Table 4, presenting the mean concentrations and standard deviations (SD) for each metal in both local and imported honey samples. All samples had detectable levels of Cd and Cr, but Pb was below the LOQ (0.076 mg/kg) in all samples. The mean Cd concentrations in local honey samples range from 0.169 (LH5) to 1.608 mg/kg (LH1). The highest concentration observed in LH1 indicates potential localized sources of contamination. Four out of 6 local samples (LH1, LH2, LH3, LH4) exceeded the maximum allowable limit (MAL) of 0.20 mg/kg. 26 In imported honey samples, Cd concentrations are more uniform, ranging from 0.230 (IH3) to 0.295 mg/kg (IH2), and all 6 samples exceed the MAL. The difference in Cd concentrations between local (0.47 ± 0.57 mg/kg) and imported (0.252 ± 0.023 mg/kg) honey approached but did not reach statistical significance (P = .068), likely due to the high variability in local samples and the limited sample size.
Toxic Metal Contents in the Studied Local and Imported Honey Samples.
Cr levels in local honey range from 0.478 (LH3) to 1.132 mg/kg (LH1), with LH1 again showing the highest contamination level. Imported honey samples exhibit Cr concentrations from 0.420 (IH2) to 0.740 mg/kg (IH4). All samples were well below the MAL of 10.0 mg/kg 2 . 6 No significant difference was observed between local and imported honey for Cr (P = .116). The Pb concentrations in all honey samples were below the LOQ (0.076 mg/kg), which is below the MAL of 0.10 mg/kg. 27 This suggests that Pb contamination is not a significant concern in these honey samples. The higher and more variable concentrations of Cd and Cr in local honey samples suggest localized environmental or anthropogenic contamination sources. The higher and more variable concentrations of Cd in local honey samples suggest localized environmental or anthropogenic contamination sources. The elevated Cd levels observed, particularly in LH1 (1.608 mg/kg), warrant detailed consideration of potential contamination sources. Cadmium can enter the beekeeping environment through multiple pathways: industrial emissions,28,29 traffic density,10,30 agricultural inputs, industrial activities, and atmospheric deposition. In contrast, the more uniform Cd levels in imported honey (range: 0.230-0.295 mg/kg) and the absence of strong correlations between Cd and other elements suggest different contamination pathways, possibly related to soil background levels in the countries of origin rather than point-source contamination. The foraging range of Apis mellifera (typically 3-5 km) means that honey composition reflects the integrated environmental quality of a substantial area around the apiary, making honey a useful bioindicator of environmental contamination. 31
Correlation of the Coefficient of Minerals and Metals in Local and Imported Samples
The Pearson correlation coefficients (r) for various minerals and metals in local and imported honey samples are presented in Figure 1. The correlation analysis helps understand the relationships among elements within honey samples and assess potential patterns of contamination or enrichment. It should be noted that the correlation analysis is based on a relatively small sample size (n = 6 per group), which increases the risk of Type II errors and may affect the stability of correlation coefficients. Therefore, these correlation patterns should be interpreted with caution and considered as preliminary indicators of elemental associations that warrant confirmation in larger-scale studies. In local honey, the strongest positive correlations were observed among Fe-Cd (r = .89), Cd-Cr (r = .81), and Fe-Cr (r = .79), indicating that these 3 elements tended to increase or decrease together across samples. Moderate positive associations were also found for Mn-Cd (r = .49) and Fe-Mn (r = .31). In contrast, Zn showed weak-to-moderate negative relationships, particularly with Cr (r = −.63), while its correlations with Fe (r = −.11) and Cd (r = −.17) were weak. Notably, Mn-Cr was essentially uncorrelated (r ≈ .008), suggesting independent variation between these 2 elements in local honey.

Correlation coefficient (r) of mineral and metal for honey samples.
Imported honey displayed a different correlation structure dominated by strong positive interrelationships among Zn, Mn, and Fe. The highest correlation in the imported dataset occurred between Zn-Mn (r = .97), followed by Fe-Zn (r = .74), Fe-Cr (r = .72), and Fe-Mn (r = .65). Correlations involving Cd were predominantly negative or near zero: Cd-Cr showed a moderate negative association (r = −.55) and Cd-Fe a weaker negative association (r = −.35), while Cd-Zn (r = −.081) and Cd-Mn (r = −.016) were negligible. Associations between Cr and Zn (r = .28) or Mn (r = .24) were weakly positive.
Taken together, the matrices indicate that local honey is characterized by a Fe-Cd-Cr cluster of strong positive correlations, whereas imported honey is characterized by a Zn-Mn-Fe cluster, with Cd tending to vary independently or inversely, particularly relative to Cr. These contrasting correlation patterns suggest that elements co-vary differently depending on honey origin, which may reflect differences in underlying co-occurrence patterns across the 2 groups. The strong positive correlations among certain elements in local honey (Fe, Cd, Cr) may indicate a common source or pathway of contamination, likely related to industrial emissions or urban pollution. The distinct patterns between local and imported honey underscore the influence of geographical factors, environmental conditions, and production practices on elemental composition.
Quality and Safety of Local and Imported Honey
Across the 12 honeys analyzed and calculated as fresh weight basis (Supplemental Table 1; 6 local, 6 imported), Fe, Zn, Mn and Cr were quantifiable in all samples, while Pb was below the method reporting limit (MRL) in every case; Cd was measurable in all samples (Table 5). The nutritional value of honey can be evaluated using essential minerals as criteria, whereas the toxic metal content is used to address safety concerns. To enable meaningful comparison between honey mineral concentrations (mg/kg) and Recommended Daily Allowances (mg/day), we assumed an average daily honey intake of 20 g (0.020 kg) per person based on national consumption patterns in Bangladesh. Using this consumption rate, the actual daily mineral intake from honey was calculated and compared to RDA values.
Central Tendency and Range of Elements in Honey (mg/kg, fresh weight).
Abbreviations: F, female; M, male; MAL, maximum allowable limit.
Note. Mean concentrations are presented on fresh weight basis (converted from dry weight using moisture content correction). For comparison with RDA values (mg/day), daily intake from honey can be calculated assuming an average consumption of 20 g/day (eg, metal intake = concentration × 0.020 kg). For Pb, values below LOQ (0.076 mg/kg) are reported as <LOQ.

Comparison of the average mineral and metal content of the studied samples with the regulatory limit.
The differences in essential mineral and toxic metal content between local and imported honey highlight the influence of regional environmental factors, soil composition, and agricultural practices on honey’s nutritional quality. 42 The significant variability in mineral and metal content between local and imported honey underscores the importance of quality control and standardization in honey production to ensure consistent nutritional benefits and safety.
Human Health Risk Assessment
To quantitatively assess the potential health risks associated with heavy metal exposure through honey consumption, we calculated Estimated Daily Intake (EDI), Target Hazard Quotient (THQ), and Hazard Index (HI) following established methodologies. 24 Table 6 presents the EDI, THQ, and HI values for toxic metals in local and imported honey samples.
Health Risk Assessment Parameters For Toxic Metals in Honey Samples.
EDI calculations assume average daily honey consumption of 20 g/person and body weight of 70 kg. For Pb, concentrations below LOQ (0.076 mg/kg) were conservatively estimated as LOQ/2 = 0.038 mg/kg for EDI calculation. 24
The health risk assessment revealed that all THQ and HI values were below 1, indicating no significant non-carcinogenic health risk from honey consumption at current contamination levels and typical consumption rates. However, the THQ for Cd in local honey (0.107) was approximately 1.6 times higher than that of imported honey (0.066), reflecting the elevated Cd concentrations in local samples. While these values remain below the threshold of concern, they suggest that cumulative exposure from multiple dietary sources could potentially approach risk levels, particularly for individuals with high honey consumption or those living in areas with elevated environmental Cd levels.
The Hazard Index (HI) values of 0.113 for local honey and 0.072 for imported honey indicate that the cumulative non-carcinogenic risk from combined metal exposure through honey consumption is well below the acceptable limit. Cadmium contributed the largest proportion to the HI (95% for local honey; 92% for imported honey), highlighting that Cd is the primary metal of concern in these honey samples.
It should be noted that this risk assessment considers only honey consumption and does not account for other dietary sources of these metals. For individuals with additional exposure routes (eg, contaminated rice, vegetables, or drinking water), the cumulative risk could be higher. Furthermore, the assessment assumes average consumption rates; individuals with higher-than-average honey intake may face greater risks.
Conclusion
This study evaluated the quality and safety of honey available in Dhaka City by analyzing the mineral and heavy metal content of both local and imported honey samples. The results revealed variability in the concentrations of essential minerals (Fe, Zn, and Mn) and toxic heavy metals (Cd, Cr, and Pb) in both local and imported honey. While all honey samples met safety standards for Pb and Cr. Cd levels exceeded the maximum allowable limit in 67% of local samples and 100% of imported samples. This finding raises concerns about potential health risks from Cd exposure through honey consumption, particularly given that Cd contributed 92% to 95% of the cumulative hazard index. Health risk assessment using Target Hazard Quotients and Hazard Index indicated no immediate non-carcinogenic risk at current consumption levels (HI <1). However, the elevated Cd concentrations warrant attention, especially for local honey, as cumulative exposure from multiple dietary sources could approach risk thresholds. The strong positive correlations among Fe, Cd, and Cr in local honey suggest common anthropogenic contamination sources, likely related to industrial emissions, traffic density, or agricultural practices in Dhaka city. The study highlights the importance of regular monitoring and stringent quality control measures in honey production and importation to ensure consumer safety. Furthermore, additional research is warranted to investigate the sources of Cadmium contamination in local honey and implement effective mitigation strategies. Overall, this study offers valuable insights into the mineral and metal content of honey in Dhaka City, contributing to a deeper understanding of the quality and safety of local and imported honey samples.
Limitations of the Study
This study has several limitations that should be acknowledged. First, the sample size (n = 12) is relatively small, which restricts the statistical power and generalizability of the findings. The results should therefore be considered preliminary, and further studies with larger sample sizes are warranted to confirm and extend these observations. Second, the sampling was limited to honey available in the Mohammadpur area of Dhaka city, while the specific botanical origins of the honey samples were not determined, which may contribute to the observed variability in mineral content. Future studies should incorporate melissopalynological analysis to correlate elemental profiles with floral sources. Third, while we identified elevated cadmium levels in some samples, the specific sources of contamination could not be definitively determined within the scope of this study. Fourth, the study relied on a single analytical technique (AAS) without cross-validation using CRMs or alternative methods such as ICP-MS. While spike recovery data demonstrate acceptable accuracy, future studies should incorporate CRM analysis and multi-technique validation to further strengthen analytical reliability. Fifth, the health risk assessment considered only honey consumption and did not account for other dietary sources or potential synergistic effects of mixed contaminant exposure. Future research should address these limitations through larger-scale studies, source apportionment analysis, and comprehensive dietary exposure assessments.
Supplemental Material
sj-docx-1-ehi-10.1177_11786302261449341 – Supplemental material for Assessing Quality and Safety: Mineral and Metal Content in Local and Imported Honey From Dhaka City
Supplemental material, sj-docx-1-ehi-10.1177_11786302261449341 for Assessing Quality and Safety: Mineral and Metal Content in Local and Imported Honey From Dhaka City by Munjuri Akter, Simana Akhter Bhuiya, Mehedi Amin, Mohammed Ariful Islam, Tazul Islam Chowdhury, Mohammed Sakhawat Hossain and Abdul Kaium in Environmental Health Insights
Footnotes
Acknowledgements
The authors acknowledge all the laboratory and office staff of the Agricultural Chemistry Department, Sher-e-Bangla Agricultural University, Dhaka-1207, for their assistance in preparing and digesting the honey samples. We acknowledge that the AI-assisted technologies ChatGPT and Co-Pilot were used in the writing process to improve the language and readability of the manuscript.
Author Contributions
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
