Skip to main content


Front. Public Health, 13 January 2023
Sec. Environmental health and Exposome
This article is part of the Research Topic The One Health Approach in the Context of Public Health View all 11 articles

Carbon sequestration potential of different forest types in Pakistan and its role in regulating services for public health

  • 1Department of Plant Sciences, Quaid-i-Azam University, Islamabad, Pakistan
  • 2Member Pakistan Academy of Sciences, Islamabad, Pakistan
  • 3Department of Botany, Government College University, Lahore, Pakistan
  • 4Audit Team, Hanmoo Convention (Oakwood Premier), Seoul, Republic of Korea
  • 5College of Hospitality and Tourism Management, Sejong University, Seoul, Republic of Korea
  • 6CBIOS (Research Center for Biosciences and Health Technologies), Universidade Lusófona de Humanidades e Tecnologias, Lisboa, Portugal

A high amount of CO2 causes numerous health effects, including headaches, restlessness, difficulty in breathing, increased heart rate, high blood pressure, asphyxia, and dizziness. This issue of increasing atmospheric CO2 can only be solved via above-ground and below-ground carbon sequestration (CS). This study was designed to determine the relationship between CS with the crown area (CA), diameter at breast height (DBH), height (H), species richness (SR), and elevation in different forest types of Pakistan with the following specific objectives: (1) to quantify the direct and indirect relationship of carbon sequestration with CA, DBH, H, and SR in various natural forest types and (2) to evaluate the effect of elevation on the trees functional traits and resultant CS. We used the linear structural equation model (SEM) for each conceptual model. Our results confirmed that the highest CS potential was recorded for dry temperate conifer forests (DTCF) i.e., 52.67%, followed by moist temperate mix forests (MTMF) and sub-tropical broad-leaved forests (STBLF). The SEM further described the carbon sequestration variation, i.e., 57, 32, 19, and 16% under the influence of CA (β = 0.90 and P-value < 0.001), H (β = 0.13 and p-value = 0.05), DBH (β = 0.07 and p-value = 0.005), and SR (β = −0.55 and p-value = 0.001), respectively. The individual direct effect of SR on carbon sequestration has been negative and significant. At the same time, the separate effect of CA, DBH, and H had a positive and significant effect on carbon sequestration. The remaining 20% of CS variations are indirectly influenced by elevation. This means that elevation affects carbon sequestration indirectly through CA, DBH, H, and SR, i.e., β = 0.133 and P-value < 0.166, followed by β = 0.531 and P-value < 0.001, β = 0.007 and P-value < 0.399, and β = −0.32 and P-value < 0.001, respectively. It is concluded that abiotic factors mainly determined carbon sequestration in forest ecosystems along with the elevation gradients in Pakistan. Quantifying the role of various forest types in carbon dioxide (CO2) reduction leads to improved air quality, which positively impacts human health. This is an imperative and novel study that links the dynamics of the biosphere and atmosphere.

1. Introduction

The forest ecosystem provides clean air through carbon sequestration for better public health. The sequestered carbon is not emitted into the atmosphere, reducing the chances of various human diseases (1). Forest ecosystems determine global climate change by removing or adding greenhouse gases such as CO2 from the atmosphere. Forests, alongside the grassland and peat swamps, collectively store more carbon than other terrestrial ecosystems (2). Forest trees store CO2 from the atmosphere in the form of wood, soil organic carbon, and other biomasses which contribute to a reduction in global warming and climate change (3). Disturbances in these ecosystems due to the overexploitation of resources cause the release of a significant amount of CO2 back into the atmosphere (4).

Carbon sequestration is essential to mitigate global climate change (5). Approximately 53% of carbon is stored in temperate and boreal regions, while the remaining 37 % is stored in tropical areas (6, 7). Half of the terrestrial carbon sink is located in forest ecosystems (8). Forest ecosystems take a large amount of CO2 from the atmosphere via photosynthesis. They also contribute a large amount of fixed carbon in the form of organic matter in the lithosphere (9). However, a small portion of integrated carbon has been stored in the belowground biomass, litter, and soil (10). According to United Nations Food and Agriculture Organization (FAO) statistics, forests contain 234 Pg of carbon in the above-ground forest, 62 Pg below ground, 41 Pg in the form of deadwood, 23 Pg in litter form, and 398 Pg in the forest soils (11). Forest ecosystems occupy the most significant part of the free-ice land surface in the entire terrestrial world. Trees are the main component of the forest ecosystem and contain the entire quantity of forest living biomass. Forests' total biomass is ~677 Petagram (PgC), and trees contribute about 80% of the world's total biomass (11).

Several biotic and abiotic factors influence carbon sequestration in forest ecosystems (1214). Trees continuously sequester CO2 from the atmosphere and store it in their different parts, i.e., stem, root leaves, and branches (15). The carbon sequestration rate depends on plant growth, individual characteristics of the tree species, the wood density, and its growing conditions. Trees store most of the atmospheric carbon because of their larger size, volume, and long-lived storage capacities, i.e., tree trunks, leaves, roots, and the soils in which they exist (16). Carbon sequestration is mainly driven through huge tree biomasses (17).

The scientific community is committed to reducing atmospheric CO2 emissions and storing them in any other form to improve the quality of the atmosphere (18). Various biological and geological mechanisms bring free carbon down from the atmosphere. However, due to the enormous rise in human population, there is a continuous degradation of the land and deforestation, which results in the irreversible loss of forest functions. Plants are the primary producers of the world's biodiversity. They capture atmospheric CO2 and convert it into glucose, the first organic molecule of life (19). Degraded forests lead to low carbon storage and poor biodiversity, which in turn cause global warming and climatic changes (20). Anthropogenic activities drastically add massive amounts of greenhouse gases, mainly CO2, into the atmosphere (16). The rapid increase in the concentration of CO2, methane, and methane dioxide is the major cause of global warming (21).

Characteristically, above-ground biomass decreases along the elevation, but precipitation influences the biomass variations alongside the elevation gradient (22). Carbon sequestration is mainly determined by the biotic determinant of trees, such as variations in diameter, height diversity, stand density, and stand basal area (20, 23). The individual tree height is not associated with diameter, but the ratio is correlated to the species' genetic nature, intraspecific competition, and abiotic drivers like precipitation, temperature, and soil types (2426). However, tree heights decline with increasing elevation, generally (27). Initially, researchers use DBH data to estimate carbon sequestration, but the precision of carbon estimation has improved using tree height data (28, 29). Proper forest administration and land improvement approaches that concurrently enhance the biodiversity and carbon stock to mitigate climate change and global warming across regional and global scales are required to bring positive changes at least at the micro-environment levels (30). However, the relationship between forest stand structure, i.e., height, diameter at breast height, crown area, and species richness along regional level elevation gradients in different forest types, have rarely been evaluated. Therefore, the current study was designed to quantify the carbon sequestration potential of various forest types in Pakistan, focusing on its biotic complexities.

1.1. Theoretical framework

The net primary productivity theory explains how forest productivity is affected by forest structure, i.e., First, the net immediate productivity drops during stand expansion due to the limited supply of water and nutrients to leaves within a forest stand. Second, how productivity changes systematically with climate change interaction and withstands developmental processes. These arguments provide predictions about variability in wood production, biomass loss rate, and carbon sequestration. The literature shows that the forest ecosystem is driven by biotic and abiotic factors such as forest age, forest origin, forest soil conditions, and geography of the region. Forest structure and growth are controlled by vegetation cover. Through aerial expansion, trees increase their potential height, producing more biomass (31, 32). The effects of abiotic and biotic factors on carbon sequestration at large regional-scale forest ecosystems have rarely been studied. Thus, the driving biotic and abiotic factors of carbon sequestration in multiple forest types are vague and need to be explored.

Stand structure is associated with differences in species composition, and there is a positive relationship between functional and taxonomic species diversity with the above-ground biomass (23, 33). The relationship between carbon sequestration with tree crown area, diameter at breast height, height, elevation, and species richness is usually positive. Still, the association and simultaneous direct and indirect impacts in different forest ecosystems remain questionable. Therefore, current studies were conducted in diverse forest ecosystems to evaluate the role of various biotic and abiotic factors on carbon sequestration at a regional level.

The most critical goal of the current work was to determine the direct and indirect role of CA, H, DBH, and SR on carbon sequestration in the diverse natural forests of Pakistan. Thereby, we also examine the effect of elevation on Carbon sequestration. We hypothesize that carbon sequestration is driven by CA, DBH, H, and SR along the elevation from sea level. We used a structural equation modeling approach to test a hypothesized causal relationship among the response variables using data from 200 plots covering sub-tropical thorn forests, sub-tropical broad-leaved forests, moist temperate mix forests, dry temperate conifer forests, and dry temperate Pinus gerardiana (Chilgoza) forest plots in Pakistan. The critical research questions are as follows. (i) How does elevation gradients, CA, DBH, H, and SR explain variation in carbon sequestration? We hypothesized that carbon sequestration, CA, DBH, H, and SR decreases along with an increase in elevation (34). (ii) How does carbon sequestration change along the elevation under the influence of explanatory variables CA, DBH, H, and SR? We also hypothesized that tree CA increases carbon sequestration. (iii) How does SR affect carbon sequestration along CA, DBH, H, and elevation? We hypothesized that SR and carbon sequestration have no significant effect on each other. (iv) How does DBH coupled with other explanatory variables, explain variations in carbon sequestration? We hypothesized that DBH has a direct impact on carbon sequestration. (v) How does H impact carbon sequestration along with the CA, DBH, SR, and elevation? We hypothesized there is no effect of H on carbon sequestration (Figure 1).


Figure 1. Abstract model explaining the direct and indirect effects of the crown area, DBH, height, species richness, and elevation on carbon sequestration. Hypothesized associations among variables are designated by –, +, or –/+, respectively.

2. Materials and methodology

2.1. Study area

This study covered a geographical area from “25°18.057” to “35°10108” north latitude and “067°11.298” to “071°56.568” east longitude in the different forests of Pakistan. It has an evaluation range between 50 and 2,700 m above sea level. The different forest types studied are the (i) sub-tropical thorn forests (STTF) of Kirthar National Park, Kirthar mountains, Sindh, (ii) sub-tropical broad-leaved forests (STBLF) of Margalla Hills, lower Himalayas, (iii) moist temperate mixed forests (MTMF) of Murree western Himalayas, (iv) dry temperate coniferous forests (DTCF) of Kumrat, Hindu Kush mountains, and (v) the dry temperate Pinus gerardiana (Chilgoza) forests (DTPGF) of Shirani, Kohe Sulaiman Mountain, Balochistan (Table 1, Figure 2).


Table 1. Studied forest types along with their elevation range, mean annual temperature and precipitation.


Figure 2. GIS-generated map of the studied forests across Pakistan.

2.2. Data collection using quantitative ecological techniques

A detailed description of tools and techniques is explained in this section.

2.2.1. Estimation of carbon sequestrated by plants annually

Carbon sequestration for each woody tree species was calculated through total DBH and total tree height for all the individuals in each plot (35).

2.2.2. Determining total tree green weight

The following algorithmic equation was applied for weighting the carbon sequestration of the woody plant species for 10 years.

Wag=0.15 × D2H    (1)

where Wag is the above-ground weight of tree species in pounds (lbs), D is the diameter of the tree trunk in inches (for trees with D>1), and H is the height of the tree in feet. The green weight is the live tree weight. First, we calculated the above-ground green weight using the following method (36).

The belowground system weight of a tree is 20% greater than the above-ground weight. We calculated the total green weight of the tree by multiplying the above-ground weight by 1.2.

Wtgw=1.2×Wag    (2)

The tree's total average mass is ~72.5% dry matter, and the moisture content of the tree is 27.5%. Hence, we calculate the tree's dry weight by multiplying the total green weight of the tree by 0.725 (37).

WDw=0.725×Wtgw    (3)

The average carbon content in the tree is generally 50% of the total tree volume (16, 37). Thus, we determined the weight of carbon in the tree by multiplying the tree's dry weight by 50% or 0.5.

Wc=0.5×WDw    (4)

The chemical composition of CO2 is one molecule of carbon and two molecules of oxygen. The carbon atomic weight is 12.001115 and the oxygen atomic weight is 15.9994. Therefore, the determination of CO2 weight in trees is the ratio of CO2 C+2 × O = 43.999915 C is 43.999915/12.001115 = 3.6663. Thus, we determined the weight of sequestered carbon dioxide in the tree by multiplying the tree's carbon weight by 3.6663, rounded to 3.67 (16, 38).

Wco2=3.67×Wc (v)    (5)

Tree crown area (CA) was calculated using the formula for the area of an ellipse (Ae):

Ae=π(0.5x)×(0.5y)    (6)

where is the crown length and (y) is the width measured perpendicularly (5, 39).

2.3. Statistical analyses

A structural equation model (SEM) was constructed based on the following conceptual hypothesis: (a) the direct effect of CA, DBH, SR, H, and elevation on the carbon sequestration (b) the direct effect of elevation on CA, DBH, H, and SR; and (c) the indirect effect of elevation, DBH, CA, H, and SR on carbon sequestration. The goodness of fit was calculated using the goodness of fit index (GFI), the chi-squared test, standardized root mean square residual (SRMSR), the comparative fit index (CFI), Akaike's Information Criterion (AIC), and the Bayesian information criterion (BIC) (40). The causal effect amongst the dependent, independent, and mediator variables was tested in the theoretical model while calculating the predicted variable's total effect, direct and indirect, on a response variable. The description of the variable used in the SEM model is represented in Table 2. A VIF test was carried out to check the multicollinearity in data. In our case, the VIF test's importance is around 2.50. Hence, we proceeded with further regression.


Table 2. Description of the independent and dependent variables.

The linear regression for each conceptual path was calculated for the complete path measurement of the SEM result, which is shown in Figure 1. However, due to the direct and indirect effects of the other response variables, the results of the bivariate associations may or may not be reliable with the SEM outcome (4144). Thus, nine series of multiple linear regression models with various response variables were used, for example, (a) the multiple effects of carbon sequestration, CA, H, DBH, and SR; (b) multiple effects of mediator and carbon sequestration; (c) multiple direct and indirect effects of carbon sequestration and response variables. The Hmisc and corrplot packages evaluated the Pearson correlation coefficient matrix.

Before statistical analysis, for normality and linearity, all continuous numerical variables were normalized and standardized and all statistical analysis was done in R 4.0.2 (45).

2.4. Empirical model

We determined the conceptual model with several empirical models. Our generalized model as

Yco2i=α+βCAi+λDBHi+δHi+SRi+θElevi+μi    (7)

where CO2 showed carbon sequestration, CA represents aerial expansion, DBH is the diameter at breast height, H indicates the tree height, and Elev is the elevation of the studied areas, where (i) indicates the number of the plots and α, β, λ, ⋎, δ, and θ are the coefficients of the response variables. In contrast, μ is an unobserved variation or error term of the model. We arranged the generalized model into six parts to find the direct, indirect, and total effects and measure all the below equations simultaneously.

Yco2i=α+βCAi+λDBHi+δHiSRi+θElevi+εi    (8)
Yco2i=α+βCAi+λDBHi+δHi+SRi+εi    (9)
Yco2i=α+βCAi+λDBHi+SRi+εi    (10)
Yco2i=α+βCAi+λDBHi+εi    (11)
Yco2i=α+βCAi+εi    (12)

3. Results

Above and belowground CS was found highest in dry temperate conifer forests (DTCF) i.e., 52.67%, followed by moist temperate mix forests (MTMF) sub-tropical broad-leaved forests (STBLF), dry temperate Pinus gerardiana (Chilgoza Forest) (DTPGF), and sub-tropical thorn forests (STTF) with 29.99, 11.66, 4.52, and 1.14% contribution, respectively (Table 3). The SEM described 80% of the CS variation i.e., 57, 32, 19, and 17% under the influence of DBH, SR, CA, and H, of the forest trees, respectively (Figure 3). CA had a strong direct effect on the CS (β = 0.90 and P-value < 0.001), followed by the effect of H, DBH and SR with the values of β = 0.13 and P-value = 0.009, β = 0.07 and P-value = 0.039, and β = −0.55 and P-value = 0.001, respectively. The individual direct effect of SR on CS had been negative and significant. At the same time, the separate effects of CA, DBH, and H were positive and significant on CS. The remaining 20% of CS variations were indirectly influenced by elevation. This means that elevation affects CS indirectly through CA, DBH, H, and SR with values of β = 0.133, followed by β = 0.531 and P-value < 0.001, and P-value < 0.166, β = 0.007 and P-value < 0.399, and β = −0.32 and P-value < 0.001, respectively (Supplementary Table 1). Whereas, the effect of the mediators (CA, DBH, H, and SR) was examined (β = 0.06 and P-value < 0.001, followed by (β = −0.28 and P-value < 0.001, β = 0.50.6 and P-value < 0.000, and β = 0.32 and P-value < 0.068), respectively, on CS (Supplementary Table 1). The total effect of all response variables jointly on CS was 44% positive and significant (Supplementary Table 1). Pearson correlation analysis shows the relative contribution of CS, elevation, CA, DBH, H, and SR (Figures 4, 5). The correlation analyses evaluated the relationship between dependent and independent variables. The analyses show that elevation had a positive correlation with CS. This means that with an increase in elevation up to a limit, carbon sequestration also increases. The relation between carbon sequestration and the biotic variables were positively correlated with each other except SR, ie., when CA, DBH, and H increase the potential for CS also increases (Figure 5).


Table 3. Primary result of the study in different forest types.


Figure 3. SEM linked elevation, CS, CA, DBH, H, and SR in Pakistan's different forests.


Figure 4. Pearson correlation among all predictor variables. CS, Carbon sequestration; DBH, Diameter at breast height; CA, Crown area; SR, Species richness.


Figure 5. Correlation between the dependent and independent variables, i.e., (A) Carbon sequestration with elevation, (B) Carbon sequestration with H, (C) Carbon sequestration with DBH, (D) Carbon sequestration with CA, (E) Carbon sequestration with SR. STTF, Sub-tropical thorn forests; STBLF, Sub-tropical broad-leaved forests; MTMF, Moist temperate mix forests; DTCF, Dry temperate conifer forests; DTPGF, Dry temperate Pinus gerardiana (Chilgoza) forests.

The direct theoretical paths are examined with each box to remove the difficulties of the regression path in the model (Supplementary Tables 1, 2). The significant standardized regression coefficients are shown in bold text and the insignificant standardized coefficients are shown in the normal text. Model fit statistic: comparative fit index =0.9000; goodness of fit index = 0.9100; standardized root mean square residual = 0.0032; chi sq. value = 10.07 and P-value = 0.122; AIC = 2215; CIF = 0.885 and GIF = 0.965, respectively in our model. Supplementary Table 1 shows the standardized direct, indirect, and total SR (b4*a1) on various elevation gradients along the regional scale. The indirect effect of CA, DBH, H, and SR significantly affected CS, and the total effect of all variables is also positive and significant (Supplementary Table 1).

The Pearson correlation coefficient measures the strength of the association between the two variables. In the first step of the Pearson correlation, we examined the relationship between two continuous variables and checked the significant effect among all the response variables through the correlation P-value. All the variables had a substantial relationship in our case, so we proceeded with further statistical analysis.

4. Discussion

This is the complete estimation of carbon sequestration in different forest types, i.e., sub-tropical thorn forest, sub-tropical broad-leaved forest, moist temperate mix forest, dry temperate conifer forest, and dry temperate Pinus gerardiana (Chilgoza) forest along regional-scale elevation gradients in Pakistan. It indicates that CA had a substantial direct effect on carbon sequestration, followed by the impact of H, DBH, elevation, and SR, respectively. We will link our findings on carbon sequestration with previous local, regional, and global studies in this section. The findings of carbon sequestration along elevation gradients with the previous study conducted at Mt Changbai, based on H, DBH, and SR, all decreased significantly with an increase in elevation. As a result, elevation negatively affects carbon sequestration (46, 47). Similar to other previous regional level studies indicate that carbon sequestration, all above-ground and below-ground living biomass were significantly and negatively related to elevation on Mt Changbai, probably due to the temperature drop along the elevation gradient, a comparative shortage of water, and low soil temperature for tree growth on high elevation (48). More importantly, this may be because the carbon dioxide did not reach the higher elevations and the trees absorbed more carbon dioxide at the lower elevations. Many researchers have reported that carbon sequestration decreased with increasing elevation (49). Growth can be restricted by water shortages, exposure, reduced temperatures, transpiration rates, deprived soil quality, and low soil temperatures (50). However, the findings vary with studies from other tropical forests, where the biomass decreases with elevation, resulting in a decline in carbon sequestration (22). Xu et al. (51) studied subtropical forests and found that the canopy density is the major biotic factor in modulating vegetation carbon sequestration in forest ecosystems (51).

Ecologists are interested in the potential functional relation between carbon sequestration and species richness (52). The relationship between CS and species richness is still uncertain and depends on other factors, i.e., elevation and physiographic factors. In the present study, we found that carbon sequestration has a significantly negative relationship with species richness. These findings are supported by other studies (53). Sharma et al. (53), from the Garwal region reported a negative correlation between SR and carbon sequestration. Potter (54) concluded that SR is essential, but not the most critical and appropriate metric in biodiversity and carbon sequestration (54). These findings show similarity with our results, i.e., carbon sequestration has a negative relationship with species richness. However, carbon sequestration was directly affected by tree height and DBH in our study. Another study strongly supports the findings of the current one, i.e., carbon density was directly correlated with tree DBH and height (55).

Similarly, Ali et al. (39), concluded that carbon sequestration in the forest is mainly driven by tree height, DBH, crown area, and stand density in the reserve forests of Pakistan (39). Nevertheless, other studies reported no correlation between carbon sequestration with species richness and tree height, suggesting that carbon sequestration has no relationship with SR in forest communities. These differences may be due to environmental conditions like precipitation, temperature, etc., which are not present in our study. Despite this, this study directly shows the effect on carbon sequestration by H, DBH, and CA. Similarly, previous studies reveal that tree crown area in complex tropical forests has a strong positive relationship to carbon sequestration across large-scale ecological gradients. These results highlight new perceptions of the tree crown area complementarity mechanism by which the species having a variety of functional traits cause high resource capture and therefore enhance the above-ground biomass (56, 57).

These findings are supported by Afzal and Akhtar (38). According to these results, the DBH and height have a significant relationship with carbon sequestration. Carbon sequestration mainly depends on the H and DBH of the trees (58). These results indicate that CA, DBH, and H, excluding SR, have directly affected carbon sequestration. Tree height and DBH are directly related to carbon sequestration (59). Ali et al. (60) conducted a study in the subtropical forest of eastern China and a skimmia superba broad leaves forest. They revealed that tree height and DBH significantly correlate with carbon sequestration in the subtropical forest (60). Similarly, Moreno et al. (61) conducted a study in southwest Spain and reported similar results to ours that tree height and DBH have significant and positive correlations with tree carbon sequestration capability (61). They concluded that tree height and DBH have a direct relationship with the above-ground biomass and biomass directly related to carbon sequestration (61). Diaz-Balteiro et al. (62) report a similar result from southwest Spain that an increase in tree height and DBH significantly increase carbon sequestration in the forest ecosystems (62). Similarly, the carbon sequestration capacity of tree height and the DBH of the terrestrial vegetation also significantly correlated with each other (63). All these studies strongly support that carbon sequestration mainly depends on tree growth i.e., CA, DBH, H, and SR.

The present study documents the relevant information on live tree CS along an elevation gradient for different forest types. The present study's findings will help us understand the CS pattern in different forest types. This study will also be helpful for researchers to further understand the relationship between CS, CA, DBH, H, and SR along an elevation gradient. Overall, the current research suggested that extensive studies are required to understand the relationship between carbon sequestration and abiotic factors, i.e. soil physicochemical properties, to understand CS's mechanism and natural driver in the wild forest ecosystem. In this study, we did not focus on the abiotic drivers of carbon sequestration, which may have a significant effect on carbon sequestration in the natural forest ecosystems of Pakistan.

5. Conclusions

The present study concluded that a higher above and belowground carbon sequestration potential is found in the dry temperate conifer forest followed by moist temperate forest, sub-tropical broad-leaved forest, dry temperate Pinus gerardiana (Chilgoza) forest, and sub-tropical thorn forest. The findings of our research show that CA directly affected carbon sequestration, followed by H, DBH, and SR. There was an insignificant and positive effect recorded of elevation on carbon sequestration. The indirect impact of elevation on carbon sequestration through CA, DBH, H, and SR was recorded positively and significantly. It is also concluded that carbon sequestration is mainly affected by elevation. Furthermore, studies are required to evaluate the relationship between carbon sequestration and multiple biotic and abiotic drivers.

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 authors.

Author contributions

SA contributed to the write-up, fieldwork, and experimental work. SK and AR provided overall supervision, and reviewed and edited the final draft. ZA and ZS helped in data analysis and initial manuscript writing. AU contributed to the field and experimental work. SY and HH provided financial support and edited the final draft. All authors contributed to the article and approved the submitted version.


This paper is a part of PhD study conducted by the first author SA. The authors indebted their deepest thanks to Mr. Shah Fahad Ali Shah, School of Business and Economics, University of Brunei Darussalam for providing help with statistical analysis. The authors also extend their appreciation to the Mr. Naeemullah Kazi, Sindh Wildlife Department Govt of Sindh Pakistan for facilitation during fieldwork.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at:


1. Jose S. Agroforestry for ecosystem services and environmental benefits: an overview. Agroforestry Syst. (2009) 76:1–10. doi: 10.1007/s10457-009-9229-7

CrossRef Full Text | Google Scholar

2. Lal R. Soil carbon sequestration impacts on global climate change and food security. Science. (2004) 304:1623–7. doi: 10.1126/science.1097396

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Zhang K, Ali A, Antonarakis A, Moghaddam M, Saatchi S, Tabatabaeenejad A, et al. The sensitivity of North American terrestrial carbon fluxes to spatial and temporal variation in soil moisture: An analysis using radar-derived estimates of root-zone soil moisture. J Geophys Res Biogeosci. (2019) 124:3208–31. doi: 10.1029/2018JG004589

CrossRef Full Text | Google Scholar

4. Daily GC, Kareiva PM, Polasky S, Ricketts TH, Tallis H, et al. Mainstreaming natural capital into decisions. Nat Capital Theory Practice Mapping Ecosyst Serv. (2011) 7:3–14. doi: 10.1093/acprof:oso/9780199588992.003.0001

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Van de Perre F, Willig MR, Presley SJ, Andemwana FB, Beeckman H, Boeckx P, et al. Reconciling biodiversity and carbon stock conservation in an Afrotropical forest landscape. Sci Adv. (2018) 4:eaar6603. doi: 10.1126/sciadv.aar6603

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sohngen B, Sedjo R. Carbon sequestration in global forests under different carbon price regimes. Energy J. (2006). doi: 10.5547/ISSN0195-6574-EJ-VolSI2006-NoSI3-6

CrossRef Full Text | Google Scholar

7. Li W, Shi Y, Zhu D, Wang W, Liu H, Li J, et al. Fine root biomass and morphology in a temperate forest are influenced more by the nitrogen treatment approach than the rate. Ecol Indic. (2021) 130:108031. doi: 10.1016/j.ecolind.2021.108031

CrossRef Full Text | Google Scholar

8. Canadell JG, Quéré CL, Raupach MR, Field CB, Buitenhuis ET, Ciais P, et al. Contributions to accelerating atmospheric CO2 growth from economic activity, carbon intensity, and efficiency of natural sinks. Proc Nat Acad Sci. (2007) 104:18866–70. doi: 10.1073/pnas.0702737104

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Kulmala M, Nieminen T, Chellapermal R, Makkonen R, Bäck J, Kerminen VM, et al. Climate Feedbacks Linking the Increasing Atmospheric Co2 Concentration, BVOC Emissions, Aerosols and Clouds in Forest Ecosystems, in Biology, Controls and Models of Tree Volatile Organic Compound Emissions. Berlin: Springer. (2013). p. 489–508. doi: 10.1007/978-94-007-6606-8_17

CrossRef Full Text | Google Scholar

10. Fan J, Zhong H, Harris W, Yu G, Wang S, Hu Z, et al. Carbon storage in the grasslands of China based on field measurements of above-and below-ground biomass. Clim Change. (2008) 86:375–96. doi: 10.1007/s10584-007-9316-6

CrossRef Full Text | Google Scholar

11. Kindermann G, Obersteiner M, Sohngen B, Sathaye J, Andrasko K, Rametsteiner E, et al. Global cost estimates of reducing carbon emissions through avoided deforestation. Proc Nat Acad Sci. (2008) 105:10302–7. doi: 10.1073/pnas.0710616105

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Phelps J, Webb EL, Adams WM. Biodiversity co-benefits of policies to reduce forest-carbon emissions. Nat Clim Chang. (2012) 2:497–503. doi: 10.1038/nclimate1462

CrossRef Full Text | Google Scholar

13. Manan F, Khan SM, Muhammad Z, Ahmad Z, Abdullah A, ur Rahman et al. Floristic composition, biological spectrum, and phytogeographic distribution of the Bin Dara Dir, in the western boundary of Pakistan. Front Glob Change. (2022) 5:1019139. doi: 10.3389/gc.2022.10191

CrossRef Full Text | Google Scholar

14. Ejaz U, Khan SM, Aqeel M, Khalid N, Sarfraz W, Naeem N, et al. Use of Parthenium hysterophorus with synthetic chelator for enhanced uptake of cadmium and lead from contaminated soils—a step toward better public health. Front Public Health. (2022) 10:1009479. doi: 10.3389/fpubh.2022.1009479

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Chavan BL, Rasal GB. Sequestered standing carbon stock in selective tree species grown in University campus at Aurangabad, Maharashtra, India. Int J Eng Sci Technol. (2010) 2:3003–7.

Google Scholar

16. Toochi E. Carbon sequestration: how much can forestry sequester CO2. Forestry Res Eng Int J. (2018) 2:148–50. doi: 10.15406/freij.2018.02.00040

CrossRef Full Text

17. Malhi YM, Brown S. Forests, carbon and global climate. Philosophical Trans Royal Soc London Series A Math Phys Eng Sci. (2002) 360:1567–91. doi: 10.1098/rsta.2002.1020

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Smith P, Ashmore MR, Black HI, Burgess PJ, Evans CD, Quine TA, et al. The role of ecosystems and their management in regulating climate, and soil, water and air quality. J Appl Ecol. (2013) 50:812–29. doi: 10.1111/1365-2664.12016

CrossRef Full Text | Google Scholar

19. Poorter L, Bongers F, Aide TM, Zambrano AMA, Balvanera P, Becknell JM, et al. Biomass resilience of Neotropical secondary forests. Nature. (2016) 530:211–4. doi: 10.1038/nature16512

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Poorter L, van der Sande MT, Thompson J, Arets EJ, Alarcón A, Álvarez-Sánchez J, et al. Diversity enhances carbon storage in tropical forests. Global Ecol Biogeogra. (2015) 24:1314–28. doi: 10.1111/geb.12364

CrossRef Full Text | Google Scholar

21. Bernstein L, Bosch P, Canziani O, Chen Z, Christ R, Riahi K. IPCC, 2007: Climate Change 2007: Synthesis Report. Geneva: IPCC. (2008).

Google Scholar

22. Girardin CA, Malhi Y, Aragao LE, Mamani M, Huaraca Huasco W, Durand L, et al. Net primary productivity allocation and cycling of carbon along a tropical forest elevational transect in the Peruvian Andes. Glob Chang Biol. (2010) 16:3176–92. doi: 10.1111/j.1365-2486.2010.02235.x

CrossRef Full Text | Google Scholar

23. Fotis P, Karkalakos S, Asteriou D. The relationship between energy demand and real GDP growth rate: The role of price asymmetries and spatial externalities within 34 countries across the globe. Energy Econ. (2017) 66:69–84. doi: 10.1016/j.eneco.2017.05.027

CrossRef Full Text | Google Scholar

24. Nagendra H, Ostrom E. Polycentric governance of multifunctional forested landscapes. Int J Commons. (2012) 6:104–33. doi: 10.18352/ijc.321

CrossRef Full Text | Google Scholar

25. Feldpausch TR, Banin L, Phillips OL, Baker TR, Lewis SL, Quesada CA, et al. Height-diameter allometry of tropical forest trees. Biogeosciences. (2011) 8:1081–106. doi: 10.5194/bg-8-1081-2011

CrossRef Full Text | Google Scholar

26. Feldpausch TR, Lloyd J, Lewis SL, Brienen RJ, Gloor M, Monteagudo Mendoza A, et al. Tree height integrated into pantropical forest biomass estimates. Biogeosciences. (2012) 9:3381–403. doi: 10.5194/bg-9-3381-2012

CrossRef Full Text | Google Scholar

27. Aiba S-I, Kitayama K. Structure, composition and species diversity in an altitude-substrate matrix of rain forest tree communities on Mount Kinabalu, Borneo. Plant Ecol. (1999) 140:139–57.

Google Scholar

28. Brown S. Measuring carbon in forests: current status and future challenges. Environ Pollut. (2002) 116:363–72. doi: 10.1016/S0269-7491(01)00212-3

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Lewis SL, Lloyd J, Sitch S, Mitchard ET, Laurance WF. Changing ecology of tropical forests: evidence and drivers. Annu Rev Ecol Evol Syst. (2009) 40:529–49. doi: 10.1146/annurev.ecolsys.39.110707.173345

CrossRef Full Text | Google Scholar

30. Villa PM, Martins SV, de Oliveira Neto SN, Rodrigues AC, Martorano LG, Monsanto LD, et al. Intensification of shifting cultivation reduces forest resilience in the northern Amazon. For Ecol Manage. (2018) 430:312–20. doi: 10.1016/j.foreco.2018.08.014

CrossRef Full Text | Google Scholar

31. Grytnes JA. Fine-scale vascular plant species richness in different alpine vegetation types: relationships with biomass and cover. J Vegetation Sci. (2000) 11:87–92. doi: 10.2307/3236779

CrossRef Full Text | Google Scholar

32. Ji X, Lee KT, Nazar LF, A. highly ordered nanostructured carbon–sulphur cathode for lithium–sulphur batteries. Nat Mater. (2009) 8:500–6. doi: 10.1038/nmat2460

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Yan Y, Jarvie S, Liu Q, Zhang Q. Effects of fragmentation on grassland plant diversity depend on the habitat specialization of species. Biol Conserv. (2022) 275:109773. doi: 10.1016/j.biocon.2022.109773

CrossRef Full Text | Google Scholar

34. Gillman LN, Wright SD, Cusens J, McBride PD, Malhi Y, Whittaker RJ, et al. Latitude, productivity and species richness. Global Ecol Biogeography. (2015) 24:107–17. doi: 10.1111/geb.12245

CrossRef Full Text | Google Scholar

35. Wright IJ, Reich PB, Westoby M, Ackerly DD, Baruch Z, Bongers F. et al. The worldwide leaf economics spectrum. Nature. (2004) 428:821–7. doi: 10.1038/nature02403

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Clark III A, Saucier JR, McNab WH. Total-tree weight, stem weight, and volume tables for hardwood species in the southeast. Georgia Forest Research Paper. (1986). p. 60.

Google Scholar

37. DeWald SJ. Heating with Wood: Producing, Harvesting and Processing Firewood. Cooperative Extension. Lincoln: Institute of Agriculture and Natural Resources (2005).

38. Afzal M, Akhtar AM. Factors affecting carbon sequestration in trees. J Agric Res. (2013) 51:61–9.

Google Scholar

39. Ali S, Khan SM, Siddiq Z, Ahmad Z, Ahmad KS, Abdullah A, et al. Carbon sequestration potential of reserve forests present in the protected Margalla Hills National Park. J King Saud Univ Sci. (2022) 34:101978. doi: 10.1016/j.jksus.2022.101978

CrossRef Full Text | Google Scholar

40. Hoyle RH. Handbook of Structural Equation Modeling. New York: Guilford Press. (2012).

Google Scholar

41. Poorter L, van der Sande MT, Arets EJ, Ascarrunz N, Enquist BJ, Finegan B, et al. Biodiversity and climate determine the functioning of Neotropical forests. Global Ecol Biogeogr. (2017) 26:1423–34. doi: 10.1111/geb.12668

CrossRef Full Text | Google Scholar

42. Ahmad Z, Khan SM, Afza R, Ullah A, Zeb SA, Issayeva KS, et al. Angiosperms distribution under the influence of microclimatic factors across a polluted ecosystem. J Hazard Mater Adv. (2023) 9:100223. doi: 10.1016/j.hazadv.2022.100223

CrossRef Full Text | Google Scholar

43. Ahmad Z, Khan SM, Page SE, Balzter H, Ullah A, Ali S, et al. Environmental sustainability and resilience in a polluted ecosystem via phytoremediation of heavy metals and plant physiological adaptations. J Clean Prod. (2023) 385:135733. doi: 10.1016/j.jclepro.2022.135733

CrossRef Full Text

44. Ullah H, Khan SM, Jahangir S, Ullah Z, Ali I, Ahmad Z, et al. Vegetation assessments under the influence of environmental variables from the Yakhtangay Hill of the Hindu-Himalayan range, North Western Pakistan. Scient Rep. (2022) 12:1–16. doi: 10.1038/s41598-022-21097-4

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Antoch J. Environment for Statistical Computing. Austria: R Foundation for Statistical Computing. (2019).

Google Scholar

46. Shu-Qing Z. Composition, structure and species diversity of plant communities along an altitudinal gradient on the northern slope of Mt. Changbai, Northeast China. Biodiversity Sci. (2004) 12:164. doi: 10.17520/biods.2004020

CrossRef Full Text | Google Scholar

47. Zhu B. Carbon stocks of main forest ecosystems in Northeast China. (MS Thesis) Beijing: Peking University. (2005).

Google Scholar

48. Luo JJ, Masson S, Behera S, Shingu S, Yamagata T. Seasonal climate predictability in a coupled OAGCM using a different approach for ensemble forecasts. J Clim. (2005) 18:4474–97. doi: 10.1175/JCLI3526.1

CrossRef Full Text | Google Scholar

49. Moser G, Hertel D, Leuschner C. Altitudinal change in LAI and stand leaf biomass in tropical montane forests: a transect study in Ecuador and a pan-tropical meta-analysis. Ecosystems. (2007) 10:924–35. doi: 10.1007/s10021-007-9063-6

CrossRef Full Text | Google Scholar

50. Way DA, Oren R. Differential responses to changes in growth temperature between trees from different functional groups and biomes: a review and synthesis of data. Tree Physiol. (2010) 30:669–88. doi: 10.1093/treephys/tpq015

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Xu L, Shi Y, Fang H, Zhou G, Xu X, Zhou Y, et al. Vegetation carbon stocks driven by canopy density and forest age in subtropical forest ecosystems. Sci Total Environ. (2018) 631:619–26. doi: 10.1016/j.scitotenv.2018.03.080

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Kirby KR, Potvin C. Variation in carbon storage among tree species: implications for the management of a small-scale carbon sink project. For Ecol Manage. (2007) 246:208–21. doi: 10.1016/j.foreco.2007.03.072

CrossRef Full Text | Google Scholar

53. Sharma CM, Mishra AK, Krishan R, Tiwari OP, Rana YS. Variation in vegetation composition, biomass production, and carbon storage in ridge top forests of high mountains of Garhwal Himalaya. J Sustain Forest. (2016) 35:119–32. doi: 10.1080/10549811.2015.1118387

CrossRef Full Text | Google Scholar

54. Potter KM, Woodall CW. Does biodiversity make a difference? Relationships between species richness, evolutionary diversity, and aboveground live tree biomass across US forests. Forest Ecol Manage. (2014) 321:117–29. doi: 10.1016/j.foreco.2013.06.026

CrossRef Full Text | Google Scholar

55. Måren IE, Sharma LN. Seeing the wood for the trees: Carbon storage and conservation in temperate forests of the Himalayas. For Ecol Manage. (2021) 487:119010. doi: 10.1016/j.foreco.2021.119010

CrossRef Full Text | Google Scholar

56. Jucker T, Bouriaud O, Coomes DA. Crown plasticity enables trees to optimize canopy packing in mixed-species forests. Funct Ecol. (2015) 29:1078–86. doi: 10.1111/1365-2435.12428

CrossRef Full Text | Google Scholar

57. Morin X, Fahse L, Scherer-Lorenzen M, Bugmann H. Tree species richness promotes productivity in temperate forests through strong complementarity between species. Ecol Lett. (2011) 14:1211–9. doi: 10.1111/j.1461-0248.2011.01691.x

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Shimamoto CY, Botosso PC, Marques MC. How much carbon is sequestered during the restoration of tropical forests? Estimates from tree species in the Brazilian Atlantic forest. Forest Ecol Manag. (2014) 329:1–9. doi: 10.1016/j.foreco.2014.06.002

CrossRef Full Text | Google Scholar

59. Arif M, Shahzad MK, Elzaki EE, Hussain A, Zhang B, Yukun C, et al. Biomass and carbon stocks estimation in Chichawatni irrigated plantation in Pakistan. Int J Agric Biol. (2017) 19:1339–49.

Google Scholar

60. Laliberté E, Zemunik G, Turner BL. Environmental filtering explains variation in plant diversity along resource gradients. Science. (2014) 345:1602–5. doi: 10.1126/science.1256330

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Maldonado-Alconada AM, Castillejo MÁ, Rey M-D, Labella-Ortega M, Tienda-Parrilla M, Hernández-Lao T, et al. System report: Iberian Dehesas. Spain. (2016).

Google Scholar

62. Diaz-Balteiro L, Alonso R, Martínez-Jaúregui M, Pardos M, et al. Selecting the best forest management alternative by aggregating ecosystem services indicators over time: a case study in central Spain. Ecol Indic. (2017) 72:322–9. doi: 10.1016/j.ecolind.2016.06.025

CrossRef Full Text | Google Scholar

63. Lal R, Smith P, Jungkunst HF, Mitsch WJ, Lehmann J, Nair PR, et al. The carbon sequestration potential of terrestrial ecosystems. J Soil Water Conserv. (2018) 73:145A−52A. doi: 10.2489/jswc.73.6.145A

CrossRef Full Text | Google Scholar

Keywords: trees functional traits, carbon sequestration, elevation, structural equation model, regional scale, air quality

Citation: Ali S, Khan SM, Ahmad Z, Siddiq Z, Ullah A, Yoo S, Han H and Raposo A (2023) Carbon sequestration potential of different forest types in Pakistan and its role in regulating services for public health. Front. Public Health 10:1064586. doi: 10.3389/fpubh.2022.1064586

Received: 08 October 2022; Accepted: 28 December 2022;
Published: 13 January 2023.

Edited by:

Mohiuddin Md. Taimur Khan, Washington State University Tri-Cities, United States

Reviewed by:

Purabi Saikia, Central University of Jharkhand, India
Keith Dana Thomsen, Lawrence Livermore National Laboratory (DOE), United States
Tahira Jatt, Shah Abdul Latif University, Pakistan

Copyright © 2023 Ali, Khan, Ahmad, Siddiq, Ullah, Yoo, Han and Raposo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Shujaul Mulk Khan, yes; Sunghoon Yoo, yes; Heesup Han, yes

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.