Tổng quan nghiên cứu
Hệ sinh thái rừng ngập mặn tại Việt Nam đóng vai trò xung yếu trong việc phòng hộ bờ biển, giảm thiểu rủi ro thiên tai và bảo tồn đa dạng sinh học cửa sông. Theo thống kê giai đoạn 1983 - 2008, diện tích rừng ngập mặn toàn quốc sụt giảm từ 252.000 ha xuống còn 156.108 ha, tương đương mức suy giảm khoảng 62% diện tích so với thời điểm năm 1943. Trước áp lực biến đổi khí hậu toàn cầu, việc phục hồi và phát triển rừng ngập mặn trở thành ưu tiên chiến lược hàng đầu của ngành lâm nghiệp. Tại tỉnh Thái Bình, với hơn 4.256 ha rừng ngập mặn hiện hữu, loài bần chua (Sonneratia caseolaris L.) giữ vị thế cây tiên phong cố định bãi bồi ven biển nhờ hệ thống rễ thở phát triển mạnh và khả năng thích nghi cao với môi trường nước lợ.
Tuy nhiên, công tác quản lý và phát triển rừng bần chua tại địa phương còn gặp nhiều rào cản kỹ thuật do thiếu hụt hệ thống cơ sở dữ liệu định lượng về động thái lập địa, cấu trúc lâm phần và mô hình sinh trưởng thực nghiệm. Nghiên cứu này được tiến hành tại xã Thụy Xuân và Thụy Hải thuộc huyện Thái Thụy, tỉnh Thái Bình — khu vực có diện tích đất ngập nước ven biển khoảng 4.404 ha trải dài trên 27 km bờ biển, nằm trong Khu dự trữ sinh quyển đồng bằng sông Hồng. Đề tài tập trung phân tích thực trạng sinh trưởng của rừng trồng thuần loài qua 5 cấp tuổi (5, 7, 9, 11 và 13 tuổi). Mục tiêu cụ thể là xác lập quy luật phân bố số cây theo đường kính thân, chiều cao vút ngọn, đường kính tán và mô hình hóa quá trình tăng trưởng sinh khối. Kết quả nghiên cứu cung cấp hệ thống luận cứ khoa học chuẩn xác, hỗ trợ nâng cao hiệu quả phòng hộ giảm sóng từ 30% đến 50% cho hệ thống đê biển quốc gia.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết áp dụng
Nghiên cứu ứng dụng nền tảng lý thuyết cấu trúc lâm phần theo quan điểm sinh thái học rừng nhiệt đới của Davis và Richard, kết hợp với các mô hình tăng trưởng sinh vật kinh điển trong toán tin lâm nghiệp. Hệ thống lý thuyết trọng tâm bao gồm:
- Lý thuyết mô hình hóa sinh trưởng rừng: Quá trình tích lũy sinh khối của cây gỗ theo thời gian được lượng hóa dưới dạng hàm toán học $y = f(A)$, trong đó các hàm phi tuyến như hàm Schumacher và hàm Gompertz được ứng dụng để mô tả động thái đường kính, chiều cao và thể tích thân cây.
- Lý thuyết phân bố cấu trúc lâm phần: Phản ánh quy luật sắp xếp cá thể trong không gian dinh dưỡng thông qua các hàm phân bố tần số thực nghiệm về đường kính ngang ngực, chiều cao vút ngọn và đường kính tán.
- Khái niệm sinh trưởng và tăng trưởng: Sinh trưởng biểu thị sự tích lũy kích thước tuyệt đối theo tuổi lâm phần, trong khi tăng trưởng biểu thị lượng biến thiên kích thước định kỳ hàng năm.
- Khái niệm độ thành thục của đất ngập mặn: Tỷ số giữa hàm lượng nước và hàm lượng hạt sét mịn trong bùn bãi bồi, chi phối trực tiếp đến khả năng đâm sâu của rễ hô hấp.
- Khái niệm rừng phòng hộ đê biển: Hệ thống quần xã thực vật ven bờ có cấu trúc phân tầng dày đặc, thực hiện chức năng cản sóng triều và giữ lắng phù sa.
Phương pháp nghiên cứu
Nghiên cứu thu thập dữ liệu sơ cấp tại thực địa thông qua phương pháp chọn mẫu điển hình ngẫu nhiên trên 5 cấp tuổi rừng (5, 7, 9, 11 và 13 tuổi). Tại mỗi cấp tuổi, 3 ô tiêu chuẩn định hình tròn diện tích 500 m² (bán kính $R = 12,62$ m) được thiết lập, tạo nên tổng dung lượng mẫu gồm 15 ô tiêu chuẩn trên toàn khu vực. Toàn bộ hoạt động đo đếm thực hiện tại thời điểm triều rút nhằm bảo đảm tính đồng nhất kỹ thuật. Các biến số điều tra gồm: đường kính ngang ngực tại vị trí 1,3 m đo bằng thước kẹp chuyên dụng, chiều cao vút ngọn đo bằng sào đo cao lâm sinh, và đường kính tán đo theo 2 hướng vuông góc chính vị.
Nhằm xác định chính xác động thái tăng trưởng thể tích thân cây qua từng năm, phương pháp giải tích thân cây đã được triển khai trên 3 cây tiêu chuẩn đại diện. Các cây mẫu sau khi hạ ngả được cưa thành các thớt đĩa tại các vị trí cách nhau 1,0 m để đếm vòng năm và tính thể tích phân đoạn theo công thức Huber. Mẫu đất được lấy tại 2 tầng phát sinh (0 - 20 cm và 20 - 40 cm) theo tiêu chuẩn TCVN 5979:2007 và TCVN 6498:1999 để xác định thành phần cơ giới, độ chua và hàm lượng dưỡng chất khoáng.
Toàn bộ chuỗi số liệu được phân tích thống kê trên phần mềm chuyên dụng. Phép chia tổ theo công thức Brooks - Caruthers được sử dụng để lập biểu phân bố thực nghiệm. Mức độ tương quan và độ tin cậy của phương trình hồi quy sinh trưởng được kiểm định nghiêm ngặt qua chỉ số $R^2$, trắc nghiệm Student ($t$-test) và trắc nghiệm Fisher ($F$-test) ở mức ý nghĩa $\alpha = 0,05$ và $\alpha = 0,01$.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Phân tích số liệu thực nghiệm trên 15 ô tiêu chuẩn cho thấy rừng trồng bần chua tại huyện Thái Thụy sinh trưởng ổn định và có sự phân hóa cấu trúc rõ rệt qua các giai đoạn tuổi:
- Động thái sinh trưởng đường kính và chiều cao: Đường kính bình quân tăng đều từ 9,50 cm ở tuổi 5 lên 10,58 cm ở tuổi 7, đạt 11,46 cm ở tuổi 9, đạt 12,32 cm ở tuổi 11 và chạm mức 13,89 cm ở tuổi 13 (mức tăng tích lũy 46,21%). Tương ứng, chiều cao vút ngọn bình quân phát triển từ 7,76 m (tuổi 5) lên 8,63 m (tuổi 9) và đạt cực đại 9,62 m (tuổi 13), thể hiện mức gia tăng chiều cao 23,97%.
- Sự biến đổi đường kính tán lá: Tán cây mở rộng dần từ 1,75 m ở tuổi 5 lên 1,98 m ở tuổi 9 và đạt 2,21 m ở tuổi 13, đóng góp trực tiếp vào khả năng giao tán sớm của lâm phần.
- Biến động mật độ lâm phần: Mật độ rừng dao động mạnh từ 1.250 cây/ha ở tuổi 5, đạt đỉnh 1.400 cây/ha ở tuổi 9 do tỷ lệ cây con sống sót cao, sau đó suy giảm nhanh xuống còn 1.060 cây/ha ở tuổi 13 (giảm 24,29% so với tuổi 9). Hiện tượng sụt giảm này phản ánh quá trình tự tỉa thưa tự nhiên khốc liệt khi tán rừng khép kín.
- Đặc tính lập địa thổ nhưỡng: Đất rừng có thành phần cơ giới sét pha thịt với tỷ lệ cấp hạt sét chiếm 34,30% đến 58,27%, cấp hạt thịt từ 22,27% đến 47,77%. Chỉ số độ chua $\text{pH}_{\text{KCl}}$ biến động trong khoảng 6,37 đến 6,80 (đất ít chua). Hàm lượng lân dễ tiêu rất phong phú (53,24 - 64,26 mg/100g đất) và kali dễ tiêu dồi dào (64,08 - 87,53 mg/100g đất), tạo điều kiện dinh dưỡng khoáng thuận lợi cho cây bần chua sinh trưởng nhanh.
Thảo luận kết quả
Cấu trúc phân bố số cây theo đường kính và chiều cao phản ánh trực quan động thái cạnh tranh không gian sống trong quần xã rừng ngập mặn. Ở tuổi 5, biểu đồ phân bố có dạng 1 đỉnh đối xứng tại giá trị 9,6 cm với hệ số biến động thấp ($C_v = 6,4%$), chứng minh lâm phần ở giai đoạn chuẩn bị giao tán có mức độ đồng đều cao, sự chênh lệch kích thước chủ yếu do chất lượng cây giống xuất vườn ban đầu. Từ tuổi 7 đến tuổi 11, các bảng phân bố thực nghiệm bắt đầu trải rộng với biên độ biến động đường kính đạt 12 - 16 cm do sự cạnh tranh ánh sáng gay gắt. Đến tuổi 13, biểu đồ đường cong phân bố lệch hẳn về phía bên phải, với nhóm cây có đường kính 13 - 16 cm chiếm tỷ trọng chi phối khoảng 45%, phản ánh sự hình thành rõ nét của các tầng tán ưu thế và bị ức chế.
Khi mô hình hóa dữ liệu tăng trưởng, hàm Schumacher thể hiện độ phù hợp vượt trội so với các dạng phương trình khác với hệ số xác định $R^2 > 0,95$, vượt qua các kiểm định $t$ và $F$ ở mức ý nghĩa cao. Đường cong sinh trưởng lý thuyết cho thấy sau tuổi 9, tốc độ tăng trưởng định kỳ về chiều cao bắt đầu chậm lại, trong khi quá trình tích lũy thể tích thân cây trên mỗi hecta diễn ra mạnh mẽ hơn. So sánh với các nghiên cứu trên rừng đước (Rhizophora apiculata) tại Nam Bộ hay keo lá tràm tại Đông Nam Bộ, loài bần chua vùng đồng bằng Bắc Bộ có chu kỳ khép tán sớm hơn khoảng 2 năm nhờ lượng phù sa bồi lắng lớn từ hệ thống cửa sông Trà Lý và sông Thái Bình (khoảng 15 đến 20 triệu tấn bùn cát mỗi năm).
D1.3 (cm) | [Tuổi 5: 9.50] -> [Tuổi 7: 10.58] -> [Tuổi 9: 11.46] -> [Tuổi 11: 12.32] -> [Tuổi 13: 13.89]
Hvn (m) | [Tuổi 5: 7.76] -> [Tuổi 7: 8.13] -> [Tuổi 9: 8.63] -> [Tuổi 11: 8.96] -> [Tuổi 13: 9.62]
N (cây/ha)| [Tuổi 5: 1250] -> [Tuổi 7: 1390] -> [Tuổi 9: 1400] -> [Tuổi 11: 1220] -> [Tuổi 13: 1060]
Đề xuất và khuyến nghị
Dựa trên các chỉ số định lượng về quy luật sinh trưởng, bốn giải pháp kỹ thuật lâm sinh then chốt được đề xuất nhằm nâng cao chất lượng rừng phòng hộ ven biển huyện Thái Thụy:
- Thực hiện tỉa thưa nuôi dưỡng lâm phần đúng thời điểm: Ban quản lý rừng phòng hộ và Chi cục Kiểm lâm Thái Bình cần tiến hành biện pháp tỉa thưa định hình ở giai đoạn 7 - 9 tuổi nhằm giảm mật độ từ 1.400 cây/ha xuống mức tối ưu khoảng 1.100 - 1.200 cây/ha. Việc can thiệp sớm giúp giảm áp lực cạnh tranh dinh dưỡng tự nhiên và hạn chế tỷ lệ cây chết đứng do thiếu sáng.
- Triển khai kỹ thuật chặt vệ sinh và mở sáng tầng tán: Định kỳ 2 - 3 năm một lần đối với các lâm phần từ 11 đến 13 tuổi, chính quyền các xã Thụy Xuân và Thụy Hải cần tổ chức chặt dọn các cá thể sâu bệnh, cây cong queo, còi cọc có đường kính dưới 8 cm. Biện pháp này giải phóng không gian sinh dưỡng, tập trung năng lượng sinh trưởng cho 100% số cây ưu thế trong quần thụ.
- Chuẩn hóa tiêu chuẩn chọn giống và kỹ thuật ươm cây con: Các cơ sở sản xuất giống lâm nghiệp trên địa bàn cần tuyển chọn hạt giống từ các cây mẹ khỏe mạnh trên 10 năm tuổi, xử lý gieo ươm trong môi trường nước lợ có độ mặn tiêu chuẩn từ 1,0% đến 1,5%. Xuất vườn cây con đạt chiều cao tối thiểu 50 cm và đường kính cổ rễ trên 0,8 cm nhằm đảm bảo tỷ lệ sống đạt trên 85% sau khi trồng ngoài bãi triều.
- Xây dựng biểu tra cứu sản lượng và số hóa dữ liệu quản lý rừng: Sở Nông nghiệp và Phát triển nông thôn tỉnh Thái Bình cần ứng dụng phương trình tương quan sinh trưởng Schumacher vào phần mềm quản lý tài nguyên rừng giai đoạn 2025 - 2030, phục vụ việc dự báo sinh khối, trữ lượng carbon và lập kế hoạch chuyển hóa rừng phòng hộ bền vững.
Đối tượng nên tham khảo luận văn
Luận văn thạc sĩ này mang lại giá trị khoa học và thực tiễn phong phú cho nhiều nhóm độc giả chuyên ngành:
- Cơ quan quản lý lâm nghiệp và Chi cục Kiểm lâm địa phương: Nắm bắt chính xác các thông số mật độ và tốc độ tăng trưởng của cây bần chua để xây dựng định mức kinh tế kỹ thuật trong các dự án trồng rừng ngập mặn chắn sóng tại dải ven biển Bắc Bộ.
- Giảng viên và nhà nghiên cứu lâm sinh học: Khai thác kho dữ liệu thực nghiệm về giải tích thân cây ngập mặn và phương pháp ứng dụng các hàm toán học phi tuyến (Schumacher, Gompertz) trong mô hình hóa cấu trúc rừng nhiệt đới.
- Học viên cao học và sinh viên ngành Quản lý tài nguyên rừng: Tiếp cận quy trình điều tra thực địa mẫu chuẩn, kỹ thuật xử lý phân bố xác suất sinh học qua công thức Brooks - Caruthers và phương pháp phân tích thổ nhưỡng ngập mặn.
- Ban quản lý các dự án biến đổi khí hậu và tổ chức bảo tồn: Vận dụng số liệu sinh trưởng chiều cao và đường kính tán để ước tính sinh khối tích lũy, tính toán khả năng hấp thụ $\text{CO}_2$ và thiết kế các công trình đai rừng phòng hộ sinh thái ven biển hiệu quả.
Câu hỏi thường gặp
Cây bần chua sinh trưởng tối ưu trong điều kiện độ mặn và nền đất như thế nào?
Loài bần chua phát triển tốt nhất ở môi trường nước lợ cửa sông có độ mặn dao động từ 1,0% đến 1,5% (trong khoảng 10‰ đến 15‰) và nhiệt độ nước từ 22°C đến 27°C. Đất ngập mặn phù hợp là dạng sét pha thịt có tỷ lệ hạt sét từ 40% đến 50%, độ chua $\text{pH}_{\text{KCl}}$ từ 6,3 đến 6,8 và giàu hàm lượng khoáng lân, kali dễ tiêu.
Tại sao mật độ rừng bần chua lại sụt giảm mạnh từ tuổi 9 đến tuổi 13?
Từ tuổi 9 trở đi, đường kính tán trung bình của cây bần chua đạt trên 1,98 m khiến tán rừng khép kín hoàn toàn. Các cá thể cây bắt đầu cạnh tranh gay gắt về ánh sáng và tầng đất dinh dưỡng. Khi không được tỉa thưa nhân tạo kịp thời, quá trình đào thải tự nhiên diễn ra làm mật độ giảm từ 1.400 cây/ha xuống còn 1.060 cây/ha ở tuổi 13.
Hàm toán học nào mô tả chính xác nhất quy luật sinh trưởng của rừng bần chua tại Thái Bình?
Hàm Schumacher là mô hình lý thuyết thích hợp nhất để biểu diễn mối tương quan giữa tuổi lâm phần với đường kính ngang ngực, chiều cao vút ngọn và thể tích thân cây. Phương trình đạt hệ số tương quan thực nghiệm $R^2 > 0,95$ và vượt qua các phép kiểm định thống kê $t$-test và $F$-test ở mức xác suất tin cậy trên 99%.
Khi nào cần can thiệp các biện pháp kỹ thuật tỉa thưa mở sáng cho rừng bần chua?
Thời điểm can thiệp lâm sinh tối ưu là khi rừng đạt từ 7 đến 9 năm tuổi. Lúc này rừng bắt đầu bước vào giai đoạn khép tán hoàn toàn. Việc chặt tỉa khoảng 15% đến 20% mật độ cây còi cọc giúp mở rộng không gian dinh dưỡng cho các cây có phẩm chất tốt phát triển tán lá tròn rộng.
Đai rừng bần chua tại huyện Thái Thụy đem lại những lợi ích bảo vệ bờ biển cụ thể nào?
Nhờ hệ thống rễ hô hấp dày đặc mọc nhô khỏi mặt bùn và tán lá rộng trung bình 2,21 m ở tuổi 13, rừng bần chua có khả năng làm giảm từ 30% đến 50% năng lượng sóng biển khi triều cường dâng cao trên 3,0 m, đồng thời giữ lắng khoảng 15 đến 20 triệu tấn phù sa bồi tụ mỗi năm dọc các cửa sông.
Kết luận
- Rừng trồng bần chua tại huyện Thái Thụy có tốc độ sinh trưởng đường kính bình quân đạt 13,89 cm và chiều cao đạt 9,62 m sau 13 năm trồng.
- Cấu trúc lâm phần chuyển đổi rõ nét từ phân bố chuẩn đồng đều ở tuổi 5 sang phân bố lệch phải phân hóa mạnh ở tuổi 13.
- Hàm tăng trưởng Schumacher phản ánh chính xác nhất quy luật phát triển sinh khối của bần chua với độ tin cậy $R^2 > 0,95$.
- Quá trình tự tỉa thưa tự nhiên làm mật độ lâm phần giảm 24,29% sau giai đoạn khép tán, đòi hỏi can thiệp tỉa thưa kỹ thuật ở tuổi 7 - 9.
- Thổ nhưỡng sét pha thịt tại địa phương chứa hàm lượng kali dồi dào ($64,08 - 87,53$ mg/100g đất), tạo điều kiện lý tưởng để nhân rộng mô hình đai rừng phòng hộ ven biển.
Nghiên cứu đã cung cấp hệ thống dữ liệu sinh trắc học hoàn chỉnh, làm cơ sở khoa học cho công tác quản lý và bảo vệ rừng ngập mặn ven biển miền Bắc. Kế hoạch tiếp theo cần tập trung số hóa dữ liệu sinh trưởng vào hệ thống bản đồ số và hoàn thiện quy trình tỉa thưa nuôi dưỡng trước năm 2026. Các cơ quan quản lý và địa phương ven biển cần áp dụng ngay các khuyến nghị kỹ thuật này để tối ưu hóa năng lực phòng hộ đê điều trước tác động của biến đổi khí hậu.