Nghiên cứu đặc tính sinh thái của cây tái sinh vên vên (anisoptera cochinchinensis pierre) trong kiểu rừng kín thường xanh và nửa rụng lá ẩm nhiệt đới ở Đồng Nai

NGHIÊN CỨU ĐẶC TÍNH SINH THÁI CỦA CÂY TÁI SINH VÊN VÊN (Anisoptera cochinchinensis Pierre) TRONG KIỂU RỪNG KÍN THƯỜNG XANH VÀ NỬA RỤNG LÁ ẨM NHIỆT ĐỚI Ở ĐỒNG NAI KS. Phạm Văn Hường – Cơ sở 2 Trường Đại học Lâm nghiệp TÓM TẮT Đông Nam Bộ nói chung và Đồng Nai nói riêng cây họ Sao – Dầu trong hệ sinh thái rừng có vai trò quan trọng, đặc biệt là những loài cây có giá trị kinh tế, sinh thái và bảo tồn. Vên vên (Anisoptera cochinchinensis) là loài cây có được những tính chất đó. Thật vậy, vên vên là loài cây đang được nghiên cứu để lựa chọn đưa vào trồng rừng để khôi phục vốn rừng. Từ trước đến nay cũng đã có một số công trình nghiên cứu về cây họ Sao – Dầu hoặc vên vên, tuy nhiên chưa làm rõ được quy luật sống, quá trình tái sinh, đặc tính sinh thái của chúng. Cho nên việc nghiên cứu đặc tính sinh thái cây tái sinh vên vên trong kiểu rừng kín thường xanh và nửa rụng lá ở Đồng Nai là việc làm hết sức có ý nghĩa. Mục đích của bài viết này là đưa ra những cơ sở dữ liệu để làm rõ đặc tính sinh thái của cây tái sinh vên vên trong kiểu rừng kín thường xanh và nửa rụng lá ẩm nhiệt đới ở Đồng Nai và là căn cứ khoa học cho việc áp dụng những biện pháp tái sinh (tự nhiên và nhân tạo) và nuôi dưỡng vên vên. Kết quả nghiên cứu cho thấy: - Độ phong phú của tái sinh vên vên ở 2 cấp tuổi đều phụ thuộc vào trạng thái rừng, độ ẩm đất và độ tàn che tán rừng. Độ ẩm đất thích hợp cho cấp tuổi 1 từ 61,8 - 82,3%, tối ưu là 73,0%; ở cấp tuổi 2 tương ứng là 62,8 - 83,9% và 73,0%. Độ tàn che tán rừng thích hợp cho cấp tuổi 1 là 0,65 – 0,85, tối ưu 0,75; ở cấp tuổi 2 từ 0,63 - 0,88, tối ưu là 0,75. - Độ phong phú của tái sinh vên vên thay đổi tùy thuộc vào sự thay đổi của 2 yếu tố độ ẩm đất và độ tàn che tán rừng. Mối liên hệ này có thể mô tả bằng những mô hình Logit Gauss 2 biến số. - Khi trạng thái rừng thay đổi, thì độ ẩm đất tầng đất mặt và độ tàn che tán rừng có ảnh hưởng khác nhau đến độ phong phú của vên vên. Những trạng thái rừng ổn đinh cao (IIIA2 và IIIA3) đảm bảo cho cây vên vên tái sinh tốt hơn so với những trạng thái rừng kém ổn định (IIB và IIIA1). Từ khoá: Cây tái sinh vên vên, độ phong phú, xác suất bắt gặp, độ ẩm đất, độ tàn che và Mô hình Logit Gauss ĐẶT VẤN ĐỀ Hiện nay một nhiệm vụ cấp thiết đang được đặt ra đối với ngành lâm nghiệp là khôi phục lại vốn rừng, nâng cao năng suất và chất lượng rừng so với tiềm năng vốn có của chúng. Nhưng muốn hoàn thành tốt nhiệm vụ đó, đòi hỏi phải có những hiểu biết đầy đủ về bản chất các quy luật sống của rừng, trước hết là các quá trình tái sinh, sự hình thành và động thái biến đổi của rừng tương ứng với những điều kiện môi trường tự nhiên khác nhau. Vì lý do đó, việc đi sâu nghiên cứu làm rõ quy luật phát sinh, sinh trưởng và phát triển của cây con; phân tích những ảnh hưởng của các điều kiện môi trường đến động thái tái sinh dưới tán rừng của cây vên vên là một việc làm cần thiết và cấp bách. Tài nguyên rừng của khu vực Đông Nam bộ nói chung, Đồng nai nói riêng với những loài cây họ Sao – Dầu chiếm giữ vai trò quan trọng trong hệ sinh thái rừng cũng như năng suất, trữ lượng rừng [2], [7], [8]. Đặc biệt là những loài cây có giá trị về kinh tế, sinh thái và bảo tồn như: vên vên, dầu song nàng, dầu con rái, sao đen,… Do vậy, trước mắt việc khôi phục lại vốn rừng, nâng cao năng suất và chất lượng rừng của khu vực cần tập trung vào việc gây trồng, phát triển và khôi phục lại những quần thụ cây họ Sao – Dầu, cụ thể là cây vên vên là rất cấp bách. Mặc dù trước đây, ở nước ta đã có một số công trình nghiên cứu về đặc tính sinh thái của một số loài cây thuộc họ Sao - Dầu ở Đông Nam Bộ ([2], [3], [4], [5], [7], [8], [10]), nhưng phạm vi và đối tượng nghiên cứu còn hạn chế, chưa đi sâu tìm hiểu quy luật sống của cây họ Sao – Dầu, đặc tính tái sinh… Đặc biệt là đối với cây vên vên, hiện tại chưa có nhiều công trình nào nghiên cứu sâu về loài. Do đó, việc kế thừa những tài liệu đã có và tiếp tục đi sâu nghiên cứu đặc tính sinh thái tái sinh tự nhiên của cây vên vên là việc làm cần thiết. Vì lý do đó, cho nên Nghiên cứu đặc tính sinh thái cây tái sinh vên vên trong kiểu rừng kín thường xanh và nửa rụng lá ẩm nhiệt đới ở Đồng Nai là rất có ý nghĩa. 1. VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU 1.1. Vật liệu, dụng cụ nghiên cứu Vật liệu sử dụng trong nghiên cứu là cây vên vên tái sinh tự nhiên trong kiểu rừng kín thường xanh và nửa rụng lá ẩm nhiệt đới. Dụng cụ sử dụng để thu thập số liệu là máy đo độ ẩm nhanh, thước dây 100m, thước kẹp kính, thước Blum-lei, dây nilon và các mẫu bảng biểu khác. 1.2. Phương pháp nghiên cứu 1.2.1. Phương pháp luận Cơ sở khoa học của phương pháp luận là dựa trên những quan niệm sau đây: (1) Sự phát sinh và phát triển cây tái sinh dưới tán rừng luôn bị kiểm soát bởi tập hợp nhiều yếu tố; trong đó một số yếu tố giữ vai trò chủ đạo, còn những yếu tố khác chỉ có vai trò thứ yếu [1], [6], [11]. Mặt khác, trong thực tiễn nhà lâm học cũng chỉ quan tâm đến những nhân tố mà con người có thể kiểm soát được thông qua những hoạt động lâm sinh. Trong số những yếu tố môi trường có ảnh hưởng lớn đến quá trình hình thành và phát triển mà con người có thể kiểm soát được, thì độ tàn che tán rừng và đặc tính lý - hóa của đất (độ ẩm đất mặt) là những yếu tố chủ yếu. Vì thế, độ phong phú của những lớp cây tái sinh tự nhiên dưới tán rừng phải được xem xét trong quan hệ với độ tàn che tán rừng và độ ẩm của tầng đất mặt. (3) Độ phong phú của những loài cây gỗ cũng thay đổi tùy theo trạng thái của quần xã hay trạng thái rừng. Ở những quần xã bị rối loạn (tương ứng với các kiểu trạng thái rừng IIA, IIB và IIIA1), độ phong phú của các loài cây gỗ có thể bị thay đổi mạnh. Nguyên nhân là do những rối loạn trong quần xã đã làm thay đổi điều kiện sống và số lượng cá thể của các loài. Ngoài ra, do trạng thái rừng khác nhau, nên mối quan hệ của các loài cây gỗ với với các yếu tố môi trường cũng thay đổi. Vì thế, khi xem xét mối quan hệ giữa độ phong phú của loài với các yếu tố môi trường thay đổi, tác giả cũng hướng vào xem xét độ bắt gặp cây tái sinh vên vên tùy thuộc vào trạng thái của quần xã hay trạng thái rừng. Những trạng thái rừng được chọn để so sánh độ phong phú của loài với các yếu tố môi trường thay đổi là IIB, IIIA1, IIIA2 và IIIA3. Như vậy, bằng cách xem xét mối quan hệ giữa độ phong phú của loài trong những quần xã khác nhau với các yếu tố môi trường thay đổi, có thể xác định được ảnh hưởng của từng nhân tố sinh thái trong các thái rừng đến độ bắt gặp cây vên vên. (4) Quan hệ giữa độ phong phú của cây vên vên với các yếu tố môi trường còn thay đổi tùy theo tuổi. Vì thế, xác định tối ưu sinh thái và tính chống chịu của các loài đối với các yếu tố môi trường ở những giai đoạn tuổi khác nhau là cần thiết [12]. (5) Khi xem xét mối quan hệ giữa độ phong phú của loài với các yếu tố môi trường, thì “Độ phong phú” của loài được sử dụng theo nghĩa hẹp, nghĩa là nó chỉ biểu thị độ bắt gặp loài cây (bắt gặp = 1, không bắt gặp = 0) trên những ô mẫu có kích thước nhất định. Theo cách thu thập số liệu về độ phong phú của loài như thế, nên cách xử lý số liệu thích hợp nhất để phân tích mối quan hệ của loài với môi trường là hồi quy logit. Kết quả tính toán hồi quy logit chỉ ra xác suất bắt gặp loài tùy theo mức độ biến đổi của biến môi trường. Theo đó, khi môi trường sống thích hợp thì xác suất bắt gặp loài sẽ có trị số cao. Ngược lại, khi môi trường sống không thích hợp thì xác suất bắt gặp loài sẽ có trị số thấp. 1.2.2. Phương pháp thu thập số liệu Trình tự đo đạc như sau: (1) Trước hết, bố trí những tuyến cắt ngang qua bốn trạng thái rừng (IIB, IIIA1, IIIA2 và IIIA3). Mỗi tuyến có bề rộng 20 m, chiều dài tùy thuộc vào trạng thái rừng. (2) Kế đến, trên mỗi tuyến cứ sau 100 m lại bố trí 1 ô mẫu với kích thước 5*5 m hay 25 m2. Dự kiến mỗi trạng thái rừng cần đo đạc 30-50 ô mẫu 25 m2. Trong mỗi ô mẫu, những dấu hiệu đo đạc bao gồm độ bắt gặp cây tái sinh vên vên và 2 biến môi trường là độ tàn che tán rừng, độ ẩm tầng đất mặt. Độ bắt gặp loài được ghi nhận bằng hai biến định danh đó là bắt gặp (mã hóa = 1) và không bắt gặp (mã hóa = 0). Để thấy rõ đặc tính sinh thái của cây tái sinh vên vên thay đổi theo tuổi, thì cây tái sinh vên vên được phân chia thành 2 nhóm tuổi: nhóm 1 là những cây tái sinh có D1.3 ≤ 10 cm và H ≤ 100 cm, nhóm 2 gồm những cây có H > 100 cm đến D1.3 ≤ 10 cm. Độ ẩm (%) của tầng đất mặt được xác định ở trung tâm ô mẫu 25 m2 bằng máy đo nhanh, còn độ tàn che tán rừng được xác định bằng mục trắc, theo dải 30*10 m. (3) Những số liệu khác cần thu thập bao gồm số liệu về khí hậu - thủy văn, đất và những hoạt động lâm sinh. Cách thức thu thập được thực hiện theo những chỉ dẫn thông thường trong lâm học. Tất cả số liệu thu thập trên ô tiêu chuẩn được ghi vào bảng ngoại nghiệp. 1.2.3. Phương pháp xử lý số liệu (1) Tính toán xác suất bắt gặp loài trong quan hệ với từng biến môi trường Phần dưới đây tính toán xác suất bắt gặp loài theo hai nhóm tuổi khác nhau - đó là nhóm cá thể có H ≤ 100 cm, nhóm cá thể có H > 100 cm đến D1.3 ≤ 10 cm. Trình tự tính toán như sau: - Trước hết, tập hợp độ bắt gặp vên vên và các biến môi trường (X1 = độ ẩm; X2 = độ tàn che tán rừng) ở cả 4 trạng thái rừng. - Tiếp đến, tính quan hệ giữa độ bắt gặp vên vên với mỗi biến môi trường. Ở đây xác suất bắt gặp loài (PX) tương ứng với một biến môi trường nhất định (Xi) được thăm dò bằng hai dạng mô hình hồi quy logit sau đây: + Mô hình sigmoid Ey = P = exp(bo + b1*Xi)/[1 + exp(bo + b1*Xi)] (1) + Mô hình logit Gauss Ey = P = exp(bo + b1*Xi + b2*Xi2)/[1 + exp(bo + b1*Xi + b2*Xi2] (2) Các tham số của mô hình 1 và 2 được ước lượng theo nguyên lý hợp lý tối đa. Để biết đường cong Logit Gauss có phù hợp hơn đường cong sigmoid hay không, thực hiện kiểm định giả thuyết (H0: b2 = 0) bằng thống kê t. Khi mô hình Logit Gauss tồn tại và b2 < 0 một cách có ý nghĩa, thì từ mô hình 2 tính những ước lượng sau đây: - Tối ưu sinh thái: U = b1/2*b2 (3) - Tính chống chịu sinh thái: T = 1/SQRT(-2*b2) (4) - Biên độ sinh thái: U ± T (5) - Xác suất lớn nhất bắt gặp loài: Pmax = exp(Y)/[1 + exp(Y)] (Với Y = exp(bo + b1*U + b2*U2)) (6) - Sau đó từ mô hình phù hợp nhất lập bảng và đồ thị mô tả xác suất bắt gặp loài tương ứng với một biến môi trường nhất định. (2) Xác định ảnh hưởng tổng hợp của các biến môi trường đến xác suất bắt gặp loài trong mỗi trạng thái rừng. Ở đây xác suất bắt gặp loài (PX) tùy theo tập hợp biến môi trường (X1 = độ ẩm; X2 = độ tàn che) được thăm dò theo những dạng mô hình hồi quy logit sau đây: + Dạng mặt phẳng với hai biến dự đoán: Loge[P/(1-P)] = b0 + b1X1 + b2X2 (7) hay P = exp(bo + b1*X1 + b2*X2)/[1 + exp(bo + b1*X1 + b2*X2)] + Dạng mặt Logit Gauss 2 biến số: Loge[P/(1-P)] = b0 + b1X1 + b2X12 + b3X2 + b4X22 (8) hay P = exp(Y)/[1+exp(Y)] (Với Y = exp(bo + b1*X1 + b2*X12 + b3*X2 + b4*X22)) Ở các mô hình 7 và 8), biến X1 = độ ẩm; X2 = độ tàn che. Các tham số của mô hình được ước lượng theo nguyên lý hợp lý tối đa. Để thấy mặt phản hồi ở mô hình 8 có giảm đều theo hướng X1 hay không, thực hiện kiểm định giả thuyết về sự tồn tại của b2 bằng thống kê t (H0: b2 ( 0; H0-: b2 < 0). Tương tự, sử dụng thống kê t để kiểm định sự tồn tại của b4 (H0: b4 ( 0; H0-: b4 < 0) nhằm xác định mặt phản hồi có giảm đều theo hướng X2 hay không. Khi các mô hình 8 tồn tại, thì những thông tin về tối ưu (U) và tính chống chịu sinh thái (T) của loài đối với biến X1, X2 được tính toán bằng việc thêm tương ứng các tham số b1 và b2, b3 và b4 của mô hình 8 vào các phương trình từ 3 đến 6. Kết quả tính toán những thông tin về tối ưu (U) và tính chống chịu sinh thái (T) ở đây có thể so sánh với kết quả tính toán từ mô hình 2. (3) Xác định sự tương tác giữa các biến giải thích trong cùng một trạng thái rừng Hai biến giải thích cho thấy ảnh hưởng tương tác với nhau nếu ảnh hưởng của biến này phụ thuộc vào giá trị của biến khác. Việc kiểm định ảnh hưởng tương tác của hai biến X1 và X2 có thể thực hiện bằng cách mở rộng phương trình (8) với biến thứ ba là tích số X1*X2. Mô hình Gauss của hàm 9 với một tích số có dạng logit như sau: loge[P/(1-P)] = b0 + b1X1 + b2X12 + b3X2 + b4X22 + b5X1X2 (9) hay P = exp(Y)/[1+exp(Y)] (Với Y = exp(bo + b1*X1 + b2*X12 + b3*X2 + b4*X22 + b5*X1*X2)) Nếu b2 + b4 0, thì phương trình 9 là một mặt phẳng với đường viền ellipsoid. Khi bề mặt là phẳng, thì tối ưu (u1, u2) cũng có thể được tính từ các hệ số của hàm 9 như sau: u1 = (b5b3 - 2b1b4)/d (10) u2 = (b5b1 - 2b3b2)/d (11) d = 4b2b4 - b52 (12) Tối ưu đối với X1 tương ứng với X2 nhất định là -(b1 + b5X2)/(2b2). Nếu b5 ( 0 thì tối ưu đối với X1 phụ thuộc vào X2, và hai biến tương tác với nhau. Để thấy rõ sự tương tác giữa X1 với X2, thực hiện kiểm định giả thuyết (H0: b5 = 0) bằng thống kê t. (4) So sánh xác suất bắt gặp loài ở 4 trạng thái rừng trong quan hệ với các biến môi trường Việc so sánh ảnh hưởng của từng biến môi trường ở bốn trạng thái rừng (IIB, IIIA1, IIIA2 và IIIA3) đến độ phong phú của cây tái sinh vên vên được thực hiện bằng cách phát triển mô hình hồi quy logit có dạng: loge(P/(1-P)) = b0 + b1X1 + b2X12 + b3*Z1 + b4*Z2 + b5*Z3 (13) hay P = exp(Y)/[1+exp(Y)] (Với Y = exp(bo + b1*X1 + b2*X12 + b3*Z1 + b4*Z2 + b5*Z3)) Trong đó X1 có thể là độ ẩm đất độ tàn che tán rừng, còn Z1, Z2 và Z3 là ba biến giả biểu thị ảnh hưởng của trạng thái rừng. Các quan sát ở trạng thái rừng IIIA3 được mã hóa tương ứng Z1 = 0, Z2 = 0 và Z3 = 0. Các quan sát ở trạng thái rừng IIB được mã hóa tương ứng Z1 = 1, Z2 = 0 và Z3 = 0. Các quan sát ở trạng thái rừng là IIIA1 được mã hóa tương ứng Z1 = 0, Z2 = 1 và Z3 = 0. Các quan sát ở trạng thái rừng IIIA2 được mã hóa tương ứng Z1 = 0, Z2 = 0 và Z3 = 1. Cách mã hóa như thế cho phép so sánh từng cặp biến phản hồi ở ba trạng thái rừng IIB, IIIA1 và IIIA2 với biến phản hồi ở trạng thái rừng IIIA3. Khi biến đổi mô hình 13, có thể thu được bốn mô hình mô tả xác suất bắt gặp loài cây gỗ quan tâm tùy thuộc vào mỗi biến môi trường trong 4 trạng thái rừng. Bốn mô hình có dạng: + Đối với trạng thái rừng IIIA3 loge(P/(1-P)) = b0 + b1X1 + b2X12 (14) hay P = exp(bo + b1*X1 + b2*X12)/[1 + exp(bo + b1*X1 + b2*X12)] + Đối với trạng thái rừng IIB loge(P/(1-P)) = (b0 + b3) + b1*X1 + b2*X12 (15) Hay P = exp(Y)/[1 + exp(Y)] (Với Y = exp((bo + b3) + b1*X1 + b2*X12) + Đối với trạng thái rừng IIIA1 loge(P/(1-P)) = (b0 + b4) + b1*X1 + b2*X12 (16) Hay P = exp(Y)/[1 + exp(Y)] (Với Y = exp((bo + b4) + b1*X1 + b2*X12) + Đối với trạng thái rừng IIIA2 loge(P/(1-P)) = (b0 + b5) + b1*X1 + b2*X12 (17) Hay P = exp(Y)/[1 + exp(Y)] (Với Y = exp((bo + b5) + b1*X1 + b2*X12) Bốn đường cong này sẽ có bốn xác suất cực đại khác nhau, nhưng tối ưu sinh thái và tính chống chịu là giống nhau. Ba hệ số b3, b4 và b5 chỉ ra sự khác biệt giữa logit lớn nhất của đường cong đối với trang thái rừng IIIA3 so với ba đường cong đối với ba trạng thái rừng IIB, IIIA1 và IIIA2. Sự khác biệt giữa các cực đại của bốn đường cong này được kiểm định bằng cách so sánh các sai tiêu chuẩn của sai lệch giữa hàm 14 với ba hàm 15, 16 và 17. Nếu trạng thái rừng không có ảnh hưởng đến độ bắt gặp loài, thì sự khác biệt là (2 với 3 độ tự do. Từ các hàm 14-17, cũng có thể tính được tối ưu sinh thái và tính chống chịu sinh thái của mỗi loài đối với từng biến môi trường. Kết quả tính toán những thông tin về tối ưu (U) và tính chống chịu sinh thái (T) ở đây cũng có thể so sánh với kết quả tính toán từ các mô hình 2, 8 và 9. Từ những kết quả tính toán, thực hiện tổng hợp thành bảng và đồ thị để thuyết minh và phân tích kết quả. 1.2.4. Công cụ xử lý số liệu Tất cả những cách thức tính toán ở mục 1.2.3 được thực hiện theo chỉ dẫn của Thái Văn Trừng (1998) [10]; Nguyễn Văn Thêm (2004) [9] và các tài liệu tham khảo khác. Công cụ tính toán là phần mềm thống kê Excel, SPSS 10.0 và Statgraphics Plus Version 3.0. 2. KẾT QUẢ NGHIÊN CỨU 2.1. Ảnh hưởng của từng yếu tố môi trường 2.1.1. Ảnh hưởng của độ ẩm đất Kết quả nghiên cứu cho thấy, độ bắt gặp cây tái sinh vên vên ở hai giai đoạn tuổi khác nhau đều phụ thuộc vào độ ẩm tầng đất mặt dưới dạng mô hình Logit Gauss như sau: Đối với cấp tuổi 1 PVeV1 = exp(Y)/(1 + exp(Y)) (2.1) (Với Y = -24,7189 + 0,6996*X1 - 0,0048*X12) Đối với cấp tuổi 2 PVeV2 = exp(Y)/(1 + exp(Y)) (2.2) (Với Y = -21,4455 + 0,6180*X1 - 0,0042*X12) Đối với toàn bộ giai đoạn tái sinh PVeV-TS = exp(Y)/(1 + exp(Y)) (2.3) (Với Y = -22,7333 + 0,6491*X1 - 0,0044*X12) Khai triển mô hình 2.1 – 2.2, có thể xác định được xác suất bắt gặp cây tái sinh vên vên ở những điều kiện độ ẩm đất khác nhau (Bảng 2.1 và 2.2). Bảng 2.1. Xác suất bắt gặp cây tái sinh vên vên trong những điều kiện độ ẩm đất khác nhau Độ ẩm đất (X1, %)  Tái sinh (D1.3 < 10 cm)    VeV-TS  VeV1  VeV2   40  0,0217  0,0121  0,0295   50  0,2178  0,1525  0,2456   60  0,5921  0,5031  0,5997   70  0,7584  0,6860  0,7474   80  0,7378  0,6440  0,7148   90  0,5114  0,3650  0,4768   100  0,1390  0,0655  0,1245   Bảng 2.2. Tối ưu, biên độ và tính chống chịu của cây tái sinh vên vên đối với độ ẩm đất Cấp tuổi  U  T  U±T  U±4T  Pmax   VeV-TS  73,0  10,6  62,4 - 83,6  30,6 - 115,5  0,7242   VeV1  73,0  10,2  61,8 - 82,3  32,2 - 113,9  0,7414   VeV2  73,0  10,9  62,8 - 83,9  29,6 - 116.5  0,7413   Phân tích số liệu của bảng 2.1 và 2.2 nhận thấy rằng: + Cây tái sinh vên vên xuất hiện trong điều kiện môi trường đất có độ ẩm từ 25% trở lên, thường bắt gặp nhất ở nơi có độ ẩm đất từ 55 - 95%. + Ở các giai đoạn tuổi khác nhau, cây tái sinh vên vên đòi hỏi nhu cầu độ ẩm đất thay đổi không lớn; trong đó ở cấp tuổi VeV1 cần độ ẩm từ 55 - 90%, đối với cấp tuổi 2 là 55 - 95%. Xu hướng chung khi điều kiện độ ẩm đất quá khô hoặc quá ẩm đều không thích hợp đối với sự xuất hiện của cây tái sinh vên vên. + Tối ưu độ ẩm tầng đất mặt đối với VeV-TS là 73,0%. Biên độ độ ẩm đất thích hợp là 62,4 - 83,6%. Phạm vi sống sót đối với độ ẩm đất tương ứng là 30,6 - 115,5%. + Mặc dù VeV-TS có thể xuất hiện trong môi trường có độ ẩm đất thay đổi từ khoảng 25% đến bão hoà nước hoàn toàn, song tối ưu và tính chống chịu lại thay đổi không rõ rệt theo tuổi. Ở cấp tuổi 1 và tuổi 2 là 73,0%. Biên độ độ ẩm đất thích hợp cho sự xuất hiện của VeV1 là 61,8 - 82,3%, không khác biệt lớn so với cấp tuổi 2 (62,8 - 83,9%). Tương tự, phạm vi chống chịu của VeV1 với độ ẩm đất từ 32,2 - 113,9%, hẹp hơn không đáng kể so với VeV2 (29,6 – 116,5%). 2.1.2. Ảnh hưởng của độ tàn che Kết quả nghiên cứu chỉ ra rằng, độ tàn che có ảnh hưởng đến xác suất bắt gặp cây tái sinh vên vên ở những cấp tuổi khác nhau. Mối quan hệ giữa xác suất bắt gặp vên vên với độ tàn che tán rừng tồn tại dưới dạng mô hình Logit Gauss như sau: Đối với cấp tuổi 1 PVeV1 = exp(Y)/(1 + exp(Y)) (2.4) (Với Y = -27,058 + 74,2638*X2 - 49,4878*X22) Đối với cấp tuổi 2 PVeV2 = exp(Y)/(1 + exp(Y)) (2.5) (Với Y = -17,794 + 49,6755*X2 - 33,0021*X22) Đối với toàn bộ giai đoạn tái sinh PVeV-TS = exp(Y)/(1 + exp(Y)) (2.6) (Với Y = -21,1918 + 58,6882*X2 - 39,1018*X22) Phân tích mô hình 2.4 – 2.6, có thể xác định được xác suất bắt gặp cây tái sinh vên vên ở những điều kiện độ tàn che tán rừng khác nhau (Bảng 2.3 và 2.4). Bảng 2.3. Xác suất bắt gặp cây tái sinh vên vên trong những điều kiện độ tàn che tán rừng khác nhau Độ tàn che tán rừng  Vên vên rái tái sinh    VeV-TS  VeV1  VeV2   0,4  0,0185  0,0051  0,0390   0,5  0,1648  0,0913  0,2303   0,6  0,4861  0,4218  0,5326   0,7  0,6748  0,6632  0,6916   0,8  0,6756  0,6639  0,6953   0,9  0,4888  0,4242  0,5454   1,0  0,1672  0,0926  0,2459   Từ số liệu của bảng 2.3 và 2.4 nhận thấy rằng: + Cây tái sinh vên vên xuất hiện trong điều kiện môi trường đất có độ tàn che tán rừng từ 0,3 trở lên, thường bắt gặp nhất ở nơi có độ tàn che tán rừng từ 0,6 - 0,9. Bảng 2.4. Tối ưu, biên độ và tính chống chịu của cây tái sinh vên vên đối với độ tàn che Cấp tuổi  (U)  (T)  (U±T)  (U±4T)  (Pmax)   VeV – TS  0,75  0,12  0,63 - 0,87  0,27 - 1,00  0,7082   VeV1  0,75  0,10  0,65 - 0,85  0,35 - 1,00  0,6906   VeV2  0,75  0,12  0,63 - 0,88  0,27 - 1,00  0,7108   + Tối ưu độ tàn che đối với VeV-TS là 0,75. Biên độ độ tàn che là 0,63 - 0,87. Phạm vi sống sót từ 0,27 - 1,00. Mặc dù VeV-TS xuất hiện nới có độ tàn che tán rừng thay đổi từ khoảng 0,3 trở lên. Nhưng VeV-TS ở cấp tuổi khác nhau có biên độ độ tàn che và phạm vi sống sót có sự khác nhau rõ rệt. Thật vậy, biên độ độ tàn che của VeV1 là 0,65 - 0,85, còn VeV2 là 0,63 - 0,88; phạm vi sống sót của VeV1 từ 0,35 - 1,00, hẹp hơn so với VeV2 là 0,27 - 1,00. Tuy nhiên, cả cấp tuổi 1 và cấp tuổi 2 đều tối ưu với độ tàn che là 0,75. Xu hướng chung ở nhưng trạng thái rừng có độ tàn che quá thấp hoặc quá cao đều không thích hợp đối với sự xuất hiện của cây tái sinh vên vên. 2.2. Ảnh hưởng tổng hợp của các yếu tố môi trường Kết quả nghiên cứu chỉ cho thấy: + Xác suất bắt gặp vên vên phụ thuộc vào yếu tố độ ẩm đất và độ tàn che theo mô hình Logit Gauss. Mô hình có dạng: PVeV-TS = exp(Y)/(1+exp(Y)) (2.7) (với Y = -30,897 + 0,4554*X1 + 42,0266*X2 - 0,0032*X12 - 27,8227*X22) + Tương tự, độ ẩm đất và độ tàn che có ảnh hưởng tương tác đến độ bắt gặp vên vên giai đoạn tái sinh. Mô hình biểu thị có dạng Logit Gauss: PVeV-TS = exp(Y)/(1+exp(Y)) (2.8) (với Y = -29,11 + 0,48*X1 + 35,81*X2 - 0,01*X12 + 0,27*X1*X2 - 36,33*X22) 2.3. Ảnh hưởng của yếu tố môi trường trong các trạng thái rừng khác nhau 2.3.1. Ảnh hưởng của độ ẩm đất Kết quả nghiên cứu cũng cho thấy, xác suất bắt gặp VeV-TS không chỉ phụ thuộc vào độ ẩm đất và còn tuỳ thuộc vào trạng thái rừng. Mô hình biểu thị độ phong phú VeV-TS, có dạng Logit Gauss như sau: PVeV-TS = exp(Y)/(1+exp(Y)) (2.9) (Với Y = -22,69 + 0,66*X1 - 0,01*X12 - 0,49*Z1 - 0,23*Z2 + 0,33*Z3) Từ mô hình 2.9, khi thay tương ứng Z1 = 0 và 1, Z2 = 0 và 1, Z3 = 0 và 1, có thể nhận được các mô hình sau đây: - Trạng thái rừng IIIA3 PVeV-TS = exp(Y)/(1+exp(Y)) (2.10) (Với Y = -22,3265 + 0,6439*X1 - 0,0044*X12) - Trạng thái rừng IIB PVeV-TS = exp(Y)/(1+exp(Y)) (2.11) (Với Y = -22,5074 + 0,6439*X1 - 0,0044*X12) - Trạng thái rừng IIIA1 PVeV-TS = exp(Y)/(1+exp(Y)) (2.12) (Với Y = -22,7126 + 0,6439*X1 - 0,0044*X12) - Trạng thái rừng IIIA2 PVeV-TS = exp(Y)/(1+exp(Y)) (2.13) (Với Y = -22,1871 + 0,6439*X1 - 0,0044*X12) Từ mô hình 2.10 – 2.13, có thể tính được xác suất bắt gặp VeV-TS trong những điều kiện độ ẩm khác nhau (Bảng 2.5) Bảng 2.5. Độ phong phú cây tái sinh vên vên tuỳ thuộc vào độ ẩm đất trong 4 trạng thái rừng Độ ẩm đất  Độ phong phú của cây VeV-TS trong 4 trạng thái rừng    IIB  IIIA1  IIIA3  IIIA3   60  0,5712  0,5203  0,6472  0,6148   70  0,7322  0,6901  0,7902  0,7661   80  0,6994  0,6546  0,7622  0,7360   90  0,4511  0,4009  0,5309  0,4961   100  0,1074  0,0893  0,1422  0,1260   Khai triển mô hình 2.10 – 2.13 và hình 2.1 cho thấy, độ ẩm tối ưu đối với VeV-TS là như nhau (71,9%), song xác suất bắt gặp cao nhất ở trong 4 trạng thái rừng lại khác nhau. Xác suất bắt gặp cao nhất tương tứng với 4 trạng thái là 0,7078; 0,5975; 0,6592 và 0,7710. So với trạng thái rừng IIIA3, độ phong phú VeV-TS trong cùng điều kiện độ ẩm (71,9%) ở trạng thái IIB, IIIA1 thấp hơn tương ứng 2,1 và 1,4 lần, còn IIIA2 cao hơn 0,8 lần (Bảng 2.6). Bảng 2.6. So sánh sai khác về độ phong phú VeV-TS tuỳ thuộc vào độ ẩm đất trong các trạng thái rừng Độ ẩm đất (%)  OR    IIIA3 - IIB  IIIA3 - IIIA1  IIIA3 - IIIA2   Cùng độ ẩm (71,9%)  2,1  1,4  0,8   Khác độ ẩm (70 và 50%)  1,6  3,0  2,1   So sánh sự khác biệt về độ phong phú VeV-TS ở 3 trạng thái IIB, IIIA1 và IIIA2 so với IIIA3 khi độ ẩm thay đổi từ 70 - 50% cho thấy, độ phong phú ở IIB, IIIA1 và IIIA2 so với IIIA3 thấp hơn tương ứng 1,6; 3,0 và 2,1 lần (Bảng 2.6) 2.3.2. Ảnh hưởng của độ tàn che Mô hình mô tả xác suất bắt gặp VeV-TS khi phụ thuộc vào độ tàn che tán rừng trong các trạng thái rừng khác nhau có dạng đường cong Logit Gauss (Hình 2.2): PVeV-TS = exp(Y)/(1+exp(Y)) (2.14) (Với Y = -23,81 + 69,21*X2 - 47,91*X22 - 0,96*Z1 - 0,58*Z2 + 0,32*Z3) Khai triển mô hình 2.14 bằng cách mã hoá các trạng thái rừng, các mô hình 2.15 – 2.18 được xây dựng. - Trạng thái rừng IIIA3 PVeV-TS = exp(Y)/(1+exp(Y)) (2.15) (Với Y = -19,9500 + 57,4232*X2 - 39,1154*X22) - Trạng thái rừng IIB PVeV-TS = exp(Y)/(1+exp(Y)) (2.16) (Với Y = -20,4455 + 57,4232*X2 - 39,1154*X22) - Trạng thái rừng IIIA1 PVeV-TS = exp(Y)/(1+exp(Y)) (2.17) (Với Y= -20,5436 + 57,4232*X2 - 39,1154*X22) - Trạng thái rừng IIIA2 PVeV-TS = exp(Y)/(1+exp(Y)) (2.18) (Với Y = -19,8422 + 57,4232*X2 - 39,1154*X22) Xác suất bắt gặp VeV-TS trong nhưng điều kiện độ tàn che tán rừng khác nhau trong các trạng thái rừng (Bảng 2.7) Bảng 2.7. Độ phong phú cây tái sinh vên vên tuỳ thuộc vào độ tàn che trong 4 trạng thái rừng Độ tàn che tán rừng  Độ phong phú của cây VeV-TS trong 4 trạng thái rừng    IIB  IIIA1  IIIA3  IIIA3   0,5  0,1805  0,1665  0,2871  0,2656   0,6  0,4817  0,4573  0,6295  0,6041   0,7  0,6420  0,6192  0,7663  0,7464   0,8  0,6128  0,5893  0,7432  0,7221   0,9  0,3898  0,3667  0,5387  0,5118   1,0  0,1055  0,0966  0,1774  0,1622   Bảng 2.8. So sánh sai khác về độ phong phú VeV-TS tuỳ thuộc vào độ tàn che trong các trạng thái rừng Độ tàn che tán rừng  OR    IIIA3 – IIB  IIIA3 – IIIA1  IIIA3 - IIIA2   Cùng độ tàn che (0,72)  0,5  1,8  0,8   Khác độ tàn che (0,7 và 0,5)  0,13  5,2  12,7   Phân tích mô hình 2.15 - 2.16 và hình 2.2 cho thấy, mặc dù tối ưu độ tàn che của VeV-TS ở trong 4 trạng thái rừng IIB, IIIA1, IIIA2 và IIIA3 là giống nhau (0,72), nhưng xác suất bắt gặp VeV-TS lại khác nhau giữa các trạng thái. Tương ứng là: 0,7657; 0,5549; 0,6465 và 0,8177. Từ mô hình 2.15 - 2.16 cho thấy, xác suất bắt gặp VeV-TS khi độ tàn che tàn rừng như nhau (0,72) giữa các trạng thái thì sai khác độ phong phú của VeV-TS giữa trạng thái IIIA3 thấp hơn 0,5 lần so với trạng thái IIB. Tương tự cao hơn 5,2 và 12,7 lần so với trạng thái rừng IIIA2 và IIIA3 (Bảng 2.8). KẾT LUẬN Từ các kết quả trên, chúng ta có thể rút ra một số kết luận sau: Độ phong phú của tái sinh vên vên ở 2 cấp tuổi đều phụ thuộc vào trạng thái rừng, độ ẩm đất và độ tàn che tán rừng. Độ ẩm đất thích hợp cho cấp tuổi 1 từ 61,8 - 82,3%, tối ưu là 73,0%; ở cấp tuổi 2 tương ứng là 62,8 - 83,9% và 73,0%. Độ tàn che tán rừng thích hợp cho cấp tuổi 1 là 0,65 – 0,85, tối ưu 0,75; ở cấp tuổi 2 từ 0,63 - 0,88, tối ưu là 0,75. Độ phong phú của tái sinh vên vên thay đổi tùy thuộc vào sự thay đổi của 2 yếu tố độ ẩm đất và độ tàn che tán rừng. Mối liên hệ này có thể mô tả bằng những mô hình Logit Gauss 2 biến số. Khi trạng thái rừng thay đổi, thì độ ẩm đất tầng đất mặt và độ tàn che tán rừng có ảnh hưởng khác nhau đến độ phong phú của vên vên. Những trạng thái rừng ổn định cao (IIIA2 và IIIA3) đảm bảo cho cây vên vên tái sinh tốt hơn so với những trạng thái rừng kém ổn định (IIB và IIIA1). TÀI LIỆU THAM KHẢO 1. George N. Baur (1976, 1979). Cơ sở sinh thái học của kinh doanh rừng mưa. Vương Tấn Nhị dịch, Nxb Khoa Học và Kỹ Thuật, Hà Nội. 2. Vũ Xuân Đề (1985). Nghiên cứu các biện pháp tổng hợp khai thác đảm bảo tái sinh rừng, cải tạo rừng và trồng rừng cây gỗ lớn gỗ quý ở miền Đông Nam Bộ. Phân Viện Lâm Nghiệp Phía Nam. 3. Bùi Đoàn, Vũ Duy Thông, (2003). Nghiên cứu gây trồng cây Vên vên (Anisoptera cochinchinensis) làm nguyên liệu gỗ dán lạng, Trung tâm NC Sinh thái và Môi trường Rừng, Viện Khoa học Lâm nghiệp Việt Nam. 4. Nguyễn Minh Đường (1985). Nghiên cứu gây trồng dầu, sao, vên vên trên các dạng đất đai trống trọc còn khả năng sản xuất gỗ lớn gỗ quý. Báo cáo khoa học 01.9.3, Phân Viện Lâm Nghiệp Phía Nam. 5. Lê Văn Mính (1985). Đặc tính sinh thái của sao, dầu, vên vên ở Đông Nam Bộ. Báo cáo khoa học 01.02.3, Phân Viện Lâm Nghiệp Phía Nam. 6. Richards, PW (1965). Rừng mưa nhiệt đới. Vương Tấn Nhị dịch, Nxb Khoa Học & Kỹ Thuật, Hà Nội, 1965. 7. Lâm Xuân Sanh (1985). Vai trò của các loài cây họ sao - dầu trong sinh thái phát sinh của các hệ sinh thái rừng ở miền Nam Việt Nam. Phân Viện Lâm Nghiệp Phía Nam. 8. Nguyễn Văn Thêm (1992). Nghiên cứu tái sinh tự nhiên của dầu song nàng (Dipterocarpus dyeri ) trong kiểu rừng kín ẩm thường xanh và nửa rụng lá ẩm nhiệt đới ở Đồng Nai, Luận án Phó tiến sĩ Khoa học Nông Nghiệp, Viện Khoa Học Lâm Nghiệp Việt Nam. 9. Nguyễn Văn Thêm (2004). Hướng dẫn sử dụng Statgraphics Plus Version 3.0 & 5.1 để xử lý thông tin trong lâm học. Nxb Nông Nghiệp, Chi nhánh Tp Hồ Chí Minh. 10. Thái Văn Trừng (1998). Những hệ sinh thái rừng nhiệt đới ở Việt Nam. Nxb Khoa Học và Kỹ Thuật, Hà Nội. 11. Kimmins, J. P., 1998. Forest ecology. Prentice – Hall, Upper Saddle River, New Jersey. 12. H. Lamprecht (1989). Silviculture in Troppics. Eschborn. (Summery) Reseach on ecological characteristics of regeneration Anisoptera cochinchinensis Pierre in Forest types closed evergreen and Semi-evergreen in Dong Nai province The South East in general and Dong Nai province in particular, Dipterocarpaceae is an important in forest ecosystems. Especially those species with economic value, ecological value, and conservation value. Anisoptera cochinchinensis is speceis that has these characteristics. Really, Anisoptera cochinchinensis is species being studied to choose to afforest to restore forest. Up to now, to have some researched the family Dipterocarpaceae or Anisoptera cochinchinensis, but not the clarify rule life, rule regeneration and ecological characteristics. So reseach on ecological characteristics of regeneration Anisoptera cochinchinensis in forest types closed evergreen and semi-evergreen in Dong Nai province is significant. The purpose of the paper is to provide the database to clarify the ecological characteristics of regeneration Anisoptera cochinchinensis in forest types closed evergreen and semi-evergreen tropical moist in Dong Nai and the scientific basis for the application of measures regeneration (natural and artificial) and nurture Anisoptera cochinchinensis. Results showed: - Abundance of regeneration Anisoptera cochinchinensis in two levels age depends on the forest type, soil humidity and forest cover. Soil humidity appropriate for level 1 age from 61.8 to 82.3%. Optimize is 73.0%; at level 2 age, respectively from 62.8 to 83.9% and 73.0%. Forest cover appropriate for level 1 age from 0.65 to 0.85. Optimize is 0.75; at age 2 from 0.63 to 0.88. Optimize is 0.75. - Abundance of regeneration Anisoptera cochinchinensis changes depending on changes in two factors: Soil humidity and forest cover. This relationship can be described by Logit Gauss model is 2 variables. - When forest type changes, the soil humidity and forest cover effect to degrees abundance different of Anisoptera cochinchinensis. The forest types are stable (IIIA2 and IIIA3), ensure for Anisoptera cochinchinensis regenerate best than the forest types unstable (IIB and IIIA1). Key words: Regeneration Anisoptera cochinchinensis, abundance, probability appear, soil humidity, forest cover, Logit Gauss model