Bài báo phân tích ảnh hưởng của các yếu tố yếu tố khí tượng (KT) đến năng suất bưởi và chuối ở Phú Thọ dựa trên tương quan và hồi qui tuyến tính nhằm xác định những biện pháp kỹ thuật phù hợp cho gia tăng năng suất. Kết quả cho thấy đối với bưởi, số ngày rét và số lần xuất hiện gió hướng đông bắc cao là những yếu tố ảnh hưởng điển hình đến giai đoạn ra hoa và quả nhỏ của bưởi. Độ ẩm cao cũng là yếu tố thuận lợi cho hoạt động của nấm và sâu bệnh hại quả ở giai đoạn bưởi chín. Đối với chuối, tốc độ gió mạnh là yếu tố ảnh hưởng khá điển hình trong tất cả giai đoạn phát triển và hình thành năng suất. Do đó, nên áp dụng những biệt pháp kỹ thuật phù hợp để giảm ảnh hưởng của các yếu tố KT này đến giai đoạn ra hoa, quả chín của bưởi và chuối như giữ ấm gốc cây, hay trồng cây chắn gió,…Bên cạnh đó, bước đầu đã xây dựng được hai mô hình ước lượng năng suất bưởi và chuối nhằm hướng tới đánh giá tác động của biến đổi khí hậu, phân vùng chuyên canh, cũng như dự báo phục vụ xây dựng kế hoạch tiêu thụ và xuất khẩu đối với bưởi và chuối
TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Bài báo khoa học Ảnh hưởng yếu tố khí tượng suất số ăn tỉnh Phú Thọ dựa tương quan hồi qui Nguyễn Đăng Mậu1*, Trịnh Hồng Dương1, Trần Thị Tâm1, Nguyễn Hữu Quyền1, Nguyễn Bình Phong2, Đặng Như Ý2, Lê Văn Phong1 Viện khoa học Khí tượng thủy văn Biến đổi khí hậu; mau.imhen@gmail.com; hoangduongktnn@gmail.com; trantam1810@gmail.com; nhquyen13@gmail.com; lephong2341999@gmail.com Đại học Tài Nguyên Môi trường Hà Nội; nbphong@hunre.edu.vn *Tác giả liên hệ: mau.imhen@gmail.com; Tel.: +84–382072468 Ban Biên tập nhận bài: 8/6/2023; Ngày phản biện xong: 19/7/2023; Ngày đăng bài: 25/8/2023 Tóm tắt: Bài báo phân tích ảnh hưởng yếu tố yếu tố khí tượng (KT) đến suất bưởi chuối Phú Thọ dựa tương quan hồi qui tuyến tính nhằm xác định biện pháp kỹ thuật phù hợp cho gia tăng suất Kết cho thấy bưởi, số ngày rét số lần xuất gió hướng đơng bắc cao yếu tố ảnh hưởng điển hình đến giai đoạn hoa nhỏ bưởi Độ ẩm cao yếu tố thuận lợi cho hoạt động nấm sâu bệnh hại giai đoạn bưởi chín Đối với chuối, tốc độ gió mạnh yếu tố ảnh hưởng điển hình tất giai đoạn phát triển hình thành suất Do đó, nên áp dụng biệt pháp kỹ thuật phù hợp để giảm ảnh hưởng yếu tố KT đến giai đoạn hoa, chín bưởi chuối giữ ấm gốc cây, hay trồng chắn gió,…Bên cạnh đó, bước đầu xây dựng hai mơ hình ước lượng suất bưởi chuối nhằm hướng tới đánh giá tác động biến đổi khí hậu, phân vùng chuyên canh, dự báo phục vụ xây dựng kế hoạch tiêu thụ xuất bưởi chuối Từ khóa: Ảnh hưởng yếu tố khí tượng; Năng suất bưởi chuối Mở đầu Nông nghiêp thường hoạt động trực tiếp bầu trời, thời tiết, khí hậu đất đai đóng vai trị quan trọng đến trồng Nơng nghiệp liên hệ phức tạp điều kiện tự nhiên, yếu tố thời tiết, khí hậu yếu tố tác động đến nông nghiệp, tác động thể qua dao động suất trồng qua năm cao hay thấp, chất lượng cao hay thấp [1] Thật vậy, việc khơng có đợt lạnh sớm vào tháng 12 tháng có tác động bất lợi đến nhu cầu tích lũy lạnh cho phân hóa mầm hoa ăn ôn đới Đối với táo, thiếu thời gian tích lũy lạnh dẫn đến hoa khơng hồn thiện, tỉ lệ đậu trái kém, chất lượng suất táo giảm [2] Nhiệt độ tăng 1,8°C, nhu cầu lạnh không đủ, thời gian hoa táo muộn hơn, kéo dài giai đoạn hoa thu hoạch dài táo [3] Ngày hoa cao điểm anh đào (Cherry) chuyển sang sớm 5,5 ngày 25 năm qua biến động nhiệt độ vào tháng tháng Tác giả [4] mối tương quan nghịch đáng kể tốc độ hoa họ cam quýt với thay đổi nhiệt độ tối đa, tối thiểu số địa điểm Iran, điều cho thấy nhiệt độ tăng làm chậm trình hoa, dẫn đến thời gian hoa kéo dài Ở Junagadh, Ấn Độ, phần lớn giống xoài hoa bị chậm lại nhiệt độ vượt nhiệt độ Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 http://tapchikttv.vn/ Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 tối ưu chúng Năng suất xoài, chuối, đu đủ hồng xiêm phát có mối tương quan nghịch với nhiệt độ lượng mưa khoảng thời gian 10 năm quận Navasari, Ấn Độ [5] Trong nghiên cứu đánh giá tác động biến đổi khí hậu sản xuất chuối toàn cầu, số kết luận rút ra; chẳng hạn khu vực có nhiệt độ trung bình ≥ 24°C, chuối sinh trưởng phát triển thuận lợi Ở Ấn Độ Argentina xu nhiệt độ cao làm giảm phát triển chuối [6] Nhiệt độ 12-13°C lượng nước sẵn có đất thấp so với trung bình làm chậm trình trao đổi chất vào thời kỳ “qua đông” Nhiệt độ cực đoan (trên 35°C) khiến 50% nho bị rụng bị thâm cháy [7] Điều thấy vải, táo bị nám nắng, nứt chúng chịu ảnh hưởng nhiệt độ cao Kích thước vải bị giảm thiếu nước, nhiệt độ tăng dẫn đến phát triển nhanh tích lũy đơn vị nhiệt nhanh [8] Tác động thời tiết khí hậu sản xuất nơng nghiệp thể nhiều mặt, ảnh hưởng đến sinh trưởng, phát triển trồng mà ảnh hưởng đến hoạt động sản xuất, hiệu suất sản xuất Ngoài ra, thiên tai hàng năm làm thiệt hại nghiêm trọng đến trồng Sự hoạt động loại côn trùng có hại, phát triển lây lan sâu bệnh có liên quan đến điều kiện khí tượng,…Vì vậy, đánh giá ảnh hưởng yếu tố KT đến suất trồng cần thiết, sở khoa học để xác định biện pháp kỹ thuật phù hợp cho gia tăng suất ăn quả, thích ứng với biến đổi khí hậu, phục vụ phát triển nông nghiệp bền vững Tuy nhiên, năm qua, cơng trình nước chủ yếu đánh giá ảnh hưởng điều kiện khí tượng nơng nghiệp đến lương thực thực phẩm [1, 9–11], chưa nhiều cơng trình nghiên cứu đánh giá cho ăn công bố [12–13] Một số nghiên cứu ảnh hưởng yếu tố KT trồng dựa sử dụng mơ hình trồng để đánh giá [9–10], dựa ngưỡng sinh học tối ưu theo giai đoạn sinh trưởng-phát triển hình thành suất trồng, đánh giá tổn thất/thiệt hại trồng điều kiện khí tượng cực đoan gây [1, 11] Bên cạnh đó, số nghiên cứu đánh giá dựa mối quan hệ yếu tố KT trồng sử dụng hồi quy tuyến tính để đánh giá tác động khí hậu tồn cầu đến đặc trưng khí hậu tiểu vùng khí hậu, tác động biến đổi khí hậu đến suất trồng [1, 9, 10–15] Các nghiên cứu thường tách suất trồng thành hai thành phần: (i) Năng suất xu (NSXT) (ii) suất thời tiết, khí hậu (NSTT) Trong NSXT tác động khoa học kỹ thuật (điều kiện giống, công nghệ, phân bón, thuốc trừ sâu, quản lý, canh tác,… ngày tốt hơn) Do đó, NSXT ngày cao điều kiện khoa học kỹ thuật ngày hồn thiện phát triển Những điều kiện khí tượng xác định cho nông nghiệp trước hết ánh sáng, nhiệt độ nước Để trồng loại cần phải bảo đảm lượng nhiệt định, hay trồng tích luỹ tổng lượng nhiệt định, phát dục chuyển trạng thái sinh trưởng, đồng thời tùy vào loại trồng khác nhau, chúng cần nhu cầu nước khác cho sinh trưởng phát triển Bên cạnh đó, yếu tố khí tượng thấp cao vượt giới hạn tối ưu gây ảnh hưởng đến suất trồng, ví dụ số ngày nắng nóng cao vào thời trổ bơng-nở hoa làm cháy phấn, dẫn đến suất trồng giảm Vì vậy, tác động yếu tố thời tiết, khí hậu đến suất trồng gọi suất thời tiết, khí hậu (NSTT) [1, 11–13, 15] Trong năm gần đây, Việt Nam ăn phát triển nhiều vùng Sơn La, Hưng Yên,…phục vụ nhu cầu tiêu thụ nước, xuất từ sản phẩm, chế phẩm ăn quả, giúp chuyển đổi cấu trồng Tỉnh Phú Thọ đẩy mạnh phát triển ăn theo hướng hiệu kinh tế cao Nhiều mặt hàng nông sản cấp mã truy xuất nguồn gốc, chuỗi giá trị nông sản bưởi Đoan Hùng, hay chè Phú Thọ Tuy nhiên chưa có cơng trình nghiên cứu đánh giá ảnh hưởng yếu tố KT đến sản xuất ăn Phú Thọ công bố Trong nghiên cứu chủ yếu khảo sát ảnh hưởng điều kiện nhiệt, mưa số yếu tố khí hậu cực đoan suất bưởi Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 chuối trung bình tỉnh Phú Thọ dựa mối tương quan hồi qui [15] Kết đánh giá nhằm xác định yếu tố KT ảnh hưởng chủ đạo vào giai đoạn nhạy cảm bưởi chuối để có biện pháp kỹ thuật phù hợp hạn chế ảnh hưởng chúng Bên cạnh đó, nhằm mục đích xác định sở khoa học hướng tới đánh giá tác động biến đổi khí hậu, phân vùng chuyên canh ăn quả, dự báo suất trồng chủ lực tỉnh Phú Thọ (bưởi chuối) phục vụ xây dựng quy hoạch sản xuất, kế hoạch tiêu thụ xuất Số liệu sử dụng phương pháp nghiên cứu 2.1 Số liệu (1) Số liệu quan trắc khí tượng trạm Minh Đài, Phú Hộ Việt Trì, thời kỳ 20102021 nghiên cứu sử dụng cho đánh giá ảnh hưởng yếu tố KT đến suất thu hoạch ăn (Hình 1) (2) Số liệu thống kê sản xuất ăn suất thu hoạch (năng suất thực thu bưởi chuối), thời kỳ 2010-2021 Năng suất tính diện tích cho sản phẩm, khơng phân biệt năm có cho sản phẩm hay khơng (diện tích đưa vào sản xuất kinh doanh), khơng tính diện tích trồng, nhỏ trắng [16–17] Số liệu sử dụng để đánh giá ảnh hưởng yếu tố KT đến suất ăn Vì niên giám thống kê đưa suất tỉnh, nên việc sử dụng số liệu trạm KT cách riêng biệt chưa phù Hình Sơ đồ trạm khí tượng Phú Thọ hợp Cũng lựa chọn trạm khí tượng gần với vùng diện tích trồng nhiều để đánh giá Tuy nhiên, để đồng đánh giá, nghiên cứu trung bình số học yếu tố KT Bảng từ số liệu KT trạm để đánh giá ảnh hưởng yếu tố KT đến suất ăn [15] Bảng Các yếu tố khí tượng sử dụng từ 2010-2021 TT 10 11 12 13 14 15 16 17 18 19 Mô tả Lượng mưa tháng (mm) Lượng mưa ngày lớn tháng (mm) Số ngày mưa tháng Số ngày có mưa 10mm tháng Số ngày có mưa 50mm tháng Số ngày có mưa 100mm tháng Nhiệt độ trung bình tháng (℃) Nhiệt độ tối cao TB tháng (℃) Nhiệt độ tối thấp TB tháng (℃) Nhiệt độ tối cao tuyệt đối tháng (℃) Nhiệt độ tối thấp tuyệt đối tháng (℃) Số nắng tháng (giờ) Số ngày rét đậm tháng Số ngày rét hại tháng Độ ẩm tối thấp trung bình tháng (%) Độ ẩm trung bình tháng (%) Tốc độ gió trung bình tháng (m/s) Trung bình tháng tốc độ gió lớn ngày (m/s) Tổng số obs xuất hướng NE tháng Ký hiệu R Rmax R0,1 R10 R50 R100 Ttb Tx Tn Txtđ Tntđ Hn RĐ RH Umin Utb fftb ffmax NE Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 TT 20 21 22 23 Mô tả Tổng số obs xuất hướng SE theo tháng Tổng số obs xuất hướng SW theo tháng Tổng số obs xuất hướng NW theo tháng Số ngày nhiệt độ đạt nắng nóng tháng (Tx ≥ 35℃) Ký hiệu SE SW NW Tx 35 2.2 Phương pháp (1) Phương pháp chuẩn hóa số liệu sử dụng để chuẩn hóa tất đặc trưng khí tượng NSTT [15]: xt -x̅ zt = (1) sx Trong zt chuẩn hóa yếu tố khí hậu; xt, x̅ , sx thành phần chuỗi, giá trị trung bình chuỗi độ lệch tiêu chuẩn đặc trưng khí tượng (2) Phương pháp tính tốn suất thời tiết khí hậu tạo nên (NSTT) xác định dựa suất thực tế suất xu (NSXT) [1, 11–13, 15] Để đánh giá ảnh hưởng yếu tố KT đến suất trồng, nghiên cứu thường tách chuỗi suất trồng trung bình tỉnh thành hai phần: (i) Năng suất xu (ii) suất thời tiết (Hình 2) Hình Sơ đồ mơ tả thành phần suất trồng [18] Có thể biểu diễn mối quan hệ ba đặc trưng suất sau [1, 11–13, 15]: YTT = YT - YXT (2) Trong đó: YTT suất thời tiết (tấn/ha); YT suất thực (tấn/ha) suất chuối bưởi; YXT suất xu (tấn/ha) - NSXT xác định dựa hàm đó, tuyến tính hàm logistic Các hàm nghiên cứu xác định dựa khảo sát ước lượng Excel với hệ số tương quan cao lựa chọn [15, 18–20] - NSTT (YTT) suất tác động thời tiết, khí hậu tạo nên [15, 18–20] (3) Phương pháp tương quan hồi qui tuyến tính với kỹ thuật bước sử dụng để xác định ảnh hưởng yếu tố khí tượng đến suất ăn dựa mối quan hệ tuyến tính (cụ thể phương pháp trình bày [21]): - Hệ số tương quan xác định: n ∑ n i=1 (xi -𝑥̅ )(yi -y̅) r= (3) 1 2 n n √ ∑i=1 (xi -𝑥̅ ) √ ∑i (y𝑖 -y̅ ) n n Trong yi NSTT xi yếu tố KT năm i (i = 1, ,12); 𝑦̅ 𝑥̅ giá trị trung bình yi hay xi , n số năm quan trắc 12 năm Độ tin cậy hệ số tương quan đánh giá dựa theo kiểm nghiệm Student: Hệ số tương quan với dung lượng mẫu (n) lớn đáng kể thỏa mãn tiêu chuẩn tương ứng với α = 5% 10% Dung lượng mẫu 12 năm bậc tự (n–2 = 10), t(α=0.05, 10) = 2,23 (t(α=0.1, 10) = 1,81), tương đương với rij= 0,58 (0,50), đạt độ tin cậy 95% (90%) Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 - Phương trình hồi qui xây dựng dựa kỹ thuật bước: (1) Tính hệ số tương quan toàn phần r1i Y Xi (i = 2…m biến), xây dựng mơ hình với X1 xác định biến sai thặng dư S(1); (2) Tính hệ số tương quan riêng r1i.2 (i = 3,…,m) chọn hệ số tương quan có giá trị lớn nhất, sau xây dựng mơ hình hồi qui với biến x1, x2 xác định biến sai thặng dư S(2); (3) So sánh giá trị biến sai thặng dư S(2) với S(1) xác định 𝑠 (2) −𝑠(1) (| 𝑠(2) | = 𝜀) Nếu ≤ 0 (giá trị tùy chọn cho trước) x2 bị loại ngược lại ≥ 0 x2 chọn cho mơ hình Qui trình tiếp tục bước để xác định biến lại Trong nghiên cứu sử dụng phần mềm thống kê SPSS (Statistical Package for the Social Sciences) để xây dựng mơ hình ước lượng suất bưởi chuối với biến phụ thuộc chuẩn sai NSTT biến độc lập đặc trưng Bảng Chất lượng phương trình đánh giá theo kiểm nghiệm Fisher [21] 2.3 Hiện trạng sản xuất ăn Phú Thọ Hình cho thấy chuối bưởi có diện tích cho thu hoạch đưa vào kinh doanh sản lượng cao tỉnh Phú Thọ Hiện nay, bưởi trồng chủ lực tỉnh việc thực chuyển đổi cấu trồng Tính đến năm 2021 diện tích cho thu hoạch bưởi 3,7 nghìn ha, sản lượng gần 32,7 nghìn tấn, tăng 800 sản lượng tăng 10,5 nghìn so với năm 2019 Trong huyện Đoan Hùng chiếm diện tích sản lượng bưởi nhiều với 1,7 nghìn khoảng 24,2 nghìn năm 2010 Từ năm 2019-2021 diện tích chuối tăng 50 sản lượng tăng 3,4 nghìn tấn, tăng chậm so với bưởi Hình Diện tích cho thu hoạch đưa vào kinh doanh (a) sản lượng (b) loại ăn tỉnh Phú Thọ giai đoạn 2019-2021 Hình thấy bưởi trồng nhiều Đoan Hùng với diện tích 1500ha với sản lượng đạt 24 nghìn năm 2020 Trong huyện khác, bưởi trồng với diện tích thấp 300ha (thấp thị xã Phú Thọ) với sản lượng khoảng 1-2 nghìn Chuối nhiều huyện Thanh Sơn với diện tích cho thu hoạch 700 với sản lượng lên đến 21,6 nghìn năm 2020, huyện Hạ Hịa Cẩm Khê trồng với diện tích 400ha, huyện khác diện tích trồng thấp 300ha Đáng ý Lâm Thao, diện tích cho thu hoạch khoảng 250ha sản lượng cao nghìn Hình Diện tích bưởi cho thu hoạch đưa vào kinh doanh sản lượng bưởi Phú Thọ năm 2020 Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 Có thể nhận thấy bưởi, chuối trồng với diện tích nhiều nhì tỉnh Phú Thọ Hai chủ lực, kế hoạch phát triển chủ lực tỉnh Phú Thọ Do nghiên cứu tập trung đánh giá ảnh hưởng yếu tố KT suất bưởi chuối tỉnh Phú Thọ Kết 3.1 Năng suất thực thu suất xu Năng suất thu hoạch (thực thu) NSXT bưởi giai đoạn 2010-2021 có xu hướng tăng Năm 2010 suất bưởi thu khoảng tấn/ha đến năm 2021 khoảng 12 tấn/ha Năng suất bưởi biến động cao đáng kể giai đoạn 2010-2014, năm gần biến động Năng suất chuối có xu hướng tăng dần, năm 2010 thu khoảng 21 tấn/ha, đến năm 2021 khoảng 25 tấn/ha, tăng tấn/ha Có thể nhận thấy suất chuối biến động, điều cho thấy chuối chịu tác động thời tiết, khí hậu (Hình 5a, 5b) Hình 5a, 5b nhận thấy suất bưởi chuối (thực thu) có biến động cao, thấp so với NSXT, biến động cao thấp cho thấy tác động thời tiết, khí hậu tạo nên Sử dụng công thức năm NSTT dương (thời tiết, khí hậu thuận lợi) năm âm (thời tiết, khí hậu bất lợi) thể Hình 5b, 5c Kết cho thấy NSTT bưởi Phú Thọ cao năm 2011, 2013, 2019 dao động khoảng từ 0,01 đến 1.06 tấn/ha, điều cho thấy năm thời tiết, khí hậu thuận lợi suất bưởi Một số năm thấp năm 2012, 2018 với giá trị khoảng từ 0,08-2,4 tấn/ha, cho thấy thời tiết, khí hậu bất lợi suất bưởi Đối với chuối, năm NSTT cao 2011, 2012, 2013, 2016 với giá trị dao động từ 0.01-0.2 tấn/ha, năm thấp năm 2010, 2015, 2017 với giá trị từ 0,05-0,4 tấn/ha (Hình 5c, 5d) Nhìn chung, ảnh hưởng yếu tố thời tiết, khí hậu đến suất bưởi rõ ràng so với chuối Hình Năng suất thu hoạch (thực thu) suất xu bưởi (a) chuối (b) Năng suất thời tiết bưởi (c) chuối (d) giai đoạn 2010-2021 3.2 Ảnh hưởng yếu tố khí tượng đến suất bưởi chuối Hệ số tương quan suất bưởi chuối với yếu tố KT theo 12 tháng Bảng 1, thời kỳ 2010-2021 thể Bảng 3, giá trị màu xanh hệ Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 số tương quan thuận, màu cam thể hệ số tương quan nghịch đạt độ tin cậy 90% Hệ số tương quan 0,50 (0,58) -0,50 (-0,58) đạt độ tin cậy 90% (95%) Có thể thấy vào tháng số ngày có lượng mưa 10mm, nhiệt độ trung bình, nhiệt độ tối cao tháng có tác động tích cực đến suất bưởi Ở Phú Thọ, tháng có nhiệt độ thấp lượng mưa ít, thời kì bưởi hoa làm quả, số ngày có lượng mưa 10mm nhiệt độ cao thuận lợi cho bưởi sinh trưởng phát triển Vào thời kì rét đậm, rét hại tần suất gió đơng bắc cao lại điều kiện bất lợi đến suất bưởi, rét đậm, rét hại làm hoa bưởi không phát triển, phát triển, dẫn đến tỷ lệ đậu thấp Tháng 4-6, nhiệt độ tối cao tuyệt đối, số ngày nắng nóng cao tốc độ gió mạnh tương quan nghịch với suất bưởi, điều cho thấy yếu tố làm giảm suất bưởi Bởi thời kỳ nhỏ, gió lớn dễ làm rụng quả, nhiệt độ cao số ngày nắng nóng nhiều, dẫn đến bị khô, cháy nám chậm phát triển, thịt phát triển chậm Tháng 10 đến tháng 12, thời kỳ bưởi chín, thời gian nắng cao thuận lợi cho suất bưởi Tương tự, điều kiện nhiệt độ thấp tuyệt đối cho thấy ảnh hưởng bất lợi đến suất bưởi, nhiên số ngày rét hại lại thuận lợi suất bưởi Điều nhiệt độ tối thấp chủ yếu có vai trò quan trọng liên quan đến khả qua đơng ăn quả, số ngày rét 150C cao có lợi cho việc kéo dài thời gian để bưởi tập trung dưỡng chất vào Nhìn chung, bưởi Phú Thọ, giai đoạn hoa, nhỏ (tháng 2-4) chịu ảnh hưởng rét đậm, rét hại gió đơng bắc Do đó, cần biệt pháp kỹ thuật phủ gốc để giữ ấm (sử dụng rơm dạ, cỏ khơ, nilong,…) Bên cạnh đó, nên trồng chắn gió đơng bắc cho vườn ăn nhằm hạn chế rụng hoa non Ở giai đoạn bưởi chín thu hoạch, độ ẩm cao dễ sinh nấm sâu bệnh hại Do đó, nên có biện pháp kỹ thuật làm giảm độ ẩm khu vực vườn ăn hun khói, trồng che bóng cản nhiệt độ hạ thấp, hay tưới nước để giảm sương muối Bảng Hệ số tương quan tuyến tính yếu tố KT đến suất bưởi Tháng/ Yếu tố Rn Rmax R0.1 R10 R50 R100 Ttb Tx Tn Txtđ Tntđ Hn RĐ RH Umin Utb fftb ffmax NE SE SW NW Tx35 10 11 12 -0,32 -0,24 -0,18 -0,32 -0,01 -0,11 -0,03 -0,02 -0,27 0,12 -0,10 -0,32 0,15 -0,11 -0,17 -0,06 -0,42 -0,14 -0,16 -0,22 -0,03 - 0,21 0,33 -0,19 0,57 0,10 0,52 0,57 0,48 0,21 0,44 0,25 -0,50 -0,61 -0,03 -0,04 -0,03 0,28 -0,57 0,07 0,23 -0,21 - 0,18 0,04 0,10 0,05 0,08 0,21 0,26 0,26 0,09 0,15 0,18 -0,49 0,18 -0,02 0,06 -0,05 -0,04 0,06 -0,12 0,37 -0,21 - 0,05 -0,04 0,09 -0,10 0,16 -0,20 -0,29 -0,08 -0,43 -0,06 -0,37 0,16 0,32 0,46 -0,05 -0,40 -0,60 -0,11 0,12 0,26 -0,71 -0,14 0,31 0,04 -0,10 -0,01 -0,01 -0,33 -0,17 -0,20 0,10 -0,47 -0,28 0,03 0,13 -0,07 -0,66 0,46 -0,30 -0,07 -0,56 -0,11 0,52 0,41 0,16 0,54 0,48 0,41 -0,24 0,02 -0,22 -0,15 -0,64 0,47 -0,02 0,05 -0,04 0,16 0,54 0,03 -0,07 0,19 0,19 -0,17 -0,15 0,07 -0,02 -0,41 -0,13 -0,22 -0,02 -0,03 -0,37 -0,11 -0,11 0,05 0,10 -0,06 -0,30 -0,11 -0,01 0,33 0,16 -0,23 -0,41 -0,31 0,10 -0,15 -0,44 -0,38 -0,20 -0,07 0,20 0,02 -0,02 -0,13 0,04 0,13 -0,05 -0,22 -0,29 -0,09 -0,34 -0,46 0,01 0,25 0,00 0,02 0,32 0,10 -0,14 -0,22 -0,02 -0,06 0,49 0,01 0,10 0,28 0,44 -0,06 -0,24 -0,17 -0,09 0,50 0,09 -0,07 0,14 0,16 0,17 0,14 0,07 -0,15 -0,37 -0,35 -0,46 0,40 -0,43 0,10 0,16 0,02 0,13 -0,07 -0,28 0,22 -0,63 0,11 -0,08 - -0,21 -0,17 -0,48 -0,08 -0,04 -0,20 -0,10 -0,01 -0,47 0,14 0,10 0,10 0,16 -0,70 -0,62 -0,07 -0,38 -0,23 -0,35 -0,16 -0,53 - -0,18 -0,15 -0,62 -0,20 0,18 -0,49 -0,10 -0,67 -0,27 -0,54 0,64 0,47 0,53 -0,70 -0,59 -0,05 0,03 0,20 -0,23 0,29 0,00 - Nhìn vào Bảng thấy, yếu tố độ ẩm tác động đến suất chuối xuyên suốt 12 tháng Chuối có đặc tính ưa ẩm, chống chịu hạn kém, nên độ ẩm cao thuận lợi cho sinh trưởng, phát triển chuối Bên cạnh đó, hệ số tương quan tốc độ gió suất chuối có tương quan âm 12 tháng, cho thấy ảnh hưởng tiêu cực đến Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 suất chuối tỉnh Phú Thọ Tác hại gió làm rách phiến, mép lá, vết rách hóa sẹo làm giảm diện tích hoạt động lá, làm giảm suất đến 20% so với khơng bị rách Gió mạnh làm gãy, đổ chuối gây thiệt hại nghiêm trọng Tháng 1-3 giai đoạn trồng chuối, vào thời kỳ này, lượng mưa 50mm, nhiệt độ cao, số nắng số ngày nắng nóng nhiều (nhiệt độ tối cao ≥ 35℃) tác động tiêu cực đến sinh trưởng phát triển chuối, giai đoạn trồng nhạy cảm, cần nhiệt độ vừa đủ mưa để thuận lợi cho việc phát triển chồi Bên cạnh đó, giai đoạn thời điểm thu hoạch chuối Phú Thọ Do đó, số ngày nắng nóng, mưa lớn nhiều làm chín nhanh dễ dập nát làm ảnh hưởng đến chất lượng Tháng lượng mưa có tương quan dương với suất chuối, cho thấy lượng mưa tăng tác động tích cực đến suất chuối, giai đoạn trưởng thành chuối, chuối cần nhu cầu nước cao Vào tháng tháng 7, nhiệt độ cao làm cho suất chuối giảm, hàm lượng nước phận chuối cao phiến rộng, số ngày nắng nóng cao làm bốc thoát cao, gây bất lợi đến suất chuối Đến tháng 10 tháng 11 nhiệt độ thấp lại tạo điều kiện thuận lợi cho chuối phát triển Nhìn chung, chuối Phú Thọ, tốc độ gió mạnh yếu tố ảnh hưởng đến suất chuối điển hình, nên có biện pháp trồng chắn gió để hạn chế ảnh hưởng đến sinh trưởng, phát triển chuối Bảng Hệ số tương quan tuyến tính yếu tố KT đến suất chuối Tháng/ Yếu tố Rn Rmax R0.1 R10 R50 R100 Ttb Tx Tn Txtđ Tntđ Hn RĐ RH Umin Utb fftb ffmax NE SE SW NW Tx35 10 11 12 -0,11 -0,27 0,34 -0,13 -0,20 -0,44 -0,61 -0,44 -0,44 -0,28 -0,29 -0,05 0,52 0,08 0,77 -0,78 -0,07 -0,11 -0,30 0,58 -0,13 - 0,34 -0,05 0,55 0,22 -0,16 -0,48 -0,62 -0,50 -0,68 0,10 -0,56 0,23 0,14 -0,08 0,55 -0,78 0,12 -0,13 -0,31 0,33 -0,14 -0,78 0,02 -0,31 0,04 0,02 -0,64 -0,25 -0,31 -0,29 -0,46 -0,12 -0,24 0,25 0,34 -0,39 0,63 -0,77 0,29 -0,39 -0,30 0,42 0,00 -0,08 0,00 0,18 -0,53 0,11 -0,18 0,30 0,22 0,29 0,00 0,13 0,09 -0,08 -0,46 0,46 -0,77 0,42 -0,15 0,19 0,41 -0,12 0,15 0,52 0,26 0,54 0,42 0,45 -0,07 -0,09 -0,29 -0,40 -0,24 -0,38 -0,01 -0,40 0,74 -0,78 0,34 -0,39 0,01 0,46 -0,29 -0,06 0,18 0,11 0,28 0,09 0,30 0,20 0,03 -0,33 -0,30 -0,56 -0,16 0,00 -0,64 0,55 -0,77 0,38 0,29 0,16 0,05 -0,21 -0,28 0,09 0,36 0,23 0,11 0,08 0,57 -0,08 -0,46 -0,37 -0,66 -0,08 -0,38 -0,07 0,78 -0,75 0,58 0,12 0,27 0,21 -0,27 -0,51 0,22 0,21 -0,35 0,18 0,22 0,02 0,26 0,21 -0,21 0,33 -0,33 0,08 -0,47 0,61 -0,78 0,28 -0,43 0,08 0,02 0,12 0,17 0,24 0,12 0,46 0,28 0,01 0,20 -0,16 -0,58 0,06 -0,16 -0,32 -0,48 0,75 -0,77 0,23 -0,10 -0,33 0,09 0,22 -0,33 -0,16 -0,16 0,04 -0,02 -0,29 -0,08 0,25 0,20 0,45 0,01 0,66 -0,08 -0,50 0,69 -0,77 0,31 -0,17 0,17 -0,03 -0,12 -0,16 -0,05 0,02 0,11 -0,02 -0,22 0,41 0,28 0,87 0,36 0,74 -0,25 -0,08 -0,51 0,62 -0,78 -0,14 -0,37 -0,15 0,18 -0,16 - -0,18 -0,04 -0,42 -0,10 0,34 -0,09 -0,06 -0,32 -0,14 -0,27 -0,17 0,20 -0,59 0,55 -0,76 0,43 0,18 0,15 0,04 0,15 - Để ước lượng suất bưởi chuối, nghiên cứu dựa phần mềm SPSS để xây dựng mơ hình hồi qui với kỹ thuật bước biến phụ thuộc chuẩn sai NSTT biến độc lập chuẩn sai yếu tố KT Kết lựa chọn hai mơ hình hồi qui thể Bảng Đánh giá chất lượng mơ hình hồi qui dựa kiểm nghiệm Fisher cho thấy hai mơ hình hồi qui ước lượng NSTT bưởi chuối dụng được, từ ước lượng suất thực thu bưởi chuối dựa công thức Kết so sánh NSTT bưởi chuối quan trắc ước lượng dựa hồi qui từ chuỗi phụ thuộc cho thấy chúng dao động đồng pha chênh lệch (Hình 5) Đối với bưởi: từ mơ hình hồi qui cho thấy số ngày nắng nóng trung bình tháng 4, độ ẩm tối thấp ngày trung bình tháng 11, tốc độ gió trung bình tháng 5, số obs xuất hướng gió tây bắc trung bình tháng ảnh hưởng bất lợi NSTT bưởi Khi đặc trưng Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 cao, NSTT bưởi giảm Mức độ ảnh hưởng bốn yếu tố NSTT bưởi mơ hình hồi qui gần tương tự (a) (b) Hình Năng suất thời tiết suất tính tốn dựa hồi qui (a) bưởi, (b) chuối giai đoạn 2010-2021 Đối với chuối: ba yếu tố KT từ mơ hình hồi qui, mức độ ảnh hưởng nhiệt độ tối thấp trung bình tháng 11 đến NSTT chuối cao (hệ số hồi qui 0,678); nhiệt độ tối thấp trung bình tháng 11 cao hơn, thuận lợi cho NSTT chuối, ngược lại Tiếp đến nhiệt độ tối cao ngày trung bình tháng số ngày có lượng mưa 50 mm tháng 3; hai đặc trưng cao, NSTT thấp ngược lại Bảng Mơ hình hồi qui ước lượng NSTT bưởi chuối STT Cây Bưởi Chuối Phương trình Y = -0.357T35_4- 0.453Umin11 0.315ffmax5 - 0.231NW5 Y = 0.678Tn11 - 0.338Tx1 0.223R50_3 Hệ số (f) Hệ số (fα=0.05) 16,526 4,12 39,965 4,07 Đánh giá f > fα : Phương trình sử dụng f > fα : Phương trình sử dụng Kết luận Nghiên cứu đánh giá ảnh hưởng yếu tố KT dựa phương pháp tương quan hồi qui tuyến tính có số nhận xét sau: - Mặc dù hệ số tương quan yếu tố KT đến suất bưởi chuối phổ biến thấp nhìn chung rét đậm, rét hại số lần xuất gió hướng bắc ảnh hưởng đáng kể đến bưởi giai đoạn hoa, non Nhiệt độ thấp, độ ẩm cao có ảnh hưởng đáng kể đến giai đoạn chín Do đó, nên áp dụng biệt pháp kỹ thuật phù hợp để hạn chế phủ gốc để giữ ấm, trồng chắn gió hạn chế rụng hoa non, hay biện pháp kỹ thuật làm giảm độ ẩm giai đoạn chín - Nhiệt độ cao số nắng, số ngày nắng nóng nhiều tác động tiêu cực đến giai đoạn chuối phát triển chồi Độ ẩm khơng khí cao năm ảnh hưởng điển hình đến suất chuối Do nên có biện pháp kỹ thuật phù hợp để giảm ảnh hưởng nhiệt độ cao vào giai đoạn đâm chồi, hạn chế độ ẩm vào giai đoạn chuối tưới nước, trồng che bóng, hun khói,… - Đã xây dựng hai mơ hình ước lượng NSTT bưởi chuối cho tỉnh Phú Thọ Các phương trình sử dụng cho ước lượng suất thực thu bưởi chuối, sử dụng cho đánh giá tác động biến đổi khí hậu đến suất bưởi chuối số liệu tương lai có đầy đủ yếu tố KT mơ hình, hướng tới xây dựng mơ hình dự báo phục vụ lập kế hoạch tiêu thụ xuất Đóng góp tác giả: Xây dựng ý tưởng nghiên cứu: N.Đ.M., T.H.D.; Lựa chọn phương pháp nghiên cứu: T.T.T., T.H.D.; Xử lý số liệu: L.V.P., Đ.N.Ý.; Tính tốn kết quả: N.B.P., Đ.N.Ý, L.V.P.; Viết thảo báo: T.H.D., T.T.T.; Chỉnh sửa báo: N.H.Q Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 10 Lời cảm ơn: Bài báo hoàn thành nhờ trợ giúp từ hai đề tài KHCN cấp “Nghiên cứu đổi cơng nghệ dự báo khí hậu nơng nghiệp Áp dụng cho dự báo điều kiện khí hậu nông nghiệp tác động đến sản xuất lúa khu vực đồng sông Hồng”, mã số TNMT.2022.02.15 Bộ “Nghiên cứu đổi công nghệ dự báo đợt rét đậm, rét hại, mưa lớn mùa đông khu vực phía Bắc Việt Nam hạn đến 10 ngày dựa quy luật hoạt động dòng xiết cận nhiệt đới”, mã số TNMT.2022.06.08 Viện Khoa học Khí tượng Thủy văn Biến đổi khí hậu chủ trì thực Lời cam đoan: Tập thể tác giả cam đoan báo cơng trình nghiên cứu tập thể tác giả, chưa công bố đâu, không chép từ nghiên cứu trước đây; khơng có tranh chấp lợi ích nhóm tác giả Tài liệu tham khảo Viết, N.V Tài ngun khí hậu nơng nghiệp Việt Nam Nhà xuất nông nghiệp 2009 Grab, S.; Craparo, A Advance of apple and pear tree full bloom dates in response to climate change in the southwestern Cape, South Africa 1973–2009 Agric For Meteorol 2011, 151(3), 406–413 El Yaacoubi, A.; El Jaouhari, N.; Bourioug, M.; El Youssfi, L.; Cherroud, S.; Bouabid, R.; Abouabdillah, A, Potential vulnerability of Moroccan apple orchard to climate change– induced phenological perturbations: effects on yields and fruit quality Int J Biometeorol 2020, 64(3), 377–387 Fitchett, J.M.; Grab, S.W.; Thompson, D.I.; Roshan, G Spatio-temporal variation in phenological response of citrus to climate change in Iran Agric For Meteorol 2014, 198, 285–293 Mali, S.C.; Swati, S.; Shrivastava, P.K Economic Evaluation of Sugarcane Based Agro forestry Systems J Tree Sci 2017, 36(1), 34–37 Calberto, G.; Staver, C.; Siles, P.; An assessment of global banana production and suitability under climate change scenarios In: Aziz E (ed) Climate change and food systems: global assessments and implications for food security and trade Rome: Food and Agriculture Organization of the United Nations (FAO) 2015, pp 262–291 Kizildeniz, T.; Pascual, I.; Irigoyen, J.J.; Morales, F Using fruit-bearing cuttings of grapevine and temperature gradient greenhouses to evaluate effects of climate change (elevated CO2 and temperature, and water deficit) on the cv red and white Tempranillo Yield and must quality in three consecutive growing seasons (2013– 2015) Agric Water Manage 2018, 202, 299–310 Nath, V.; Kumar, G.; Pandey, S.D.; Pandey, S Impact of climate change on tropical fruit production systems and its mitigation strategies Proceeding of the Climate change and agriculture in India: Impact and adaptation Cham: Springer 2019, pp 129–146 Dương, T.H cs Nghiên cứu ứng dụng mơ hình ORYZA2000 để đánh giá rủi ro khí hậu nơng nghiệp đề xuất phương án quản lý sản xuất lúa thích hợp phục vụ phát triển nơng nghiệp bền vững đồng sông Hồng Đề tài cấp Cơ sở Viện KHKTTV BĐKH 2014 10 Giang, N.T Nghiên cứu ứng dụng mơ hình động thái để xác định công thức luân canh trồng vùng phù sa đồng sông Cửu Long Luận án tiến nông nghiệp, Trường Đại học Nông nghiệp 2012 11 Hà, N.T cs Nghiên cứu dự báo suất ngô, đậu tương, lạc xây dựng công nghệ giám sát KTNNở Việt Nam Đề tài cấp Bộ 2008 12 Giai, N.S.; Quyền, N.H Về khả sử dụng số liệu quan trắc vải ngành nông nghiệp nghiên cứu yếu tố KT hình thành suất Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 2015 Tạp chí Khí tượng Thủy văn 2023, 752, 1-11; doi:10.36335/VNJHM.2023(752).1-11 11 13 Quyền, N.H Nghiên cứu đánh giá phân vùng khí hậu nơng nghiệp phục vụ phát triển ngành trồng trọt tỉnh Sơn La Luận án tiến sĩ khí tượng khí hậu học, Viện Khoa học Khí tượng Thuỷ văn Biến đổi khí hậu 2021, tr 155 14 Hường, C.T.T Nghiên cứu tác động biến đổi khí hậu tồn cầu đến số cực trị khí hậu tượng khí hậu cực đoan Việt Nam Luận án Tiến sĩ Khoa học trái đất 2015, tr 158 15 Han, X.; Chang, L.; Wang, N.; Kong, W.; Wang, C Effects of Meteorological Factors on Apple Yield Based on Multilinear Regression Analysis: A Case Study of Yantai Area, China Atmosphere 2023, 14, 183 16 Trực tuyến: https://www.gso.gov.vn/du-lieu-dac-ta/2019/12/htcttk-cap-tinh-nangsuat-mot-so-loai-cay-trong-chu-yeu/ 17 Trực tuyến: https://www.gso.gov.vn/du-lieu-dac-ta/2019/03/nong-nghiep/ 18 Hung, N.T.; Quyen, N.H Effect of the change of climate indicators on agricultural yields in Son La province VN J Sci Technol 2017, 6, 756–766 19 Irwin, S.H Revisiting South American Corn and Soybean Yield Trends and Risks Farmdoc Daily 2020, 10, 212 20 Wang, Y.; Yang, T.; Zhao, R.; Jiang, K.; Ma, X A New Method for Determining Critical Irrigation Period for Large Regions Based on Precipitation-Meteorological Yield Integral Regression Relationship—A Case Study of Winter Wheat in Shaanxi Province, China Appl Sci 2019, 9(23), 5181 21 Tân, P.V Phương pháp thống kê khí hậu Nhà xuất Đại học Quốc gia Hà Nội 2005 Effects of memotological factors on yield of some fruit trees in Phu Tho province basis on correlation and regression Nguyen Dang Mau1*, Trinh Hoang Duong1, Tran Thi Tam1, Nguyen Huu Quyen1, Nguyen Binh Phong2, Dang Nhu Y 2, Le Van Phong1 Vietnam Institute of Meteorology, Hydrology and Climate Change; mau.imhen@gmail.com; hoangduongktnn@gmail.com; trantam1810@gmail.com; nhquyen13@gmail.com; lephong2341999@gmail.com Hanoi University of Natural Resources and Environment; nbphong@hunre.edu.vn Abstract: The article analyzes the influence of meteorological factors (MF) on the yield of pomelo and banana in Phu Tho based on linear correlation and regression to determine the appropriate technical measures to increase yield of them The results showed that for pomelos, the high number of cold days and northeasterly winds were typical factors affecting the flowering stage and the small fruit High humidity is also a favorable factor for the activity of fungi and pests at the ripe grapefruit stage For bananas, strong wind speed is a typical influencing factor in all stages of growth and yield formation Therefore, it should be applied to reduce the influence of these factors on the flowering and ripening stage of pomelo and banana such as keeping the stump warm, or planting windbreaks, etc In addition, we have been built two models to estimate the yield of pomelo and banana, aiming to assess the impact of climate change, zoning specialized crops, productivity forecast serving the development of consumption and export plans for pomelos and bananas Keyword: Influence of meteorological factors; Yields of pomelo and banana trees