Phytosociology and Dendrochronological Study of Central Karakoram National Park, Northern Areas (Gilgit-Baltistan), Pakistan By ALAMDAR HUSSAIN
Under the supervision of Professor Dr. Moinuddin Ahmed
A thesis submitted to the Graduate Research Management Council in partial fulfillment for the degree of doctor of Philosophy
2013
Laboratory of Dendrochronology and Plant Ecology, Department of Botany, Federal Urdu University of Arts Science and Technology, Karachi-75300, Pakistan
In the best name of Allah, Almighty, the most Beneficent, and the most Merciful
“And use split the earth in fragments and produce there corn, and grapes and nutrition plants, and olives and dates and enclosed gardens dense with lofty trees, and fruits and fodder, for use and convenience to and your cattle”
(Al-Quran, Part 30, Surah Abas 80, Ayat 25-32)
I
Federal Urdu University of Arts, Sciences & Technology GulshanEIqbal Karachi Department Of Botany
CERTIFICATE
Certified that Mr. Alamdar Hussain was enrolled for the program M.Phil leading to Ph.D. He has passed the required course work successfully and fulfilled all the criteria of Higher Education Commission (HEC) for the degree of doctorate. The dissertation titled “Phytosociology and Dendrochronological Study of Central Karakoram National Park, Northern Areas (Gilgit Baltistan), Pakistan” submitted by him is satisfactory and confide for the partial requirement of the award of Doctor of Philosophy degree.
‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐ ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐ Chairperson Supervisor Department of Botany Prof. Dr. Moinuddin Ahmed (Foreign Professor) II
III
DEDICATION
This thesis is lovingly dedicated to my respective parents, brothers,
sisters, relatives and honorable teachers who have been my constant
source of inspiration. They have given me the drive and discipline to
tackle any task with enthusiasm and determination. Without their love
and support this thesis would not have been made possible. IV
ACKNOWLEDGEMENTS
First of all my greatest thanks to Almighty ALLAH for bestowing upon me the courage to face the complexities of life and complete this thesis successfully. I would like to gratefully acknowledge the enthusiastic supervision of Professor Dr. Moinuddin Ahmed (Foreign Professor) for the unforgettable supervision, useful comments, remarks and encouragement through out the learning process of this thesis. I have no word to explain his generosity and I proudly accept him as an ideal person for my whole life. Furthermore I would like to thank Professor Dr. S. Shahid Shaukat (Eminent Scholar) for introducing me to the statistical analysis as well for the support on multivariate analysis and valuable comments and suggestions on the manuscript. I am grateful to my best colleague and Research associate Mr. Muhammad Akbar, who has supported me throughout this work, both on the field and laboratory. I would be remiss if I did not mention Dr. Muhammad Faheem Siddiqui, Dr. Muhammad Uzair Khan, Dr. Kanwal Nazim, Dr. Muhammad Wahab and Dr. Nasrullah Khan, for their constant moral support as well as for useful information for the completion of the thesis. I wish to avail myself this opportunity to express a sense of gratitude and love to my friend and senior Dr. Toqeer Ahmed Rao for his generosity, support, valuable suggestion and encouragement during the completion of the thesis. I take this opportunity to express my gratitude to my father Gulam Jan and elder brother Dildar Hussain, who have been instrumental in the successful completion of this thesis. I can’t say thank them enough for their tremendous financial and moral support. I feel motivated and encouraged every time and without their encouragement and guidance this thesis would not have materialized and the reason why I am here. The guidance and support received from all the teaching and non-teaching members of Botany Department who contributed to this thesis in various ways. I am grateful for their constant support and help. I cannot forget to pay my thank to Laboratory staff members Mr. Azhar Ali Kazmi, Mr. Mirza Waqar Ali Baig, Mr. Abdul Basit and Mr. Muhammad Fahad for their support and love. I am very thankful to Mr. Ahsan (GIS laboratory supervisor) for his support and cooperation in mapping and other facilities. V
My thanks and appreciations also go to my class fellows Mr. Faisal Hussain, Hina Zafar, Fariha Naz, Farzana Usman, Muhammad Usama Zafar, Jawaria Sultana, Shaheena Arshad, Tuba Tahir and people who have willingly helped me out with their abilities for their support and love during the completion of this thesis. It is my pleasure to acknowledge my elder sister Miss. Shaheen Abbas (Assistant Professor) and brother Muhammad Ilyas for their encouragement support and love. I wish to express my sincere gratitude to Dr. Sher Wali (Assistant Professor KIU) for his help in identification of plant specimens. I take an immense pleasure in thanking Mr. Muhammad Ismail (Forest Conservator Gilgit-region), Mr. Syed Yasir Abbas Rizvi (Ecologist, CKNP), Mr. Babar Khan (WWF Head, Gilgit) for their guidance and support during field work. I am deeply indebted to Seerat Hussain, my cousin Ashiq Hussain, Arif Hussain, Sajid Ali, Kamran Haider, Muzahir Hussain with family, for their love hospitality and support in the field. Last but not the least I cannot express my gratitude to my mother in words, whose unconditional love and prayers have been my greatest strength. Without any doubt the most astonishingly patience women I have ever seen. The constant love and support of my younger sister Sajida Batool is also sincerely acknowledged. I cannot forget my brothers Imdad Hussain, Sakhawat Hussain and Sajid Hussain for their unconditional support and love throughout my life.
Alamdar Hussain VI
TABLE OF CONTENTS
Certificate……………………………………………………………………………… I Thesis similarity index result………………………………………………………….. II Dedication……………………………………………………………………………... III Acknowledgments…………………………………………………………………….. IV Table of Contents……………………………………………………………………… VI List of tables…………………………………………………………………………… XV List of figures………………………………………………………………………….. XVII Abstract in Urdu……………………………………………………………………….. 1 Abstract………………………………………………………………………………... 3 General introduction…………………………………………………………………... 5 Aims and objectives…………………………………………………………………… 6 CHAPTER 1 DESCRIPTION OF THE STUDY AREA 1.1-Geography of Gilgit-Baltistan…………………………………………………….. 7 1.2-Karakorum Mountainous Range…………………………………………………... 8 1.3-Geology of Gilgit- Baltistan………………………………………………………. 9 1.4-Economic structure of Gilgit-Baltistan……………………………………………. 10 1.5-Central Karakoram National Park (CKNP)……………………………………….. 11 1.6-Management of Central Karakoram National Park……………………………….. 11 1.7-Climate of Gilgit-Baltistan………………………………………………………... 12 1.8-Ecological Zonation……………………………………………………………….. 15 1.8.2-Montane Dry Temperate Coniferous forests……………………………………. 15 1.8.3- Montane Dry Temperate Broadleaved forests………………………………….. 15 1.8.4- Sub-Alpine forests……………………………………………………………… 15 1.8.5-Northern Dry Scrub……………………………………………………………... 16 1.9-Land used in Gilgit-Baltistan……………………………………………………… 17 1.10-Medicinal Plants…………………………………………………………………. 17 1.11-Wildlife…………………………………………………………………………... 20 1.12-Water Resources…………………………………………………………………. 20 1.13-Forest Cover in Gilgit-Baltistan…………………………………………………. 21 1.13.1-Protected forests………………………………………………………… 21 VII 1.13.2-Private Forests…………………………………………………………... 22 1.14-Polices and legislation…………………………………………………………… 23 1.15-Problems and issues……………………………………………………………… 23 PART-1 ECOLOGY CHAPTER 2 REVIWE OF LITERATURE 2.1-Introduction……………………………………………………………………….. 26 2.2-Literature review …………………………………………………………………. 26 CHAPTER 3 STRUCTURE AND FUTURE TREND OF THE VEGETATION OF CKNP 3.1- Introduction………………………………………………………………………. 36 3.2-Materials and Methods……………………………………………………………. 38 3.2.1-Size class structure……………………………………………………………… 38 3.2.2-Weibull distribution…………………………………………………………….. 38 3.3-Results…………………………………………………………………………….. 40 3.3.1-Vegetation description and size class structure of each stand……………….. 40 A.FORESTED AREA 3.3.1.1-Stand No 1-Bagrot……………………………………………………… 40 3.3.1.2-Stand No 2-Haramosh…………………………………………………... 40 3.3.1.3-Stand No 3-Hopar………………………………………………………. 41 3.3.1.4-Stand No 4-Stak -1……………………………………………………… 41 3.3.1.5-Stand No 5-Stak -2……………………………………………………… 41 3.3.1.6-Stand No 6-Rakaposhi-1………………………………………………... 42 3.3.1.7-Stand No 6-Rakaposhi-1………………………………………………... 42 3.3.1.8-Stand No 8-Rakaposhi-3……………………………………………….. 42 3.3.1.9-Stand No 9 Rakaposhi- 4……………………………………………….. 43 B.NON-FORESTED AREA
3.3.1.10-Stand No 10-Bagrot…………………………………………………… 43 3.3.1.11-Stand No 11-Hopar……………………………………………………. 44 3.3.1.12-Stand No 12-Stak 1……………………………………………………. 44 3.3.1.13-Stand No 13-Stak-2……………………………………………………. 44 3.3.1.14-Stand No 14-Stak 3……………………………………………………. 45 VIII 3.3.1.15-Stand No 15-Thally 1………………………………………………….. 45 3.3.1.16-Stand No 16-Thally 2………………………………………………….. 46 3.3.1.17-Stand No 17-Kowardo………………………………………………… 46
3.3.1.18-Stand No 18-Arandu-1………………………………………………… 46 3.3.1.19-Stand No 19-Arandu- 2……………………………………………… 47 3.3.1.20-Stand No 20-Arandu- 3……………………………………………… 47 3.3.1.21-Stand No 21 Shigar 1………………………………………………….. 48 3.3.1.22-Stand No 22-Shigar 2………………………………………………….. 48 3.3.1.23-Stand No 23-Shimshal- 1-1…………………………………………… 48 3.3.1.24-Stand No 24-Shimshal- 1-2…………………………………………… 49 3.3.1.25-Stand No 25-Shimshal- 2-1…………………………………………… 49 3.3.1.26-Stand No 26-Shimshal- 2-2…………………………………………… 50 3.3.1.27-Stand No 27-Braldu -1-1……………………………………………… 50 3.3.1.28-Stand No 28-Braldu- 1-2……………………………………………… 51 3.3.1.29-Stand No 29-Braldu- 2-1……………………………………………… 51 3.3.1.30-Stand No 30-Braldu- 2-2……………………………………………… 52 3.3.1.31-Stand No 31-Chungo- 1……………………………………………….. 52 3.3.1.32-Stand No 32-Chungo- 2……………………………………………….. 53 3.3.2-Overall diameter distribution of the tree species…………………………….. 65 3.3.3-The Weibull function………………………………………………………… 67 3.3.4-Density ha-1 of trees and shrubs……………………………………………… 69 3.3.5-Basal area m2ha-1 of trees and shrubs………………………………………... 70 3.3.6-Correlation of density ha-1 with basal area m2 ha-1 and topographic factors with density ha-1 72 3.3.7-Correlation of topographic factors with basal area (m2 ha-1)………………… 80 3.4-Discussion and conclusion………………………………………………………… 85 3.4.1-Size class structure…………………………………………………………… 85 3.4.2-Overall diameter distribution of dominant species…………………………... 88 3.4.3-The Weibull function………………………………………………………… 89
IX CHAPTER 4 VEGETATION COMMUNITY ANALYSIS 4.1-Introduction ...... 93 4.2-Materials and methods…………………………………………………………….. 94 4.3-Results……………………………………………………………………………... 96 4.3.1-Forest Community and Pure Stands…………………………………………. 96 4.3.2-Shrubs and Herbs Communities……………………………………………... 96 4.3.1.1-Picea-Pinus wallichiana community………………………………………. 96 4.3.1.2-Juniperus excelsa Pure Stand……………………………………………… 97 4.3.1.3-Picea smithiana Pure Stand……………………………………………….. 98 4.3.1.4-Pinus wallichiana pure stand………………………………………………. 98 4.3.2.1-Rosa-Hippophae Community……………………………………………… 100 4.3.2.2-Hippophae-Ribes alpestre Community……………………………………. 100 4.3.2.3-Rosa-Ribes orientale Community………………………………………….. 101 4.3.2.4-Rosa-Berberis lycium Community………………………………………… 101 4.3.2.5-Hippophae-Tamarix indica Community…………………………………… 101 4.3.2.6-Berberis lycium-Tamarix indica Community……………………………… 102
MULTIVARIATE ANALYSIS CHAPTER 5 ORDINATION AND CLASSIDFICATION OF VEGETATION 5.1-Introduction………………………………………………………………………... 111 5.2-Purpose of the study ...... 112 5.3.-Materials and methods……………………………………………………………. 112 5.3.1-Data analysis…………………………………………………………………. 112 5.3.2-Soil analysis………………………………………………………………….. 113 5.3.3-Environmental variables……………………………………………………... 113 5.3.4-Ward’s clustering Method…………………………………………………… 113 5.3.5-DCA Ordination……………………………………………………………… 114 5.4-Results…………………………………………………………………………….. 114 5.4.1-Classification………………………………………………………………… 114 5.4.2-Ward’s cluster Analysis (Forested and non forested vegetation)……………. 114 5.4.2.1-Group I………………………………………………………………….. 114 X 5.4.2.2-Group II………………………………………………………………… 115 5.4.2.3-Group III………………………………………………………………... 116 5.4.2.4-Group IV………………………………………………………………... 116 5.4.3-Ward’s clustering Analysis (Understorey vegetation)……………………….. 123 5.4.3.1-Group I………………………………………………………………...... 123 5.4.3.2-Group II………………………………………………………………… 123 5.4.3.3-Group III………………………………………………………………... 124 5.4.4.4-Group IV………………………………………………………………... 124 5.4.4.5-Group V………………………………………………………………… 125 5.4.4.6Group VI………………………………………………………………… 126 5.4.5-DCA Ordination……………………………………………………………… 130 5.4.6-Ordination of forested and non-forested vegetation…………………………. 130 5.4.7-Correlation of ordination axes with environmental variables and soil nutrients (Forested and non forested vegetation)……………………………………… 134 5.4.8- Ordination of Understorey vegetation………………………………………. 136 5.4.9-Correlation of ordination axes with environmental variables and soil nutrients (Understorey vegetation)……………………………………………………. 140 5.4.10-Univariate analysis of variance (Forested and non-forested vegetation)…… 142 5.4.11-Univariate analysis of variance (Understorey vegetation)………………….. 145 5.4.12– Discussion and conclusion………………………………………………… 148 CHAPTER 6 VEGETATION-ENVIRONMENTS RELATIONSHIP 6.1-Introduction ...... 153 6.2-Objectives of the study……………………………………………………………. 154 6.3-Materials and methods……………………………………………………………. 154 6.3.1-Classification………………………………………………………………… 154 6.3.2-PCA Ordination……………………………………………………………… 155 6.4-Results……………………………………………………………………………... 155 6.4.1-Classification of forested and non-forested vegetation on the basis of 155 environmental variables……………………………………………………………….. 6.4.1.1-Group-I…………………………………………………………………. 155 6.4.1.2-Group-II………………………………………………………………… 156 6.4.1.3-Group-III………………………………………………………………... 156 6.4.1.4-Group-IV………………………………………………………………... 157 XI 6.4.2-Classification of Understorey vegetation on the basis of environmental 163 variables……………………………………………………………………………….. 6.4.2.1-Group-I…………………………………………………………………. 163 6.4.2.2-Group II A………………………………………………………………. 163 6.4.2.3-Group II B………………………………………………………………. 164 6.4.2.4-Group III………………………………………………………………... 165 6.4.2.5-Group IV A……………………………………………………………... 165 6.4.2.6-Group IV B……………………………………………………………... 166 6.4.3-PCA Ordination……………………………………………………………… 170 6.4.3.1-Ordination of forested and non-forested vegetation…………………..... 170 6.4.3.2-Ordination of Understorey vegetation………………………………….. 174 6.4.4-Correlation of ordination axes with environmental variables, edaphic factors and soil nutrients (Forested and non-forested vegetation)…………………………….. 178 6.4.5-Correlation of ordination axes with environmental variables, edaphic factors and soil nutrients (Understorey vegetation)…………………………………………… 180 6.4.6-Univariate analysis of variance (Forested and non-forested vegetation)…….. 182 6.4.7-Univariate analysis of variance (Understorey vegetation)…………………… 185 6.4.8-Physico-chemecal status of CKNP…………………………………………... 188
6.4.8.1-Elevation………………………………………………………………... 188 6.4.8.2-Slope……………………………………………………………………. 188 6.4.8.3-Conductivity……………………………………………………………. 188 6.4.8.4-Salanity…………………………………………………………………. 188 6.4.8.5-pH………………………………………………………………………. 189 6.4.8.6-MWHC…………………………………………………………………. 189 6.4.8.7-TDS……………………………………………………………………... 189 6.4.8.8-Organic matter………………………………………………………….. 189 6.4.8.9-Nitrogen………………………………………………………………… 189 6.4.8.10-Phosphorus…………………………………………………………….. 190 6.4.8.11-Potassium……………………………………………………………… 190 6.4.8.12-Calcium………………………………………………………………... 190 6.4.8.13-Magnesium……………………………………………………………. 190 6.4.8.14-Cobalt………………………………………………………………….. 190 6.4.8.15-Manganese…………………………………………………………….. 191 XII 6.4.8.16-Zinc……………………………………………………………………. 191 6.4.8.17-Iron…………………………………………………………………….. 191 6.4.8.18-Sulphur………………………………………………………………… 191 6.5-Discussion and conclusion………………………………………………….. 196 PART II DENDROCHRONOLOGY CHAPTER 7 AN INTRODUCTION TO DENDROCHRONOLOGY 7.1-Introduction ...... 201 7.2-Dendroclimatology………………………………………………………………... 201 7.3-Brief history of dendrochronology ...... 202 7.4-Principles of dendrochronology ...... 204 7.4.1-Uniformitarian ………………………………………………………………. 204 7.4.2-Limiting Factors……………………………………………………………… 205 7.4.3-Modeling growth-environmental relationship……………………………….. 205 7.4.4-Ecological Amplitude………………………………………………………... 205 7.4.5-Site Selection………………………………………………………………… 205 7.4.6-Crossdating…………………………………………………………………... 206 7.4.7-Replication…………………………………………………………………… 206 7.4.8-Sensitivity……………………………………………………………………. 206 7.4.9-Standardization………………………………………………………………. 206 7.4.10-Calibration and verification………………………………………………… 206 7.4.11-The domain of climate……………………………………………………… 207 7.5-Chracteristic features of Picea smithina tree rings…………………………….. 209 CHAPTER 8 LITERATURE REVIEW 8.1-Literature review………………………………………………………………….. 210 CHAPTER 9 MATERIALS AND METHODS 9.1-Site description……………………………………………………………………. 214 9.2-Field methods……………………………………………………………………… 214 9.3-Mouting and crossdating…………………………………………………………... 215 9.4-Measurement using Velmex ...... 216 9.5-Age and growth rates……………………………………………………………… 216 9.6-Chronology preparation…………………………………………………………… 216 XIII 9.6.1-COFECHA………………………………………………………………………. 216 9.6.2-ARSTAN ...... 217 9.6.3-Standard chronology…………………………………………………………….. 218 9.6.4-Residual chronology……………………………………………………………. 218 9.6.5-ARSTAN chronology…………………………………………………………… 218 9.7-Growth-climate response………………………………………………………….. 219 CHAPTER 10 AGE AND GROWTH RATES 10.1-Introduction………………………………………………………………………. 221 10.2-Materials and Methods…………………………………………………………... 222 10.3-Results……………………………………………………………………………. 222 10.3.1-Age and growth rates of seedlings………………………………………….. 222 10.3.2-Growth rate of past seedling………………………………………………... 225 10.3.3-Age and growth rate of tress………………………………………………... 227 10.4-Discussion and conclusion……………………………………………………. 230 10.4.1-Age and growth rates of seedlings………………………………………….. 230 10.4.2-Growth rate of past seedling……………………………………………….. 231 10.4.3-Age and growth rate of tress………………………………………………... 232 CHAPTER 11 CHRONOLOGY DEVELOPMENT 11.1-Introduction………………………………………………………………………. 234 11.2-Materials and Methods…………………………………………………………... 235 11.3-Results…………………………………………………………………………… 235 11.3.1-COFECHA statistics………………………………………………………... 235 11.3.2-Raw chronology……………………………………………………………. 237 11.3.3-Residual chronology………………………………………………………... 238 11.3.4-Standard chronology……………………………………………………….. 239 11.3.5-ARSTAN chronology………………………………………………………. 240 11.3.6-RBAR, EPS and sample depth……………………………………………… 242
XIV CHAPTER 12 GROWTH-CLIMATE RESPONSE 12.1-Introduction……………………………………………………………………… 246 12.2-Materials and Methods………………………………………………………….. 248 12.3-Results……………………………………………………………………………. 249 12.3.1-Correlation coefficients of residual chronology vs local climate ………….. 249 12.3.2-Response coefficients of residual chronology vs local climate…………….. 250 12.3.3-Correlation coefficients of residual chronology vs grid……………………. 251 12.3.4-Response coefficients of residual chronology vs grid……………………… 252 12.3.5-Correlation coefficients of standard chronology vs local climate………….. 253 12.3.6-Response coefficients of standard chronology vs local climate……………. 254 12.3.7-Correlation coefficients of standard chronology vs grid…………………… 255 12.3.8-Response coefficients of standard chronology vs grid……………………... 256 CHAPTER 13 GENERAL DISCUSSION AND CONCLUSION……………….. 263 REFRENCES………………………………………………………………………… 268 APPENDICES………………………………………………………………………... 298 PUBLICATIONS……………………………………………………………………. 320
XV LIST OF TABLES
Table 1.1 Land distribution in Gilgit Baltistan………………………………………... 17 Table 1.2 Some important medicinal gymnospermic plant species from CKNP……... 18 Table 1.3 Protected Forests cover in Gilgit-Baltistan………………………………… 21 Table 1.4 Trees volume in private forests in Chilas, Darel and Tangir………………. 22 Table 1.5 Overall status of private forests (ha)………………………………………. 22 Table 1. 6 Environmental characteristics of study sites of CKNP……………………. 25 Table 3.1 The Weibull function parameters of three dominant conifer species of 67 CKNP…………………………………………………………………………………. Table 3.2 Mean values of density ha-1, basal area m2 ha-1 and IVI…………………… 71 Table 3.3 Correlation of stand density ha-1 with stand basal area m2 ha-1, slope and 73 elevation with density ha-1…………………………………………………………….. Table 3.4 Correlation of species density ha-1 with species basal area m2 ha-1, slope 73 and elevation with density ha-1………………………………………………………… Table 3.5 Correlation of topographic (slope and elevation) factors with basal area m2 80 ha-1……………………………………………………………………………………... Table 4. 1 Communities of CKNP with IVI, absolute values and topographic range... 103 Table 4.2 Phytosociological attributes and absolute values of Forest, Bushes and Herbs from CKNP……………………………………………………………………... 104 Table 5.1 Four groups and one isolated stand obtained from Ward’s cluster analysis of forested and non forested species from 32 stands based on density ha-1 and environmental variable (elevation, slope)……………………………………………... 120 Table 5.2 showing understorey mean frequency of forested and non-forested vegetation …………………………………………………………………………….. 121 Table 5.3 Mean values of environmental variables (topographic and edaphic) and soil nutrients based on forested and non forested groups derived from Ward’s cluster analysis using 32stands of CKNP. (Mean ± SE)……………………………………… 122 Table 5.4 Showing means of groups of circular plot species (understorey vegetation) on the basis of frequency and environmental variables using Ward’s cluster analysis.. 128 Table 5.5 Mean values of environmental variables (topographic and edaphic) and soil nutrients based on circular plot groups derived from Ward’s cluster analysis using 32stands of CKNP. (Mean ± SE)……………………………………………….. 129 Table 5.6 Relationship (correlation coefficients) of environmental variables (topographic variables and edaphic variables) and soil nutrients with 3 DCA ordination axes obtained by forested and non forested vegetation data based on density ha-1…………………………………………………………………………….. 135 Table 5.7 Relationship (correlation coefficients) of environmental variables (topographic variables and edaphic variables) and soil nutrients with 3 DCA ordination axes obtained by understorey vegetation data based on frequency………... 141 Table 5.8 Analysis of variance of individual environmental variables (topographic and edaphic) and soil nutrients. Four groups were derived by Ward's cluster analysis using forested and non forested vegetation data of 32 stands of CKNP Gilgit- Baltistan, Pakistan……………………………………………………………………... 143 Table 5.9 Analysis of variance of individual environmental variables (topographic and edaphic and) and soil nutrients. Six groups were derived by Ward's cluster XVI analysis using understorey vegetation data of 32 stands of CKNP Gilgit-Baltistan, Pakistan………………………………………………………………………………... 146 Table 6.1 Four groups isolated from Ward’s cluster analysis of forested and non forested species from 32 stands based on environmental data and density ha-1 (elevation, slope), edaphic factors and soil nutrients………………………………….. 160 Table 6.2 Showing means of groups of forested and non forested species frequency on the basis of environmental variables (topographic and edaphic) and soil nutrients using Ward’s cluster analysis………………………………………………………….. 161 Table 6.3 Mean values ± SE of environmental variables (topographic, and edaphic) and soil nutrients based on forested and non forested groups derived from Ward’s cluster analysis using 32stands of CKNP. (Mean ± SE)………………………………. 162 Table 6.4 Showing means of groups of circular plot species (understorey vegetation) frequency on the basis of environmental variables and soil nutrients using Ward’s cluster analysis………………………………………………………………………… 168 Table 6.5 Mean values ± SE of environmental variables (topographic, and edaphic) and soil nutrients based on circular plot groups derived from Ward’s cluster analysis using 32stands of CKNP. (Mean ± SE)……………………………………………….. 169 Table 6.6 Relationship (correlation coefficients) of environmental variables (topographic variables and edaphic variables) and soil nutrients with 3 DCA ordination axes obtained by forested and non forested vegetation data based on density ha-1…………………………………………………………………………….. 179 Table 6.7 Relationship (correlation coefficients) of environmental variables (topographic variables and edaphic variables) and soil nutrients with 3 DCA ordination axes obtained by understorey vegetation data based on frequency………... 181 Table 6.8 Analysis of variance of individual environmental variables (topographic and edaphic) and soil nutrients of 32 stands from CKNP Gilgit-Baltistan, Pakistan…. 183 Table 6.9 Analysis of variance of individual environmental variables (topographic and edaphic and soil nutrients. Six groups were derived by Ward's cluster analysis using understorey vegetation data of 32 stands of CKNP Gilgit-Baltistan, Pakistan………………………………………………………………………………... 186 Table 6.10 Physical and chemical properties of the vegetation of CKNP…………... 195 Table 11.1 COFECHA Statistics of Picea smithiana from Stak Valley……………… 236 Table 11.2 Descriptive statistics of Raw, Standard, ARSTAN and Residual chorology of Picea smithiana…………………………………………………………. 241 Table 11.3 EPS, SNR and Rbar values……………………………………………….. 242 Table 12.1 Summary of various correlation and response function analysis using different chronologies with local climate and grid data. Only significant elements are shown………………………………………………………………………………….. 257
XVII LIST OF FIGURES
Fig. 1.1 Mean monthly temperature (Co) and precipitation (mm) of Gilgit and Skardu (1972-2011)…………………………………………………………………………… 14 Fig. 1.2 Map of the study area with sampling Locations……………………………… 24 Fig. 3.1 dbh size classes of forested and non forested vegetation of CKNP………….. 64 Fig.3.2 Overall dbh size class structure of dominant conifer tree species on the basis of mean density ha-1……………………………………………………………. 66 Fig 3.3 Generalized Weibull distribution models of dominant coniferous species…… 68 Fig. 3.4 Correlation of stand density ha-1 with stand basal area m2 ha-1 , slope and elevation with density ha-1…………………………………………………………….. 75 Fig. 3.6 Correlation of forested and non-forested species density ha-1 with species basal area m2 ha-1, slope and elevation with species density ha-1……………………... 79 Fig. 3.7 Correlation of forested and non-forested topographic factors (slope and elevation) with stand and species basal area m2 ha-1…………………………………... 84 Fig.5.2 Dendrogram resulting from cluster analysis based on frequency of understorey vegetation………………………………………………………………… 127 Fig.5.3 DCA ordination axes 1 and 2 of forested and non forested vegetation data based on density ha-1. The groups derived from Ward’s cluster analysis are not superimposed on 2-D ordination axes...... 132 Fig.5.4 DCA ordination axes 1 and 3 of forested and non forested vegetation data based on density ha-1. The groups derived from Ward’s cluster analysis are not superimposed on 2-D ordination axes…………………………………………………. 132 Fig.5.5 DCA ordination axes 2 and 3 of forested and non forested vegetation data based on density ha-1. The groups derived from Ward’s cluster analysis are superimposed on 2-D ordination axes…………………………………………………. 133 Fig.5.6 DCA among axes 1 and 2 of understorey vegetation data based on frequency. 137 Fig.5.7 DCA among axes 1 and 3 of understorey vegetation data based on frequency. 138 Fig.5.8 DCA among axes 2 and 3 of understorey vegetation data based on frequency (Not superimposed)...... 139 Fig.6.1 Dendrogram, based on Information level and Euclidean distance of the 32 stands of forested and non-forested environmental and vegetation data are presenting four groups…………………………………………………………………………….. 159 Fig.6.2 Dendrogram, based on Information level and Euclidean distance of the 32 stands of understorey vegetation………………………………………………………. 167 Fig.6.5 DCA ordination axes 2 and 3 of forested and non forested vegetation data based on density ha-1………………………………………………………………….. 173 Fig.6.6 DCA among axes 1 and 2 of understorey vegetation data based on frequency. 175 Fig.6.7 DCA among axes 1 and 3 of understorey vegetation data based on frequency. 176 Fig.6.8 DCA among axes 2 and 3 understorey vegetation data based on frequency….. 177 Fig.6.9 Whisker box plots of Physico-chemical factors of CKNP……………………. 194 Fig.9.1 Map of study area, circle shows sampling site (Stak valley of CKNP)………. 220 Fig.10.1 dbh vs age histograms analysis of Picea smithiana seedlings……………….. 223 XVIII Fig.10.2 dbh vs age regression analysis of Picea smithiana seedlings………………... 223 Fig.10.3 Age vs growth rates histogram analysis of Picea smithiana seedlings……… 224 Fig. 10.4 Age vs growth rates regression analysis of Picea smithiana seedlings……... 224 Fig.10.5 Growth rates of seedlings in various time periods…………………………... 226 Fig.10.6 dbh vs age histograms analysis of Picea smithiana…………………………. 227 Fig. 10.7 dbh vs age regression analysis of Picea smithiana...... 227 Fig. 10.8 Age vs growth rates histograms analysis of Picea smithiana……………… 228 Fig.10.9 Age vs growth rates regression analysis of Picea smithiana...... 228 Fig.11.1 Raw chronology plot of Picea smithiana……………………………………. 237 Fig.11. 2 Residual chronology plot of Picea smithiana………………………………. 238 Fig.11.3 Standard chronology plot of Picea smithiana……………………………….. 239 Fig.11.4 ARSTAN chronology plot of Picea smithiana……………………………… 240 Fig. 11.5 Graphs of rbar , EPS (Expressed population signal) and sample depth of Picea smithiana from Stak valley……………………………………………………... 243 Fig.12.1 Correlation coefficients of residual chronology vs local climate……………. 249 Fig.12.2 Response coefficients of residual chronology vs local climate……………… 250 Fig.12.3 Correlation coefficients of residual chronology vs grid……………………... 251 Fig.12.4 Response coefficients of residual chronology vs grid……………………….. 252 Fig.12.5 Correlation coefficients of standard chronology vs local climate…………… 253 Fig.12.6 Response coefficients of standard chronology vs local climate……………... 254 Fig.12.7 Correlation coefficients of standard chronology vs grid…………………….. 255 Fig.12.8 Response coefficients of standard chronology vs grid………………………. 256
1
2
Abstract
ABSTRACT This study focuses on quantitative community description, structure, multivariate analyses, vegetation-environment relationships and dendrochronological potential of Central Karakoram National Park Gilgit-Baltistan, Pakistan. Thirty two stands of forested and non-forested vegetation were sampled using point centered quarter method (PCQ) for trees while 3×5m quadrat were employed for non-forested vegetation sampling respectively. On the basis of phytosociological analysis, one mixed community, two pure conifer stands, a juniper forest and 6 non-forested communities were recognized. Advanced multivariate techniques including Ward’s agglomerative clustering method were used to seek group structure. For the forest vegetation groups were based on density ha-1 of tree species while the typification of understorey vegetation was performed on the basis of species frequency (%). DCA ordination was used to examine the distribution pattern of vegetation, to seek trends and gradients in vegetation. The groups derived from cluster analysis were superimposed on the ordination plane. Vegetation-environment relationship was investigated employing environmental variables including topographic factors, edaphic factors and soil nutrients using PCA ordination. It is observed that elevation and slope are significantly correlated with vegetation. Among the edaphic factors conductivity and soil pH were significantly correlated with vegetation while among the soil nutrients Nitrogen, Potassium, Calcium, Cobalt and Iron attains significant relation with vegetation. An attempt was also made to investigate structure and future trends of the vegetation of CKNP. In the forested vegetation Picea smithiana, Pinus wallichiana and Juniperus excelsa while in non-forested vegetation Rosa webbiana, Hippophae rhamnoides, Berberis lycium, Ribes orientale, Ribes alpestre and Tamarix indica exhibited high densities ha-1. Some gaps were observed in different size classes of these species which indicates that the size structures of these species were deteriorating with the passage of time. It may be due to overgrazing, illegal cutting, logging, soil erosion and competition. Most of the stands did not show the ideal situation and no inverse J- shaped curve was formed. Seven stands showed the positive skweness distribution, 5 stands attained flat distribution, 4 stands normal distribution,3 stands distributed in rectangular manner, 3 stands gave bimodal shape, 3 stands attained unimodal while the
3
Abstract
remaining stands distributed with U-shaped and leptokurtic shape. It is also noticed that these species have lesser recruitment which is alarming as the vegetation may completely disappear with the passage of time. Size class distribution of Picea smithiana, Pinus wallichiana and Juniperus excelsa was performed using Weibull probability function. It was observed that Weibull model gave a good fit for all three tree species examined. Age and growth rates of Picea smithiana trees and seedlings from Stak valley were investigated. It was observed that 90 cm dbh tree may attain 400 years age while growth rate ranges 3.9 to 17.4 year/cm in trees. In the seedling, the maximum age of seedling was 126 years while growth rates ranged 2 to 15.8 year/cm. Dendrochronological potential of Picea smithiana was also investigated using 22 corssdated wood samples out of 32 samples whereas remaining cores were rejected due to complacent ring. The results of master series showed that Picea smithiana from Stak valley covered highest age of 330 (1680-2009) years whereas average length of this forest was 156 years. The maximum length of two or more portions attained 230 years (1780-2009). The EPS, SNR and Rbar values were observed as 0.94, 18.06 and 0.68 respectively. An attempt was also made to check the relationship between growth and climate of Picea smithiana form Stak valley of CKNP. The results show that, in case of temperature, previous November and previous December were observed to be significantly positively correlated with tree growth. June, July temperature was also seen to be significantly positively correlated with tree growth. However, April temperature was negatively correlated with tree growth in both correlation and response function analysis. While in case of precipitation, tree ring indices showed significant positive relationship with April precipitation. The results show the suitability of Picea smithiana with respect to dendroclimatic potential. With extended sampling chronology can be extended to series more than five hundred years.
4
General introduction
GENERAL INTRODUCTION
My research objective covers two aspects of the study, first relates with the phytosociological studies while the second part is about the dendrochronological potential of tree species from CKNP. I discussed the work of these two studies separately in my review of literature. Central Karakorum National Park is located in Gilgit-Baltistan area which lies in the dry temperate zone of Pakistan. The altitude ranges from 2000m to 8000m with the coordinates at North 35° to 36.5 while at East 74°to 77 °. Due to the unique variety of flora and fauna, it is decaled as national park in 1993. The park is rich of beneficial plants s including medicinal and economically important species. In the study area three gymnospermic tree species included Picea smithiana, Pinus wallichiana and Juniperus excelsa while in non-forested species Rosa webbiana, Hippophae rhamnoides, Berberis lycium, Tamarix indica, Juniperus communis, Ribes orientale and Ribes alpestre were the common and dominant shrub species in the park. These species are great economic, environmental and medicinal importance. In this thesis an attempt was made to investigate the ecology of this forested and non-forested vegetation and also present the dendrochronological potential of conifer tree species of the national park. The first part of my thesis covers quantitative description of vegetation, soil analysis, multivariate analysis, future trend and structure of the vegetation while in the second part an endeavor was made to investigate age and growth rates, develop the chronology and growth- climate response of gymnospermic trees of the Stak valley of CKNP. Due to the importance and commercial value of above mentioned forests different NGOs conducted the qualitative study in this region. Therefore current studies present the first detailed quantitative data of this national park. This study may be useful for the management of forests, ecology and to understand vegetation environment relationship of this area. In addition dendrochronological investigation may provide growth climate response of the species.
5
Aims and objectives
AIMS AND OBJECTIVES The main purpose of my work is as under:
¾ To perform quantitative sampling in gymnospermic forests as well as herbs and shrubs species of the National Park. ¾ Quantitative description of the vegetation. ¾ To present the population structure and its future trend. ¾ To present classification ordination in order to expose the underlying group structure and measure gradient that controls the vegetation dynamics. ¾ Physico-chemical analysis of soil from the forested and non-forested vegetation of CKNP ¾ To estimate age and growth rates of conifer species, using Dendrochronological technique. ¾ To develop the standardized chronology of conifer species from (CKNP). ¾ To explore the tree growth-climate response.
6
Chapter 1: Description of the study area
CHAPTER-1 DESCRIPTION OF THE STUDY AREA 1.1-Geography of Gilgit-Baltistan Gilgit–Baltistan formerly known as the Northern Areas was renamed by national assembly on 29 August, 2009 as “Gilgit-Baltistan Empowerment and Self-Governance Order 2009”. Geographically it borders with Afghanistan in west, Xinjiang area of China in north, India in the east and Khyber Pakhtunkhwa province of Pakistan in the west , Afghanistan's Wakhan corridor to the north, China to the east and north- east, Azad Kashmir to the south-west and in the south-east Jammu and Kashmir . The territory of Gilgit-Baltistan consists of seven districts namely Gilgit, Ghizir, Diamar, Hunza Nagar, Astore, Skardu and Ghanche. Gilgit–Baltistan covers an area of 72,971 km² and is highly mountainous. It has an estimated population approximately 2 millions (Census 1998). Its administrative center is the Gilgit city. There are different languages are spoken in different regions like Shina, Balti, Brushushki, Khowar and Wakhi but the dominant language is Shina which is spoken in all seven districts of Gilgit-Baltistan while Urdu and English are also spoken. This area is unique in its culture, tourism, natural resources and topography that undeniably make the area glorious. The area is rich of flora and fauna and its alpine pastures are full of economical, medicinal and aromatic plants. Mountains of Gilgit-Baltistan are covered with snow which is the prime watershed of Indus River, providing the fresh water to the whole country. The culture and costumes of the Gilgit-Baltistan is more or less similar in all seven districts which includes food, music, dressing, folk dance, and traditional sports. The region of Gilgit-Baltistan is famous for hospitality which attracts the tourists. The tourists, nature lovers, and adventurers also attracted due to the presence of K-2 world’s 2nd highest mountain peak. After north and south Pole, the area holds the largest amount of snow. The glaciers are the main source of fresh water for drinking, agriculture and electric power generation. The region is also attractive due to the junction of world's three famed mountain ranges Karakorum, Hindukush and Himalaya. The unique diversity of flora, fauna, watersheds, tourism, pastures, glaciers, lakes, culture, minerals
7
Chapter 1: Description of the study area and mines attracts the researchers, visitors, tourists, trekkers, adventurers from all over the world. The areas is accessible both by bus and airplane from Islamabad to Gilgit- Baltistan 1.2-Karakoram Mountain Range The word “Karakoram” originated from the Turkish word “Kara-koram” which means “black rock”. The Karakoram mountain range begins from the Chilnji, south east of Wakhjir Pass which extends to the Ladakh beyond the Siachen glaciers to the extreme ,it touch the Tibet in China. The Karakoram Range is 300 miles long and 60 miles wide which originated from the Tethys Sea 50 million years ago due to the collision between Indian and Asian tectonic plates. The results of this collision, Indian plate penetrated the edge under the Asian plates. The evidence of that tectonic moment still felt in Nanga Parbat with 7mm Richter magnitude scale of frequent earthquakes and landslide. The positive effect of this collision is that hot springs in this range erupting which provides the dynamic geology to the Karakoram Range. There are seven famous glaciers are located in Karakoram Range namely, Batura, Hisper, Panamah, Baltoro, Siachen, Rimo and Saser which are collectively called Mustagh. The altitude of the mostly mountains are above 6000m while K-2 (Karakorum II,8611m) is the highest peak of this region along with 133 peaks above 7000m and hundreds more above 6000m which include Rakapohi, Hidden peak, Broad peak, Gasherbrum 1 to IV, Kiangshish, Masherbrum, Spatnik,. The slope of Karakoram Range is very steep with gigantic glaciers. Forty glaciers are found in the Karakorum range including famous glaciers which are Siachen, Hisper-Hopar, Batura, Baltoro, Biafo, Balafond, Godwin Austin Gondgoro, Chogholunggma, Chogolunsa, Chomaig, Gasheburum, Sapolago, Bozgil, Broad peak Glacier, Weyin, Wirjerab, and Khurdopin.
8
Chapter 1: Description of the study area
1.3-Geology of Gilgit- Baltistan
The geological history of the massive mountains of Gilgit-Baltistan is very ancient, some worlds highest peak group including Rakaposhi, Ultar, Diran, Broadpeak, Muztagh towers, Gasherbrum, Mashabrum, Baltoro, Trango Towers, Batura, Saltoro Kangri and many others peaks are found in Gilgit-baltistan (Trench, 1992). According to Stein (1987) these mountains covered with snow and glaciers and approximately 12 % of the regions formed by the glaciers which including Hispar 62 km, Siachen 62 km, Biafo – 61 km , Baltoro – 58 km , Batura – 58 km ,Gasherbrun – 38 km , Chungo Lungma – 38 km , Passu – 32 km , Nabandi – 32 km , Baraldu – 30 km ,Rupal 29 km and Snow lake sim glacier – 20 km , along with hundreds of other glaciers. The region contains mostly metamorphic rocks, due to the tectonic activity or magmatism in past, resulting in the creation of precious gemstones. The area is also a centre of tectonic activity where earthquakes and land sliding are common. Atabad Lake is the recent disaster of these activities which occurred in 4 January 2010. Other stories of geological activities i.e disasters of the earthquake, lightening flesh, floods, land sliding are frequently reported from the different regions of Pakistan in the past. The mountain ranges are rich for gem stones which including aquamarine, tourmaline, topaz, quartz etc. There are many types of gemstone mines in Gilgit-Baltistan around of 30 types are reported from this region. The famous sites for mining are Haramosh, Nagar, Shigar, Khaplu and Rundu. These valleys are the boundaries of CKNP where thousand of miners worked in above 2000 mines. These
9
Chapter 1: Description of the study area mines are the public property however, in some places localities formed a society which decided to work in mining. The gem stone of the region is famous nationally and internationally due to the commence of divers tourist. Although these gemstones benefit for the localities and the economy of the country but there are some health, environmental and safety issues are also associated. It is scientific fact that silicosis and carbon monoxide affect the lungs which are the lethal agents for human beings. On the other hand wildlife of the area is directly or indirectly affected by mining activates. There are several problems are safety cannot be denied because the miners are not well trained and many causalities are reported from the region. 1.4-Economic structure of Gilgit-Baltistan The economy of Gilgit Baltistan is endurance of extensive and diverse. The area depends upon agriculture, forestry, tourism, gem stones and mineral resources. The unique beauty of this area attracts the thousand of national and international tourists. The region is rich for precious gem stones and minerals which including Nickel, Cobalt, Copper, Lead, Tin, Bismuth, Mica, Quartz, Zircon, Coal and Actinolite that are famous for their good quality and economic importance directly or indirectly polished and converted into jewelry form. These gemstones are introduced and supplied to the national and international markets. Some mineral resources of Iron, Silver, Gold, Zinc, Marble, Granite, Sulpher, Calcite, Fluorite, Lime Stone, Arsenic, Spinal, Garnet, Epidote, Topaz, Moon Stone, Pargasite, Tourmaline, Aquamarine, Pyrite and feldspar are also reported from the region by Hussainabadi, (2003). Some cultural goods, fruits, vegetables are also the source of income. The other source of the high income is hotel industry which rapidly increases due to the scenic beauty of Gilgit-Baltistan convinced the national and foreign visitors, researchers and adventurers. Although deforestation is a bad practice but still going on, in various localities. The wood of Gilgit-Baltistan forest is very important both commercially and economically (Hussain, 2010).The Local people cut the forests for the purpose of fuel, shelter and to earn money from the market. It is concluded that the major sources of income in this region are tourism, mines and minerals, agriculture, vegetables, fruit farms, forestry and hotel industry.
10
Chapter 1: Description of the study area
1.5-Central Karakoram National Park (CKNP) Central Karakoram National Park (CKNP) is located in Gilgit-Baltistan of Pakistan. It falls into the administrative boundaries of district Gilgit, Skardu and Ghanche. It is one of the 24 national parks of Pakistan. Because of its unique and diverse habitat of flora and fauna, it was declared as National Park in the year 1993.It is the biggest National park in Pakistan, covering an area of 10,000 km2 .The CKNP extends from 35°N to 36.5°N latitude and from 74°E to 77°E longitude. The elevation ranges from 2000m-8000m. Sixty peaks over 7,000 m, and ten of the world highest and most famous mountains, including Gasherbrum, Broad Peak, and Masherbrum are located within the park' boundary and the world’s 2nd highest peak K-2 is also situated in this park. According to the 1998 census there are 350 settlements (permanent and non- permanent) and total population of these settlements is approximately 211,000. The fresh water source is mainly glaciers. The largest glaciers i.e Siachen (75km long), Baltoro (57km), and Hispur-Biafo (122km) are originating within the park. The main valleys included in the park are Bagrot, Haramosh, Hopar-Hisper, Shimshal I, Shimshal II, Braldu I, Braldu II, Hushe, Thalley, Stak-Tormik, Shigar, Nar, Kowardo, Chungo, and Basha (Arandu). 1.6-Management of Central Karakoram National Park National Parks in Gilgit-Baltistan are the responsibility of the Forest Department falling under the Directorate of Parks and Wildlife. The Park is governed by the 1975 Northern Areas Wildlife Preservation Act. Though declared National Park in 1993 there is neither physical indication nor is public awareness on the ground to suggest the existence of CKNP. Currently CKNP Directorate is developed and director of CKNP is the responsible for all matters collaboration with different institution involved in and around the park includes the Aga Khan Rural Support Program (AKRSP), IUCN, and the World Wildlife Fund. However there are a number of smaller NGOs working in or adjacent to the park with a range of conservation projects.
11
Chapter 1: Description of the study area
1.7-Climate of Gilgit-Baltistan Due to the influence of higher mountains, rain shadows in some laces and high precipitation in other regions present. The Karakorum and Hindukush create barriers between monsoon lands of South-Asia to south and west deserts of the Central-Asia to their north. Therefore in the summer season this area is influence by the residue of the monsoon system coming from south and during the winter and spring it is influenced by westerly depressions originating in the Mediterranean and Caspian Sea. The climate is varies between lowlands and valleys. The valleys are dry with the annual precipitation around 200 mm but totals can go up as high as 600 mm at elevations of 13,000 ft (Kreutzmann, 2006, 2000a and 2000b). In the month of July and August the temperature of Gilgit-Baltistan is maximum around 20-35 Co in the valleys while in cold arid regions in January with an average of -10 to 00 Co. Archer (2010) found that both Gilgit and Skardu annual rainfall has increases for all seasons. The annual one day maximum rainfall has risen from 12 to 30 mm at Skardu and from 9 to 28mm at Gilgit. Archer (2003 and 2004) suggested that the changes in precipitation at high elevation can be interfered by increment of ice cores and snowfall.
12
Chapter 1: Description of the study area
Monthly Mean Pricipitation (mm) of Gilgit (1972-2011) 30.0 25.0 20.0 (mm) 15.0 10.0 5.0 0.0 Precipitation Jan Feb Mar Apr May Jun Jul Agu Sep Oct Nov Dec Months
Monthly Mean Maximum Temperature of Gilgit (1972-2011) 40.0 )
30.0 C
(
20.0
10.0
0.0 Temperature Jan Feb Mar Apr May Jun Jul Agu Sep Oct Nov Dec Months
Mean Monthly Minimum Temperature of Gilgit (1972-2011)
20
15 )
(c 10 5 0 Temperature ‐5 JanFebMar Apr May Jun Jul Agu Sep Oct Nov Dec Months
13
Chapter 1: Description of the study area
Monthly Mean Pricipitation (mm) of Skardu (1972-2011)
50 40 30 (mm)
20 10 0 Pricipitation Jan Feb Mar Apr May Jun Jul Agu Sep Oct Nov Dec Months
Mean Monthly Maximum Temperature of Skardu (1972-2011) 40
30 )
c
(
20
10
0
Temperature JanFebMar Apr May Jun Jul Agu Sep Oct Nov Dec Months
Mean Monthly Minimum Temperature of Skardu (1972-2011) 20 ) 10 (c
0
‐10
Temperature Jan Feb Mar Apr May Jun Jul Agu Sep Oct Nov Dec Months
Fig. 1.1 Mean monthly temperature (Co) and precipitation (mm) of Gilgit and Skardu (1972-2011)
Data source: Pakistan Metrological Department
14
Chapter 1: Description of the study area
1.8-Ecological Zonation According to Rao and Marwat (2003) on the basis of ecological Zonation five main types of forests are existed in Northern Areas i.e Montane Sub-Tropical Scrub, Montane Dry Temperate Coniferous, Montane Dry Temperate Broadleaved, Sub-Alpine and Northern Dry Scrub. 1.8.1- Montane Sub-Tropical Scrub forests These forests are found between the elevation 750 to 3,900m a.s.l. These forests divided into Montane sub-tropical scrub comprises of Capparis spinosa, Pisticia, Artemisia, Saccharum, Dodonaea, Berberis, Rosa moschata and Daphne oloides. This area covered along the Indus River up to Raikot and Bunji. 1.8.2-Montane Dry Temperate Coniferous forests This Zonation is consisted by dry deodar (Cedrus deodara), blue pine (Pinus wallichiana), fir (Abies spectabilis), spruce (Picea smithiana), chilghoza (Pinus gerardiana) and juniper (Juniperus spp.). These important species are found in district Diamer, some parts of districts Gilgit, Skardu and two villages in Ghizer district. 1.8.3- Montane Dry Temperate Broadleaved forests These forests contain the species oak (Quercus ilex), ash (Fraxinus spp.), poplar (Populus), willow (Salix) and Artemisia. 1.8.4- Sub-Alpine forests This zone consists of snowfall up to 3 m/year with slight rainfall. In this zone plant species included birch, willow, and juniper, Ephedra, Vibernum, Andropogon, Berberis, Lonicera and Ribes. 1.8.5-Northern Dry Scrub This is the zonation of scattered scrub vegetation. In this zone, important and common species are included sea buckthorn (Hippophae rhamnoides) Willow (Salix) and stunted Juniper trees also grow on the hill sides.
15
Chapter 1: Description of the study area
Mix forest at Bagrot valley Picea smithiana forest at Haramosh Valley
Shrub communities at Thally Valley A mix-forest view of Stak valley
1.9-Land used in Gilgit-Baltistan As we discussed that Gilgit-Baltistan is covered of high mountains which are rich of glaciers therefore melting of glaciers would affect the land use patterns around the glacier, e.g., agriculture, forestry and cultivated areas and also affect the habitat of plants and animals. The classification of land use in Gilgit-Baltistan is as under:
16
Chapter 1: Description of the study area
Table 1.1 Land distributions in Gilgit Baltistan
Type Area (ha) Percentage (%) Mountain/ Lakes/ Rivers/ 4,810 66 Glaciers Forest 646 09 Rangeland 1,646 23 Cultivated Area 58 1 Cultivable waste 90 1 Total 7,250 100 Source: Forest department of Gilgit-Baltistan
17
Chapter 1: Description of the study area
1.10-Medicinal Plants Gilgit-Baltistan is rich of a wide variety of plants which have high medicinal and economic importance. Wali and Khatoon (2007) reported some important gymnospermic tree species are medicinally used by the local communities (Table 1.2).
Table 1.2 Some important medicinal gymnospermic plant species from CKNP
S.No Name of species Local Name Purpose of usage 1 Juniperus communis L. Mitthary Kidney stone, Urine problem, Leucorrhea and Tuberculosis 2 Juniperus excelsa. M.Bieb Cheleh Kidney stone, Urine problem, Weakness of urinary bladder 3 Pinus wallichiana A.B.Jackson Cheenh Resin is used against wound 4 Pinus gerardiana Wall.ex. Lamb Cheenh Resin is used against wound 5 Picea smithiana (Wall). Bois. Kachul Heart problems 6 J.turkestanica Komarov Cheenh Kidney stone, Urine problem,Weakness of urinary bladder 7 Ephedra gerardiana Wall.ex Stapf Soom Asthma, Cough and other respiratory problems 8 E.intermedia Schrenk & Meyer Shaay soom Wound, Asthma, Rheumatism and Gouts Source: Wali and Khatoon (2007).
Some other plant species are also practiced medicinally such as Betula utilis, Berberis lycium, Hippophae rhamnoides, Ribes alpester, Ribes orientale, Rosa webbiana, Rubus irritans, Spiraea canescens etc.
18
Chapter 1: Description of the study area
Some dominant shrubs and herbs of CKNP
Ribes orientale Rosa webbiana
Hippophae rhamnoides Spiraea canescens
Thymus linearis Taraxacum sp.
Artemisia brevifolium Anaphalis virgata
19
Chapter 1: Description of the study area
1.11-Wildlife According to Virk et al., (2003) the Gilgit-Baltistan has one of the most diverse avi fauna of the mountain region of the world. But little information is available on the distribution, status, diversity and ecology of many of these bird species. The most comprehensive account of the avifauna of Pakistan comes from Robert (1991). In spite of ruthless hunting wildlife in Gilgit-Baltistan supports a wide variety of endangered species of mammals and birds like Marco Polo sheep, blue sheep, markhor, black bear, brown bear, chakor and ram chakor. Fifty-four mammal species are reported from Gilgit-Baltistan. These species consist of one shrew, 10 bats, 18 carnivores, 6 artiodactyls, 3 lagomorphs, and 16 rodents. Upper Hunza and the triangle between Indus and Astore rivers are considered the “Hot Spot” for large animal diversity including snow leopard, Marco Polo sheep, Himalayan brown bear, black bear, musk deer, flare horned markhor, Laddakh urial, blue sheep, and Himalayan lynx. Other species like wolf (Canis lupis), Urial (Ovis orientalis) and Marmot (Marmota caudata and Marmota bobock) also reported in different areas of Gilgit-Baltistan. 1.12-Water Resources As we know that glaciers are asset of the country and the mountains of the Gilgit- Baltistan are covered by glaciers and snow. These glaciers and the snow are the main source of water for Gilgit-Baltistan and Pakistan. Irrigation is also essential in this region because of meager precipitation in the valleys. The peoples of Gilgit-Baltistan used water from the water nallahs which fed by glaciers and snow melt water for their drinking and irrigation purposes. The water is contaminated due to the interruption of human and animal feces. Different NGOs and GBPWD are working on to purify the water. Inspite of these precautions there is still an immense need to improve drinking water quality in Gilgit-Baltistan.
20
Chapter 1: Description of the study area
1.13-Forest Cover in Gilgit-Baltistan Forest in Gilgit-Baltistan is divided into two categories i.e protected and private forest. 1.13.1-Protected forests Those forests which are ownership to the Gilgit-Baltistan Government have been declared as “protected forests” under the Pakistan forest Act (1927). The other category of forests is “private forests” which are owned by the local communities. These forests are also legally enclosed under the Gilgit Private Forests Regulations (1970). According to the Act (1976) protected forests are the property of the government. However local communities may have used these forests for grazing and fuel wood. These protected forests are found in Gilgit, Astore and Baltistan and regularized under the Northern Areas Forest Rules 1983. According to Ali (2004) the total area of protected forests in Gilgit-Baltistan is 64,512 ha with the dominant species including spruce (Picea smithiana), silver fir (Abies pindrow), blue pine (Pinus wallichiana), Juniper (Juniperus excels). The total shrubs cover is 381,200 ha.
Table 1.3 Protected Forests cover in Gilgit-Baltistan.
District Area (ha) Forest Types Gilgit and Nagar 17,028 Montane Dry Temperate and Sub- alpine Ghizer (Punial) 7,740 Montane Dry Temperate and Sub- alpine Astore 30,960 Montane Dry Temperate and Sub- alpine Baltistan 9,288 Montane Dry Temperate and Sub- alpine Total 65016
Source: Rao (2003).
1.13.2-Private Forests Private forests are the property of local people which are regularized under the Gilgit Private Forests Regulations, 1970 and the rules notified in 1975. In an agreement with the Accession Deed of 1952 with the government of Pakistan the tribal communities of Chilas, Darel and Tangir in Diamer district own the private forests of the Gilgit-Baltistan. The local communities agreed that the management of these forest
21
Chapter 1: Description of the study area by the Northern Areas Forest Department (NAFD). This forest department has been helping in harvesting of private forest, working schemes and to resolve the problems. The Local communities have a 100% ownership and the government gets 50% royalties from the revenue (Ali, 2004).
Table 1.4 Trees volume in private forests in Chilas, Darel and Tangir
Period Total trees (million) Total volume (m3) Volume to be harvested (m3) 2002-2013 5.51 9,169,924 1,281,507
Table 1.5 Overall status of private forests (ha)
Govt. protected forests in Gilgit, Skardu, Ghizer & Astore 65,000 Private forests (govt. managed) in Chilas, Darel and Tangir 2,17,000 Private owned farmland forestry in seven districts Gilgit-Baltistan 3,25,000 Total 6,06,000
Source of Tables 2 and 3: Rao (2003).
22
Chapter 1: Description of the study area
1.14-Polices and legislation The first forest policy of Pakistan was declared in 1955 and updated in 1992, 1975, 1980 and 1988 as part of the National Agricultural Policy. In 1991 National forest policy revised and integrated formal policy on the federal level (Ahmed and Mahmood, 1998). The purposes of these polices was to protect and the productivity of the forest. In 2001 a draft has been prepared by the National Forest Policy which is also applicable for Northern Areas. In this report it is stated that use of natural resources with community participation was allowed and recommended that timber harvesting can be used for poverty alleviation.
1.15-Problems and Issues Although there are the legislation and policies to protect the forests, however some of the following problems and issues are common in the forest areas: ¾ Lack of community collaboration with government departments. ¾ Insufficient funds for the management and protection of forests. ¾ Illegal cutting and overgrazing. ¾ Lack of awareness in local communities. ¾ Inadequate research approach. ¾ Lack of skilled staff or professional in the forest department. ¾ No alternative jobs for forest people. ¾ No alternative fuel except wood. ¾ Poverty. ¾ Lack of education.
23
Chapter 1: Description of the study area
Fig. 1.2 Map of the study area with sampling Locations
24
Chapter 1: Description of the study area
Table 1. 6 Environmental characteristics of study sites of CKNP
Stands Main Location Lat.(N) Long.(E) Ele.(m) Asp. Slo. CN Forest areas 1 Bagrot 36.02918 74.60156 3130 E 45° Mod. 2 Haramosh 35.88388 74.88433 3296 E/S 53° Ope 3 Hoper 36.14278 74.94410 3486 E 49° Clo 4 Stak 1 35.75901 75.05340 3344 E 35° Mod 5 Stak 2 35.77398 75.04300 3600 E 20° Mod 6 Rakaposhi 1 36.12485 74.94081 3444 N 70° Mod 7 Rakaposhi 2 36.18000 74.66000 3263 N 59° Mod 8 Rakaposhi 3 36.04021 74.54186 3188 N 64° Mod 9 Rakaposhi 4 36.15703 74.92910 3512 N/E 70° Mod Shrubs/Herbs area 10 Bagrot 36.03400 74.57735 27774 N Pla. Ope 11 Hopar 36.16258 74.84320 3353 N/E 30° Ope 12 Stak 01 35.73915 75.10943 2949 E/N 35° Ope 13 Stak 02 35.74459 75.05935 2782 E/S 20° Ope 14 Stak 03 35.74053 75.05651 2742 E Pla Ope 15 Thallay 01 35.17268 76.33680 3300 E 20° Ope 16 Thallay 02 35.17575 76.33440 3500 E/N 25° Ope 17 Kowardo 35.40611 75.60833 3559 E 50° Ope 18 Arandu 1 35.93333 75.70375 2790 S/W 30° Ope 19 Arandu 2 35.83333 75.73858 2815 S 25° Ope 20 Arandu 3 35.79565 75.73868 2875 S/W 35° Ope 21 Shigar 1 35.68970 75.88908 2527 N/E 40° Ope 22 Shigar 2 35.72278 75.79650 2444 E 35° Ope 23 Shimshal 1-1 36.73260 75.53743 3047 E/S Pla Ope 24 Shimshal 1-2 36.73423 75.54813 3065 E/S Pla Ope 25 Shimshal 2-1 36.73228 75.55151 3076 E/S Pla Ope 26 Shimshal 2-2 36.72853 75.55553 3097 E/S Pla Ope 27 Braldu 1-1 35.67020 75.76706 2895 E 25° Ope 28 Braldu 1-2 35.70150 75.75516 2910 E 20° Ope 29 Braldu 2-1 35.70516 75.75353 2948 E/S 35° Ope 30 Braldu 2-2 35.71848 75.85377 3055 E 30° Ope 31 Chungo 1 35.81715 75.68707 3010 N 40° Ope 32 Chungo 2 35.81808 75.92186 3109 N/E 35° Ope Lat = Latitude, Long = Longitude, Ele = Elevation, Asp = Aspect, Slo = Degree of slope, CN = Canopy, Mod = Moderate, Ope = Open, Clo = Closed, CS = Covered surface, Pla = Plain.
25
PARTI ECOLOGY
Chapter 2: Review of literature
CHAPTER- 2 REVIEW OF LITERATURE 2.1-Introduction This chapter focuses on the brief literature review from the different parts of Pakistan regarding observational studies, quantitative studies, multivariate analysis, and vegetation relationship with environmental variables and physico-chemical analysis. 2.2-Literature review Many researchers conducted the observational studied from different sites of Pakistan (Chaudhri, 1952, 1953, 1960; Khan, 1955 & 1960). Hussain (1960) carried out the vegetation survey of Ayub National park. Stewart (1961) worked on the flora of Deosai plains. Jafri (1962) classify the vegetation of Bolan Pass. Khan and Repp (1961) described the Riverian forest and irrigated plantation of Pakistan while Nasir and Webster (1965) worked on the vegetation of Hushe valley, Baltistan. Similarly Rafi (1965 and 1973) presented the vegetation of Baluchistan. Champion et al., (1965) conducted a detail survey of different parts of Pakistan, they also described the vegetation of different climatic zones (i.e subtropical, moist temperate, dry temperate and sub-alpine) of Northern Pakistan. Rep and Beg (1966) studied the Juniper forest of Ziarat Balochistan. Shaukat and Qadir (1969) quantitatively describe the vegetation of calcareous hills around Karachi and they described the lithosere succession. Hussain (1969) studied the relationship between vegetation and edaphic factors from Wah Garden Cambellpur District Rawalpendi. Hussain and Qadri (1970) investigated the relationship between the plants species growth and quality of the physical and chemical properties from the surrounding areas of Karachi. They found correlation of plants density with soil texture and CaCo3.Shaukat and Qadir (1971) presented Multivariate analysis of the vegetation of calcareous hill around Karachi, and stated that potential continuity in vegetation with the aid of an indirect gradient analysis. Shaukat and Hussain (1972) analyzed the vegetation of Khadeji-fall area and described the hydrosere and lithosere succession and also studied quantitatively four stands representing different aspects of the hills. Their results revealed that the vegetation of the hills is principally constituted by perennial xerophytic shrubs like Commiphora
26
Chapter 2: Review of literature
wightii, Grewia tenax, Euphorbia coducifolia and Acacia senegal. The distribution of vegetation types was correlated with obvious physiographic factors. Malik et al., (1973) studied in the forests of Malakand Division and observed the nutrients condition of representative soil profiles, minerals compositions Khan and Ahmed (1976) quantitatively studied the medicinal plants in Rawalpindi north, Rawalpindi south and Muree forest division and described the ethnobotanical importance of these plants. Ahmed and Qadir (1976) carried out quantitative phytosociological survey along the way of Gilgit to Gopis, Yasin and Phunder. They sampled 46 different locations and recognized ten communities on the basis on physiognomy, floristic composition and importance value index. The attributes, maturity and homogeneity of every stand were also studied. First stand ordination of Skardu Northern area was presented by Ahmed (1976). Shaukat et al., (1976) presented Phytosociological study of Gadap area, Southern Sindh. Twenty two stands in Gadap area of Southern Sindh were sampled quantitatively and soil samples were analyzed physico-chemically. They studied vegetation composition and structure, dominant group of vegetation, diversity relations of leading dominant groups, objective classification and correlation of edaphic variables with the vegetation types. Ahmed (1973) and Ahmed et al., (1978) presented multivariate approaches to the analysis of the vegetation-environmental complex of Gharo, Dhabeji and Manghopir industrial areas. They studied the vegetation characteristics, vegetation ordination, relationships between environmental variables and stand composition, species ordination, the behavior of species along environmental gradients. Chaghtai et al., (1983a, b) described the ecology of a dry stream bed in Peshawar and phytosociology of the Muslim graveyards of Kohat division respectively. According to them the variety of habitats depends on the availability of moisture and the extent of biotic disturbance. The phytosociology of the Muslim graveyards of Kohat explained that the numbers of species
in a stand were controlled by the sand and CaCo3 proportions of the soil. Ahmed (1986) presented vegetation of some foothills of Himalayan range in Pakistan. Six communities have been recognized on the basis of species dominance at 17 locations near road side on the great Silk Road from Gilgit to Passu. Chagthai et al., (1987) gave a detail description
of soil factors i.e TDS, Sodium, Potassium and Calcium, NO3, PO4 and CaCo3 with
27
Chapter 2: Review of literature
vegetation. According to them CaCo3 is an important factor to control the supremacy of species. Tareen and Qadir (1987) reported sixteen plant communities on the basis of important index value in the plains of Quetta district and also they correlate soil with type of vegetation and concluded that total coverage as well as species diversity tended to be relatively high in protected areas and graveyards. Ahmed (1988a, b) studied plant communities of some northern temperate forests of Pakistan. Chaghtai et al., (1988) presented ecology of an upland forest near Noshehra, NWFP. They observed that the lower valley slope was dominated by arboreal vegetation, the middle by tall shrubs and the top exposed by grasses. Kayani et al., (1988) also studied the relationship between soil properties and vegetation types. Ahmed et al., (1989) described the natural regeneration potential of Juniperus excelsa Balochistan. They also studied the regenerating seedlings from 60 mature stands on Juniper track ranged from zero to 219 ha-1 with a mean of 52 ha-1. Seedling density and basal area were significantly correlated (P<0.001) while tree basal area and seedling density were also significantly correlated (P<0.05), indicating that seedlings are sciophytic and are found under the shade of trees. Height average seedling density and basal area were recorded from west facing slope, in addition, future trend of the seedling population suggesting that Juniper forest are not deteriorating. Shaukat & Uddin (1989) investigated the vegetation-environmental data set of Gadap area using Canonical Correlation Analysis (CCA). They also provided three new alternative models that could be used when singular value decomposition algorithm of CCA failed to provide any worthwhile results. Qadir et al., (1989) presented phytosociology of woodland communities of Hazarganji National Park Quetta. According to them the soil of the communities was coarse-textured, calcareous and non –saline. Ahmed et al., (1990) described population structure and dynamics of Juniperus excelsa in Balouchistan. Ahmed (1991) et al., described vegetation structure and dynamics of Pinus gerardiana forests in Balochistan and also noticed that diameter distributions within the stands were mostly skewed and unimodal with gaps appearing in large size classes. Cross sections of 16 trees were also used to determine age and growth rate. Moreover, adequate recruitment of Pinus seedlings was observed. The average growth rate was estimated 0.08 cm/year radial growth. However, trees on high elevations
28
Chapter 2: Review of literature
and cooler slopes grow faster. They analyzed soil variables which showed no correlation with density, basal area or importance values. They also suggested that the present degraded stage of the forests in the study area is of anthropogenic origin. Tareen and Qadir (1991) presented Phytosociology of the Hills of Quetta District. Fifty seven plant communities were recognized. Out of forty five communities were grouped into 8 steppe type viz. Twenty eight species were found as indicator species of specific soil condition. Hussain et al., (1993) described phytosociology of the vanishing tropical dry deciduous forests in district Swabi. According to him the variation in the dominant species was due to the edaphic and biotic disturbance. He suggested that the existing vegetation might further change due to underground seepage of water from nearby Tarbela dam. Rasool (1994) described the status and management polices of the protected areas included forested vegetation in the Northern Areas of Pakistan. Hussain et al., (1994) presented phytosociology of the vanishing tropical deciduous forests in district Swabi. The study deals with the multivariate analysis of vegetation of Swabi District. He examined the physico-chemical analysis of soil with the help of polar ordination. Soil pH, CaCo3 and P2O5 were found as the controlling factors in the distribution of vegetation. Hussain and Mustafa (1995) reported ecological studies of plants in relation to animal, found in Nasirabad valley Hunza Pakistan. Malik et al., (1994) investigated the vegetation of Samani and Dhirkot Hills and stated that heavy illegal logging, soil erosion and overgrazing has led to degradation of forests into scrubs. Ashraf, (1995) focused on the phytosociology of the vegetation in Pir- Chinasi hills and he recognized ten different communities. Among the tree species Pinus wallichiana and Abies pindrow were found to be the dominant species, while understorey species consisted of Viburnum grandiflorum, Indigofera, Elsholtzia, Sorbaria and Sibbaldia. Pinus wallichiana appeared to be the only regenerating species while Abies pindrow, Picea smithiana and Juglans regia did not regenerate. Malik and Zandiyal, (1996) also reported ten plant communities in Machyara Hills and declared that Cedrus deodara, Pinus wallichiana, Abies pindrow and Bistorta were the most common species of the vegetation. Rasool (1998) had provided a detail account of the economically important plants of northern areas. Alpine deserts have little values as grazing lands due to the absence of forage and difficult topography. Alpine pastures were subjected to heavy
29
Chapter 2: Review of literature
grazing during summer. Siddiqui et al., (1999) investigated the climate change impact and adaptation approaches of forest ecosystem in Pakistan. Tareen and Qadir (2000) carried out the soil organic matter, MWHC, CaCo3, conductivity, Bicarbonate, Chloride, Calcium, Magnesium, Sodium and Potassium from the plains Harnai, Sinjawi to Duki region of Pakistan. Ahmed and Khattak (2001) studied the vegetation of Islamabad and they recognized 17 plant communities on the basis of highest importance value. Acacia modesta was the dominant tree followed by Broussonetia papyrifera and Dalbergia sissoo respectively. Jaffari et al., (2003) studied the relationship of chemical and physical properties of soil with vegetation type. They used PCA and CCA methods of ordination to determine the soil factors included acidity (pH), texture, electrical conductivity, and - lime, Calcium, Magnesium, Sodium, Potassium, Cl , CO3 and HCO3Gilani (2003) surveyed the Astore area to provide information on the conservation of plant diversity. Malik and Malik, (2004) recognized seven plant communities and fifty eight plant species from Kotli hills and also reported the deforestation and overgrazing were common in the study area. Ali et al., (2005) described in detailed about deforestation and its causes from the forest of Basho valley Skardu. Shaukat et al., (2005) compared the techniques CA, DCA and CCA (canonical correspondence analysis) using a field data set from Lower Sindh and found that of the three techniques CCA was most useful in exposing the underlying correlation structure between vegetation and environment. Hussain et al., (2006) also recognized 28 different plant communities of trees, shrubs and herbs vegetation. Ahmed et al., (2006) presented phytosociology and structure of Himalayan forest from different climatic zones of Pakistan. A quantitative phytosociological survey was conducted in 184 sampling stands in various climatic zones of Himalayan forest of Pakistan. Based on floristic composition and importance value index, 24 different communities and 4 monospecific forests were recognized. The quantitative description and their population structure were also presented. Most of the communities showed similar floristic composition however they were different in quantitative values. Parveen and Hussain (2007) conducted a study for the plant biodiversity and phytosociological attributes of the Gorakh hill and they provide inclusive inventory of the area. Shahbaz et al., (2007) gave a detailed description about
30
Chapter 2: Review of literature
the forest policies and implementation on forests of Pakistan. Malik (2007) described the phytosociological attributes of different plant communities of Pir Chinasi Hills Azad Jammu Kashmir. However, Dasti et al., (2007) also investigated vegetation composition and multivariate analysis on the Pothwar Plateau and recognized five plant communities the basis of cluster analysis. He also noticed that the application of the classification result may be used to ordination in terms of topography, redistribution of rain water, the nature of bedrock and soil depth. The relationship between vegetation and soil characteristics were evaluated in the Cholistan desert by Arshad et al., (2008). Noureen et al., (2008) used the Calligonum polygonoides species to improve the fertility of soil from the Cholistan desert, Pakistan. They observed that association between these two parameters are common and also concluded that ecological characteristics are the responsible for the distribution of plants seem to be salinity, organic matter and ionic concentration. Wazir et al., (2008) studied the multivariate analysis of vegetation of Chapursan valley an alpine Meadow of Pakistan .They performed classification and ordination of the species on the basis of cluster analysis and recognized 5 vegetation types and were discussed in terms of topography and edaphic heterogeneity. Perveen et al., (2008) described plant biodiversity and phytosociological attributes of Dureji (Khirthar range) .An Inventory of plant species of Dureji game reserve was prepared on the basis of filed trips conducted in different parts of the year particularly in winter, they studied phonological status of each species i.e. flowering and fruiting condition, species diversity, phytosociological attributes, and some ecological parameters such as density, relative density, cover, relative cover, frequency and relative frequency were investigated. Wahab et al., (2008) described the Phytosociology and dynamics and age and growth rates relation of some forests from Afghanistan. Abbas et al., (2009) also investigated the vegetation of Gery Goral Range in Azad Kahmir with respect to the phytosociology. Karim et al., (2009) studied the effect of canopy cover on the organic and inorganic content of soil from Cholistan desert Ahmed et al., (2009) described the vegetation structure of Olea ferruginea Royle forests of lower Dir district of Pakistan and recognized ten different communities with similar floristic composition but different quantitative values. Correlation was checked among density/ basal area, elevation/
31
Chapter 2: Review of literature
density, and elevation/ basal area. They concluded that as the forest showed no recruitment since last 10-15 years; therefore no future trends could be predicted for these forests. The evaluation of ecological aspect of roadside vegetation around Havalian city using multivariate techniques (DCA and CCA) was given by Sheikh et al., (2009) . Jabeen and Ahmed (2009) used multivariate techniques to explore the relationship between vegetation and environment from Ayubia National Park Rawalpindi and they observed that Pb, Cr and Cd were not correlated with the vegetation whereas Zn showed a signification relation with the vegetation. They recognized 5 major communities and 63 plant species and also observed that copper, Zinc and Lead concentrations are the most important factors influencing the vegetation of the study area. Ajaib et al., (2009) analyzed the soil and community description from Goharabad valley District Diamer, Gilgit-Baltistan, Pakistan and they demonstrated that the changing of community composition is due to the deforestation, overgrazing, human influence, soil erosion. Ahmed et al., (2010) studied the floristic composition and communities of deodar forest from Himalayan range of Pakistan. Saima et al.,(2010) described the floristic composition during moon soon in Ayubia National Park, Abbotabad. Khan et al., (2010) described the phytosociology, structure and chemical analysis of soils in Quercus baloot Griff, forest from district Chitral and concluded that Monotheca buxifolia and Quercus baloot showed good regeneration potential but the associated broad leaved species were at the risk of elimination. Ahmed et al., (2010) described the status of vegetation analysis in Pakistan. They divided into five categories i.e observational, quantitative phytosociology, multivariate analysis of ordination, population structure and advanced multivariate techniques. They also gave a description about the trend of this analysis. Siddiqui et al., (2010a, b) described the vegetation of moist temperate forests of Northern areas of Pakistan using Ward’s cluster analysis, TWINSPAN and ordination (DCA and PCA).Ahmed et al., (2011) described the multivariate approach to evaluate the structure and dynamics of Cedrus deodara in Hindu Kush and Himalayan range of Pakistan and they used 13 tree species in the cluster analysis and found six groups while 46 understorey species were reported. They also studied the environmental characteristics using DCA ordination such
32
Chapter 2: Review of literature
as elevation, slope, aspect and canopy and other edaphic factors such as soil compaction, MWHC, salinity, pH and conductivity. Mashwani et al., (2011) classified the vegetation of Saif-ul-Mulook lake Western Himalaya and also described the floristic composition of vegetation using multivariate techniques of TWINSPAN and DCA. Ahmed and Ann (2011) explored the vegetation dynamics and community description of Aubia National Park, Rawalpindi using CCA ordination. They found a significant relation among species abundance, soil moisture and pH. Shaheen et al.,(2011) described in detail, structural diversity, vegetation dynamics and anthropogenic impacts of district Bagh ,Kashmir forests. They observed 0.58 to 1.96 Simpson’s diversity, 1.49 to 1.37 Menhinick‘s diversity while evenness was recorded 0.23 to 0.66. They also stated that this range of diversity was similar with the forests of Himalayan Range of Pakistan but the density, basal area and seedling count were very low. Nimatullah et al., (2011) studied the edaphic factors of soil from 87 different sites of D.I.Khan Division. Khan (2011) gave a detail note about vegetation relationship with soil from forest of Chitral using multivariate analysis and also studied the community structure. On the other hand Wahab (2011) studied the forest of Dir District regarding the status of soil nutrients, physical and chemical properties of forest soil, forests description. Siddiqui (2011) investigated the moist temperate forest of the Pakistan and provide a database for the better management of forests. Ahmed et al., (2011) investigated the soil characters and soil nutrients availability of an open scrub rangeland in the sub-mountainous Himalayan tract of Pakistan. Siddiqui et al., (2011) described communities of moist temperate areas of Pakistan. Khan et al., (2011) described the regeneration potential of Monotheca buxifolia and also investigated the structure and diversity of lower Dir districts forests using multivariate analysis. Akbar et al., (2010, 2011) studied the phytosociology and structure of Skardu district. Hussain et al., (2010, 2011) described the phytosociology and structure of a few sites from Central Karakoram National Park. Ali et al.,(2012) analyzed the ecological ranking in the districts of Pakistan. Khalid et al., also (2012) conducted a study to evaluate the importance of soil nutrients, physical and chemical properties from the Chakwal district. Khan (2012) analyzed the community of Quercus baloot forest from district dir forests at the mean elevation of 1524m to 1753m.He also analyzed the relationship
33
Chapter 2: Review of literature
among edaphic factors, topographic factors and soil nutrients including calcium, magnesium, sodium, potassium and nitrogen. He observed a significant relation among vegetation, potassium, conductivity, TDS, MWHC, salinity and pH. Sikandar and Pandit (2012) phytosociology studied the forest trees and seedling of Kashmir and gave a brief description about the status of forests. Irshad et al., (2012) described the protection impact of the Himalyan and moist temperate forest of Galayat, Pakistan. On the other side Fahad and Bano (2012) quantitatively described the endangered species from the two site of Gilgit and also explore the medicinal important of species. Khan et al., (2012) described the vegetation dynamics, diversity and change of climate from the vegetation of Western Himalaya. Ilyas et al., (2012) described the vegetation composition and pressure to the montane temperate forest from Swat valley. The common growth of plants species recorded 43% perennial, 23% herbs, 16% shrubs and 15% trees. They also concluded that the vegetation of Swat valley is deteriorating with the passage of time due to the deforestation. Khan and Hussain (2013) presented ordination and classification of Karak district Khyber Pakhtunkhwa using multivariate techniques of HCA and DCA. Khan et al., (2013) used the multivariate techniques of cluster analysis and DCA ordination to analyze the vegetation the vegetation and environment relationship from the forests of dir district Hindukush range of Pakistan. They found a significant relation among vegetation, elevation and slope at the probability level P<0.05. They also observed weak correlation among the vegetation, many soil variables and edaphic factors and concluded that the weak correlation due to the anthropogenic disturbance. Siddiqui et al.,(2013) found 8 conifer tree communities with 3 monospecific stands from 41 locations of moist temperate areas of Himalayan and Hindukush range of Pakistan. They also observed the relationship between vegetation and environment and gave a detail description. Central Karakorum National Park is an important national park due to its unique flora and fauna. Initially some researchers and NGO conducted observational studies but with the passage of time quantitative parameters were introduce. In 2009 WWF studied over all vegetation of CKNP, using remote sensing and satellite images techniques; however no quantitative phytosociological work, advance multivariate analysis and population dynamics of vegetation from this park. Therefore present study carried out the
34
Chapter 2: Review of literature
vegetation analysis of study area using quantitative techniques. Also make an attempt to evaluate the population structure of the vegetation. However no quantitative phytosociological work, advance multivariate analysis and dynamics of gymnosperm trees were presented, so far from CKN Park.
35
Chapter 3: Structure and future trend of the vegetation of CKNP
CHAPTER-3 STRUCTURE AND FUTURE TREND OF THE VEGETATION OF CKNP
3.1- Introduction This chapter reveals the structure, present status and future trend of the vegetation of Central Karakoram National Park. The distributions of vegetation, density, basal area, canopy, cover are often used in bio-mass distribution in forest communities (Goff and Zedler, 1968). Spatial forest structure is an important parameter to determine the habitat and diversity of the species (Pommerning, 2002). Structure, composition and the function are the important factors of the forest (Gairola et al., 2008). Different researchers stated that using of diameter distribution transversely a range of diameter class is an aspect of stand structure (Koop et al., 1994). The structure and future trend of forests would be better by selected cutting. (Uuttera et al., 2000). According to FAO (2009) the forest cover of Pakistan is 2 % while 4 % mentioned by forest service which include all trees planted in gardens ,cities, along rivers, canals and agricultural. Ahmed (2008) reported that most of the forest cover found in Northern Areas of Pakistan. These forests are deteriorating with the passage of time due to poor management, less research and anthropogenic disturbance. Forest management can be characterized by silviculture practices which maintain the age class of forests (Matthews, 1999; Schütz, 2001). The cutting of these fortes should be selective with extensive practice of silviculture which could be changed in different times, depending on the need of forest owners and market situation. Additionally apply the system of proportion among small, medium and large dbh size tress. Secondly the possibility to cut trees from different species, aspect of stands (Brang, 2001). Single tree selection silviculture maintains the stand growth and maintains the different size classes and a reverse J-shaped diameter distribution (Matthews 1991). The structure of the forest has been studied at different levels through different communities (Ogden and Ahmed, 1987; Kimmins, 1987). The structure of vegetation determined by the presence of species, quantitative
36
Chapter 3: Structure and future trend of the vegetation of CKNP
relation between the species, species distribution and interaction between them (Boncina, 2000), or even succession (Begon et al., 1990; Cook, 1996).Ahmed (1986) used this method for the vegetation of some foothills of Himalayan range of Pakistan. Ahmed et al. (1990 a, b) also described the status and population structure of Juniperus excelsa in Baluchistan. Ahmed et al., (1991) worked on the vegetation structure and dynamics of Pinus gerardiana forest of Baluchistan. Malik (2005) conducted a comparative study with special reference to range conditions on the vegetation of Ganga Chotti and Bedori Hills District Bagh of Azad Jammu Kashmir. Wahab (2008) et al., described the phytosociology and dynamics of some forests of Afghanistan. Ahmed et al., (2009) described the vegetation structure of Olea ferruginea Royle forests of lower Dir district of Pakistan. Khan et al., (2010) described the phytosociology, structure and chemical analysis of soils in Quercus baloot Griff, forest from district Chitral .Akbar et al., (2010, 2011) also studied the phytosociology and structure of Skardu district while Hussain et al., (2010,2011) studied the quantitative study from the valleys of Central Karakoram National Park, Gilgit. Other researchers Siddique (2011, Khan (2011) and Wahab (2011) studied the structure, present status and future trend of different species from regions of the country. Beside the above studies no such type of study has been conducted to explore the present status, future trend and diameter distribution in terms of probability distribution of vegetation from the Central Karakoram National Park .This study is anticipated to play an important role in modeling diameter distribution and elsewhere in Pakistan. This type of modeling is widely used in forestry. Diameter distribution models are used to obtain an estimate of tree size distribution. This predicted distribution is needed for the further computation of forest volume characteristics and effective management of forests and also expected to be helpful in preserving and conserving the flora of the Park. Objectives of the study ¾ To determined the present status and future trend of the vegetation of Central Karakoram National Park. ¾ To explore the diameter distribution of the forests from CKNP
37
Chapter 3: Structure and future trend of the vegetation of CKNP
3.2-Materials and Methods 3.2.1-Size class structure For size class structure, interval of classes in forested areas was 10cm while in non-forested areas, class interval was 50cm. Ten classes of dbh in forested while 13 classes in non-forested vegetation were exhibited. Diameter at breast height of trees and cover of the shrubs were arranged and counted in each class and covert into density ha-1. Class density of each stand was plotted through MS Excel bar chart. Dbh size classes showed in x- axis while density ha -1 showed in vertical bars. Structure of forested vegetation is divided into three classes i.e. young classes (10-30cm), middle classes (40-, 70cm) large classes (80-100cm) following Ahmed, 1984; Khan, 2011; Wahab, 2011 and Siddique, 2011. Overall diameter of size class is also taken following by Ahmed, 1984; Khan, 2011; Wahab, 2011 and Siddiqui, 2011. 3.2.2-Weibull distribution Many techniques used to explore the modeling diameter distribution. Many distribution functions, such as normal, gamma, Johnson’s SB, Gram-Charlier, beta and Weibull, have been used in describing diameter distributions for forest stands (Von Gadow, 1984; Borders et al. , 1987 ) .Several researchers have reported Weibull function as the most suitable one to portray the diameter distributions (Bailey and Dell 1973, Von Gadow, 1984; Borders et al. , 1987 ). Weibull function is introduced by Bailey and Dell (1973) to describe the distribution of forest. The reputation of the Weibull function is based on its relative simplicity and flexibility Bailey and Dell (1973). Kinerson et al., (1974) fitted a nonlinear least squares model to cumulative crown class frequency data. The model gave a satisfactory fit to their data. The diameter distribution model can be used to obtain the distribution of trees into diameter classes. ( Hyink and Moser, 1983 ). Diameter distribution is also used in the possible outcome of disturbances in the forests. (Hett and Louks, 1976: Denslow, 1995; Baker et al., 2005; Cooms and Allen, 2007). This function is also helpful to explore the structure development of forest (Goff and West, 1975; Poorter et al., 1996 and Zenner,2005).
38
Chapter 3: Structure and future trend of the vegetation of CKNP
For diameter distribution and histograms CumFreq program was used with option of Weibull function (Bailey and Dell 1973).
Three parameter Weibull probability distribution function is given by