Nghiên cứu thực trạng gây trồng loài cây bần chua (sonneratia caseolaris (l ) engler) và cây trang (kandelia obovata sheue, liu yong) ở vùng ven biển bắc bộ

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC LÂM NGHIỆP HÀ ĐÌNH LONG NGHIÊN CỨU THỰC TRẠNG GÂY TRỒNG LOÀI CÂY BẦN CHUA (Sonneratia caseolaris (L.) Engler) VÀ CÂY TRANG (Kandelia obovata Sheue, Liu & Yong) Ở VÙNG VEN BIỂN BẮC BỘ CHUYÊN NGÀNH: LÂM HỌC MÃ SỐ: 8620201 LUẬN VĂN THẠC SĨ LÂM HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS TS ĐỖ ANH TUÂN TS LÊ VĂN THÀNH Hà Nội, 2022 i CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc LỜI CAM ĐOAN Tôi xin cam đoan, cơng trình nghiên cứu tơi, thực sở kết nghiên cứu đề tài cấp Quốc gia: “Nghiên cứu tuyển chọn số lồi ngập mặn thích hợp phục vụ trồng rừng phịng hộ ứng phó với biến đổi khí hậu vùng Bắc Bộ Tây Nam Bộ” TS Lê Văn Thành làm chủ nhiệm đề tài, làm cộng tác viên Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình nghiên cứu khác Nếu nội dung nghiên cứu trùng lặp với cơng trình nghiên cứu cơng bố, tơi xin hồn tồn chịu trách nhiệm tuân thủ kết luận đánh giá luận văn Hội đồng khoa học Hà Nội, ngày 19 tháng năm 2022 NGƯỜI CAM ĐOAN Hà Đình Long ii LỜI CẢM ƠN Để hồn thành Chương trình đào tạo sau Đại học Trường Đại học Lâm nghiệp, thực đề tài “Nghiên cứu thực trạng gây trồng loài Bần chua (Sonneratia caseolaris (L.) Engler) Trang (Kandelia obovata Sheue, Liu & Yong) vùng ven biển Bắc bộ” Trước tiên, xin gửi lời cảm ơn sâu sắc đến TS Lê Văn Thành PGS TS Đỗ Anh Tuân, người thầy tận tình bảo, hướng dẫn truyền đạt kiến thức quý báu, giúp tơi hồn thiện luận văn Tơi xin chân thành cảm ơn thầy, cô Ban giám hiệu Trường Đại học Lâm nghiệp; phòng đào tạo sau đại học; Các thầy, cô khoa Lâm Học dạy cho kiến thức, kinh nghiệm quý báu Tôi xin trân trọng cảm ơn ban lãnh đạo Viện Nghiên cứu Sinh Thái Môi Trường rừng toàn thể cán đồng nghiệp quan bảo, giúp đỡ tơi q trình hồn thành luận văn Cuối cùng, muốn gửi lời cảm ơn đến Bố Mẹ Vợ tôi, người thân gia đình, bạn bè tơi – Những người ủng hộ, tin tưởng động viên tơi suốt q trình học tập vừa qua Trong q trình hồn thành luận văn, thân có nhiều cố gắng kinh nghiệm nghiên cứu chưa nhiều, hạn chế mặt thời gian nên luận văn khơng thể tránh khỏi có khiếm khuyết định Tôi mong nhận bảo, đóng góp thầy giáo, chuyên gia đồng nghiệp để luận văn hồn thiện Tơi xin chân thành cảm ơn! Hà Nội,ngày 19 tháng năm 2022 TÁC GIẢ Hà Đình Long iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC TỪ VIẾT TẮT v DANH MỤC CÁC BẢNG, BIỂU vi DANH MỤC CÁC HÌNH ẢNH viii ĐẶT VẤN ĐỀ Chương TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Trên giới 1.1.1 Vai trò RNM phòng hộ ven biển ứng phó với BĐKH 1.1.2 Nghiên cứu loài Bần chua 1.1.3 Nghiên cứu loài Trang 1.2 Trong nước 1.2.1 Vai trò RNM phịng hộ ven biển ứng phó với BĐKH 1.2.2 Nghiên cứu loài Bần chua 12 1.2.3 Nghiên cứu loài Trang 17 Chương ĐỐI TƯỢNG, PHẠM VI, MỤC TIÊU, NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 20 2.1 Đối tượng, phạm vi nghiên cứu 20 2.1.1 Đối tượng nghiên cứu 20 2.1.2 Phạm vi nghiên cứu 20 2.2 Mục tiêu nghiên cứu 20 2.2.1 Mục tiêu tổng quát 20 2.2.2 Mục tiêu cụ thể 20 2.3 Nội dung phương pháp nghiên cứu 20 2.3.1 Nội dung nghiên cứu 20 2.3.2 Phương pháp nghiên cứu 21 iv Chương KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 31 3.1 Điều tra đánh giá trạng RNM phòng hộ Bần chua Trang vùng ven biển Bắc Bộ, lựa chọn lâm phần tuyển chọn trội/mẹ 31 3.1.1 Diện tích Phân bố rừng trồng lồi Bần chua Trang khu vực nghiên cứu 31 3.1.2 Hiện trạng sinh trưởng phát triển RNM Bần chua Trang khu vực nghiên cứu 33 3.1.3 Chọn lâm phần trội/mẹ 49 3.2 Điều kiện lập địa nơi trồng/phân bố Bần chua Trang 70 3.2.1 Điều kiện lập địa nơi trồng Bần chua 70 3.2.2 Điều kiện lập địa nơi trồng Trang 80 3.3 Bổ sung kỹ thuật gây trồng Bần chua Trang 90 3.3.1 Kỹ thuật nhân giống, trồng Bần chua 90 3.3.2 Kỹ thuật nhân giống, trồng Trang 91 KẾT LUẬN, TỒN TẠI VÀ KIẾN NGHỊ 94 TÀI LIỆU THAM KHẢO 97 PHỤ LỤC v DANH MỤC TỪ VIẾT TẮT TT Ký hiệu Giải thích BĐKH Biến đổi khí hậu KVNC Khu vực nghiên cứu RNM Rừng ngập mặn OTC Ô tiêu chuẩn CV Hệ số biến thiên Hvn Chiều cao vút D00 Đường kính gốc Hdc Chiều cao cành Dt Đường kính tán 10 UBND Ủy ban nhân dân 11 TCVN Tiêu chuẩn Việt Nam vi DANH MỤC CÁC BẢNG, BIỂU Bảng 2.1 Phân cấp mức độ sâu bệnh hại lâm phần 24 Bảng 3.1 Diện tích nơi trồng Bần chua Trang 31 Bảng 3.2 Sinh trưởng Bần chua phường Yên Giang, thị xã Quảng Yên, tỉnh Quảng Ninh 34 Bảng 3.3 Sinh trưởng Bần chua xã Đại Hợp, huyện Kiến Thuỵ, TP Hải Phòng 35 Bảng 3.4 Sinh trưởng Bần chua xã Đơng Hồng, huyện Tiền Hải, Thái Bình 37 Bảng 3.5 Sinh trưởng Bần chua xã Nam Điền, xã Nghĩa Hải, huyện Nghĩa Hưng, tỉnh Nam Định 38 Bảng 3.6 Sinh trưởng Bần chua xã Kim Đông, huyện Kim Sơn, tỉnh Ninh Bình 40 Bảng 3.7 Sinh trưởng Trang xã Quảng Phong, huyện Hải Hà, tỉnh Quảng Ninh 42 Bảng 3.8 Sinh trưởng Trang xã Đại Hợp, huyện Kiến Thụy, TP Hải Phòng 43 Bảng 3.9 Sinh trưởng Trang xã Thụy Xuân, huyện Thái Thụy, tỉnh Thái Bình 45 Bảng 3.10 Sinh trưởng Trang xã Nghĩa Hùng, huyện Nghĩa Hưng, tỉnh Nam Định 46 Bảng 3.11 Sinh trưởng Trang xã Kim Chung, Kim Sơn, Ninh Bình 48 Bảng 3.12 Sinh trưởng Bần chua địa điểm điều tra tỉnh 49 Bảng 3.13 Sinh trưởng Trang địa điểm điều tra tỉnh 51 Bảng 3.14 Các tiêu sinh trưởng trội Bần chua trồng năm 2008 xã Đơng Hồng, huyện Tiền Hải, tỉnh Thái Bình 53 Bảng 3.15 Các tiêu sinh trưởng trội Bần chua trồng năm 2008 xã Kim Đơng, huyện Kim Sơn, tỉnh Ninh Bình 58 vii Bảng 3.16 Các tiêu sinh trưởng trội Trang trồng năm 2007 xã Thuỵ Xuân, huyện Thái Thuỵ, tỉnh Thái Bình 62 Bảng 3.17 Các tiêu sinh trưởng trội Trang trồng năm 2007 xã Nghĩa Hùng, huyện Nghĩa Hưng, tỉnh Nam Định 66 Bảng 3.18 Kết điều tra điều kiện gây trồng Bần chua 70 Bảng 3.19 Kết phân tích tiêu lý hóa tính 20 mẫu đất nơi trồng Bần chua 77 Bảng 3.20 Kết điều tra điều kiện gây trồng Trang 80 Bảng 3.21 Kết phân tích tiêu lý hóa tính 20 mẫu đất nơi trồng Trang 87 viii DANH MỤC CÁC HÌNH ẢNH Hình 3.1 Lập OTC đo đếm Bần chua phường Yên Giang – TX Quảng Yên – Tỉnh Quảng Ninh 33 Hình 3.2 Đo đếm tái sinh Ơ dạng phường Yên Giang – TX Quảng Yên – Tỉnh Quảng Ninh 33 Hình 3.3 Rừng bần chua Quảng Yên 35 Hình 3.4 Thảm thực vật rừng Bần chua Quảng Yên 35 Hình 3.5 Rừng Bần chua xã Đại Hợp, huyện Kiến Thuỵ, TP Hải Phịng 36 Hình 3.6 Rừng bần chua xã Đơng Hồng, huyện Tiền Hải, Thái Bình 38 Hình 3.7 Rừng trồng bần chua xã Nghĩa Hải, Nghĩa Hưng, Nam Định 39 Hình 3.8 Rừng trồng Bần chua xã Kim Chung, Kim Sơn, Ninh Bình 41 Hình 3.9 Lập OTC, Ơ dạng đo đếm tái sinh rừng Trang 41 Hình 3.10 Rừng trồng Trang xã Quảng Phong, huyện Hải Hà, Quảng Ninh 43 Hình 3.11 Rừng trồng Trang Kiến Thuỵ - Hải Phịng 44 Hình 3.12 Rừng trồng 2009 Thuỵ Xuân – Thái Thuỵ 46 Hình 3.13 Rừng Trang trồng năm 2004 xã Nghĩa Hùng, huyện Nghĩa Hưng, tỉnh Nam Định 47 Hình 14 Rừng trồng Trang Kim Sơn – Ninh Bình 49 Hình 3.15 Hình ảnh chọn trội Bần chua xã Đơng Hồng, huyện Tiền Hải, tỉnh Thái Bình 57 Hình 3.16 Cây trội Bần chua xã Kim Đơng, huyện Kim Sơn, 61 Hình 3.17 Cây trội Trang xã Thuỵ Xuân, huyện Thái Thuỵ, tỉnh Thái Bình 65 Hình 3.18 Hình ảnh trội Trang xã Nghĩa Hùng, huyện Nghĩa Hưng, tỉnh Nam Định 69 Hình 3.19 Độ lún bàn chân khu vực điều tra 76 Hình 3.20 Khoan mẫu lấy mẫu đất OTC 89 Hình 3.21 Vườn ươm Bần chua xã Đơng Hồng, huyện Tiền Hải, tỉnh Thái Bình 90 ĐẶT VẤN ĐỀ Việt Nam với bờ biển dài 3260 km, chạy dọc từ Móng Cái đến Hà Tiên, hàng năm chịu nhiều thiên tai đặc biệt bão lũ nước đứng đầu giới dễ bị tổn thương BĐKH gây (Chương trình Phát triển Liên Hợp Quốc, 2007) BĐKH thực làm cho thiên tai, đặc biệt bão, lũ hạn hán ngày khốc liệt; riêng giai đoạn 2000 - 2010, số người thiệt mạng thiên tai 5.045 người, số người tích 617 người, tổng thiệt hại khoảng 91.000 tỷ đồng (INDC Việt Nam, 2015) Một giải pháp phòng chống thiên tai, ứng phó với BĐKH gây cho vùng ven biển phát triển bền vững hệ sinh thái RNM Ngoài tác dụng cố định phù sa, lấn biển, chống xói lở bờ biển, RNM cịn làm giảm tốc độ gió, giảm lượng sóng, hạn chế tác hại sóng biển vùng bị bão lốc xốy, giảm 50% tổn hại thiên tai gây so với vùng khơng có RNM bảo vệ (Kathiresan, 2012) Vai trò RNM chống lại mối nguy hiểm ven biển nhận ý từ sau hậu đợt sóng thần năm 2004 Ấn Độ Dương, nhiều nghiên cứu tốn học RNM có khả làm giảm cường độ sóng; lượng sóng hấp thụ phụ thuộc phần lớn vào mật độ rừng, đường kính thân cây, rễ cây, đường kính tán tầng (World Bank, 2016) RNM với lớp tán dày với thân, cành tạo thành lớp rào vật liệu mềm giảm sức công phá sóng biển (Phan Nguyên Hồng, 2007) Luật Lâm nghiệp năm 2017 khẳng định, RNM tác nhân quan trọng hạn chế tác hại BĐKH, giảm nhẹ thiệt hại nước biển dâng, phòng hộ bảo vệ bờ biển khỏi sóng biển xâm lấn… Đánh giá vai trò RNM, năm gần đây, Chính phủ Việt Nam đưa sách để phục hồi phát triển hệ sinh thái RNM ven biển nhằm phịng chống thiên tai, ứng phó với BĐKH Để thúc đẩy nhanh việc phát triển rừng ven biển có RNM, ngày 23/8/2016 Chính phủ ban hành Nghị định 119/2016/NĐ-CP, Nghị định đưa số sách đặc biệt cho phát triển bền vững hệ sinh thái 104 58 Chen, Y ; Liao, B ; Huang, Chen, Y ; Lin, W and Li, M , 2011 Studies on the quantitative attenuatation effect on artificial mangrove Kanelia obovata and Sonneratia apelatat on waves Oceanologia et Limnologia Sinica, Vol 42, No 6, pp 764-770 59 Chen, Y ; Liao, B , Zheng, S ; Li, M and Song, X , 2004 Dynamics and species-diversities of artificial Sonneratia apetala, Sonneratia caseoloaris and Kandelia candel communities The Journal of Applied Ecology, Vol 15, Issue 6, pp 924 – 928 60 Cheng, S ; Tam, N F Y; Li, R ; Shen, X ; Niu, Z ; Chai, M and Qiu, G Y , 2017 Temporal variations in physiological responses of Kandelia obovata seedlings exposed to multiple heavy metals Marine Pollution Bulletine, Vol 124, Issue 2, pp 1089-1095 61 Christensen, B and Wium-Andersen, S , 1977 Seasonal growth of mangrove trees in Southern Thailand I The phenology of Rhizophora apiculata Bl Aquatic Botany, Vol 3, pp 281- 286 62 Clarke, A and Johns, L , 2002 Fisheries Guidelines for Mangrove nurseries: Construction, propagation and planting Fish Habitat Guideline FHG 004 Department of Primary Industries Queensland Fishereis Service, March 2002 63 Clough, B ; Tan, D T ; Phuong, D X and Buu, D C , 2000 Canopy leaf area index and litter fall in stands of the mangrove Rhizophora apiculata of different age in the Mekong Delta, Vietnam Aquatic Botany, Vol 66, Issue 4, pp 311 – 320 64 Cousins, J and Saenger, P 2002 Developing a protocol for in vitro propagation of the grey mangrove – Avicennia marina Refereed proceedings of the 7th meeting, Armidale, NSW, 20 – 23 January 2002, University of New England, Armidale, NSW 105 65 Duke, N C , 1988 Phenologies and litter fall of two mangrove trees, Sonneratia alba Sm And S caseolaris (L ) Engl , and their Putative hybrid, S x gulngai N C Duke Australian Journal of Botany, Vol 36, Issue 4, pp 473 – 482 66 Drexler, J Z , 2001 Maximum longevities of Rhizophora apiculata and R mucronata propagules Pacific Science, Vol 55, Issue 1, pp 17 – 22 67 Guan, G ; Wang, Y ; Cheng, H ; Jiang, Z and Fei, J , 2015 Physiological and biochemical response to drought stress in the leaves of Aegicerascorniculatum and Kandelia obovata Ecotoxicology, Vol 24, Issue – 8, pp 1688-1676 68 Hawlader, N I , 1999 Forest Resources Management Project Mid-term review, Ministry of Environment and Forest, Dhaka, p 103 69 Hideki, H , Oliva, S and Atsuo I, 1998 Nursery manual for mangrove species The development of sustainable mangrove management project Ministry of Forest and Estate crops, Indonesia and Japan International Cooperation Agency, 50 pages 70 IUCN, 2010 The IUCN red list of threatened species 2010 71 Judy, Z D , 2001 Maximum longevities of Rhizophora apiculata and R mucronata propagules Pacific Science, Vol 55, Issue 1, pp 17 – 22 72 Kathiresan, K , 2012 Importance Ecosystems of Asia: Status, Challenges and Management Strategies Springer New York, New York, pp 233 – 256 73 Kathiresan, K and Thangam, T S 1990 A note on the effects of salinity and pH on growth of Rhizophora seedlings Indian Forester, Vol 116, Issue 3, pp, 243 – 244 74 Kathiresan, K and Rajendran, N , 2002 Growth of mangrove (Rhizophora apiculata) seedlings as influenced by GA3, linght and salinity International Journal of Tropical biolohy and Biology and conservation, Vol 2, pp 525 – 536 106 75 Kathiresan, K and Ravikumar, S , 1995 Vegetative propagation through air- layering in two species of mangroves Aquatic Botany, Vol 50, pp 107 – 110 76 Kathiresan, K ; Salmo, III S G ; Fernando, E S ; Peras, J R ; Sukardjo, S ; Miyagi, T ; Ellison, J ; Koedam, N E ; Wang, Y ; Primavera, J ; Jin, E O ; Wan-Hong, Y J and Ngoc Nam, V , 2010 Sonneratia caseolaris 77 Kesari, V ; Krishnamachari, A and Rangan, L , 2008 Effect of auxins on adventitious rooting from stem cuttings of candidate plus tree Pongamia pinnata (L ), a potential biodiesel plant Trees, Vol 23, Issue 3, pp 597 – 604 78 Kesava, K R , Venkaiah, K and Bramramba, B , 1994 Rooting of stem cuttings of Avicennia officinalis Linn and Avicennia alba BL : A tool for afforestation of blanks in mangrove forests The Indian Forester A Pioneer monthly journal in forestry research and education, Vol 120, Issue 79 Khan, N I and Kabir, E , 2017 Ecology of Kandelia obovata (S , L ) Yong: A Fast-growting mangrove in Okinawa, Japan Chapter 18 in book: Participatory mangrove management in a changing climate Doi: 10 1007/978-4-431-56481-2_18 80 Kumar, M A , 2010 Conservation of Aegiceras corniculatum (L ) Blanco (River mangrove, Khalsi): A new approach of vegetative propagation through hypocotylar juvenile stem cuttings International Journal of Biodiversity and Conservation, Vol 2, pp 162 – 172 81 Lewis, R R , 2005 Ecological engineering for successful management and restoration of mangrove forests Ecol Eng 24, 403–418 82 Li, Y ; Zheng, D Z and Liao, B W cộng sự, 1995 Thử nghiệm dẫn giống gieo ươm Bần chua Thông luận Khoa học Kỹ thuật Lâm nghiệp, Trang 21 - 22 (Dịch từ tài liệu tiếng Trung) 107 83 Li, Y ; Zheng, D Z and Chen, H X cộng sự, 1998 Bước đầu nghiên cứu dẫn giống thực vật Ngập mặn Bần chua Tạp chí Khoa học Lâm nghiệp, Số 11 (1), Trang 39 - 44 (Dịch từ tài liệu tiếng Trung) 84 Liang, S ; Liang, M ; Yuan-ling, W U ; Qiang, XIE ; Qi-jie, ZAN and Wang, Y , 2005 Study on the population structure of Sonneratia caseolaris at Futian Mangrove Reserve of Shenzhen Journal of Guangxi Normal University(Natural Science), 2005-02 85 Liao, B , Zheng, D and Li, Y , 1997 Seed germination conditions of Sonneratia carseolaris of mangrove Journal of Central South Forestry University, Vol 86 Lim, Z K , Ngoh, G H P , Goh, M M and Loh, T Y K , 2012 Investigating the effect of soil pH on the germination of Avicennia alba seedlings Project report Little Green Dot Student Research Grant Nature Society (Singapore) 87 Mark, D S ; Susan, R ; Carmen, L ; Imèn, M ; Lynne, Z H ; Christine, C S and Michael, W B , 2014 The role of ecosystems in coastal protection: Adapting to climate change and coastal hazards Ocean & Coastal Management, Vol 90, March 2014, pp 50 – 57 88 Mazda; Wolanski; King; Save, A and Magi, 1997 Drag fore due to vegetation in mangrove swamps Mangroves and salt marshes, Vol 1, pp 237 – 135 89 Nandy, P ; Alam, M J and Haider, M R , 2004 Establishment of mangrove seed production area for Sonneratia apetala Journal of Tropical Forest Science, Vol 16, Issue 3, pp 363-368 90 Primvera, J H , 2004 Mangrove greenbelts and utilization of mangrove ecosystem TERRAPUB, pp 204 – 220 91 Ravishankar, T and Ramasubramanian, R , 2004 Manual on Mangrove nursery raising techniques M S Swaminathan Research Foundation Chennai, India 108 92 Samuel, C S and Jane, G S , 1984 The Mangrove Ecosystem: Research Methods Paris: Unesco, 1984: 1-251 93 Sheue, C ; Liu, H and Yong, J W H , 2003 Kandelia obovata (Rhizophoraceae), a new mangrove species from Eastern Asia Taxon, Vol 52, Issue 2, pp 287 – 294 94 Siddiqi, N A et al,1993 Mangrove nurseries in Bangladesh International Society for Mangrove ecosystem, Okinawa, Japan, occasional papers, Vol 1, pp 97 – 110 95 Silva, K H W L and Amarashighe, M D , 2010 Vegetative propagation of some selected mangrove species from Negombo estuary, Sri Lanka Sri Lanka Journal of Aquatic Sciences, Vol 15, pp 25 – 38 96 Silva, K H W L and Amarasinghe, M D , 2013 Vegetative propagation of some selected mangrove species from Negombo estuary, Sri Lanka Sri Lanka Journal of Aquatic Sciences, Vol 15 97 Somboon, J R P , 1990 Palynological study of mangrove and marine sediments of the Gulf of Thailand Journal of Southeast Asian Earth Science, Vol 4, Issue 2, pp 85-97 98 Sufang, C ; Renchao, Z ; Yelin, H ; Meng, Z ; Guili, Y ; Cairong, Z and Suhua, S , 2011 Transcriptome sequencing of a highly salt tolerant mangrove species Sonneratia alba using Illumina platform Marine Genomics, Vol 4, Issue 2, pp 129 – 136 99 Thatoi, H ; Mishra, P K ; Mohanty, J R and Acharjyo, L N , 2001 Vegetative propagation by stem cutting with auxins of four mangrove (and associate) species of Bhitarkanika, India Journal of Tropical Forest Science, Vol 13, No 1, pp 223-227 100 Trench, C A and Webber, M K , 2011 Nursery propagation of Jamaican coastal forest species International Symposium on Tropical Horiculture 101 UNEP, 2005 After the tsunami Rapid environment assessment UNEP Nairobi, Keynia 109 102 William, N ; Curtis, E W and Serigio, F , 2016 Bottom sediments affect Sonneratia mangrove forests in the prograding Mekong delta, Vietnam Estuarine, Coastal and Shelf Sciences, Vol 177, Vol 5, pp 60 – 70 103 World Bank Group, 2016 Managing coasts with natural solutions: Guidelines for measuring and valuing the coastal protection services of mangroves and coral reefs Wealth accounting and the valuation of ecosystem services partnership (WAVES), World Bank, Washington, DC 104 Youssef, T , 1997 Approaches in mangrove planting: some options for Darwin mangroves International Marine Biological Workshop, pp 321333 105 Zan, Q.J.; Wang, B.S.; Wang, Y.J and Li, 2003 Ecological assessment on the introduced Sonneratia caseolaris and S apetala at the Mangrove Forest of Shenzhen Bay, China Acta Botanica Sinica, Vol 45, Issue 5, pp 544 – 551 106 Zheng, J.; Wang, J W.; Chen, Q X cộng 2010 Thử nghiệm di chuyển số loài Ngập mặn lên trồng vùng ven biển phía nam tỉnh Chiết Giang Tạp chí Viện Lâm học Tây Nam (Dịch từ tài liệu tiếng Trung Quốc) 107 Zeng, X.; Luzhen, C.; FungYee, T.N.; Huang, J.; Xu, H and Lin, G., 2008 Seedling emergence and dispersal pattern of the introduced Sonneratia caseolaris in Shenzhen Bay, China Biodiversity Science, Vol.16, Issue 3, pp 236 – 244 Nguồn Website 108 Nparks Flora & Fauna Web: https://florafaunaweb.nparks.gov.sg/ 109 Plant Resources of South-East Asia: https://uses.plantnetproject.org/en/ 110 Useful Tropical Plants Website: http://tropical.theferns.info/ 111 WWF Hong Kong: http://assets.wwf.org.hk/ PHỤ LỤC Phụ lục 01: Kết so sánh lâm phần trồng Bần chua tỉnh NPar Tests Descriptive Statistics N Mean Std Deviation Minimum Maximum Duong kinh goc (cm) 2865 7.3714 2.84844 60 26.10 Chieu cao vut ngon (m) 2865 3.2860 95727 25 5.60 Chieu cao duoi canh (m) 2865 38.7054 60.97953 15 340.00 Duong kinh tan (m) 2865 7760 39290 00 3.00 Goc phan canh 2865 31.8073 9.02324 3.00 75.00 Lam phan Trang 2865 3.31 1.402 Kruskal-Wallis Test Ranks Lam phan Trang Duong kinh goc (cm) 239 576.77 LP Ban chua Hai Phong 927 854.67 LP Ban chua Thai Binh 277 2349.36 LP Ban chua Nam Dinh 555 1800.49 LP Ban chua Ninh Binh 867 1759.37 239 121.35 LP Ban chua Hai Phong 927 1734.93 LP Ban chua Thai Binh 277 1354.87 LP Ban chua Nam Dinh 555 1530.00 LP Ban chua Ninh Binh 867 1434.62 239 133.04 LP Ban chua Hai Phong 927 893.00 LP Ban chua Thai Binh 277 1429.72 LP Ban chua Nam Dinh 555 1740.00 LP Ban chua Ninh Binh 867 2173.26 2865 LP Ban chua Quang Ninh 239 1116.50 LP Ban chua Hai Phong 927 1395.82 LP Ban chua Thai Binh 277 2253.75 LP Ban chua Nam Dinh 555 1760.32 LP Ban chua Ninh Binh 867 1088.24 Total Goc phan canh 2865 LP Ban chua Quang Ninh Total Duong kinh tan (m) 2865 LP Ban chua Quang Ninh Total Chieu cao duoi canh (m) Mean Rank LP Ban chua Quang Ninh Total Chieu cao vut ngon (m) N LP Ban chua Quang Ninh 2865 239 2484.48 LP Ban chua Hai Phong 927 1163.18 LP Ban chua Thai Binh 277 1757.88 LP Ban chua Nam Dinh 555 1479.61 LP Ban chua Ninh Binh 867 1298.01 Total 2865 Test Statisticsa,b Duong kinh Chieu cao vut Chieu cao duoi Duong kinh Goc phan goc (cm) ngon (m) canh (m) tan (m) canh Kruskal-Wallis H 1296.110 736.902 1757.466 552.046 627.175 4 4 000 000 000 000 000 Df Asymp Sig a Kruskal Wallis Test b Grouping Variable: Lam phan Trang Homogeneous Subsets Duong kinh goc (cm) Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Ban chua Quang Ninh 239 LP Ban chua Hai Phong 927 LP Ban chua Ninh Binh 867 8.1543 LP Ban chua Nam Dinh 555 8.3641 LP Ban chua Thai Binh 277 Sig 4.7678 5.6339 10.9931 1.000 1.000 170 1.000 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed Chieu cao vut ngon (m) Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Ban chua Quang Ninh 239 LP Ban chua Nam Dinh 555 3.4086 LP Ban chua Ninh Binh 867 3.4699 LP Ban chua Thai Binh 277 LP Ban chua Hai Phong 927 9254 3.4699 3.5224 3.5224 3.5786 Sig 1.000 162 231 199 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed Chieu cao duoi canh (m) Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Ban chua Quang Ninh 239 4203 LP Ban chua Hai Phong 927 1.8778 LP Ban chua Thai Binh 277 LP Ban chua Ninh Binh 867 71.7382 LP Ban chua Nam Dinh 555 74.5912 19.6921 Sig .675 1.000 411 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed Duong kinh tan (m) Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Ban chua Ninh Binh 867 6122 LP Ban chua Quang Ninh 239 6217 LP Ban chua Hai Phong 927 LP Ban chua Nam Dinh 555 LP Ban chua Thai Binh 277 7113 8932 1.4032 Sig .665 1.000 1.000 1.000 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed Goc phan canh Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Ban chua Hai Phong 927 LP Ban chua Ninh Binh 867 LP Ban chua Nam Dinh 555 LP Ban chua Thai Binh 277 LP Ban chua Quang Ninh 239 Sig 28.5739 29.6251 32.4955 36.1011 45.6904 1.000 1.000 1.000 1.000 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed 1.000 Phụ lục: Kết so sánh lâm phần trồng Trang tỉnh NPar Tests Descriptive Statistics N Mean Std Deviation Minimum Maximum Duong kinh goc (cm) 2865 7.3714 2.84844 60 26.10 Chieu cao vut ngon (m) 2865 3.2860 95727 25 5.60 Chieu cao duoi canh (m) 2865 1.2412 81349 10 4.00 Duong kinh tan (m) 2865 7760 39290 00 3.00 Goc phan canh 2865 31.8073 9.02324 3.00 75.00 Lam phan Trang 2865 3.31 1.402 Kruskal-Wallis Test Ranks Lam phan Trang Duong kinh goc (cm) 239 576.77 LP Trang Hai Phong 927 854.67 LP Trang Thai Binh 277 2349.36 LP Trang Nam Dinh 555 1800.49 LP Trang Ninh Binh 867 1759.37 239 121.35 LP Trang Hai Phong 927 1734.93 LP Trang Thai Binh 277 1354.87 LP Trang Nam Dinh 555 1530.00 LP Trang Ninh Binh 867 1434.62 239 570.47 LP Trang Hai Phong 927 2074.67 LP Trang Thai Binh 277 1082.00 LP Trang Nam Dinh 555 1729.76 LP Trang Ninh Binh 867 906.86 2865 LP Trang Quang Ninh 239 1116.50 LP Trang Hai Phong 927 1395.82 LP Trang Thai Binh 277 2253.75 LP Trang Nam Dinh 555 1760.32 LP Trang Ninh Binh 867 1088.24 Total Goc phan canh 2865 LP Trang Quang Ninh Total Duong kinh tan (m) 2865 LP Trang Quang Ninh Total Chieu cao duoi canh (m) Mean Rank LP Trang Quang Ninh Total Chieu cao vut ngon (m) N LP Trang Quang Ninh 2865 239 2484.48 LP Trang Hai Phong 927 1163.18 LP Trang Thai Binh 277 1757.88 LP Trang Nam Dinh 555 1479.61 LP Trang Ninh Binh 867 1298.01 Total 2865 Test Statisticsa,b Kruskal-Wallis H Duong kinh goc Chieu cao vut Chieu cao duoi Duong kinh tan (cm) ngon (m) canh (m) (m) Goc phan canh 1296.110 736.902 1291.797 552.046 627.175 4 4 000 000 000 000 000 df Asymp Sig a Kruskal Wallis Test b Grouping Variable: Lam phan Trang Homogeneous Subsets Duong kinh goc (cm) Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Trang Quang Ninh 239 LP Trang Hai Phong 927 LP Trang Ninh Binh 867 8.1543 LP Trang Nam Dinh 555 8.3641 LP Trang Thai Binh 277 Sig 4.7678 5.6339 10.9931 1.000 1.000 170 1.000 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed Chieu cao vut ngon (m) Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Trang Quang Ninh 239 LP Trang Nam Dinh 555 9254 3.4086 LP Trang Ninh Binh 867 LP Trang Thai Binh 277 LP Trang Hai Phong 927 3.4699 3.4699 3.5224 3.5224 3.5786 Sig 1.000 162 231 199 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed Chieu cao duoi canh (m) Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Trang Quang Ninh 239 LP Trang Ninh Binh 867 LP Trang Thai Binh 277 LP Trang Nam Dinh 555 LP Trang Hai Phong 927 4203 7174 9106 1.5149 1.8778 Sig 1.000 1.000 1.000 1.000 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed Duong kinh tan (m) Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Trang Ninh Binh 867 6122 LP Trang Quang Ninh 239 6217 LP Trang Hai Phong 927 LP Trang Nam Dinh 555 LP Trang Thai Binh 277 7113 8932 1.4032 Sig .665 1.000 1.000 1.000 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed Goc phan canh 1.000 Duncana,b Subset for alpha = 0.05 Lam phan Trang N LP Trang Hai Phong 927 LP Trang Ninh Binh 867 LP Trang Nam Dinh 555 LP Trang Thai Binh 277 LP Trang Quang Ninh 239 28.5739 Sig 29.6251 32.4955 36.1011 45.6904 1.000 1.000 1.000 1.000 Means for groups in homogeneous subsets are displayed a Uses Harmonic Mean Sample Size = 422.720 b The group sizes are unequal The harmonic mean of the group sizes is used Type I error levels are not guaranteed 1.000