Ảnh hưởng của khí hậu đến sinh trưởng thông 3 lá

Thông ba lá (Pinus keysia Royle ex Gordon) là loài cây tự nhiên của khu hệ thực vật núi vừa và cao (Thái Văn Trừng, 1998). Nó phân bố tự nhiên ở Ấn Độ, Miến Điện, Thái Lan, Trung Quốc, Lào và Việt Nam…. Ở Việt Nam, Thông ba lá phân bố ở khu vực Tây Nguyên, Hà Giang, Yên Bái và Lào Cai... Do Thông ba lá cho gỗ có chất lượng tốt và giá trị thương phẩm cao, nhu cầu thị trường lớn, đồng thời nó là loài cây dễ trồng, nên hiện nay Thông ba lá đã được trồng rộng rãi ở Tây Nguyên (Nguyễn Ngọc Lung, 1988, 1999). Tại Lâm Đồng, Thông ba lá đã được đưa vào trồng rừng từ năm 1978 trên đất đã mất rừng tự nhiên từ 5 10 năm. Mục tiêu của trồng rừng Thông ba lá là sản xuất gỗ với năng suất cao và chất lượng tốt để đáp ứng nhu cầu về gỗ làm nhà cửa, đồ mộc gia dụng (bàn, ghế…) và nguyên liệu giấy; đồng thời góp phần bảo vệ đất và nguồn nước, tạo cảnh quan đẹp dùng vào mục đích tham quan – du lịch (Nguyễn Ngọc Lung, 1988). Để đạt được mục tiêu đề ra, nhận thấy bên cạnh việc chọn lựa lập địa thích hợp, rừng Thông ba lá cần phải được nuôi dưỡng theo một chương trình lâm sinh khoa học. Nhưng muốn xây dựng được một chương trình lâm sinh khoa học, trước hết cần phải có những hiểu biết tốt về đặc tính sinh thái học của Thông ba lá. Trước đây đã có một số công trình nghiên cứu về rừng Thông ba lá ở Lâm Đồng, trong đó đáng kể nhất là những nghiên cứu về sinh trưởng và năng suất rừng của Nguyễn Ngọc Lung (1988, 1999). Ngoài ra là nghiên cứu thử nghiệm các phương thức khai thác tái sinh (Phó Đức Đỉnh, 1995); nghiên cứu sinh trưởng, năng suất và sinh khối rừng Thông ba lá (Lê Hồng Phúc, 1995). Một số nghiên cứu gần đây cũng đã hướng vào xem xét ảnh hưởng của khí hậu đến sinh trưởng của Thông ba lá (Phạm Trọng Nhân, 2001; Nguyễn Văn Thêm, 2003, 2004). Tuy vậy, những nghiên cứu này mới chỉ xem xét ảnh hưởng của khí hậu tới sinh trưởng của rừng trồng Thông ba lá tại khu vực Đà Lạt. Trong khi đó, Lâm Đồng là nơi cây Thông ba lá phân bố tự nhiên và được trồng trên không gian rất rộng, trải dài từ độ cao từ 500 m đến 1.500 m. Mặt khác, cho đến nay khoa học và thực tiễn vẫn còn thiếu nhiều thông tin về vai trò của các yếu tố khí hậu và phi khí hậu (địa hình, đất, hoạt động của con người…) đối với sinh trưởng và phát triển của rừng Thông ba lá. Xuất phát từ đó, đề tài “Ảnh hưởng của một số yếu tố khí hậu đến sinh trưởng của Thông ba lá (Pinus kesyia Royle ex Gordon) ở tỉnh Lâm Đồng” đã được đặt ra.2. Mục tiêu nghiên cứu 2.1. Mục tiêu chungXác định phản ứng của Thông ba lá đối với những yếu tố khí hậu ở ba khu vực Bảo Lộc, Di Linh và Đà Lạt để làm cơ sở khoa học cho kinh doanh rừng Thông ba lá.2.2. Mục tiêu cụ thểTừ mục tiêu chung, đề tài xác định 4 mục tiêu cụ thể sau đây:(1) Xây dựng ba chuỗi chỉ số bề rộng vòng năm chuẩn hóa của Thông ba lá ở ba khu vực Bảo Lộc, Di Linh và Đà Lạt, tỉnh Lâm Đồng.(2) Xác định điều kiện khí hậu của những tháng có ảnh hưởng rõ rệt đến tăng trưởng bề rộng vòng năm của Thông ba lá ở ba khu vực Bảo Lộc, Di Linh và Đà Lạt.(3) Xây dựng những mô hình để mô tả mối quan hệ giữa tăng trưởng bề rộng vòng năm của Thông ba lá với những yếu tố khí hậu.(4) Phân cấp điều kiện khí hậu thuận lợi cho sinh trưởng của Thông ba lá dựa trên những yếu tố khí hậu có ảnh hưởng rõ rệt đến sinh trưởng của Thông ba lá.2.3. Phạm vi nghiên cứuPhạm vi nghiên cứu của đề tài là quần thể Thông ba lá. Địa điểm nghiên cứu tại ba khu vực Bảo Lộc, Di Linh và Đà Lạt. Những yếu tố khí hậu được xem xét bao gồm nhiệt độ không khí, lượng mưa, độ ẩm không khí và số giờ nắng của 12 tháng trong năm. Những phản ứng của Thông ba lá đối với những yếu tố khí hậu được đánh giá thông qua biến động bề rộng vòng năm trên thân cây.Nội dung nghiên cứu tập trung làm rõ những khác biệt về khí hậu của ba khu vực Bảo Lộc, Di Linh và Đà Lạt; đặc điểm vòng năm và chỉ số vòng năm của Thông ba lá; mối liên hệ giữa tăng trưởng bề rộng vòng năm của Thông ba lá với nhiệt độ không khí, lượng mưa, độ ẩm không khí và số giờ nắng của những tháng và mùa trong năm. Từ kết quả nghiên cứu, đề xuất những mô hình dự đoán sinh trưởng của Thông ba lá và bảng phân cấp điều kiện khí hậu thuận lợi và khó khăn đối với sinh trưởng của Thông ba lá ở ba khu vực Bảo Lộc, Di Linh và Đà Lạt tỉnh Lâm Đồng.2.4. Ý nghĩa của đề tài Những kết quả nghiên cứu của đề tài luận văn này mang lại những ý nghĩa sau đây:(1) Về lý luận, đề tài góp phần làm rõ mối quan hệ giữa sinh trưởng của Thông ba lá với những yếu tố khí hậu. (2) Về thực tiễn, đề tài luận văn cung cấp những mô hình để dự đoán ảnh hưởng của khí hậu đến sinh trưởng của Thông ba lá. Từ đó có thể xây dựng những phương thức quản lý bền vững rừng Thông ba lá.

MỞ ĐẦU Đặt vấn đề Thông ba (Pinus keysia Royle ex Gordon) loài tự nhiên khu hệ thực vật núi vừa cao (Thái Văn Trừng, 1998) Nó phân bố tự nhiên Ấn Độ, Miến Điện, Thái Lan, Trung Quốc, Lào Việt Nam… Ở Việt Nam, Thông ba phân bố khu vực Tây Nguyên, Hà Giang, Yên Bái Lào Cai Do Thông ba cho gỗ có chất lượng tốt giá trị thương phẩm cao, nhu cầu thị trường lớn, đồng thời loài dễ trồng, nên Thông ba trồng rộng rãi Tây Nguyên (Nguyễn Ngọc Lung, 1988, 1999) Tại Lâm Đồng, Thông ba đưa vào trồng rừng từ năm 1978 đất rừng tự nhiên từ - 10 năm Mục tiêu trồng rừng Thông ba sản xuất gỗ với suất cao chất lượng tốt để đáp ứng nhu cầu gỗ làm nhà cửa, đồ mộc gia dụng (bàn, ghế…) nguyên liệu giấy; đồng thời góp phần bảo vệ đất nguồn nước, tạo cảnh quan đẹp dùng vào mục đích tham quan – du lịch (Nguyễn Ngọc Lung, 1988) Để đạt mục tiêu đề ra, nhận thấy bên cạnh việc chọn lựa lập địa thích hợp, rừng Thông ba cần phải nuôi dưỡng theo chương trình lâm sinh khoa học Nhưng muốn xây dựng chương trình lâm sinh khoa học, trước hết cần phải có hiểu biết tốt đặc tính sinh thái học Thông ba Trước có số công trình nghiên cứu rừng Thông ba Lâm Đồng, đáng kể nghiên cứu sinh trưởng suất rừng Nguyễn Ngọc Lung (1988, 1999) Ngoài nghiên cứu thử nghiệm phương thức khai thác - tái sinh (Phó Đức Đỉnh, 1995); nghiên cứu sinh trưởng, suất sinh khối rừng Thông ba (Lê Hồng Phúc, 1995) Một số nghiên cứu gần hướng vào xem xét ảnh hưởng khí hậu đến sinh trưởng Thông ba (Phạm Trọng Nhân, 2001; Nguyễn Văn Thêm, 2003, 2004) Tuy vậy, nghiên cứu xem xét ảnh hưởng khí hậu tới sinh trưởng rừng trồng Thông ba khu vực Đà Lạt Trong đó, Lâm Đồng nơi Thông ba phân bố tự nhiên trồng không gian rộng, trải dài từ độ cao từ 500 m đến 1.500 m Mặt khác, khoa học thực tiễn thiếu nhiều thông tin vai trò yếu tố khí hậu phi khí hậu (địa hình, đất, hoạt động người…) sinh trưởng phát triển rừng Thông ba Xuất phát từ đó, đề tài “Ảnh hưởng số yếu tố khí hậu đến sinh trưởng Thông ba (Pinus kesyia Royle ex Gordon) tỉnh Lâm Đồng” đặt Mục tiêu nghiên cứu 2.1 Mục tiêu chung Xác định phản ứng Thông ba yếu tố khí hậu ba khu vực Bảo Lộc, Di Linh Đà Lạt để làm sở khoa học cho kinh doanh rừng Thông ba 2.2 Mục tiêu cụ thể Từ mục tiêu chung, đề tài xác định mục tiêu cụ thể sau đây: (1) Xây dựng ba chuỗi số bề rộng vòng năm chuẩn hóa Thông ba ba khu vực Bảo Lộc, Di Linh Đà Lạt, tỉnh Lâm Đồng (2) Xác định điều kiện khí hậu tháng có ảnh hưởng rõ rệt đến tăng trưởng bề rộng vòng năm Thông ba ba khu vực Bảo Lộc, Di Linh Đà Lạt (3) Xây dựng mô hình để mô tả mối quan hệ tăng trưởng bề rộng vòng năm Thông ba với yếu tố khí hậu (4) Phân cấp điều kiện khí hậu thuận lợi cho sinh trưởng Thông ba dựa yếu tố khí hậu có ảnh hưởng rõ rệt đến sinh trưởng Thông ba 2.3 Phạm vi nghiên cứu Phạm vi nghiên cứu đề tài quần thể Thông ba Địa điểm nghiên cứu ba khu vực Bảo Lộc, Di Linh Đà Lạt Những yếu tố khí hậu xem xét bao gồm nhiệt độ không khí, lượng mưa, độ ẩm không khí số nắng 12 tháng năm Những phản ứng Thông ba yếu tố khí hậu đánh giá thông qua biến động bề rộng vòng năm thân Nội dung nghiên cứu tập trung làm rõ khác biệt khí hậu ba khu vực Bảo Lộc, Di Linh Đà Lạt; đặc điểm vòng năm số vòng năm Thông ba lá; mối liên hệ tăng trưởng bề rộng vòng năm Thông ba với nhiệt độ không khí, lượng mưa, độ ẩm không khí số nắng tháng mùa năm Từ kết nghiên cứu, đề xuất mô hình dự đoán sinh trưởng Thông ba bảng phân cấp điều kiện khí hậu thuận lợi khó khăn sinh trưởng Thông ba ba khu vực Bảo Lộc, Di Linh Đà Lạt tỉnh Lâm Đồng 2.4 Ý nghĩa đề tài Những kết nghiên cứu đề tài luận văn mang lại ý nghĩa sau đây: (1) Về lý luận, đề tài góp phần làm rõ mối quan hệ sinh trưởng Thông ba với yếu tố khí hậu (2) Về thực tiễn, đề tài luận văn cung cấp mô hình để dự đoán ảnh hưởng khí hậu đến sinh trưởng Thông ba Từ xây dựng phương thức quản lý bền vững rừng Thông ba Chương TỔNG QUAN 1.1 Khái quát Khí hậu – thực vật Khí hậu thực vật (Dendroclimatology) khoa học khôi phục lại khí hậu khứ cách sử dụng lớp vòng năm thân gỗ (Douglass, 1936; Douglass, 1937; Fritts, 1971; Bitvinskas, 1974) Khí hậu thực vật phân môn khoa học niên đại hay tuổi thọ gỗ Tiếp đầu ngữ Dendro xuất phát từ tiếng Hylạp, Dendron, có nghĩa gỗ Từ chronology tên ngành khoa học nghiên cứu thời gian xuất kiện đặc biệt Bởi cháy rừng, sâu bệnh hại rừng, khai thác rừng, hạn hán, lũ lụt ảnh hưởng đến cấu trúc, màu sắc bề rộng vòng năm, nên biến động bề rộng vòng năm dấu hiệu để xác định năm xuất kiện đặc biệt (Douglass, 1936; Douglass, 1937; Koerber And Wickman, 1970; Fritts, 1971; Fritts, 2001; Bitvinskas, 1974 ; Cook, 1985; Vaganov, 1996; Vương Văn Quỳnh Trần Tuyết Hằng, 1996; Nguyễn Văn Thêm, 2001, 2003; Phạm Trọng Nhân, 2001) Theo Bitvinskas (1974) Fritts (1971), từ bề rộng vòng năm truy tìm xác thời gian xuất kiện đặc biệt (khí hậu, lửa, sâu hại…) vì, tăng trưởng bề rộng vòng năm ấn định khí hậu, địa hình - đất, lửa, sâu hại… Những năm có khí hậu thuận lợi không thuận lợi (mưa nhiều hạn hán…) loài gỗ ghi lại thay đổi bề rộng màu sắc vòng năm Chính thế, người ta xem gỗ “nhà biên niên sử” tự nhiên Theo Bitvinskas (1974), Fritts (1971, 2001) Kozlowski (1971), vòng năm ghi lại xác tượng thời tiết năm mà chúng hình thành, nên đặc tính vòng năm sử dụng để truy tìm biến động khí hậu xuất định kỳ hay theo chu kỳ định Ngoài ra, phương pháp giúp dự đoán biến đổi khí hậu tương lai 1.2 Lịch sử nghiên cứu khí hậu – thực vật Những nghiên cứu mối quan hệ vòng năm gỗ với khí hậu tìm thấy công trình nghiên cứu Douglass (1936, 1937) Theo Bitvinksas (1974), phát triển vòng năm thân gỗ nguồn thông tin kiện tự nhiên xảy thời gian hình thành chúng Vì thế, thông qua đặc điểm vòng năm gỗ, xác định lượng mưa, gió, tuyết, lửa rừng hoạt động núi lửa cách hàng trăm năm Bằng phương pháp khí hậu thực vật, Kohler (1949) Kozlowski (1971) xác lập mối quan hệ chặt chẽ tăng trưởng bề rộng vòng năm gỗ với hoạt động mặt trời, khôi phục dự báo biến động trình tự nhiên lũ lũt hạn hán Phương pháp sinh khí hậu học sử dụng để nghiên cứu động thái nguồn nước, chế độ thủy văn, qui luật biến đổi khí hậu, mùa sai tái sinh rừng, suất diễn rừng, sâu bệnh hại rừng, lửa rừng ảnh hưởng người tới rừng (Douglass,1936); Schulman and Bryson, 1965; Bitvinskas, 1974) Nhiều nhà nghiên cứu cho rằng, việc xây dựng chuỗi biến động vòng năm thời gian dài thang chuẩn biến động vòng năm vùng địa lý riêng biệt, phân tích ảnh hưởng yếu tố sinh thái, đặc biệt hoạt động mặt trời, đến sinh trưởng suất rừng (Douglass, 1936; Douglass, 1937; Eklund, 1957) Phương pháp khí hậu thực vật áp dụng để phân tích ảnh hưởng khí hậu đến sinh trưởng gỗ Khi nghiên cứu tương quan nhiệt độ lượng mưa với biến động bề rộng vòng năm loài Pinus longaeva, Oberhuber (2002) nhận thấy bề rộng vòng năm nhỏ ảnh hưởng nhiệt độ thấp Fritts (1972) phát thấy sinh trưởng loài Picea glauca dọc theo kênh đào dòng suối phụ thuộc vào lượng mưa hàng năm Vào năm khô hạn, tăng trưởng vòng năm nhiều so với năm có lượng mưa lớn (dẫn theo Nguyễn Văn Thêm, 2010) Khi sử dụng phương pháp sinh khí hậu để phân tích biến động vòng năm Pinus sylvestris Varônhezơ (Nga), Vương Văn Quỳnh (1990) (dẫn theo Nguyễn Văn Thêm, 2010) rằng, thuộc cấp sinh trưởng khác có phản ứng không giống với khí hậu Ở lâm phần non, tăng trưởng Pinus sylvestris phụ thuộc chặt chẽ vào khí hậu Hoạt động mặt trời ảnh hưởng trực tiếp gián tiếp đến tăng trưởng Pinus sylvestris Sinh trưởng thuộc cấp sinh trưởng phụ thuộc vào hoạt động mặt trời Theo Fritt (1980) (dẫn theo Nguyễn Văn Thêm, 2010), sinh trưởng hai loài Abies lasiocarpa Pseudotsuga menziesli có mối liên hệ chặt chẽ với nhiệt độ lượng mưa Chỉ số tăng trưởng bề rộng vòng năm loài Pseudotsuga menziesli có mối quan hệ tuyến tính dương với lượng mưa từ tháng năm trước đến tháng 1, 2, tháng năm sau Ngược lại, số tăng trưởng bề rộng vòng năm loài Abies lasiocarpa có quan hệ tuyến tính dương với lượng mưa tháng 11 12 năm trước tháng 2, năm sau Lượng mưa lớn giúp cho loài Abies lasiocarpa tăng trưởng thời gian dài từ tháng 11 đến tháng Nghiên cứu Fritt Mayer cho thấy số tăng trưởng hai loài có tương quan dương với nhiệt độ tháng (tháng cuối mùa tăng trưởng) Bằng phương pháp hồi qui tuyến tính, Schulman Bryson (1965) dự đoán vòng năm loài Quercus rubra đạt tối đa lượng nước bốc tháng thấp, tổng lượng mưa tháng tháng cao, nhiệt độ bình quân tháng năm trước thấp lượng nước bốc tháng năm trước cao (dẫn theo Viện điều tra quy hoạch rừng, 1995) Phạm Trọng Nhân (2001) Nguyễn Văn Thêm (2003, 2004) ứng dụng phương pháp khí hậu thực vật để nghiên cứu phản ứng Thông ba với khí hậu khu vực Đà Lạt Kết nghiên cứu rằng, tăng trưởng bề rộng vòng năm Thông ba có quan hệ tuyến tính âm chặt chẽ với nhiệt độ không khí trung bình tháng 2, 3, 4, 9, 10 tập hợp tháng - Sự gia tăng số nắng tháng 2, tập hợp tháng từ tháng đến tháng 10 có ảnh hưởng xấu đến tăng trưởng bề rộng vòng năm Thông ba Biến động số độ ẩm không khí hàng tháng năm có ảnh hưởng không rõ rệt đến biến động số tăng trưởng bề rộng vòng năm Thông ba Sự gia tăng số thuỷ nhiệt tháng 2, 10 – 12 kéo theo suy giảm số tăng trưởng bề rộng vòng năm Thông ba Ngược lại, gia tăng số thủy nhiệt tháng – tháng lại có khuynh hướng kéo theo nâng cao số tăng trưởng bề rộng vòng năm Thông ba Biến động số tăng trưởng bề rộng vòng năm Thông ba phụ thuộc chặt chẽ vào biến động tổ hợp số nhiệt độ tháng 2, số lượng mưa tháng số nắng tháng Biến động số tăng trưởng bề rộng vòng năm Thông ba có mối quan hệ chặt chẽ với biến động tổ hợp số nhiệt độ tháng 9, số lượng mưa tháng số nắng tháng Giữa biến động số tăng trưởng bề rộng vòng năm Thông ba tổ hợp số nhiệt độ, số lượng mưa số nắng tháng 2, tồn mối quan hệ chặt chẽ 1.3 Thảo luận chung Từ tổng quan phương pháp khí hậu thực vật nghiên cứu quan hệ Thông ba với khí hậu, nhận thấy cần thảo luận thêm số vấn đề sau đây: (1) Phương pháp niên đại thực vật khí hậu thực vật phát triển mạnh Mỹ nhiều nước châu Âu khác Chúng cho phép xác định phản ứng gỗ với biến động khí hậu thông qua biến động bề rộng vòng năm theo năm lịch khác Vì thế, đề tài luận văn thạc sĩ áp dụng kiến thức hai môn khoa học niên đại thực vật khí hậu thực vật để phân tích phản ứng Thông ba với khí hậu (2) Khi nghiên cứu ảnh hưởng khí hậu đến sinh trưởng gỗ lâm phần, người ta thường tập trung xem xét ảnh hưởng nhiệt độ không khí, lượng mưa, độ ẩm không khí số nắng Sở dĩ bốn yếu tố khí hậu quan tâm nhiều chúng yếu tố định trình vật lý xảy khí trình sống gỗ Ngoài ra, chúng yếu tố dễ đo đạc… Vì thế, đề tài luận văn xem xét ảnh hưởng nhiệt độ không khí, lượng mưa, độ ẩm không khí số nắng 12 tháng mùa năm đến sinh trưởng Thông ba (3) Cho đến nay, phạm vi khu vực Đà Lạt, chưa có công trình nghiên cứu ảnh hưởng khí hậu đến tăng trưởng bề rộng vòng năm Thông ba khu vực khác tỉnh Lâm Đồng Nhận thấy rằng, sống điều kiện môi trường có điều kiện khí hậu phi khí hậu khác nhau, phản ứng Thông ba với khí hậu có biểu khác Vì thế, nghiên cứu ảnh hưởng khí hậu đến tăng trưởng bề rộng vòng năm Thông ba khu vực khác cần thiết Kết nghiên cứu cho phép: (a) đánh giá xác vai trò yếu tố khí hậu sinh trưởng phát triển Thông ba lá, (b) xây dựng mô hình dự đoán sinh trưởng Thông ba lá, (c) phân cấp điều kiện khí hậu thuận lợi khó khăn sinh trưởng Thông ba lá, (d) xây dựng phân chia rừng Thông ba theo điều kiện khí hậu khác nhau… Chương NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Đối tượng nghiên cứu Đối tượng nghiên cứu rừng Thông ba tự nhiên rừng trồng Những quần thụ phân bố độ cao từ 500 đến 1.500 m so với mặt biển Địa điểm nghiên cứu thực Bảo Lộc, Di Linh Đà Lạt Thời gian nghiên cứu từ tháng 4/2011 đến tháng 10/2011 Tại Bảo Lộc, mẫu gỗ thu thập rừng trồng Thông ba 32 tuổi Rừng Thông ba có mật độ trung bình 400 cây/ha; chiều cao bình quân 19,7 m; đường kính ngang ngực trung bình 32 cm; độ tàn che 0,5; tiết diện ngang 33,12 m2/ha; trữ lượng 331,16 m3/ha Hình 2.1 Bản đồ địa hình khu vực Bảo Lộc Tọa độ tâm ô thu thập mẫu gỗ (hệ VN 2000) 504.260 m kinh Đông; 1.289.096 m vĩ Bắc Địa hình đồi núi, bị chia cắt mạnh; độ cao tuyệt đối 530 m so với mặt biển; độ dốc trung bình từ 100 - 150 Đất feralit vàng đỏ phát triển đá Bazan Tại Di Linh, mẫu gỗ thu thập rừng Thông ba tự nhiên loại Mật độ trung bình 190 cây/ha; chiều cao bình quân 20,5 m; đường kính ngang ngực trung bình 45 cm; độ tàn che 0,4; tiết diện ngang 30,2 m2/ha; trữ lượng 310,98 m3/ha Tọa độ tâm ô thu thập mẫu gỗ (hệ VN 2000) 539.927 m kinh Đông; Hình 2.2 Bản đồ địa hình khu vực Di Linh 1.267.768 m vĩ Bắc Địa hình đồi núi với độ cao tuyệt đối 1.020 m so với mặt biển; độ dốc trung bình từ 20250 Đất feralit nâu đỏ phát triển đá granit Tại Đà Lạt, mẫu gỗ thu thập rừng Thông ba tự nhiên loại Mật độ trung bình 260 cây/ha đến; chiều cao bình quân 22,0 m; đường kính ngang ngực trung bình 52 cm; độ tàn che 0,6; tiết diện ngang 54,6 m2/ha; trữ lượng 594,6 m3/ha Tọa độ tâm ô thu thập mẫu gỗ (hệ VN 2000) 584.520 m kinh Đông; 1.312.736 m vĩ Hình 2.3 Bản đồ địa hình khu vực Đà Lạt Bắc Độ cao tuyệt đối 1.540 m so với mặt biển; độ dốc trung bình từ 220 Đất feralit vàng xám phát triển đá granit 10 mô hình dự đoán dựa nhiều nhân tố dự đoán, chi phí dự đoán cao Mặc dù vậy, có nhiều yếu tố khí hậu ảnh hưởng đến sinh trưởng Thông ba lá, đồng thời ảnh hưởng thay đổi tùy theo thời gian năm, nên việc phân cấp điều kiện khí hậu thuận lợi hay khó khăn cho sinh trưởng Thông ba xây dựng dựa mô hình biểu thị ảnh hưởng tổng hợp nhiều yếu tố khí hậu Xuất phát từ đó, yếu tố khí hậu (yếu tố dự đoán) có vai trò lớn tăng trưởng bề rộng vòng năm Thông ba đánh giá xếp hạng theo cấp; cấp – xấu, cấp – xấu, cấp – bình thường, cấp – tốt cấp - tốt Ảnh hưởng tổng hợp nhiều yếu tố khí hậu lên tăng trưởng bề rộng vòng năm Thông ba đánh giá theo tổng số cấp 3.6.1 Phân cấp điều kiện khí hậu Bảo Lộc Đối với khu vực Bảo Lộc, điều kiện khí hậu thuận lợi khó khăn cho sinh trưởng Thông ba dự đoán dựa mô hình 3.11 Mô hình dự đoán có dạng: Kd = 4,58121 – 4,09568*T3 + 0,489292*N3 R2 = 54,62%; Se = 0,1201; P = 0,0018 Phân tích số liệu khí hậu Bảo Lộc cho thấy trị trung bình, lớn nhất, nhỏ nhất, phạm vi biến động T3 N3 ghi lại Bảng 3.27 Bảng 3.27 Đặc trưng thống kê nhiệt độ không khí số nắng tháng khu vực Bảo Lộc tỉnh Lâm Đồng Thống kê (1) Trung bình Nhỏ Lớn Phạm vi biến động Nhiệt độ tháng (T3) (2) 22,3 21,0 23,2 2,2 Số nắng tháng (N3) (3) 222,3 153,0 265,6 112,6 Từ số liệu bảng 3.27 cho thấy phân chia nhiệt độ số nắng tháng thành cấp, khoảng cách cấp nhiệt độ số nắng tương 77 ứng 0,40C 22,5 Theo đó, từ mô hình 3.11 bảng 3.27 phân chia cấp thời tiết thuận lợi khó khăn cho tăng trưởng bề rộng vòng năm Thông ba Bảo Lộc (Bảng 3.28) Bảng 3.28 Phân cấp mức độ thuận lợi nhiệt độ không khí số nắng tháng tăng trưởng vòng năm Thông ba Bảo Lộc Phân cấp yếu tố khí hậu: T3 (0C) N3 (giờ) (1) (2) (3) < 21,7 > 256,5 21,7-22,1 233,5-256,5 22,1-22,5 210,5-233,5 22,5-22,9 187,5-210,5 > 22,9 < 187,5 TT Cấp sinh trưởng Cấp thời tiết (4) tốt tốt trung bình xấu xấu (5) Khi so sánh biến động số tăng trưởng bề rộng vòng năm Thông ba (Kd) với biến động tổng số cấp thời tiết tháng (X) Bảo Lộc từ năm 1988-2006 (Phụ lục 53), nhận thấy chúng có tương đồng cao (73,7%) Bằng thuật toán thống kê nhận thấy Kd Thông ba Bảo Lộc với X tồn quan hệ tuyến tính dương chặt chẽ theo dạng (Hình 3.47): Kd = 0,641701 + 0,0574334*X (3.34) r = 0,679; R2 = 46,15%; Se = ±0,1269; P = 0,0014 Nói chung, năm có tổng số cấp thời tiết tổng hợp tính theo nhiệt độ không khí số nắng tháng Bảo Lộc từ trở lên dẫn đến gia tăng bề rộng vòng năm Thông ba Ngược lại, số cấp thời tiết nhỏ 6, tăng trưởng bề rộng vòng năm Thông ba giảm Từ mối liên hệ dự đoán mức thuận lợi thời tiết tăng trưởng Thông ba Bảo Lộc 78 cách theo dõi nhiệt độ không khí số nắng tháng Sau thực đánh giá mức độ thuận lợi thời tiết tăng trưởng Thông ba Bảo Lộc theo cấp Bảng 3.28 Chỉ số Kd Tổng cấp thời tiết tháng (X) Năm Hình 3.47 Biến động số tăng trưởng bề rộng vòng năm Thông ba (Kd) tổng số cấp thời tiết tháng Bảo Lộc (X) 3.6.2 Phân cấp điều kiện khí hậu Di Linh Đối với khu vực Di Linh, điều kiện khí hậu thuận lợi khó khăn cho sinh trưởng Thông ba dự đoán dựa mô hình 3.21 Mô hình dự đoán có dạng: Kd = 1,70533 – 0,84033*T2 – 1,90828*T3 + 2,0349*R5 R2 = 80,09%; Se = ±0,052; P < 0,001 Phân tích số liệu khí hậu Di Linh cho thấy, trị trung bình, lớn nhất, nhỏ phạm vi biến động T2, T3 R5 ghi lại Bảng 3.29 Từ số liệu Bảng 3.29 cho thấy phân chia nhiệt độ không khí tháng và độ ẩm không khí tháng thành cấp, khoảng cách cấp nhiệt độ độ ẩm không khí tương ứng 0,4 0C 2% Theo đó, từ mô hình 3.21 Bảng 79 3.29, phân chia cấp thời tiết thuận lợi khó khăn cho tăng trưởng bề rộng vòng năm Thông ba Di Linh (Bảng 3.30) Bảng 3.29 Đặc trưng thống kê nhiệt độ không khí tháng 2, độ ẩm không khí tháng khu vực Di Linh tỉnh Lâm Đồng Thống kê (1) Trung bình Nhỏ Lớn Phạm vi biến động Nhiệt độ không khí (T,0C): tháng tháng (2) (3) 20,4 21,5 19,4 20,1 21,5 22,2 2,1 2,1 Độ ẩm không khí (R%): tháng (4) 82 76 86 10 Bảng 3.30 Phân cấp mức độ thuận lợi nhiệt độ không khí tháng 2, độ ẩm không khí tháng tăng trưởng vòng năm Thông ba Di Linh TT (1) Phân cấp yếu tố khí hậu: T3 (0C) T2 (0C) R5 (2) (3) (4) < 19,8 < 20,9 > 84 19,8-20,2 20,9-21,3 82-84 20,2-20,6 21,3-21,7 80-82 20,6-21,0 21,7-21,9 78-80 > 21,0 > 21,9 < 78 Cấp sinh trưởng Cấp thời tiết (5) tốt tốt trung bình xấu xấu (6) Những tính toán cho thấy số Kd Thông ba với biến động tổng số cấp thời tiết (X) Di Linh từ năm 1987 - 2006 tồn tương đồng cao (90,0%) Bằng thuật toán thống kê nhận thấy Kd Thông ba Di Linh với tổng số cấp thời tiết tổng hợp (X) tồn quan hệ tuyến tính dương chặt chẽ theo dạng (Hình 3.48): 80 Kd = 0,737533 + 0,02837*X (3.35) r = 0,788; R2 = 62,1%; Se = ±0,0676; P < 0,001 Nói chung, năm có tổng số cấp thời tiết tổng hợp tính theo nhiệt độ không khí tháng và độ ẩm không khí tháng Di Linh từ trở lên dẫn đến gia tăng bề rộng vòng năm Thông ba Ngược lại, số cấp thời tiết nhỏ 9, tăng trưởng bề rộng vòng năm Thông ba giảm Từ mối liên hệ này, dự đoán mức thuận lợi thời tiết tăng trưởng Thông ba Di Linh cách theo dõi nhiệt độ không khí tháng và độ ẩm không khí tháng Sau đánh giá mức độ thuận lợi thời tiết tăng trưởng Thông ba Di Linh dựa theo số liệu Bảng 3.29 Chỉ số Kd Tổng cấp thời tiết (X) Năm Hình 3.48 Yếu tố động số tăng trưởng bề rộng vòng năm Thông ba (Kd) tổng số cấp thời tiết tổng hợp Di Linh (X) 3.6.3 Phân cấp điều kiện khí hậu Đà Lạt Điều kiện khí hậu thuận lợi khó khăn cho sinh trưởng Thông ba Đà Lạt dự đoán dựa theo mô hình 3.32 Mô hình dự đoán có dạng: Kd = 4,410 - 1,009*T1 - 1,4*T6 - 0,152*M10 - 0,857*R12 R2 = 65,18%; Se = ±0,0747; P < 0,001 Phân tích số liệu khí hậu Đà Lạt cho thấy trị trung bình, lớn nhất, nhỏ phạm vi biến động T1, T6, M10 R12 dẫn Bảng 3.31 81 Từ số liệu Bảng 3.31 cho thấy phân chia T T6, M10 R12 thành cấp, khoảng cách cấp nhiệt độ, lượng mưa độ ẩm không khí tương ứng 0,40C, 77 mm 5% Theo đó, từ mô hình 3.30 số liệu Bảng 3.31 phân chia cấp thời tiết thuận lợi khó khăn cho tăng trưởng bề rộng vòng năm Thông ba khu vực Đà Lạt (Bảng 3.32) Bảng 3.31 Đặc trưng thống kê nhiệt độ không khí tháng 6, lượng mưa tháng 10 độ ẩm không khí tháng 12 khu vực Đà Lạt Nhiệt độ (T,0C): Mưa (mm): Độ ẩm (%): tháng tháng tháng 10 tháng 12 (1) (2) (3) (4) (5) Trung bình 15,8 19,0 249,6 80 Nhỏ 14,7 18,4 42,0 66 Lớn 16,8 20,2 425,7 92 Phạm vi biến động 2,1 1,7 383,7 26 Bảng 3.32 Phân cấp mức độ thuận lợi nhiệt độ tháng 6, mưa tháng 10 Thống kê độ ẩm không khí tháng 12 tăng trưởng vòng năm Thông ba Đà Lạt TT (1) Phân cấp yếu tố khí hậu: T6 (0C) R12 T1 ( C) M10 (2) (3) (4) (5) Cấp sinh trưởng Cấp thời tiết 87,5 (6) (7) tốt tốt trung bình xấu xấu Những tính toán cho thấy số Kd Thông ba với biến động tổng số cấp thời tiết tổng hợp (X) khu vực Đà Lạt từ năm 1987-2006 (Phụ lục 55) tồn tương đồng cao (90,0%) Chỉ số Kd với X tồn mối quan hệ tuyến tính dương chặt chẽ theo dạng (Phụ lục 57; Hình 3.49): Kd = 0,63258 + 0,02994*X (3.36) r = 0,707; R2 = 50,05%; Se = ±0,0839; P < 0,001 Nói chung, năm có tổng số cấp thời tiết tổng hợp tính theo nhiệt độ không khí tháng 6, lượng mưa tháng 10 độ ẩm không khí tháng 12 khu 82 vực Đà Lạt từ 12 trở lên dẫn đến gia tăng bề rộng vòng năm Thông ba Ngược lại, số cấp thời tiết nhỏ 12, tăng trưởng vòng năm Thông ba giảm Từ mối liên hệ này, dự đoán mức thuận lợi thời tiết tăng trưởng Thông ba khu vực Đà Lạt cách theo dõi nhiệt độ không khí tháng 6, lượng mưa tháng 10 độ ẩm không khí tháng 12 Sau đánh giá mức độ thuận lợi thời tiết tăng trưởng Thông ba dựa theo số liệu Bảng 3.32 Chỉ số Kd Tổng cấp thời tiết (X) Năm Hình 3.49 Yếu tố động số tăng trưởng vòng năm Thông ba (Kd) tổng số cấp thời tiết tổng hợp (X) khu vực Đà Lạt 3.7 Thảo luận chung kết nghiên cứu (1) Chuỗi niên đại vòng năm Thông ba thu thập 23 vòng năm khu vực Bảo Lộc (1986 - 2009), 39 vòng năm khu vực Di Linh (1972 - 2010) 58 vòng năm khu vực Đà Lạt (1953 - 2009) Kết nghiên cứu bề rộng vòng năm Thông ba có biến động lớn không theo tuổi, mà theo khu vực lấy mẫu Bề rộng vòng năm Thông ba có tượng tự tương quan tính nhạy cảm cao Tất tượng giải thích mẫu có tuổi nguồn gốc khác (tự nhiên nhân tạo), mọc nơi có điều kiện khí hậu phi khí hậu (địa hình – đất, hướng dốc, quần xã thực vật…) không giống (2) Chỉ số tăng trưởng bề rộng vòng năm Thông ba có biến động lớn; biến động mạnh xảy Thông ba phân bố khu vực 83 Di Linh Ngoài ra, số tăng trưởng bề rộng vòng năm Thông ba có tượng tự tương quan tính nhạy cảm cao Hiện tượng số tăng trưởng bề rộng vòng năm có biến động mạnh tính nhạy cảm cao giải thích khác biệt nguồn gốc mẫu biến động điều kiện môi trường sống Nghiên cứu cho thấy mọc ba khu vực (Bảo Lộc, Di Linh Đà Lạt) có khác khí hậu yếu tố môi trường khác (địa hình – đất, hướng dốc, quần xã thực vật…), nên biến động ba chuỗi số tăng trưởng bề rộng vòng năm Thông ba không tương đồng với Hiện tượng tương đồng biến động vòng năm xảy Thông ba mọc khu vục gần kề So với khu vực Bảo Lộc Đà Lạt, bề rộng vòng năm lẫn số tăng trưởng bề rộng vòng năm Thông ba Di Linh có biến động mạnh tính nhạy cảm cao Điều giải thích yếu tố khí hậu khu vực Di Linh thay đổi mạnh so khu vực Bảo Lộc Đà Lạt (3) Kết nghiên cứu cho thấy biến động số tăng trưởng bề rộng vòng năm Thông ba ba khu vực Bảo Lộc, Di Linh Đà Lạt tồn mối liện hệ âm với biến động nhiệt độ không khí trung bình tháng năm Điều chứng tỏ Thông ba cần chế độ nhiệt thấp Tuy vậy, ảnh hưởng nhiệt độ đến tăng trưởng Thông ba thay đổi tùy theo địa phương Khi địa hình nâng cao dần từ 500m (so với mặt nước biển) Bảo Lộc đến 1000m Di Linh 1.500m Đà Lạt, ảnh hưởng nhiệt độ không khí đến tăng trưởng vòng năm Thông ba thay đổi Tại khu vực Bảo Lộc, Thông ba có quan hệ âm chặt chẽ với nhiệt độ không khí tháng Trái lại, khu vực Di Linh Thông ba có quan hệ âm chặt chẽ với nhiệt độ không khí tháng 3, khu vực Đà Lạt tháng Nói chung, nâng cao nhiệt độ không khí vào mùa khô từ tháng đến có ảnh hưởng xấu đến tăng trưởng vòng năm Thông ba (4) Kết nghiên cứu biến động lượng mưa hàng tháng năm có ảnh hưởng đến tăng trưởng bề rộng vòng năm Thông ba Tuy vậy, quan hệ Thông ba với lượng mưa thay đổi rõ rệt tùy theo nơi Tại khu vực Bảo Lộc, gia tăng lượng mưa hàng tháng năm có ảnh 84 hưởng tốt đến tăng trưởng bề rộng vòng năm Thông ba Trái lại, mọc khu vực Di Linh Đà Lạt, Thông ba đòi hỏi lượng mưa cao từ tháng đến tháng 6, mưa nhỏ từ tháng đến tháng 12 Hiện tượng xảy giải thích Thông ba đòi hỏi nhiệt độ không khí thấp Nhiệt độ trung bình hàng tháng năm Bảo Lộc (21,90C) cao khu vực Di Linh (21,30C) Đà Lạt (17,90C) Vì thế, mưa lớn khu vực Bảo Lộc có tác dụng làm giảm chế độ nhiệt Trái lại, mưa lớn vào tháng cuối năm, nhiệt độ Di Linh Đà Lạt hạ xuống thấp Điều làm giảm tính chống chịu Thông ba với nhiệt độ (5) Kết nghiên cứu chứng tỏ độ ẩm không khí thay đổi có ảnh hưởng đến tăng trưởng vòng năm Thông ba Tại khu vực Bảo Lộc Di Linh, gia tăng độ ẩm không khí từ tháng 1- có ảnh hưởng tốt đến tăng trưởng bề rộng vòng năm Thông ba Nhưng độ ẩm cao vào tháng đến tháng 12 lại có ảnh hưởng xấu đến tăng trưởng bề rộng vòng năm Thông ba Tại khu vực Đà Lạt, độ ẩm không khí cao vào tháng năm có khuynh hướng làm giảm tăng trưởng bề rộng vòng năm Thông ba (6) Kết nghiên cứu nắng nhiều có ảnh hưởng xấu đến tăng trưởng bề rộng vòng năm Thông ba Tuy vậy, mối liên hệ chặt chẽ số tăng trưởng bề rộng vòng năm Thông ba với số nắng biểu rõ vào tháng Bảo Lộc tháng Di Linh (7) Kết nghiên cứu cho thấy Thông ba đòi hòi chế độ thủy nhiệt khác tùy theo địa phương Tại khu vực Bảo Lộc Đà Lạt, hệ số thủy nhiệt cao có ảnh hưởng tốt đến Thông ba Trái lại, phân bố khu vực Di Linh, Thông ba đòi hỏi hệ số thủy nhiệt cao vào tháng mùa khô (8) Bằng phương pháp phân tích hàm phản hồi bước, nhận thấy tăng trưởng bề rộng vòng năm Thông ba khu vực Bảo Lộc tồn mối quan hệ chặt chẽ với nhiệt độ không khí số nắng tháng Thông ba khu vực Di Linh phản hồi rõ rệt với nhiệt độ không khí tháng 2, độ ẩm không khí tháng Tương tự, Thông ba khu vực Đà Lạt tồn mối quan hệ chặt chẽ với 85 nhiệt độ không khí tháng 6, lượng mưa tháng 10 độ ẩm không khí tháng 12 Khi nghiên cứu ảnh hưởng khí hậu đến tăng trưởng bề rộng vòng năm rừng trồng Thông ba 22 tuổi Đà Lạt, Phạm Trọng Nhân (2001) Nguyễn Văn Thêm (2003, 2004) xác nhận nhiệt độ cao nắng nhiều có ảnh hưởng xấu đến Thông ba Phản ứng Thông ba biểu rõ rệt với nhiệt độ không khí trung bình tháng 9; nắng tháng 2, 9; lượng mưa tháng 9-10 Hai tác giả cho biến động độ ẩm không khí hàng tháng năm có ảnh hưởng không rõ rệt đến biến động số tăng trưởng bề rộng vòng năm Thông ba Nói chung, phản ứng Thông ba biểu rõ rệt với thay đổi yếu tố khí hậu Mức độ thay đổi Thông ba tùy thuộc vào vị trí địa lý nguồn gốc rừng Tuy vậy, từ nghiên cứu đề tài luận văn tác giả khác, nhận thấy nhiệt độ không khí cao từ tháng đến tháng 3, mưa lớn vào mùa mưa (tháng - 10) tháng có nắng nhiều gây ảnh hưởng xấu tăng trưởng bề rộng vòng năm Thông ba (9) Từ mối quan hệ chặt chẽ số tăng trưởng bề rộng vòng năm Thông ba với số khí hậu, phân chia cấp điều kiện khí hậu thuận lợi khó khăn cho tăng trưởng bề rộng vòng năm Thông ba Tại khu vực Bảo Lộc, năm có tổng số cấp thời tiết tổng hợp tính theo nhiệt độ không khí số nắng tháng từ trở lên dẫn đến gia tăng bề rộng vòng năm Thông ba Tại khu vực Di Linh, năm có tổng số cấp thời tiết tổng hợp tính theo nhiệt độ không khí tháng 2, độ ẩm không khí tháng từ trở lên dẫn đến gia tăng bề rộng vòng năm Thông ba Tại khu vực Đà Lạt, năm có tổng số cấp thời tiết tổng hợp tính theo nhiệt độ không khí tháng 6, lượng mưa tháng 10 độ ẩm không khí tháng 12 từ 12 trở lên dẫn đến gia tăng bề rộng vòng năm Thông ba 3.8 Một số đề xuất 3.8.1 Dự đoán tăng trưởng Thông ba 86 Từ kết qủa nghiên cứu cho thấy số tăng trưởng bề rộng vòng năm Thông ba có tượng tự tương quan với Ngoài ra, số tăng trưởng bề rộng vòng năm Thông ba có mối liện hệ chặt chẽ với yếu tố khí hậu Đó hai để dự đoán tăng trưởng Thông ba 3.8.1.1 Dự đoán tăng trưởng Thông ba dựa theo tượng tự tương quan số tăng trưởng bề rộng vòng năm Kết nghiên cứu số tăng trưởng bề rộng vòng năm Thông ba năm (Kd H) năm trước (KdT) tồn mối quan hệ chặt chẽ với theo dạng: + Đối với Thông ba khu vực Bảo Lộc Kd(H) = 1,54467 – 0,55550*Kd(T) (3.37) r = -0,501; R2 = 25,13%; Se = ±0,1269; P = 0,0243 + Đối với Thông ba khu vực Di Linh Kd(H) = 1,55269 – 0,56688*Kd(T) (3.38) r = -0,549; R2 = 30,13%; Se = ±0,0942; P = 0,0149 + Đối với Thông ba khu vực Đà Lạt Kd(H) = 1,58524 – 0,598296*Kd(T) (3.39) r = -0,586; R2 = 34,40%; Se = ±0,0981; P = 0,0016 Từ mô hình 3.37 - 3.39 việc dự đoán số tăng trưởng bề rộng vòng năm Thông ba thực sau: + Trước hết, xác đinh vòng năm trở lên Thông ba khoan tăng trưởng Những vòng năm xác định gần với năm + Kế đến, từ chuỗi bề rộng vòng năm, tính số tăng trưởng bề rộng vòng năm di động năm + Tiếp theo, tính số tăng trưởng bề rộng vòng năm năm sau cách thay số tăng trưởng bề rộng vòng năm năm vào mô hình 3.37 - 3.39 3.8.1.2 Dự đoán tăng trưởng Thông ba dựa theo yếu tố khí hậu 87 Tăng trưởng Thông ba dự đoán dựa mô hình phản hồi Thông ba yếu tố khí hậu Để đạt mục đích này, tác giả đề xuất sử dụng mô hình dự đoán sau đây: + Đối với khu vục Bảo Lộc TT Mô hình dự đoán tăng trưởng Thông ba Kd = 5,8242 – 4,84744*T3 Kd = -1,3911 + 2,3703*R3 Kd = -0,01299 + 0,98410*N3 Kd = 4,58121 – 4,09568*T3 + 0,489292*N3 Công thức (3.5) (3.8) (3.9) (3.11) + Đối với khu vực Di Linh TT Mô hình dự đoán tăng trưởng Thông ba Kd = 5,27874 – 1,17466*T2 – 3,11467*T3 Kd = 3,14745 – 2,15853*T2 Kd = 5,06576 – 4,07576*T3 Kd = 1,70533 – 0,84033*T2 – 1,90828*T3 + 2,0349*R5 Công thức (3.13) (3.14) (3.15) (3.19) + Đối với khu vực Đà Lạt TT Mô hình dự đoán tăng trưởng Thông ba Kd = 0,95183 + 0,04084*M1 Kd = 2,20404 – 1,20884*R12 Kd = 0,95193 + 0,04073*K1 Kd = 4,410 – 1,009*T1 – 1,4*T6 – 0,152*M10 – 0,857*R12 Công thức (3.24) (3.26) (3.28) (3.30) Trong thực tế, việc dự đoán số tăng trưởng bề rộng vòng năm Thông ba thực theo bước Bước Thu thập yếu tố khí hậu tương ứng với mô hình dự đoán ba khu vực Sau tính số hậu (nhiệt độ, mưa, ẩm) theo phương pháp bình quân di động năm khu vực Bảo Lộc năm khu vực Di Linh Đà Lạt Bước Xác định Kd Thông ba cách thay số khí hậu vào mô hình tương ứng Nếu Kd lớn 1, yếu tố hậu (nhiệt độ, mưa, ẩm) tháng thuận lợi cho tăng trưởng bề rộng vòng năm Thông ba Ngược lại, Kd nhỏ 1, yếu tố hậu (nhiệt 88 độ, mưa, ẩm) tháng không thuận lợi cho tăng trưởng bề rộng vòng năm Thông ba 3.8.1.3 Dự đoán điều kiện khí hậu thuận lợi cho tăng trưởng Thông ba Để dự đoán điều kiện khí hậu thuận lợi hay khó khăn cho tăng trưởng bề rộng vòng năm Thông ba ba khu vực Bảo Lộc, Di Linh Đà Lạt, tác giả đề xuất sử dụng tương ứng Bảng 3.28, 3.30 3.32 Theo đó, trình tự bước dự báo điều kiện khí hậu thuận lợi hay khó khăn cho tăng trưởng bề rộng vòng năm Thông ba sau: + Trước hết, thu thập yếu tố khí hậu tương ứng với bảng phân cấp điều kiện khí hậu thuận lợi hay khó khăn cho tăng trưởng bề rộng vòng năm Thông ba khu vực + Kế đến, từ yếu tố khí hậu, tra Bảng 3.28, 3.30 3.32 để xác định số cấp thời tiết tương ứng với yếu tố khí hậu + Tiếp theo, tính tổng số cấp thời tiết tương ứng với ba khu vực Bảo Lộc, Di Linh Đà Lạt + Sau cùng, thay tổng số cấp thời tiết tổng hợp vào ba mô hình 3.32, 3.33 3.34 tương ứng với ba khu vực Bảo Lộc, Di Linh Đà Lạt Nếu số Kd lớn 1, tổng số cấp thời tiết tổng hợp năm thuận lợi cho tăng trưởng Thông ba Ngược lại, số Kd nhỏ 1, tổng số cấp thời tiết tổng hợp năm không thuận lợi cho tăng trưởng Thông ba 89 KẾT LUẬN VÀ ĐỀ NGHỊ Kết luận Từ kết qủa nghiên cứu đến kết luận sau đây: (1) Tại khu vực Bảo Lộc có tổng nhiệt độ năm 7.988 0C, nhiệt độ không khí trung bình năm 21,90C; lượng mưa trung bình năm 2.973,0 mm; độ ẩm trung bình năm 86%; số nắng trung bình tháng 169,4 giờ; hệ số thủy nhiệt trung bình năm 3,67 Tại khu vực Di Linh có tổng nhiệt độ năm 7.766 0C, nhiệt độ không khí trung bình năm 21,3 0C; lượng mưa trung bình năm 1.410,0 mm; độ ẩm trung bình năm 81%; số nắng trung bình tháng 199,7 giờ; hệ số thủy nhiệt trung bình năm 1,75 Tại khu vực Đà Lạt có tổng nhiệt độ năm 6.5370C, nhiệt độ không khí trung bình năm 17,9 0C; lượng mưa trung bình năm 1.805,0 mm; độ ẩm trung bình năm 84%; số nắng trung bình tháng 182,9 giờ; hệ số thủy nhiệt trung bình năm 2,67 Theo phân loại khí hậu Thái Văn Trừng (1998), khí hậu Bảo Lộc thuộc cấp I (mưa ẩm), tháng khô, hạn kiệt Khí hậu Di Linh 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) tháng kiệt 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 kiệt (2) Chỉ số tăng trưởng bề rộng vòng năm Thông ba có tượng tự tương quan tính nhạy cảm cao (3) Phản ứng bề rộng vòng năm Thông ba khí hậu thay đổi tùy theo nơi Tại khu vực Bảo Lộc, bề rộng vòng năm Thông ba có phản ứng rõ rệt nhiệt độ không khí số nắng tháng Tại khu vực Di Linh, nhiệt độ không khí tháng 2, độ ẩm không khí tháng có ảnh hưởng rõ rệt đến tăng trưởng bề rộng vòng năm Thông ba Tại khu vực Đà Lạt, bề rộng vòng năm Thông ba phản ứng rõ rệt với nhiệt độ không khí tháng 6, lượng mưa tháng 10 độ ẩm không khí tháng 12 90 (4) Tăng trưởng bề rộng vòng năm Thông ba dự đoán dựa theo mối quan hệ với yếu tố khí hậu tổng số cấp thời tiết tổng hợp Thông ba Bảo Lộc, Di Linh Đà Lạt sinh trưởng tốt tổng số cấp thời tiết tổng hợp tương ứng lớn 6, 12 Đề nghị Đề tài xác định mối liên hệ tăng trưởng Thông ba với yếu tố khí hậu Tuy vậy, tác giả nhận thấy nhiều vấn đề chưa làm sáng tỏ Vì thế, tác giả kiến nghị quan tâm đến mối quan hệ Thông ba với khí hậu cần tiếp tục nghiên cứu số vấn đề sau đây: (1) Xác định thời kỳ bắt đầu kết thúc mùa sinh trưởng năm Thông ba Đây khoa học để giải thích mối liên hệ tăng trưởng bề rộng vòng năm Thông ba với yếu tố môi trường (khí hậu, địa hình – đất…) (2) Tiếp tục nghiên cứu phản ứng Thông ba với khí hậu yếu tố môi trường khác (địa hình-đất, quần xã thực vật, lửa, côn trùng…) khu vực khác tỉnh Lâm Đồng Những thông tin cần thiết để xây dựng ngân hàng vòng năm đồ phân bố Thông ba dựa theo phản ứng với khí hậu 91 [...]... gồm: (1) Đặc điểm khí hậu Bảo Lộc, Di Linh và Đà Lạt (2) Đặc điểm bề rộng vòng năm của Thông ba lá (3) Phản ứng của Thông ba lá với khí hậu Bảo Lộc (4) Phản ứng của Thông ba lá với khí hậu Di Linh (5) Phản ứng của Thông ba lá với khí hậu Đà Lạt (6) Phân cấp điều kiện khí hậu đối với sinh trưởng của Thông ba lá (7) Một số đề xuất 2 .3 Phương pháp nghiên cứu 2 .3. 1 Cơ sở khoa học Cơ sở của phương pháp luận... thập theo lý lịch rừng 2 .3. 2 .3 Xác định ảnh hưởng của khí hậu đến sinh trưởng của Thông ba lá Để xác định ảnh hưởng của khí hậu đến sinh trưởng của Thông ba lá, đã chọn 9 cây mẫu để thu thập vòng năm (khu vực Bảo Lộc 3 cây, khu vực Di Linh 3 cây và khu vực Đà Lạt 3 cây) Đối với mỗi cây mẫu, đã thu thập mẫu gỗ bằng khoan tăng trưởng; trong đó mỗi cây được khoan 2 mẫu ở vị trí D 1 .3 (cm) theo 2 hướng lên... vật và hoạt động của con người) Nói khác đi, sinh trưởng của Thông ba lá chịu ảnh hưởng tổng hợp của nhiều yếu tố môi trường bên ngoài (khí hậu, địa hình – đất, tác động của sinh vật và con người) và bên trong quần thể (kết cấu và cấu trúc quần thể, sự tương tác giữa cá cá thể, tiểu khí hậu ) Vì thế, phân tích ảnh hưởng tổng hợp của các yếu tố môi trường đến sinh trưởng của Thông ba lá cần được đặt... sống 3. 2.2 Đặc điểm vòng năm của Thông ba lá ở khu vực Di Linh 3. 2.2.1 Tăng trưởng bề rộng vòng năm của Thông ba lá ở Di Linh Số vòng năm của Thông ba lá tại Di Linh được thu thập là 39 năm, tương ứng với năm lịch từ 1972 đến 2010 Từ số liệu vòng năm thực nghiệm, đã xây dựng mô hình mô tả quy luật yếu tố đổi Zr theo A dưới dạng (Bảng 3. 10; Hình 3. 9): 35 Zr = 9 ,31 1 83* exp(-0, 033 34*A) – 0 ,39 268 (3. 2)... 4,40 2,89 8 ,32 5,85 4,08 2,81 4 7,47 3, 38 3, 18 2,50 8, 03 5,64 3, 93 2,71 5 6,90 5,44 2,99 2,55 7,76 5,45 3, 79 2,60 6 4 ,32 10,79 4,20 2,15 7,49 5,25 3, 65 2,51 7 7,51 9,01 3, 88 2,74 7, 23 5,07 3, 52 2,41 8 7,90 6,49 2,60 2,41 6,98 4,89 3, 39 2 ,32 9 6,28 6,55 2,99 2, 13 6,74 4,72 3, 27 2, 23 + Bề rộng vòng năm trung bình thực tế 39 năm là 4, 73 ± 0,41 mm; dao động từ 2, 13- 10,96 mm; biến động rất mạnh ( 53, 85%) Chuỗi... 25,2 2 63, 0 7,4 18,1 29,5 3 112,2 39 ,3 76,5 76,0 4 216,1 84,0 166,2 155,4 5 255,2 158,0 210 207,7 6 32 7,7 194,9 196,8 239 ,8 7 39 7,4 206,4 219,1 274 ,3 8 525,0 235 ,0 2 53, 4 33 7,8 9 4 03, 5 194,8 274,1 290,8 10 35 2,4 170,4 249,6 257,5 11 185,0 79,7 99,2 121 ,3 12 77,8 29,7 34 ,3 47 ,3 247,8 117,5 150,4 171,9 2.9 73, 0 1.410,0 1.804,0 2.0 63, 0 Trung bình tháng Trung bình cả năm Nguồn: Tính toán từ số liệu của các... (2) (3) 23 23 4,19 4,20 1,65 1 ,38 0 ,34 0 ,30 1,69 1,97 7,78 6,44 39 ,41 32 ,81 0,7 93 1,00 Thống kê (1) + Số vòng năm (n) + ∆r (mm) +S + Sx + Zrmin (mm) + Zrmax (mm) + V% + R+ 33 + msx 0,200 0,054 3. 2.1.2 Đặc điểm chỉ số vòng năm của Thông ba lá ở Bảo Lộc Từ chuỗi bề rộng năm thực tế của Thông ba lá ở Bảo Lộc (Bảng 3. 6) đã tính được chỉ số tăng trưởng bề rộng vòng năm theo phương pháp bình quân di động 3. .. 230 ,67495*exp(-0,000889*A) -224, 034 48 (3. 1) R2 = 69,4%; P < 0,01 Bảng 3. 6 Chuỗi bề rộng vòng năm của Thông ba lá ở khu vực Bảo Lộc Lý thuyết Thực nghiệm Zr(mm) Năm 1980 1990 2000 1980 1990 2000 0 1 2 3 4 5,15 5,90 5 ,38 6,42 7,78 3, 66 2,74 3, 26 2,94 2,97 5,62 5,41 5,21 5,01 4,80 3, 58 3, 38 3, 18 2,98 2,78 5 4,52 1,69 4,60 2,58 6 6,74 4,05 1,92 6,44 4,40 2 ,37 7 5,87 3, 84 3, 18 6, 23 4,19 2,17 8 4,86 2,74 2,50 6, 03 3,99 1,97 9... biểu đồ, phân tích quy luật biến đổi của các yếu tố khí hậu ở ba khu vực Bảo Lộc, Di Linh và Đà Lạt Để loại bỏ khuynh hướng biến đổi theo thời gian, những yếu tố khí hậu trung bình tháng cũng được chuyển thành chỉ số tương tự như chỉ số tăng trưởng bề rộng vòng năm 2 .3. 3.6 Xác định phản ứng của Thông ba lá với khí hậu Phản ứng của Thông ba lá với khí hậu được xác định thông qua mối quan hệ giữa chỉ số... vòng năm của Thông ba lá ở khu vực Bảo Lộc 3. 2.1.1 Tăng trưởng bề rộng vòng năm của Thông ba lá ở Bảo Lộc Số vòng năm của Thông ba lá tại Bảo Lộc được thu thập là 23 năm, tương ứng với năm lịch từ 1986 đến 2009 Kết quả nghiên cứu cho thấy quá trình yếu tố đổi bề rộng vòng năm của Thông ba lá (Zr, mm) theo tuổi cây (A, năm lịch) có thể mô tả bằng hàm số mũ dưới dạng (Bảng 3. 6; Hình 3. 7): Zr = 230 ,67495*exp(-0,000889*A)