Đánh giá giá trị kinh tế và tiềm năng giảm phát thải khí nhà kính CO2 của rừng ngập mặn Rú Chá, tỉnh Thừa Thiên Huế

Trần Thị Tú, Trần Hiếu Quang

Viện Tài nguyên và Môi trường- Đại học Huế

E-mail: tutrancreb@yahoo.com

TÓM TẮT

Kết quả đã xác định được 27 loài TVNM thuộc 26 chi, 22 họ của 2 ngành Dương xỉ (Polypodiophyta) và Ngọc lan (Magnoliophyta). TVNM ở Rú Chá có rất nhiều giá trị sử dụng như cây cho gỗ, củi đốt, làm thuốc, thực phẩm… Rú Chá bị tác động mạnh bởi bốn nguyên nhân: sản xuất nông nghiệp, xây dựng cơ sở hạ tầng, san lấp đất rừng và hoạt động khác. Giá trị kinh tế mà rừng ngập mặn Rú Chá ước tính khoảng 1,24 tỷ đồng/năm, trong đó 86,4% giá trị sử dụng trực tiếp, 8,1% giá trị sử dụng gián tiếp và 5,5% giá trị phi sử dụng. Tổng lượng carbon tích lũy ước tính trung bình đạt 50,2 tấn/ha; tiềm năng giảm phát thải khí nhà kính trung bình đạt 184,0 tấn CO₂/ha, tương ứng với 39,86 triệu đồng/ha/năm. Như vậy, tổng giá trị giảm phát thải khí nhà kính CO₂ của Rú Chá đạt 231,6 triệu đồng/năm.

Từ khóa: Đa dạng loài, giá trị kinh tế, khí nhà kính, rừng ngập mặn Rú Chá.

ABSTRACT

The results have identified 27 species mangrove flora of 26 genera, 22 families, 2 phylums included Polypodiophyta and Magnoliophyta. Mangrove flora at Ru Cha have many valuable uses, such as timber, firewood, medicinal, food, etc. There were four causes to effect strongly Ru Cha mangroves such as agriculture, infrastructure construction, land reclamation and other activities. The total economic value of Ru Cha mangroves was estimated 1.24 billion VNĐ/year, including 86.4% of direct value, 8.1% of indirect value and 5.5% of non-using value. In Ru Cha mangroves, total cumulative carbon reached 50.2 tons/ha; the potential to reduce greenhouse gas emissions was 184.0 tons CO₂/ha which corresponds to 39.86 million VNĐ/ha/year. Therefore, the total value of GHG emission reduction of Ru Cha reached 231.6 million VNĐ/year.

Keywords: Economic value, greenhouse gas (GHG), Ru Cha mangroves, species diversity.

I. ĐẶT VẤN ĐỀ

Rừng ngập mặn Rú Chá thuộc địa phận quản lý hành chính của thôn Thuận Hòa, xã Hương Phong, thị xã Hương Trà, tỉnh Thừa Thiên Huế. Tổng diện tích Rú Chá hiện còn khoảng 5,81 ha. Do diện tích Rú Chá hiện còn ít, mặt khác hoạt động nông nghiệp và nuôi trồng thủy sản là hai hoạt động sản xuất chính ở địa phương nên chính quyền cũng chưa quan tâm nhiều về diện tích đất sản xuất lâm nghiệp, trong đó có diện tích rừng ngập mặn Rú Chá.

Rú Chá là một trong những hệ sinh thái rừng ngập mặn (HST RNM) còn lại ở khu vực đầm phá Tam Giang- Cầu Hai, bên cạnh thảm thực vật ngập mặn ở cửa sông Bù Lu- Cảnh Dương, thôn Tân Mỹ, xã Phú Tân, huyện Phú Vang và ở đầm Lập An, huyện Phú Lộc. RNM Rú Chá có chức năng như một vùng đệm sinh thái giữa đất liền và đầm phá. Ngoài ra, đây còn là bãi giống lý tưởng cho nhiều loài thủy sinh như các loài cá, loài giáp xác… Bên cạnh đó, Rú Chá là nơi phục vụ cho công tác nghiên cứu khoa học, thử nghiệm các giải pháp bảo tồn rừng ngập mặn nhằm chọn lựa giải pháp tối ưu nhất cho các khu vực tương đồng trong khu vực. Đồng thời, Rú Chá nằm ven theo đầm phá Tam Giang- Cầu Hai được nhiều người biết đến như là một điểm du lịch sinh thái hấp dẫn trong tương lai.

Mục tiêu của nghiên cứu này là đánh giá hiện trạng, cập nhật thông tin về sự đa dạng loài thực vật ngập mặn (TVNM), ước tính các giá trị kinh tế và khả năng giảm phát thải khí nhà kính của rừng ngập mặn Rú Chá làm cơ sở cho việc quản lý, bảo tồn và phát triển nguồn tài nguyên RNM nơi đây.

II. ĐỐI TƯỢNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU

2.1. Đối tượng nghiên cứu

Thảm TVNM ở Rú Chá, xã Hương Phong, thị xã Hương Trà, tỉnh Thừa Thiên Huế.

2.2. Phương pháp nghiên cứu

2.2.1. Tổng hợp tài liệu

Tiến hành thu thập các số liệu, thông tin liên quan đến TVNM ở Rú Chá, xã Hương Phong, thị xã Hương Trà, tỉnh Thừa Thiên Huế.

2.2.2. Phương pháp nghiên cứu có sự tham gia của người dân địa phương (PRA)

Phỏng vấn trực tiếp, thảo luận nhóm, gửi phiếu điều tra thông tin, điều tra theo tuyến để thu thập các thông tin về hiện trạng TVNM ở Rú Chá, về các hoạt động kinh tế- xã hội và tác động của nó đến TVNM. Hoạt động điều tra tiến hành phỏng vấn người dân theo tiêu chí là người lớn tuổi, đã sinh sống ở Rú Chá ít nhất từ năm 1985 và cán bộ chính quyền xã Hương Phong với số lượng là 38 phiếu điều tra trong tổng số 117 hộ có đời sống liên quan đến Rú Chá (chiếm 32,5%).

2.2.3. Khảo sát thực địa

Tiến hành điều tra thành phần loài thực vật theo tuyến nghiên cứu, lập 5 ô tiêu chuẩn (ÔTC) kích thước 100m² (10m*10m), dùng để điều tra cây tầng cao có D_1,3 ≥ 5cm. Trong mỗi ÔTC, lập ra 5 ô dạng bản diện tích 4m² với kích thước (2m*2m) để điều tra cây bụi, thảm mục và vật rơi rụng (VRR); trong đó 4 ô ở 4 góc, 1 ô ở trung tâm ÔTC. Tổng cộng có 25 ô dạng bản (2m*2m). Điều tra theo 3 tuyến như sau: tuyến thứ 1 đi xuyên qua và vòng xung quanh Rú chính và lên trên Cồn Miếu với 3 ÔTC (A, B và C); tuyến thứ 2 đi xuyên qua và vòng quanh Rú dưới với 1 ÔTC D; tuyến thứ 3 đi xuyên qua và vòng quanh Rú trên với 1 ÔTC E. Sử dụng máy định vị vệ tinh GPS Map 78S (hãng GARMIN, Đài Loan) để xác định tọa độ các khu vực có TVNM.

Hình 1. Vị trí lấy mẫu thực vật (a) và đất (b) ở Rú Chá

2.2.4. Phương pháp phân tích và đánh giá nhanh trữ lượng carbon và hấp thụ CO₂

Lượng carbon (C) tích lũy trong RNM được ước tính tổng hợp từ các thành phần, gồm C tích lũy trong thực vật (cây tầng cao, cây bụi, thảm mục + VRR) và C trong đất. Do vậy, nghiên cứu tập trung áp dụng tổng hợp các phương pháp “đánh giá nhanh” để định lượng tương đối lượng C hiện tại tích lũy trong các lâm phần theo hướng dẫn của IPCC- Eggleston H.S., et al. (2006) [2],  Ketterings Q.M. và nnk (2001) [7].

2.2.5. Đánh giá giá trị kinh tế

Đánh giá giá trị kinh tế bằng cách điều tra, phỏng vấn người dân và thống kê số liệu từ chính quyền địa phương; từ đó ước tính ra các giá trị mà RNM mang lại. Giá trị kinh tế của RNM bao gồm: (1) Giá trị sử dụng trực tiếp (gỗ, củi đốt, khai thác và nuôi trồng thủy sản, khai thác chim nước, dược liệu, làm cảnh, thức ăn…) được xác định thông qua giá cả thị trưởng. (2) Giá trị sử dụng gián tiếp (điều hòa vi khí hậu, điều tiết nước ngầm, cung cấp nơi ở, chất dinh dưỡng, xử lý ô nhiễm, khả năng tích lũy carbon và hấp thụ CO₂…) được xác định bằng chi phí thay thế. Trong phạm vi bài báo này, chúng tôi tập trung xác định giá trị sử dụng gián tiếp về khả năng tích lũy carbon và hấp thụ CO₂ của RNM Rú Chá. (3) Giá trị sử dụng gián tiếp (phòng chống thiên tai, cản sức gió, chống xói mòn, bảo vệ đất …) được đánh giá thông qua chi phí thiệt hại tránh được; (4) giá trị chọn lựa, tồn tại trong việc sẵn lòng chi trả của người dân cho việc xây dựng quỹ bảo vệ, bảo tồn RNM Rú Chá thông qua đánh giá ngẫu nhiên.

2.2.6. Xử lý số liệu

Sử dụng phần mềm Excel để phân tích số liệu điều tra về kinh tế- xã hội; thống kê thành phần loài và tính toán giá trị kinh tế, tích lũy carbon và giảm phát thải khí nhà kính CO₂ của RNM.

III. KẾT QUẢ VÀ THẢO LUẬN

3.1. Hiện trạng và đa dạng thành phần loài TVNM ở Rú Chá

Nhóm khảo sát tiến hành điều tra theo tuyến để điều tra thành phần loài, tra cứu xác định tên khoa học các loài thực vật, sắp xếp các loài theo các đơn vị phân loại; đồng thời, đối chiếu với các tài liệu nghiên cứu trước đây như Phan Nguyên Hồng (1999) [5], Hoàng Công Tín và nnk (2012) [8], Dự án IMOLA II (2010) [10], Trần Thị Tú và nnk (2013) [9] (Hình 2).

Hình 2. So sánh sự đa dạng sinh học ở các khu vực phân bố TVNM ở Thừa Thiên Huế

Kết quả nghiên cứu đã xác định được 27 loài TVNM ở Rú Chá, ngành Dương Xỉ (Polypodiophyta) chỉ có 1 họ và 1 loài (chiếm 4,5% tổng số họ và 3,7% tổng số loài); ngành Ngọc Lan (Magnoliophyta) chiếm ưu thế với 21 họ, 25 chi và 26 loài (chiếm 95,5% tổng số họ và 96,3% tổng số loài). Trong ngành Ngọc Lan, số lượng các taxon trong lớp Ngọc Lan (Magnoliopsida) chiến ưu thế với 17 họ, 21 chi và 21 loài; lớp Loa kèn (Liliopsida) ít hơn với 4 họ, 4 chi và 5 loài. Trong số 27 loài TVNM ở Rú Chá, có 10 loài cây ngập mặn chính thức (MS) (chiếm 37% tổng số loài) và 17 loài cây tham gia ngập mặn (MAS) (chiếm 63% tổng số loài). Theo Phan Nguyên Hồng (1999) [5] thì ở Việt Nam có 34 loài cây ngập mặn thực sự và trên 40 loài cây tham gia vào RNM. Như vậy, số loài cây ngập mặn thực sự ở Rú Chá chỉ chiếm 29,4% tổng số loài cây ngập mặn thực sự ở Việt Nam. Các loài cây tham gia và di cư vào RNM Rú Chá tương đối nhiều, có 17 loài thuộc 16 chi, 12 họ có đại diện ở nơi đây, chiếm tới 42,5% tổng số loài của RNM. Các loài cây này thường gặp ở những nơi tiếp giáp với RNM, gần mép nước, ven bờ.

Dựa theo tiêu chuẩn đánh giá và quan trắc RNM của Phạm Nhật và nnk (2003) [6] và FAO (2007) [3], tính đa dạng của RNM thấp khi rừng có từ 1- 3 loài, trên 10 loài tính đa dạng của rừng cao, thì TVNM ở Rú Chá có độ đa dạng loài cao. Trong đó, Rú Chá có 5 loài TVNM chủ yếu (Giá (Excoecaria agallocha L.), Quao nước (Dolichandrone spathacea (L.f.) Schum.), Tra hoa vàng (Hibicus tiliaceus L.), Ngọc nữ biển (Clerodendrum inerme (L.) Gaertn.) và Ô rô trắng (Acanthus ebracteatus (L.) Vahl.)), thuộc 5 họ khác nhau đã có từ trước, trong đó có 3 loài thân gỗ, hai loài cây bụi. Ở đây, Giá (Excoecaria agallocha L.) được coi là loài cây ngập mặn tiên phong trong quá trình lấn biển. Ngoài ra, ba loài TVNM được trồng từ năm 2003 là Đước vòi, Vẹt dù và Sú.

Bảng 1. Diện tích phân bố TVNM Rú Chá      

Khu vực	Diện tích (ha)	Tỷ lệ (%)
Rú trên	0,68	11,7
Rú chính	3,56	61,3
Rú dưới	1,10	18,9
Cồn Miếu	0,47	8,1
Tổng cộng	5,81	100,0

Bảng 2. Diện tích Rú Chá bị mất do tác động của một số hoạt động chủ yếu

TT	Hoạt động ảnh hưởng đến TVNM	Diện tích bị mất (ha)
1	Hoạt động phát triển sản xuất nông nghiệp	6,00
	Nuôi trồng thủy sản	6,00
2	Hoạt động xây dựng cơ sở hạ tầng	2,65
	Xây dựng nhà ở	0,19
	Xây dựng nghĩa địa, mồ mả	2,46
3	San lấp đất rừng ngập mặn để xây dựng các công trình dân dụng (đường xá, đê kè,…)	0,90
4	Các hoạt động dân sinh khác	Không xác định được
	Tổng cộng	9,55

Qua phân tích dữ liệu bản đồ và điều tra khảo sát, TVNM ở Rú Chá phân bố chủ yếu ở 4 khu vực (Bảng 1). Diện tích TVNM hiện nay còn tập trung ở Rú chính, Rú dưới và Rú trên (chiếm 91,9%); một phần nhỏ ở Cồn Miếu (chiếm 8,1%). Tuy nhiên, tuyến đê ngăn mặn và vùng ao nuôi trồng thủy sản đã làm chia cắt sự kết nối liên tục giữa các rú với nhau. Các khu vực còn lại, thảm TVNM có diện tích quá nhỏ và phân bố rời rạc do bị chia cắt nên không xác định được diện tích hoặc bị chuyển sang mục đích sử dụng khác (xây dựng mồ mả) (Bảng 2).

3.2. Đánh giá giá trị kinh tế của RNM Rú Chá

3.2.1. Giá trị sử dụng trực tiếp của TVNM Rú Chá

Hình 3. Giá trị sử dụng của các loài TVNM Rú Chá

(* Một số loài có thể cho nhiều giá trị sử dụng khác nhau)

Bảng 3. Tỷ lệ khai thác và giá trị sử dụng trực tiếp từ RNM Rú Chá

Giá trị sử dụng dựa theo các tài liệu của Võ Văn Chi (1997) và Phạm Hoàng Hộ (2001). Với 27 loài có giá trị sử dụng chiếm 100% tổng số loài thực vật ngập mặn ở đây. Công dụng của các loài TVNM được trình bày ở Hình 3. Giá trị sử dụng trực tiếp của rừng ngập mặn Rú Chá ước tính khoảng 1.070.850.406 đồng/năm (Bảng 3).

3.2.2. Giá trị sử dụng gián tiếp

Tính đa dạng về thành phần loài, nhất là đa dạng về di truyền tạo cho sinh vật của Rú Chá sống ổn định trong môi trường thường xuyên biến động của bãi triều lầy, đồng thời giúp cho chúng tham gia vào các bậc dinh dưỡng khác nhau của hệ thống các mắt xích thức ăn, nhằm khai thác tối đa nguồn năng lượng và vật chất dưới dạng sản phẩm sơ cấp được phức hợp từ các CNM tạo ra trong quá trình quang hợp. Rú Chá không chỉ hình thành nên năng suất cao dưới dạng cây rừng mà hằng năm còn cung cấp một sản lượng vật rơi rụng khá lớn để làm giàu cho đất rừng và vùng cửa biển Thuận An. Kết quả nghiên cứu cho thấy lượng vật rơi rụng hàng năm ở RNM Rú Chá ước tính khoảng 0,017 tấn/ha/năm, tương ứng với khả năng hấp thụ 0,063 tấn CO₂/ha/năm.

Kết quả nghiên cứu ở Bảng 4, Hình 4 và Hình 5 cho thấy: (1) Sinh khối thực vật ở RNM Rú Chá đạt là 14,8 tấn/ha. (2) Lượng carbon tích lũy trong thực vật là 6,9 tấn/ha. Lượng CO₂ hấp thụ trong thực vật là 25,4 tấn/ha. Lượng carbon tích lũy trong đất trung bình đạt 43,2 tấn/ha. Tổng lượng CO₂ tích lũy trong đất là 158,5 tấn/ha. Như vậy, sự biến động hàm lượng carbon tích lũy trong đất tùy thuộc vào tính chất đất, độ chặt của đất và mức độ phân hủy xác bã thực vật. Khi hàm lượng carbon trong đất tăng thì đồng nghĩa với khả năng hấp thụ CO₂ cũng cao.

Bảng 4. Mật độ và trữ lượng gỗ cây tầng cao ở Rú Chá

Vị trí	Ô tiêu chuẩn	Mật độ, N/ha (cây/ha)	Trữ lượng gỗ, M (m3/ha)	Đường kính TB, DBHtb (cm)	Chiều cao TB, Htb (m)	Số cành/ cây	Đường kính tán, CD (m)	Độ tàn che (của 2 loài), (%)
Rú Chính	A	900	186,3	17,3	9,2	2,5	5,0	89,0
	B	1.000	160,3	14,3	8,0	4,5	4,5	85,0
	C	2.700	159,0	9,8	5,8	5,0	4,0	86,0
Rú Dưới	D	2.700	99,9	8,8	7,0	3,0	3,5	83,5
Rú Trên	E	2.500	89,2	8,5	5,7	6,0	2,5	81,0
Trung bình		1.960	138,9	11,7	7,1	4,2	3,9	84,9

Kết quả phân tích cho thấy, cây thân gỗ tầng cao có độ tuổi càng lớn thì khả năng tích lũy carbon cũng như hấp thụ CO₂ trong các bộ phận thân cao cao hơn nhiều so với cây thân gỗ có độ tuổi nhỏ và còn non. Bên cạnh đó, vùng đất ngập nước Rú Chá vẫn chịu ảnh hưởng của thủy triều và nước ngập, nên lượng thảm mục và VRR thường bị rửa trôi, sinh khối không nhiều như ở rừng trồng và rừng phòng hộ ở núi cao, hoặc vùng ít ngập nước. Như vậy, tổng lượng carbon tích lũy trong lâm phần chủ yếu từ nguồn carbon đất (chiếm 86,2%), phần còn lại là từ carbon thực vật (chiếm 13,8%) (Hình 4 và Hình 5). Số liệu về trữ lượng carbon thực vật phản ánh rõ quy luật giữa hiện trạng cấu trúc loài, mật độ và đường kính thân cây với lượng carbon tích lũy trong thành phần thực vật.

Hình 4. Tổng C và CO­2 trong đất

Hình 5. Tổng C và CO2 tích lũy ở Rú Chá

Bảng 5. Giá trị tích lũy carbon và giá trị kinh tế của RNM Rú Chá năm 2015

Bảng 6. Giá trị du lịch sinh thái rừng ngập mặn Rú Chá

Số người tham quan TB (người/tuần)	Thời gian du lịch từ tháng 4 – 7	Chi phí chi trả du lịch (đồng/ tuần) năm 2015	Giá trị du lịch sinh thái (đồng/năm)
25	4	150.000	60.000.000

Giá trị tích lũy carbon do TVNM mang lại ước tính khoảng 5.506.795 đồng/ha/năm và của toàn bộ rừng ngập mặn Rú Chá với diện tích 5,81 ha là khoảng 39.862.118 đồng/ha/năm,. Tổng giá trị tích lũy carbon của RNM Rú Chá đạt 231.598.908 đồng/năm. (Bảng 5). Ước tính giá trị du lịch sinh thái rừng ngập mặn Rú Chá thu được là 60.000.000 đồng/năm (Bảng 6).

3.2.3. Giá trị phi sử dụng

Bảng 7. Mức sẵn lòng chi trả của người dân cho quỹ 1 và quỹ 2

Bảng 8. Nguồn tài trợ từ các chương trình dự án trong và ngoài nước

Bảng 9. Tổng giá trị kinh tế của rừng ngập mặn Rú Chá

Giá trị phi sử dụng bao gồm giá trị chọn lựa (Value 1), giá trị để lại (Value 2) và giá trị tồn tại (A). Ước tính giá trị phi sử dụng của rừng ngập mặn Rú Chá là 58.864.620 đồng/năm (Bảng 7 và Bảng 8). Như vậy, tổng giá trị kinh tế mà rừng ngập mặn Rú Chá mang lại bao gồm giá trị sử dụng trực tiếp, giá trị sử dụng gián tiếp và giá trị phi sử dụng ước tính khoảng 1.238.721.618 đồng/năm (Bảng 9).

KẾT LUẬN

Kết quả nghiên cứu hiện trạng thành phần loài TVNM ở Rú Chá đã xác định được 27 loài thuộc 26 chi, 22 họ của 2 ngành Dương xỉ (Polypodiophyta) và Ngọc lan (Magnoliophyta). Diện tích RNM Rú Chá bị tác động mạnh bởi bốn nguyên nhân: sản xuất nông nghiệp, xây dựng cơ sở hạ tầng, san lấp đất RNM và hoạt động khác. Rú Chá có nhiều giá trị kinh tế và giá trị sử dụng khác nhau. Tổng lượng carbon tích lũy ước tính trung bình đạt 50,2 tấn/ha; tiềm năng giảm phát thải khí nhà kính trung bình của Rú Chá đạt 184,0 tấn CO₂/ha, tương ứng với 39,86 triệu đồng/ha/năm. Như vậy tổng giá trị giảm phát thải khí nhà kính CO₂ của Rú Chá đạt 231.598.908 đồng/năm.

 

TÀI LIỆU THAM KHẢO

• Võ Văn Chi (2007), Sách tra cứu tên cây cỏ Việt Nam, Nxb. Giáo dục, Hà Nội.
• Eggleston H.S., Buendia L., Miwa K., Ngara T., (eds) (2006), IPCC Guidelines for National Greenhouse Gas Inventories, Prepared by the Natinal Greenhouse Gas Inventories Programme, Published in IGES, Japan.
• FAO and Wetlands International (2007), Mangrove Guidebook for Southeast Asia, Printed by Dharmasarn Co. Ltd.
• Phạm Hoàng Hộ (2001), Cây cỏ Việt Nam, tập I, II, III. Nxb. Trẻ, Tp. Hồ Chí Minh.
• Phan Nguyên Hồng (1999), Rừng ngập mặn Việt Nam, Nxb. Nông nghiệp, Hà Nội.
• Phạm Nhật và nnk (2003), Sổ tay Hướng dẫn điều tra và giám sát đa dạng sinh học, Nxb. Giao thông vận tải, 315-331.
• Ketterings Q.M., Coe R., Van Noordwijk M., Reducing uncertainty in the use of allometric biomass equations for predicting above ground tree biomass in mixed secondary forests, Forest Ecology and Management Journal, (120), (2001), 199-209.
• Hoàng Công Tín, Mai Văn Phô, 2012. Thành phần loài và đặc điểm phân bố của thực vật ngập mặn ở Thừa Thiên Huế, Tạp chí Khoa học Lâm nghiệp, số 1, 2085-2094.
• Trần Thị Tú, Phạm Ngọc Dũng, Trần Hiếu Quang, 2013. Đánh giá những tác động của hoạt động kinh tế- xã hội đến thực vật ngập mặn ở đầm Lập An, huyện Phú Lộc, tỉnh Thừa Thiên Huế, Tuyển tập các công trình khoa học Viện Môi trường và Phát triển Bền vững (1995-2012), Nxb. Khoa học và Kỹ thuật, Hà Nội, 421-429, ISBN 978-604-670069-2.
• UBND tỉnh Thừa Thiên Huế (2010), Báo cáo tổng kết Nghiên cứu rừng ngập mặn Rú Chá, Hương Phong, Thừa Thiên Huế, Dự án Quản lý tổng hợp hoạt động đầm phá- Dự án IMOLA II, Huế.

………………………………………………………………………………….

Nguồn: Trích dẫn từ kỷ yếu Hội nghị môi trường toàn quốc lần thứ IV, Bộ tài nguyên và Môi trường, Hà Nội, 29/09/2015

…………………………………………………………………………………..