ẢNH HƯỞNG của KHÍ hậu đến TĂNG TRƯỞNG của DU SAM, BẠCH TÙNG và ĐỈNH TÙNG ở KHU vực đà lạt và đức TRỌNG TỈNH lâm ĐỒNG TT

1 MỞ ĐẦU Đặt vấn đề Du sam (Keteleeria evelyniana Masters), Bạch tùng (Dacrycarpus imbricatus (Blume) de Laub) Đỉnh tùng (Cephalotaxus mannii Hook f.) loài gỗ lớn, quý, có giá trị cao kinh tế, sinh thái môi trường Tại Lâm Đồng, ba loài gỗ phân bố tự nhiên kiểu rừng hỗn giao rộng kim, ẩm, nhiệt đới núi thấp (Rkh) thuộc huyện Lạc Dương, Đơn Dương, Di Linh, Đức Trọng Bidoup Núi Bà Hiện ba loài gỗ xếp vào nhóm IIA theo Nghị định số 32/2006/NĐ-CP Thủ tướng phủ cấp V Sách đỏ Việt Nam (Nguyễn Đức Tố Lưu Philip Ian Thomas, 2004) Xác định đặc tính sinh thái học kỹ thuật trồng rừng từ loài gỗ lớn có giá trị cao khoa học kinh tế nhiệm vụ quan trọng lâm học Mặc khác, quản lý bảo vệ rừng cần có kiến thức sinh trưởng lồi gỗ rừng Những mơ hình dự đốn tăng trưởng lồi gỗ xây dựng dựa mối quan hệ chúng với yếu tố khí hậu Tuy vậy, để đạt mục đích này, lâm học cần phải hiểu rõ mối quan hệ tăng trưởng loài gỗ với yếu tố khí hậu Tăng trưởng bề rộng vòng năm gỗ dấu hiệu rõ nét ảnh hưởng tuổi, khí hậu yếu tố mơi trường khác (địa hình, đất, sinh vật, người…) Vì thế, nhiều nhà nghiên cứu giới (Fritts, 1965, 1976; Bitvinskas, 1974; Cook, 1987) sử dụng phương pháp niên đại thực vật (Dendrochronology) khí hậu thực vật (Dendroclimatology) để xác định niên đại lớp vòng năm ảnh hưởng yếu tố khí hậu đến sinh trưởng loài gỗ Ở Việt Nam, hai phương pháp áp dụng để xác định ảnh hưởng khí hậu đến sinh trưởng Thông ba (Pinus keysiya ex Gordon) (Phạm Trọng Nhân, 2003; Phạm Trọng Nhân ctv, 2011), Du sam (Nguyễn Văn Nhẫn, 2012) Bạch tùng (Nguyễn Văn Thêm, 2012) Tuy vậy, nghiên cứu chưa làm sáng tỏ vai trò yếu tố khí hậu yếu tố mơi trường khác sinh trưởng Du sam, Bạch tùng Đỉnh tùng Những nghiên cứu chưa xây dựng mơ hình dự đốn tăng trưởng ba loài gỗ dựa biến động yếu tố khí hậu Vì thế, lâm học chưa hiểu rõ mối quan hệ sinh trưởng Du sam, Bạch tùng Đỉnh tùng với yếu tố khí hậu Xuất phát từ vấn đề đặt đây, nghiên cứu sử dụng phương pháp niên đại thực vật khí hậu thực vật để xác định thời gian hình thành vòng năm phân tích phản ứng Du sam, Bạch tùng Đỉnh tùng với biến động yếu tố khí hậu khu vực Đức Trọng Đà Lạt thuộc tỉnh Lâm Đồng Từ mối liên hệ này, xây dựng mơ hình dự đốn sinh trưởng ba loài gỗ dựa theo yếu tố khí hậu Mục tiêu nghiên cứu 2.1 Mục tiêu tổng quát Xác định mối quan hệ tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng với biến động yếu tố khí hậu 2.2 Mục tiêu cụ thể (1) Phân tích đặc trưng khí hậu để làm sở cho việc phân cấp chế độ khô ẩm khu vực nghiên cứu (2) Xác định tháng mà điều kiện khí hậu có ảnh hưởng rõ rệt đến tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng (3) Phân tích ảnh hưởng khí hậu đến tăng trưởng Du sam, Bạch tùng Đỉnh tùng điều kiện môi trường khác (4) Xây dựng mơ hình dự đốn tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng theo yếu tố khí hậu Phạm vi nghiên cứu Phạm vi nghiên cứu đề tài biến động bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng mối quan hệ với biến động yếu tố khí hậu Địa điểm nghiên cứu đặt khu vực Đức Trọng Đà Lạt thuộc tỉnh Lâm Đồng Nghiên cứu phân tích phản ứng tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng biến động yếu tố khí hậu: nhiệt độ khơng khí trung bình tháng (T, 0C); tổng lượng mưa tháng (M, mm); độ ẩm khơng khí trung bình tháng (Rh%); tổng số nắng tháng (N, giờ); tổng lượng nước bốc tháng (P, mm) hệ số thủy nhiệt tháng (K) Ý nghĩa đề tài Về lý luận, nghiên cứu cung cấp khoa học để xác định đặc tính sinh thái Du sam, Bạch tùng Đỉnh tùng Về thực tiễn, nghiên cứu cung cấp hàm dự đoán tăng trưởng số khoa học để xây dựng biện pháp lâm sinh bảo tồn Du sam, Bạch tùng Đỉnh tùng Những đóng góp luận án Luận án xây dựng ba chuỗi số bề rộng vòng năm chuẩn hóa Du sam, Bạch tùng Đỉnh tùng khu vực Đức Trọng Đà Lạt thuộc tỉnh Lâm Đồng Luận án ba yếu tố khí hậu kiểm sốt chặt chẽ tăng trưởng bề rộng vòng năm Du sam nhiệt độ khơng khí vào tháng 4, lượng mưa vào tháng lượng nước bốc vào tháng 10 Sự nâng cao nhiệt độ khơng khí vào tháng có ảnh hưởng xấu tăng trưởng bề rộng vòng năm Du sam Trái lại, nâng cao lượng mưa vào tháng lượng nước bốc vào tháng 10 điều kiện tốt tăng trưởng bề rộng vòng năm Du sam Ba yếu tố khí hậu kiểm sốt chặt chẽ tăng trưởng bề rộng vòng năm Bạch tùng lượng mưa vào tháng 11, số nắng vào tháng tháng Sự nâng cao ba yếu tố dẫn đến ảnh hưởng xấu tăng trưởng bề rộng vòng năm Bạch tùng Ba yếu tố khí hậu kiểm sốt chặt chẽ tăng trưởng bề rộng vòng năm Đỉnh tùng lượng mưa vào tháng 1, nhiệt độ khơng khí trung bình vào tháng – số nắng vào tháng 11 Sự nâng cao lượng mưa vào tháng nhiệt độ khơng khí trung bình vào tháng – điều kiện xấu tăng trưởng bề rộng vòng năm Đỉnh tùng Trái lại, nắng nhiều vào tháng 11 điều kiện tốt tăng trưởng bề rộng vòng năm Đỉnh tùng Luận án độ cao địa hình có ảnh hưởng đến phản ứng Du sam, Bạch tùng Đỉnh tùng với khí hậu Khi sống độ cao 1.400 - 1.600 m so với mặt nước biển, ba loài gỗ phản ứng với thay đổi yếu tố khí hậu rõ rệt so với độ cao 1.400 m Luận án xây dựng ba mô hình ước lượng tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng dựa biến động yếu tố khí hậu Bố cục luận án bao gồm phần mở đầu, chương phần kết luận Chương 1: Tổng quan Chương 2: Nội dung phương pháp nghiên cứu Chương 3: Kết nghiên cứu thảo luận Kết luận đề nghị Luận án bao gồm 175 trang; 62 bảng; 18 hình đồ thị; 38 phụ lục Luận án tham khảo 62 tài liệu nước nước Chƣơng TỔNG QUAN Đề tài tổng quan phương pháp phân tích khí hậu thực vật Tổng quan tóm tắt từ 59 tài liệu tham khảo Dưới thảo luận chung Khoa học niên đại thực vật khí hậu thực vật phát triển mạnh mẽ từ nửa đầu kỷ XX Hai phương pháp sử dụng để phân tích mối quan hệ tăng trưởng bề rộng vòng năm gỗ với khí hậu yếu tố mơi trường khác (địa hình, đất, QXTV, động đất, lũ lụt, lửa rừng, sâu bệnh ) Thông qua mối liên hệ này, nhà sinh thái học không xác định đặc tính sinh thái học lồi gỗ, mà mức độ ảnh hưởng tốt hay xấu yếu tố khí hậu Vì thế, đề tài áp dụng phương pháp niên đại thực vật khí hậu thực vật để phân tích vai trò yếu tố khí hậu tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng khu vực Đức trọng Đà Lạt tuộc tỉnh Lâm Đồng Cho đến khoa học thực tiễn Việt Nam thiếu thơng tin ảnh hưởng địa hình, điều kiện QXTV đất trống đến mối quan hệ sinh trưởng Du sam, Bạch tùng Đỉnh tùng với khí hậu Mặt khác, chưa có mơ hình dự đốn tăng trưởng ba loài gỗ dựa yếu tố khí hậu 5 Vì thế, đề tài phát triển chuỗi số bề rộng vòng năm chuẩn hóa sử dụng chúng để xác định phản ứng Du sam, Bạch tùng Đỉnh tùng với biến động yếu tố khí hậu Đây sở khoa học để xác định đặc tính sinh thái học dự đoán sinh trưởng Du sam, Bạch tùng Đỉnh tùng Khí hậu cấu thành nhiều yếu tố khác Đời sống gỗ rừng phụ thuộc lớn vào hai trình nhiệt ẩm Vì thế, đề tài phân tích phản ứng Du sam, Bạch tùng Đỉnh tùng với yếu tố khí hậu: T, M, Rh, N, P K Ở khu vực phía Nam nước ta, số liệu khí hậu đo đạc từ đầu thập niên 1980 trở lại Vì thế, đề tài phân tích phản ứng tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng yếu tố khí hậu khu vực nghiên cứu từ năm 1980 đến năm 2014 Chƣơng ĐỐI TƢỢNG, NỘI DUNG, PHƢƠNG PHÁP NGHIÊN CỨU 2.1 Đối tƣợng địa điểm nghiên cứu Đối tượng nghiên cứu quần thể Du sam, Bạch tùng Đỉnh tùng Ba loài gỗ mọc tự nhiên kiểu Rkh Nghiên cứu thực khu vực Đức Trọng Đà Lạt thuộc tỉnh Lâm Đồng Tọa độ địa lý: 108022’13’’ - 108030’24’’ kinh độ Đông; 11046’15’’ 11055’27’’ vĩ độ Bắc Độ cao địa hình từ 1.000 – 1.700 m so với mặt nước biển; độ dốc 250 Đất vàng đỏ phát triển đá granite Khu vực nghiên cứu có khí hậu ơn hòa núi cao Nhiệt độ khơng khí trung bình hàng năm 180C, cao 18,90C thấp 16,90C Lượng mưa trung bình năm 1.823 mm, cao 2.357 mm, thấp 1.354 mm Độ ẩm khơng khí trung bình 84%, cao 88%, thấp 80% 6 2.2 Nội dung nghiên cứu (1) Đặc điểm khí hậu khu vực nghiên cứu (2) Đặc điểm bề rộng vòng năm số bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng (3) Phản ứng Du sam, Bạch tùng Đỉnh tùng biến động yếu tố khí hậu (4) Xác định vai trò yếu tố khí hậu tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng (5) Ảnh hưởng độ cao địa hình đến mối quan hệ tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng với yếu tố khí hậu (6) Phản ứng Du sam quần xã thực vật đất trống với biến động yếu tố khí hậu (7) Dự đoán tăng trưởng Du sam, Bạch tùng Đỉnh tùng dựa theo yếu tố khí hậu 2.3 Phƣơng pháp nghiên cứu 2.3.1 Phƣơng pháp luận Phương pháp luận đề tài phương pháp mô tả phương pháp mơ hình hóa Phương pháp mơ tả sử dụng để mô tả đặc trưng quần thụ, bề rộng vòng năm số bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng Phương pháp mơ hình hóa mơ hình tốn áp dụng để phân tích mối quan hệ biến động bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng với biến động yếu tố khí hậu Các mơ hình xây dựng dựa giả định bề rộng vòng năm (Zr, mm) gỗ hàm yếu tố khí hậu yếu tố mơi trường khác, nghĩa Zr = f(khí hậu; địa hình, đất, lửa rừng, sâu bệnh, QXTV) Từ quan điểm đây, hướng tiếp cận đề tài việc phát triển ba chuỗi số bề rộng vòng năm chuẩn hóa Du sam, Bạch tùng Đỉnh tùng Tiếp đến xây dựng ma trận tương quan để phân tích khuynh hướng cường độ mối quan hệ tăng trưởng bề rộng vòng năm ba lồi gỗ với yếu tố khí hậu Tiếp theo sử dụng phương pháp phân tích hàm phản hồi tuyến tính đa biến bước để xác định biến khí hậu đóng vai trò lớn tăng trưởng bề rộng vòng năm ba lồi gỗ Sau xây dựng hàm dự đốn tăng trưởng bề rộng vòng năm ba lồi gỗ dựa theo yếu tố khí hậu có ảnh hưởng lớn Nghiên cứu dựa ba giả thuyết (1) Tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng thay đổi theo tuổi Giả thuyết kiểm định cách so sánh bề rộng vòng năm ba lồi gỗ mẫu có kích thước khác (2) Những yếu tố khí hậu tháng năm ảnh hưởng không giống đến tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng Giả thuyết kiểm định phân tích ma trận tương quan tăng trưởng bề rộng vòng năm ba lồi gỗ với biến động yếu tố khí hậu tháng khác Vai trò yếu tố khí hậu tăng trưởng bề rộng vòng năm ba lồi gỗ phân tích hàm phản hồi tuyến tính đa biến (3) Khí hậu ảnh hưởng đến tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng lớn so với yếu tố môi trường khác (địa hình điều kiện QXTV) Giả thuyết kiểm định phân tích so sánh phản ứng tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng với biến động yếu tố khí hậu cấp địa hình điều kiện QXTV khác 2.3.2 Phƣơng pháp thu thập số liệu Số liệu vòng năm Du sam, Bạch tùng Đỉnh tùng thu thập quần thụ thuộc kiểu Rkh Chuỗi bề rộng vòng năm Du sam thu thập quần thụ có mật độ 665 cây/ha, tiết diện ngang 51,7 m2/ha, trữ lượng gỗ 481,8 m3 Chuỗi bề rộng vòng năm Bạch tùng thu thập quần thụ có mật độ 820 cây/ha, tiết diện ngang 36,8 m2/ha, trữ lượng gỗ 296,9 m3 Chuỗi bề rộng vòng năm Đỉnh tùng thu thập quần thụ có mật độ 790 cây/ha, tiết diện ngang 37,5 m2/ha, trữ lượng gỗ 260,5 m3 Những quần thụ phân bố ba cấp độ cao 1.000 - 1.200 m, 1.200 – 1.400 m 1.400 – 1.600 m so với mặt nước biển Vị trí mẫu nằm cách Trạm khí tượng thủy văn Đà Lạt 20 km Để loại bỏ yếu tố đất, mẫu thu thập loại đất vàng đỏ phát triển đá granite Những mẫu Du sam Bạch tùng chọn để thu thập vòng năm có D > 80 cm, Đỉnh tùng có D > 50 cm Thời gian Du sam sống đất trống khoảng 34 năm từ 1980 – 2014 Vì thế, đề tài phân tích ảnh hưởng điều kiện QXTV đất trống đến mối quan hệ tăng trưởng bề rộng vòng năm Du sam với khí hậu Chuỗi bề rộng vòng năm thân thu thập vị trí cách mặt đất 120 – 130 cm khoan tăng trưởng Pressler Số mẫu thu thập Du sam, Bạch tùng Đỉnh tùng Ngoài ra, chuỗi bề rộng vòng năm Du sam mọc đất trống thu thập mẫu có kích thước D = 200 cm H = 28 m Tổng số mẫu ba loài 25 Những yếu tố khí hậu (T, M, Rh, N, P, K) 12 tháng năm thu thập từ 1980 đến 2014 Trạm khí tượng thủy văn Đà Lạt 2.3.3 Phƣơng pháp xử lý số liệu Các chuỗi bề rộng vòng năm mẫu chuyển thành chuỗi số bề rộng vòng năm (Kd) Sau chuỗi số bề rộng vòng năm mẫu tính trung bình để nhận chuỗi số bề rộng vòng năm chuẩn hóa (Kd) Chuỗi số Kd chuẩn hóa Du sam xây dựng từ mẫu với tổng số 192 năm (1822 - 2013) Chuỗi số Kd chuẩn hóa Bạch tùng xây dựng từ mẫu với tổng số 201 năm (1813 - 2013) Chuỗi số Kd chuẩn hóa Đỉnh tùng xây dựng từ mẫu với tổng số 127 năm (1887 - 2013) Sau xác định tín hiệu khí hậu (mSx), tín hiệu quần thể (EPS) tỷ lệ tín hiệu khí hậu với tín hiệu nhiễu loạn (SNR) chuỗi số bề rộng vòng năm Mối quan hệ tăng trưởng bề rộng vòng năm ba lồi gỗ với yếu tố khí hậu tháng năm phân tích ma trận tương quan Vai trò yếu tố khí hậu tăng trưởng bề rộng vòng năm ba lồi gỗ phân tích hàm phản hồi tuyến tính đa biến bước nhằm loại bỏ tượng cộng tuyến tính biến khí hậu Các hàm phản hồi xây dựng yếu tố khí hậu có ảnh hưởng rõ rệt (P < 0,05) đến tăng trưởng bề rộng vòng năm ba lồi gỗ Vai trò yếu tố khí hậu đánh giá thơng qua độ lớn hệ số hồi quy chuẩn hóa (HSHQ) Những yếu tố khí hậu mà HSHQ theo thứ tự giảm dần cho biết mức độ giảm dần vai trò chúng tăng trưởng bề rộng vòng năm ba lồi gỗ Sau yếu tố khí hậu đóng vai trò lớn sử dụng để xây dựng hàm dự đoán tăng trưởng bề rộng vòng năm ba lồi gỗ Ảnh hưởng độ cao địa hình, QXTV đất trống đến phản ứng tăng trưởng bề rộng vòng năm ba loài gỗ với yếu tố khí hậu phân tích so sánh ma trận tương quan Nếu ba loài gỗ mọc điều kiện môi trường khác (độ cao địa hình, QXTV đất trống) phản ứng giống với thay đổi khí hậu, khí hậu yếu tố kiểm sốt mạnh tăng trưởng bề rộng vòng năm ba lồi gỗ Trái lại, ba loài gỗ mọc độ cao khác nhau, QXTV đất trống phản ứng không giống với thay đổi khí hậu, ba yếu tố chi phối đến mối quan hệ ba loài gỗ với biến động yếu tố khí hậu Cơng cụ tính tốn bảng tính Excel, phần mềm thống kê Statgraphics Plus version 4.0 SPSS 10.0 Bảng tính Excel sử dụng để tập hợp số liệu vẽ biểu đồ đồ thị Hai phần mềm thống kê Statgraphics Plus version 4.0 SPSS 10.0 sử dụng để tính tốn thống kê mơ tả phân tích mối quan hệ tăng trưởng bề rộng vòng năm với khí hậu Chƣơng KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 3.1 Đặc điểm khí hậu khu vực Đà Lạt Khu vực Đà Lạt có khí hậu ơn hòa Tổng nhiệt độ năm khoảng 6.561 C; bình quân tháng năm 18,0C Lượng mưa bình quân năm 1.800 mm Độ ẩm khơng khí bình qn năm 84% Hệ số thủy nhiệt bình qn năm 2,68 Những năm có nhiệt độ cao lượng mưa thấp Theo hệ thống phân loại chế độ khô ẩm Việt Nam Thái văn Trừng (1999), khí hậu Đà Lạt thuộc cấp II (hơi ẩm): tháng khô (tháng 12 năm trước đến tháng năm sau), tháng hạn (tháng 2) khơng có tháng kiệt Nói chung, lượng mưa năm khu vực Đà Lạt đủ cung 10 cấp nước cho thực vật Tuy vậy, tình trạng thiếu hụt nước thực vật xảy tháng mùa khô 3.2 Đặc điểm bề rộng vòng năm số bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng 3.2.1 Đặc điểm bề rộng vòng năm Đặc điểm bề rộng vòng năm (Zr, mm) Du sam, Bạch tùng Đỉnh tùng ghi lại Bảng Nói chung, bề rộng vòng năm thay đổi tùy theo lồi gỗ tuổi mẫu Vì thế, phân tích ảnh hưởng khí hậu đến tăng trưởng bề rộng vòng năm ba lồi gỗ này, chuỗi bề rộng vòng năm cần phải chuẩn hóa để loại bỏ yếu tố tuổi Bảng Đặc trưng thống kê bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng Loài (1) Du sam Bạch tùng Đỉnh tùng N (năm) (2) 222 205 176 Zr (mm) (3) 2,12 1,87 1,56 ±S ±Se CV% ZrMin ZrMax (4) (5) 1,11 0,09 0,82 0,06 0,53 0,04 (6) 53,4 41,4 34,2 (7) 0,45 0,62 0,51 (8) 6,67 5,45 3,82 3.2.2 Chỉ số bề rộng vòng năm Du sam, Bạch tùng, Đỉnh tùng Phân tích số bề rộng vòng năm trung bình (Kd) Du sam, Bạch tùng Đỉnh tùng cho thấy, chuỗi số Kd Du sam, Bạch tùng Đỉnh tùng tương ứng 192 năm, 201 năm 127 năm Chỉ số Kd trung bình ba loài gỗ 1,0 Biên độ dao động số Kd loài Bạch tùng lớn (Kd = 0,65 – 1,41), thấp loài Du sam (Kd = 0,78 – 1,21) Hệ số biến động số Kd cao loài Bạch tùng (11,0%), thấp loài Du sam (9,0%) Cả ba lồi gỗ có tính nhạy cảm cao biến động yếu tố khí hậu; Bạch tùng (mSx = 0,151) Đỉnh 11 tùng (mSx = 0,150) có tính nhạy cảm cao so với Du sam (mSx = 0,130) Biến động số bề rộng vòng năm cá thể loài tương đồng với (rm = 0,689 với P < 0,001 Du sam; rm = 0,599 với P < 0,001 Bạch tùng; rm = 0,635 với P < 0,001 Đỉnh tùng) Điều chứng tỏ cá thể loài phản ứng tương tự biến động môi trường Tín hiệu quần thể nhận giá trị cao (EPS = 0,801 Bạch tùng; 0,814 Du sam 0,839 Đỉnh tùng) Điều chứng tỏ gỗ quần thể loài phản ứng tương đồng biến động môi trường Những tín hiệu khí hậu (SNR) biểu rõ rệt lớp vòng năm; tín hiệu khí hậu Đỉnh tùng (5,2) biểu rõ rệt so với Du sam (4,4) Bạch tùng (4,0) 3.3 Phản ứng Du sam, Bạch tùng Đỉnh tùng biến động yếu tố khí hậu 3.3.1 Phản ứng Du sam yếu tố khí hậu Những phân tích thống kê cho thấy, số Kd biểu mối quan hệ rõ rệt với T1 (r = -0,381; P = 0,029), T2 (r = -0,364; P = 0,037), T4 (r = -0,601; P = 0,000), T6 (r = -0,495; P = 0,003), T7 (r = -0,601; P < 0,001), T10 (r = -0,459; P = 0,007), T1-4 (r = -0,506; P = 0,003) T5-10 (r = -0,609; P < 0,001) Chỉ số Kd biểu mối quan hệ rõ rệt với M (r = 0,438; P = 0,011), M4 (r = 0,345; P = 0,049), M5 (r = 0,410; P = 0,018), M10 (r = -0,351; P = 0,045), M12 (r = -0,383; P = 0,028), M1-4 (r = 0,433; P = 0,012) M11-12 (r = -0,378; P = 0,030) Chỉ số Kd tồn mối quan hệ rõ rệt với Rh4 (r = 0,514; P = 0,002) Mối quan hệ số Kd với số nắng biểu tương đối rõ rệt vào tháng (r = 0,330; P = 0,061) tháng 11 (r = -0,308; P = 0,081) Chỉ số Kd biểu mối quan hệ rõ rệt với P4 (r = -0,540; P = 0,001), P10 (r = 0,637; P < 0,001), P12 (r = 0,409; P = 0,018) P11-12 (r = 0,402; P = 0,020) Chỉ số Kd biểu mối quan hệ rõ rệt với K4 (r = 0,375; P = 0,032), K1-4 (r = 0,338; P < 0,035) K11-3 (r = 0,374; P = 0,032) 12 3.3.2 Phản ứng Bạch tùng yếu tố khí hậu Kết nghiên cứu cho thấy, số Kd với biểu mối quan hệ rõ rệt với T6 (r = -0,468; P = 0,000), T7 (r = -0,357; P = 0,042), T10 (r = -0,412; P = 0,017) T5-10 (r = -0,404; P = 0,020) Chỉ số Kd biểu mối quan hệ rõ rệt với M3 (r = 0,364; P = 0,037), M6 (r = 0,327; P = 0,063), M11 (r = -0,685; P = 0,000), M1-4 (r = 0,381; P = 0,029), M11-12 (r = -0,617; P = 0,001) M11-3 (r = -0,395; P = 0,023) Chỉ số Kd tồn mối quan hệ rõ rệt với Rh5 (r = -0,371; P = 0,034), Rh7 (r = 0,361; P = 0,039), Rh11 (r = -0,353; P = 0,044), Rh11-12 (r = -0,422; P = 0,014) Rh11-3 (r = -0,373; P = 0,032) Chỉ số Kd biểu mối quan hệ rõ rệt với N1 (r = -0,562; P = 0,001), N4 (r = -0,541; P = 0,001), N7 (r = -0,325; P = 0,065); N9 (r = 0,331; P = 0,060); N11 (r = 0,441; P = 0,010) N1-4 (r = -0,496; P = 0,003) Chỉ số Kd biểu mối quan hệ rõ rệt với P3 (r = -0,387; P = 0,026), P4 (r = -0,326; P = 0,064), P9 (r = 0,489; P = 0,004), 10 (r = 0,332; P = 0,059), P11 (r = 0,513; P = 0,002), P1-4 (r = 0,379; P = 0,030) P11-12 (r = 0,416; P = 0,016) Chỉ số Kd có mối quan hệ rõ rệt với K3 (r = 0,362; P = 0,029), K6 (r = 0,349; P = 0,047), K11 (r = -0,684; P = 0,000), K1-4 (r = 0,382; P = 0,028), K11-12 (r = 0,618; P = 0,000), K11-3 (r = -0,391; P = 0,025) 3.3.3 Phản ứng Đỉnh tùng yếu tố khí hậu Những phân tích thống kê cho thấy, số Kd biểu mối quan hệ rõ rệt với T1 (r = -0,600; P < 0,001), T2 (r = -0,410; P = 0,018), T3 (r = -0,440; P = 0,010), T4 (r = -0,457; P = 0,008), T5 (r = -0,503; P = 0,003), T1-4 (r = -0,631; P < 0,001) T11-3 (r = -0,423; P < 0,014) Chỉ số Kd biểu mối quan hệ rõ rệt với M1 (r = -0,436; P = 0,011), M3 (r = 0,364; P = 0,037) M8 (r = 0,446; P = 0,009) Chỉ số Kd biểu mối quan hệ rõ rệt với Rh7 (r = 0,315; P = 0,074) Rh11-12 (r = 0,376; P = 0,031) Chỉ số Kd biểu mối quan hệ rõ rệt với N1 (r = -0,372; P = 0,033), N3 (r = =0,517; P = 0,003), N6 (r = -0,413; P = 0,017), N11 (r = 0,546; P = 0,001), N1-4 (r = -0,343; P = 0,013), N5-10 (r = -0,350; P = 0,046) N11-12 (r = 0,360; P = 0,040) Sự nâng cao lượng nước bốc vào tháng có ảnh hưởng tốt đến tăng trưởng bề rộng vòng năm Đỉnh tùng Trái lại, nâng cao lượng nước bốc vào tháng – 12 lại dẫn đến ảnh hưởng xấu tăng trưởng 13 bề rộng vòng năm Đỉnh tùng Sự gia tăng hệ số thủy nhiệt vào tháng 2, 4, 5, 6, 10, 12, - - 10 có ảnh hưởng tốt đến tăng trưởng bề rộng vòng năm Đỉnh tùng Trái lại, gia tăng hệ số thủy nhiệt vào tháng 1, 3, 7, 8, 9, 11, 11 - 12 tháng 11 năm trước đến tháng năm sau điều kiện xấu tăng trưởng bề rộng vòng năm Đỉnh tùng 3.4 Xác định vai trò yếu tố khí hậu tăng trƣởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng 3.4.1 Vai trò yếu tố khí hậu tăng trƣởng bề rộng vòng năm Du sam Kết phân tích hồi quy tương quan đa tuyến tính bước cho thấy, số Kd với yếu tố T1, T2, T4, T6, T7, T10, T1-4 T5-10 tồn mối quan hệ chặt chẽ (R2 = 52,0%) Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Du sam phụ thuộc nhiều vào yếu tố T4 (HSHQ = |-0,2575|); T10 (HSHQ = |-0,2395|) T7 (HSHQ = |-0,2282|); thấp T6 (HSHQ = |-0,0722| Mặt khác, tăng trưởng bề rộng vòng năm Du sam phụ thuộc vào T5-10 (HSHQ = |-0,4333|) lớn so với T1-4 (HSHQ = |-0,2025|) Nói chung, nâng cao nhiệt độ khơng khí có ảnh hưởng xấu đến tăng trưởng bề rộng vòng năm Du sam Tăng trưởng bề rộng vòng năm Du sam tồn mối quan hệ chặt chẽ (R2 = 50,5%) với yếu tố M3, M4, M5, M10, M12, M1-4 M11-12 Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Du sam phụ thuộc nhiều vào M5 (HSHQ = |0,4280|); M10 (HSHQ = |0,2268|) M4 (HSHQ = |0,2179|); thấp M12 (HSHQ = |-0,0089|) Nói chung, nâng cao lượng mưa từ tháng - có ảnh hưởng tốt tăng trưởng bề rộng vòng năm Du sam Trái lại, nâng cao tổng lượng mưa tháng 10 - 12 điều kiện xấu tăng trưởng bề rộng vòng năm Du sam Chỉ số Kd tồn mối quan hệ chặt chẽ (R2 = 49,3%) với yếu tố P4, P10, P12 Trong mối quan hệ này, yếu tố P10 (HSHQ = |0,3531|) đóng vai trò lớn so với P4 (HSHQ = |-0,3198|) P12 (HSHQ = |0,2334|) Ngồi ra, tăng trưởng bề rộng vòng năm Du sam có quan hệ với P1-4 P11-12 dạng hàm (9) - (10); P11-12 14 (HSHQ = |0,2798|) đóng vai trò lớn so với P1-4 (HSHQ = |0,2332|) Nói chung, nâng cao lượng nước bốc vào mùa khô (tháng - 4) điều kiện xấu tăng trưởng bề rộng vòng năm Du sam Trái lại, lượng nước bốc nhiều vào cuối mùa mưa (tháng 10 - 12) điều kiện tốt tăng trưởng bề rộng vòng năm Du sam Kết phân tích hàm phản hồi đa biến bước (Bảng 2) cho thấy, số Kd Du sam có quan hệ chặt chẽ (R2 = 65,2%) với yếu tố khí hậu (T4, T10, M5, M10, P4 P10) dạng hàm (1) (2) Sai lệch hàm (1) 3,7% Trong mối quan hệ này, ba yếu tố kiểm soát lớn tăng trưởng bề rộng vòng năm Du sam T4 (HSHQ = |-0,3321|), M5 (HSHQ = |0,2939|) P10 (HSHQ = |0,2255|) Bảng Phân tích ảnh hưởng tổng hợp nhiệt độ khơng khí, lượng mưa lượng nước bốc tăng trưởng bề rộng vòng năm Du sam TT (1) Hệ số hồi quy Chưa chuẩn Chuẩn hóa hóa (2) (3) (4) Hằng số 3,3524 T4 -1,7823 -0,3321 T10 -0,7091 -0,1833 M5 0,0743 0,2939 M10 -0,0392 -0,1494 P4 -0,0395 -0,1192 P10 0,1424 0,2255 Hàm (1) (2) Tham số R2 ±Se (5) (6) 66,5 0,051 MAP E (7) 3,7 3.4.2 Vai trò yếu tố khí hậu tăng trƣởng bề rộng vòng năm Bạch tùng Kết phân tích hồi quy tương quan đa tuyến tính bước cho thấy, số Kd Bạch tùng tồn mối quan hệ chặt chẽ (R = 24,9%) với yếu tố T6, T7 T10 Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc nhiều vào T6 (HSHQ = 15 |-0,2730|) T10 (HSHQ = |-0,2166|); thấp T7 (HSHQ = |0,0902|) Nói chung, nâng cao nhiệt độ khơng khí dẫn đến suy giảm tăng trưởng bề rộng vòng năm Bạch tùng Tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc chặt chẽ (R = 51,6%) vào yếu tố mưa (M3, M11, M1-4, M11-12 M11-3) Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc lớn vào M11 (HSHQ = |-0,6522|); M3 (HSHQ = |0,1841|) M1-4 (HSHQ = |0,1786|); thấp M11-3 (HSHQ = |0,0502|) Nói chung, nâng cao lượng mưa từ tháng 11 đến tháng 12 điều kiện bất lợi tăng trưởng bề rộng vòng năm Bạch tùng Trái lại, mưa lớn từ tháng – điều kiện tốt tăng trưởng bề rộng vòng năm Bạch tùng Biến động số Kd Bạch tùng phụ thuộc tương đối chặt chẽ (R = 36,0%) vào yếu tố độ ẩm không khí (Rh5, Rh7, Rh11, Rh11-12 Rh11-3) Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc lớn vào biến động Rh5 (HSHQ = |-0,3085|); R11-3 (HSHQ = |0,1975|) Rh11 (HSHQ = |-0,1971|); thấp Rh7 (HSHQ = |0,1152|) Nói chung, nâng cao độ ẩm khơng khí tháng 5, 11, 11 – 12 tháng 11 năm trước đến tháng năm sau dẫn đến suy giảm tăng trưởng bề rộng vòng năm Bạch tùng Lượng nước bốc nhiều vào tháng điều kiện thuận lợi tăng trưởng bề rộng vòng năm Bạch tùng Giữa số Kd yếu tố nắng (N1, N4, N7, N9, N11 N1-4) tồn mối quan hệ chặt chẽ (R2 = 48,5%) Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc lớn vào biến động N1 (HSHQ = |-0,2898|); N4 (HSHQ = |-0,2732|) N9 (HSHQ = |0,1792|); thấp N11 (HSHQ = |0,1471|) Nói chung, nâng cao số nắng vào tháng 1, dẫn đến suy giảm tăng trưởng bề rộng vòng năm Bạch tùng Trái lại, nắng nhiều vào tháng 11 điều kiện thuận lợi tăng trưởng bề rộng vòng năm Bạch tùng Chỉ số Kd tồn mối quan hệ chặt chẽ (R2 = 37,2%) với yếu tố lượng nước bốc (P3, P4, P9, P10, P11) Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc lớn vào biến động P11 (HSHQ = |0,3101|); P9 (HSHQ = |0,2807|) P3 (HSHQ = |-0,1398|); thấp P10 (HSHQ = |0,0053|) Mặt khác, số Kd tồn mối quan hệ tương đối chặt chẽ (R = 21,0%) với 16 P1-4 P11-12; yếu tố P11-12 (HSHQ = |0,2883|) đóng vai trò lớn so với P1-4 (HSHQ = |-0,2325|) Nói chung, nâng cao lượng nước bốc vào cuối mùa mưa (tháng 11 – 12) điều kiện tốt tăng trưởng bề rộng vòng năm Bạch tùng Trái lại, lượng nước bốc nhiều từ tháng – yếu tố bất lợi tăng trưởng bề rộng vòng năm Bạch tùng Nói chung, tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc lớn vào biến động yếu tố: T6, M11, Rh5, N1, N4, P9 P11 Bằng phân tích hàm phản hồi đa biến bước (Bảng 3) cho thấy, biến động số Kd Bạch tùng có quan hệ chặt chẽ (R = 65,1%) với yếu tố khí hậu dạng hàm (3) (4) Sai lệch hàm (3) 4,2% Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc lớn vào biến động M11 (HSHQ = |0,3363|); N4 (HSHQ = |-0,2659|) N1 (HSHQ = |-0,2318|); thấp P11 (HSHQ = |-0,1213|) Điều chứng tỏ ba yếu tố kiểm sốt lớn tăng trưởng bề rộng vòng năm Bạch tùng lượng mưa tháng 11, số nắng tháng tháng Bảng Phân tích ảnh hưởng tổng hợp nhiệt độ khơng khí, lượng mưa, độ ẩm khơng khí, số nắng lượng nước bốc tăng trưởng bề rộng vòng năm Bạch tùng TT (1) Hệ số hồi quy Tham số Chưa chuẩn hóa Chuẩn hóa (2) (3) (4) Hằng số 2,5995 T6 -0,7321 -0,1298 M11 -0,0402 -0,3363 Rh5 -0,4321 -0,1317 N1 -0,2111 -0,2318 N4 -0,2480 -0,2659 P9 0,1369 0,2027 P11 -0,0737 -0,1213 Hàm (3) (4) MAP E (5) (6) (7) 65,1 0,058 4,2 R2 ±Se 17 3.4.3 Vai trò yếu tố khí hậu tăng trƣởng bề rộng vòng năm Đỉnh tùng Kết phân tích hồi quy tương quan đa tuyến tính bước cho thấy, số Kd Đỉnh tùng tồn mối quan hệ chặt chẽ (R = 40,6%) với yếu tố nhiệt độ khống khí (T1, T2, T4, T5) Mặt khác, số Kd với yếu tố T1-4 T11-3 tồn mối quan hệ chặt chẽ (R2 = 35,8%) Trong mối quan hệ này, yếu tố T1 (HSHQ = |-0,3601|) ảnh hưởng đến tăng trưởng bề rộng vòng năm Đỉnh tùng lớn so với T2 (HSHQ = |-0,0918|), T4 (HSHQ = |-0,1492|) T5 (HSHQ = |0,1735|) Mặt khác, yếu tố T1-4 (HSHQ = |-0,5665|) đóng vai trò lớn so với T11-3 (HSHQ = |-0,0107|) Từ phân tích cho thấy, gia tăng nhiệt độ khơng khí vào mùa khô từ tháng – điều kiện xấu tăng trưởng bề rộng vòng năm Đỉnh tùng Tăng trưởng bề rộng vòng năm Đỉnh tùng có quan hệ chặt chẽ (R = 47,6%) với yếu tố mưa (M1, M3, M6 M8) Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Đỉnh tùng phụ thuộc lớn vào M1 (HSHQ = |-0,4412|); M8 (HSHQ = |0,3237|) M3 (HSHQ = |0,2888|); thấp M6 (HSHQ = |0,1110|)Nói chung, nâng cao lượng mưa vào tháng dẫn đến suy giảm tăng trưởng bề rộng vòng năm Đỉnh tùng Trái lại, mưa lớn vào tháng 3, có ảnh hưởng tốt tăng trưởng bề rộng vòng năm Đỉnh tùng Tăng trưởng bề rộng vòng năm Đỉnh tùng có mối quan hệ chặt chẽ (R2 = 46,8%) với yếu tố nắng (N1, N3, N6, N11) Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Đỉnh tùng phụ thuộc lớn vào N3 (HSHQ = |-0,3585|); N11 (HSHQ = |0,2943|) N6 (HSHQ = |-0,2201|); thấp N1 (HSHQ = |0,0857|) Mặt khác, tăng trưởng bề rộng vòng năm Đỉnh tùng tồn mối quan hệ chặt chẽ (R2 = 30,2%) với N1-4, N5-10, N11-12; yếu tố N1-4 (HSHQ = |-0,3137|) đóng vai trò lớn so với N5-10 (HSHQ = |-0,2751|) N1112 (HSHQ = |0,1939|) Từ phân tích cho thấy, gia tăng số nắng vào mùa khô từ tháng - mùa mưa từ tháng – 10 điều kiện bất lợi tăng trưởng bề rộng vòng năm Đỉnh tùng Trái lại, nắng nhiều vào tháng 11 12 yếu tố thuận lợi tăng trưởng bề rộng vòng năm Đỉnh tùng Từ phân tích cho thấy, tăng trưởng bề rộng vòng năm Đỉnh tùng bị kiểm sốt chặt chẽ yếu tố khí hậu: T1, T1-4, 18 M1, M8, N3 N11 Bằng phân tích hồi quy tương quan đa biến bước (Bảng 4) cho thấy, biến động số Kd có quan hệ chặt chẽ (R2 = 63,7%) với yếu tố khí hậu dạng hàm (5) (6) Sai lệch hàm (5) 4,6% Trong mối quan hệ này, tăng trưởng bề rộng vòng năm Đỉnh tùng phụ thuộc lớn vào biến động M1 (HSHQ = |-0,2890|); T1-4 (HSHQ = |-0,2698|) N11 (HSHQ = |0,2400|); thấp M8 (HSHQ = |0,1059|) Nói chung, ba yếu tố kiểm sốt lớn tăng trưởng bề rộng vòng năm Đỉnh tùng M 1, T1-4 N11 Sự nâng cao lượng mưa tháng nhiệt độ không khí tháng - dẫn đến ảnh hưởng xấu tăng trưởng bề rộng vòng năm Đỉnh tùng Trái lại, mưa nhiều vào tháng nắng nhiều vào tháng 11 yếu tố thuận lợi tăng trưởng bề rộng vòng năm Đỉnh tùng Bảng Phân tích ảnh hưởng tổng hợp nhiệt độ khơng khí, lượng mưa số nắng đến tăng trưởng bề rộng vòng năm Đỉnh tùng Hệ số hồi quy Chưa chuẩn hóa Chuẩn hóa (1) (2) (3) (4) Hằng số 3,1575 T1 -0,4903 -0,1230 T1-4 -1,5660 -0,2698 M1 -0,0246 -0,2890 M8 0,0316 0,1059 N3 -0,2206 -0,1931 N11 0,1150 0,2400 Hàm (5) (6) TT Tham số R2 ±Se MAPE (5) 63,7 (6) 0,066 (7) 4,6 3.5 Ảnh hƣởng độ cao địa hình đến mối quan hệ tăng trƣởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng với yếu tố khí hậu Ở độ cao 1.200 m so với mặt nước biển, tăng trưởng bề rộng vòng năm Du sam phụ thuộc chặt chẽ vào biến động yếu tố: T4, T7, T11-12, M4 M5 Ở độ cao 1.400 m so với mặt nước biển, trưởng bề rộng vòng năm Du sam phụ thuộc chặt chẽ vào biến động yếu tố: 19 T1, T6 T10, M1, M9 M11-12 Đối với toàn khu vực nghiên cứu, tăng trưởng bề rộng vòng năm Du sam phụ thuộc nhiều vào yếu tố: T4, T7, T5-10, M4, M5 M10 Như vậy, độ cao địa hình chi phối đến mối quan hệ tăng trưởng bề rộng vòng năm Du sam với nhiệt độ khơng khí lượng mưa Phản ứng Du sam với nhiệt độ khơng khí độ cao 1.400 m so với mặt nước biển biểu rõ rệt so với độ cao 1.200 m Đối với Bạch tùng, sống độ cao 1.400 m so với mặt nước biển, tăng trưởng bề rộng vòng năm bị kiểm soát chặt chẽ yếu tố: T9, T10, M6, M11, M11-12 M111-3 Trái lại, độ cao 1.600 m so với mặt nước biển, tăng trưởng bề rộng vòng năm Bạch tùng bị kiểm sốt yếu tố: T6, T7, T5-10, M3 M10 Đối với toàn khu vực nghiên cứu, tăng trưởng bề rộng vòng năm Bạch tùng phụ thuộc nhiều vào yếu tố: T6, T7, T10, M3 M11 Như vậy, độ cao địa hình chi phối đến mối quan hệ tăng trưởng bề rộng vòng năm Bạch tùng với nhiệt độ khơng khí lượng mưa Phản ứng Bạch tùng với nhiệt độ không khí độ cao 1.600 m so với mặt nước biển biểu rõ rệt so với độ cao 1.400 m Đối với Đỉnh tùng, mọc độ cao 1.400 m so với mặt nước biển, tăng trưởng bề rộng vòng năm bị kiểm sốt chặt chẽ yếu tố: T11, T11-3, M3 M10 Trái lại, sống độ cao 1.600 m so với mặt nước biển, tăng trưởng bề rộng vòng năm Đỉnh tùng bị kiểm soát chặt chẽ yếu tố: T2, T5, T1-4, M1 M8 Đối với toàn khu vực nghiên cứu, tăng trưởng bề rộng vòng năm Đỉnh tùng phụ thuộc lớn vào yếu tố: T1 T1-4, M1, M3 M8 Nói chung, độ cao địa hình chi phối đến mối quan hệ tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng với nhiệt độ không khí lượng mưa Nhiệt độ khơng khí gia tăng làm giảm tăng trưởng bề rộng vòng năm ba lồi gỗ Nói cách khác, ba lồi gỗ đòi hỏi chế độ nhiệt thấp Vì lý này, nên chúng phân bố vùng núi cao Ngoài ra, lượng mưa vào cuối mùa khô đến đầu mùa mưa cuối mùa mưa ảnh hưởng lớn đến tăng trưởng bề rộng vòng năm ba loài gỗ Sự khác biệt phản ứng Du sam, Bạch tùng Đỉnh tùng mọc độ cao địa hình khác với biến động yếu tố khí hậu giải thích tính nhạy cảm chúng với thiếu hụt nước cục 20 3.6 Phản ứng Du sam quần xã thực vật đất trống với biến động yếu tố khí hậu Kết nghiên cứu cho thấy, phản ứng Du sam với nhiệt độ khơng khí lượng mưa thay đổi tùy theo môi trường sống Khi mọc QXTV, tăng trưởng bề rộng vòng năm Du sam bị kiểm soát chặt chẽ yếu tố: T4, T10, M5 M10 Trái lại, mọc đất trống, tăng trưởng bề rộng vòng năm Du sam bị kiểm soát chặt chẽ yếu tố: T5, T12, M3 M8 Sự khác biệt phản ứng Du sam QXTV đất trống với biến động yếu tố khí hậu giải thích tính nhạy cảm lồi gỗ thay đổi điều kiện ánh sáng 3.7 Dự đoán tăng trƣởng Du sam, Bạch tùng Đỉnh tùng dựa theo yếu tố khí hậu Những phân tích hồi quy tương quan cho thấy, hàm dự đoán số Kd Du sam, Bạch tùng Đỉnh tùng dựa theo ba yếu tố khí hậu có dạng tương ứng hàm (7) – (9) KdDS = 2,56065 – 1,89668*T4 + 0,0679015*M5 + 0,267313*P10 (7) R = 59,2%; ±Se = 0,053; MAE = 0,041; MAPE = 4,2% KdBT = 1,50133 – 0,225009*N1 – 0,223064*N4 – 0,0526812*M11 (8) R2 = 57,0%; ±Se = 0,061; MAE = 0,043; MAPE = 4,5% KdDT = 3,21597 – 2,33996*T1-4 – 0,0250939*M1 + 0,152877*N11 (9) R2 = 55,4%; ±Se = 0,069; MAE = 0,049; MAPE = 5,0% Điều kiện thời tiết thuận lợi không thuận lợi tăng trưởng bề rộng vòng năm ba lồi gỗ phân chia thành cấp: tốt, bình thường xấu Ba cấp điều kiện thời tiết tương ứng với ba cấp sinh trưởng tốt (Kd > 1,05), trung bình (Kd = 0,95 – 1,05) (Kd < 0,95) Theo đó, năm có T4 < 18,50C, M5 > 325 mm P10 > 70 mm điều kiện tốt tăng trưởng bề rộng vòng năm Du sam Trái lại, năm có T4 > 19,50C, M5 < 125 mm P10 < 40 mm điều kiện xấu tăng trưởng bề rộng vòng năm Du sam Điều kiện tốt cho tăng trưởng bề rộng vòng năm Bạch tùng N < 200 giờ, N4 < 170 M11 < 15 mm Trái lại, năm có N1 > 270 giờ, N4 > 230 M11 > 180 mm điều kiện xấu tăng trưởng bề 21 rộng vòng năm Bạch tùng Điều kiện tốt cho tăng trưởng bề rộng vòng năm Đỉnh tùng M1 < mm, T1-4 < 170C N11 > 220 Trái lại, năm có M1 > 20 mm, T1-4 > 180C N11 < 125 điều kiện xấu tăng trưởng bề rộng vòng năm Đỉnh tùng Để dễ dàng cho việc đánh giá điều kiện thời tiết thuận lợi hay khó khăn tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng, ba biến khí hậu chuyển thành số khí hậu tổng hợp (X) Điểm số số X tích số ba biến khí hậu Theo đó, điểm số số X Du sam tính theo cơng thức: X DS = (T4*M5*P10)/10^3 Điểm số số X Bạch tùng tính theo cơng thức: XBT = (N1*N4*M11)/10^6 Điểm số số X Đỉnh tùng tính theo cơng thức: XBT = (T1-4*M1*N11)/10^3 Những phân tích thống kê cho thấy, hàm dự đoán số Kd = f(X) Du sam, Bạch tùng Đỉnh tùng tương ứng có dạng hàm (10) - (12) KdDS = (0,693802 + 0,0579755*ln(XDS))^2 r2 = 99,8%; ±Se = 0,001; P < 0,001 KdBT = exp(0,111185 – 0,0541944*sqrt(XBT)) r2 = 98,7%; ±Se = 0,010; P < 0,001 KdDT = 1/(0,916023 + 0,00368734*XDT) r2 = 88,5%; ±Se = 0,037; P < 0,001 (10) (11) (12) Chỉ số Kd Du sam, Bạch tùng Đỉnh tùng dự đoán cách thay điểm số ba số khí hậu vào ba hàm (11) – (12) Kết phân tích thống kê cho thấy, Du sam sinh trưởng tốt, trung bình xấu tương ứng với ba cấp điểm số số khí hậu X DS > 300, 125 – 300 < 125 Bạch tùng sinh trưởng tốt, trung bình xấu tương ứng với ba cấp điểm số số khí hậu XBT < 2, – > Đỉnh tùng sinh trưởng tốt, trung bình xấu tương ứng với ba cấp điểm số số khí hậu XDT < 10, 10 - 40 > 40 22 KẾT LUẬN VÀ ĐỀ NGHỊ Kết luận Khu vực Đà Lạt có tổng nhiệt độ khơng khí năm 6.5610C Nhiệt độ khơng khí trung bình năm 18C Lượng mưa trung bình năm 1.800 mm Độ ẩm khơng khí trung bình tháng 84% Số nắng trung bình tháng 181 Hệ số thủy nhiệt trung bình năm 2,68 Theo phân cấp chế độ khô ẩm Thái văn Trừng (1999), chế độ khô ẩm khu vực Đà Lạt thuộc cấp II (Hơi ẩm) Bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng có biến động lớn theo tuổi Chỉ số tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng có biến động nhỏ so với bề rộng vòng năm Biến động số tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng tương tự Những tín hiệu khí hậu biểu rõ rệt lớp vòng năm Du sam, Bạch tùng Đỉnh tùng Ba yếu tố khí hậu kiểm sốt mạnh tăng trưởng bề rộng vòng năm Du sam nhiệt độ khơng khí vào tháng 4, lượng mưa vào tháng lượng nước bốc vào tháng 10 Sự nâng cao nhiệt độ không khí vào tháng có ảnh hưởng xấu tăng trưởng bề rộng vòng năm Du sam Trái lại, nâng cao lượng mưa vào tháng lượng nước bốc vào tháng 10 điều kiện tốt tăng trưởng bề rộng vòng năm Du sam Ba yếu tố khí hậu kiểm sốt mạnh tăng trưởng bề rộng vòng năm Bạch tùng lượng mưa vào tháng 11, số nắng vào tháng tháng Sự nâng cao ba yếu tố dẫn đến ảnh hưởng xấu tăng trưởng bề rộng vòng năm Bạch tùng Ba yếu tố khí hậu kiểm sốt mạnh tăng trưởng bề rộng vòng năm Đỉnh tùng lượng mưa vào tháng 1, nhiệt độ khơng khí trung bình vào tháng – số nắng vào tháng 11 Sự nâng cao lượng mưa vào tháng nhiệt độ khơng khí trung bình vào tháng – điều kiện xấu tăng trưởng bề rộng vòng năm Đỉnh tùng Trái lại, nắng nhiều vào tháng 11 điều kiện tốt tăng trưởng bề rộng vòng năm Đỉnh tùng Độ cao địa hình có ảnh hưởng đến phản ứng Du sam, Bạch tùng Đỉnh tùng với khí hậu Khi sống độ cao 1.400 - 1.600 m so với 23 mặt nước biển, ba loài gỗ phản ứng với thay đổi yếu tố khí hậu rõ rệt so với độ cao 1.400 m Phản ứng tăng trưởng bề rộng vòng năm Du sam với thay đổi yếu tố khí hậu phụ thuộc vào điều kiện môi trường khác Trong điều kiện QXTV, tăng trưởng bề rộng vòng năm Du sam bị kiểm sốt chặt chẽ nhiệt độ khơng khí vào tháng 4, 10, lượng mưa vào tháng 10 Trái lại, sống đất trống, tăng trưởng bề rộng vòng năm Du sam bị kiểm sốt chặt chẽ nhiệt độ khơng khí vào tháng tháng 12, lượng mưa vào tháng Tăng trưởng bề rộng vòng năm Du sam dự đoán theo biến động nhiệt độ khơng khí vào tháng 4, lượng mưa vào tháng lượng nước bốc vào tháng 10 Tăng trưởng bề rộng vòng năm Bạch tùng dự đoán theo lượng mưa vào tháng 11, số nắng vào tháng Tăng trưởng bề rộng vòng năm Đỉnh tùng dự đốn theo lượng mưa vào tháng 1, nhiệt độ khơng khí trung bình vào tháng – số nắng vào tháng 11 Cấp điều kiện thời tiết thuận lợi khó khăn tăng trưởng ba lồi gỗ dự đốn theo điểm số số khí hậu tổng hợp Đề nghị Tác giả kiến nghị quan nghiên cứu quan tâm đến đặc tính sinh thái Du sam, Bạch tùng Đỉnh tùng cần tiếp tục nghiên cứu vấn đề (1) Xây dựng chuỗi niên đại vòng năm chuẩn hóa Du sam, Bạch tùng Đỉnh tùng khơng phạm vi tỉnh Lâm Đồng, mà khu vực khác nước (2) Xác định thời kỳ bắt đầu kết thúc trình sinh trưởng (mùa sinh trưởng) năm tình trạng vật hậu Du sam, Bạch tùng Đỉnh (3) Phân tích phản ứng Du sam, Bạch tùng Đỉnh biến động yếu tố mơi trường khác (địa hình, đất, QXTV, hoạt động lâm sinh…) 24 NHỮNG CƠNG TRÌNH ĐÃ CƠNG BỐ VỀ LUẬN ÁN Nguyễn Văn Nhẫn, 2017 Ảnh hưởng khí hậu đến tăng trưởng đường kính Bạch tùng (Dacrycarpus imbricatus (Blume) de Laub) khu vực Đà Lạt Đức Trọng thuộc tỉnh Lâm Đồng Tạp chí Khoa học lâm nghiệp; số 1-2017, trang 51 – 59 Nguyễn Văn Nhẫn, 2017 Ảnh hưởng nhân tố khí hậu đến tăng trưởng đường kính Đỉnh tùng (Cephalotaxus mannii Hook f.) khu vực Đà Lạt Đức Trọng thuộc tỉnh Lâm Đồng Tạp chí Nơng nghiệp Phát triển nông thôn; Số 7/2017, trang 127 – 133 Nguyễn Văn Nhẫn Nguyễn Văn Thêm, 2017 Phản ứng tăng trưởng đường kính Du sam (Keteleeria evelyniana Masters) thay đổi nhân tố khí hậu khu vực Đà Lạt Đức Trọng thuộc tỉnh Lâm Đồng Tạp chí KHKT Nơng lâm nghiệp, Trường Đại học Nơng Lâm Tp Hồ Chí Minh; Số 4/2017, trang 74 – 80  ... tích ảnh hưởng khí hậu đến tăng trưởng Du sam, Bạch tùng Đỉnh tùng điều kiện mơi trường khác (4) Xây dựng mơ hình dự đốn tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng theo yếu tố khí hậu. .. chuẩn hóa Du sam, Bạch tùng Đỉnh tùng khu vực Đức Trọng Đà Lạt thuộc tỉnh Lâm Đồng Luận án ba yếu tố khí hậu kiểm sốt chặt chẽ tăng trưởng bề rộng vòng năm Du sam nhiệt độ khơng khí vào tháng... Phản ứng Du sam, Bạch tùng Đỉnh tùng biến động yếu tố khí hậu (4) Xác định vai trò yếu tố khí hậu tăng trưởng bề rộng vòng năm Du sam, Bạch tùng Đỉnh tùng (5) Ảnh hưởng độ cao địa hình đến mối