Luận án nghiên cứu thống kê ảnh hưởng của biến đổi khí hậu đến năng suất lúa ở miền bắc việt nam

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ QUỐC DÂN - NGUYỄN ĐĂNG KHOA NGHIÊN CỨU THỐNG KÊ ẢNH HƯỞNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN NĂNG SUẤT LÚA Ở MIỀN BẮC VIỆT NAM LUẬN ÁN TIẾN SĨ NGÀNH KINH TẾ HỌC HÀ NỘI – 2023 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ QUỐC DÂN - NGUYỄN ĐĂNG KHOA NGHIÊN CỨU THỐNG KÊ ẢNH HƯỞNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN NĂNG SUẤT LÚA Ở MIỀN BẮC VIỆT NAM Chuyên ngành: Kinh tế học (Thống kê kinh tế) Mã số: 9310101 LUẬN ÁN TIẾN SĨ Người hướng dẫn khoa học: PGS.TS TRẦN THỊ KIM THU HÀ NỘI – 2023 i LỜI CAM ĐOAN Tôi cam kết danh dự cá nhân nghiên cứu tự thực khơng vi phạm quy định liêm học thuật nghiên cứu khoa học Trường Đại học Kinh tế Quốc dân Hà Nội, ngày tháng năm 2023 Nghiên cứu sinh Nguyễn Đăng Khoa ii MỤC LỤC LỜI CAM ĐOAN i MỤC LỤC ii DANH MỤC CHỮ VIẾT TẮT iv DANH MỤC BẢNG v DANH MỤC BIỂU ĐỒ vi DANH MỤC HÌNH viii PHẦN MỞ ĐẦU CHƯƠNG 1: CƠ SỞ LÝ LUẬN VÀ TỔNG QUAN NGHIÊN CỨU VỀ TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN SẢN XUẤT LÚA 1.1 Cơ sở lý luận biến đổi khí hậu 1.1.1 Khí hậu biến đổi khí hậu 1.1.2 Nguyên nhân biến đổi khí hậu 14 1.1.3 Đo lường khí hậu đánh giá biến đổi khí hậu 18 1.2 Hoạt động sản xuất lúa suất lúa 25 1.2.1 Khái niệm chung hoạt động trồng trọt 25 1.2.2 Hoạt động sản xuất lúa (trồng lúa) 28 1.2.3 Năng suất lúa 31 1.3 Các kênh tác động biến đổi khí hậu đến suất lúa 35 1.4 Tổng quan nghiên cứu tác động biến đổi khí hậu đến suất lúa 39 1.4.1 Tổng quan nghiên cứu tác động BĐKH đến kết hoạt động sản xuất nông nghiệp giới 40 1.4.2 Tổng quan nghiên cứu tác động BĐKH đến sản xuất lúa Việt Nam 47 1.5 Khung nghiên cứu tác động BĐKH đến suất lúa miền Bắc Việt Nam 54 CHƯƠNG 2: PHƯƠNG PHÁP NGHIÊN CỨU 56 2.1 Xây dựng mô hình kinh tế lượng phân tích ảnh hưởng BĐKH đến suất lúa 56 2.1.1 Tiếp cận xây dựng mơ hình phân tích 56 2.1.2 Mô hình kinh tế lượng phân tích ảnh hưởng BĐKH đến suất lúa 57 2.1.3 Giải thích biến đưa vào mơ hình 58 2.1.4 Phương pháp ước lượng mơ hình dự đốn thống kê 64 2.2 Nguồn liệu cách thức thu thập liệu 72 2.2.1 Nguồn liệu 72 2.2.2 Cách thức thu thập, tổng hợp liệu 73 iii CHƯƠNG 3: THỰC TRẠNG BIẾN ĐỔI KHÍ HẬU VÀ SẢN XUẤT LÚA Ở KHU VỰC MIỀN BẮC VIỆT NAM GIAI ĐOẠN 2001 – 2020 76 3.1 Tổng quan khu vực Miền Bắc Việt Nam 76 3.1.1 Đặc điểm điều kiện tự nhiên khí hậu 76 3.1.2 Đặc điểm tình hình phát triển kinh tế xã hội miền Bắc 79 3.2 Thực trạng sản xuất lúa miền Bắc giai đoạn 2001 - 2020 82 3.2.1 Diện tích trồng lúa 82 3.2.2 Sản lượng suất lúa hàng năm 84 3.2.3 Đặc điểm canh tác lúa địa phương miền Bắc 87 3.2.4 Tình hình sản xuất lúa hai mùa vụ miền Bắc 89 3.3 Biến đổi khí hậu khu vực miền Bắc giai đoạn 2001 - 2020 92 3.3.1 Biến động nhiệt độ 92 3.3.2 Lượng mưa 96 3.3.3 Biến động thiên tai miền Bắc nước ta giai đoạn 2001 – 2020 99 3.3.4 Biến động tượng thời tiết cực đoan 101 3.4 Thực trạng sách thích ứng BĐKH sản xuất lúa 105 CHƯƠNG 4: PHÂN TÍCH THỐNG KÊ TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN NĂNG SUẤT LÚA Ở MIỀN BẮC GIAI ĐOẠN 2001 – 2020 109 4.1 Đặc điểm nguồn số liệu định hướng phân tích 109 4.1.1 Đặc điểm nguồn số liệu sử dụng cho phân tích 109 4.1.2 Định hướng phân tích 110 4.2 Thống kê mô tả biến mơ hình phân tích 111 4.2.1 Mô tả đặc trưng phân phối biến sử dụng mơ hình 111 4.2.2 Mô tả mối quan hệ tương quan biến khí hậu suất lúa 117 4.3 Phân tích thống kê ảnh hưởng BĐKH đến suất lúa Miền Bắc Việt Nam giai đoạn 2001 – 2020 119 4.3.1 Tác động BĐKH đến suất lúa vụ Đông Xuân 120 4.3.2 Tác động BĐKH đến suất lúa vụ Mùa 124 4.4 Dự đoán suất lúa miền Bắc với liệu năm 2021 127 4.4.1 Dự đoán suất lúa Đông Xuân miền Bắc với liệu năm 2021 128 4.4.2 Dự đoán suất lúa Mùa tỉnh miền Bắc năm 2021 129 KẾT LUẬN 133 DANH MỤC CÁC CƠNG TRÌNH NGHIÊN CỨU CỦA TÁC GIẢ LIÊN QUAN ĐẾN ĐỀ TÀI LUẬN ÁN 139 DANH MỤC TÀI LIỆU THAM KHẢO 140 PHỤ LỤC 155 iv DANH MỤC CHỮ VIẾT TẮT Chữ viết tắt Giải nghĩa Tiếng Việt BĐKH Biến đổi khí hậu BQT Bình qn trượt ĐBSCL Đồng sơng Cửu Long ĐBSH Đồng sông Hồng NN&PTNN TD-MNPB TN&MT Ttb Nông nghiệp phát triển nông thôn Trung du miền Núi phía Bắc Tài ngun mơi trường Nhiệt độ trung bình Tiếng Anh IPCC International Panel on Climate Change (Ủy ban liên phủ biến đổi khí hậu) WB World Bank (Ngân hàng giới) UN United Nations (Liên hiệp quốc) UNDP United Nations Development Programme (Chương trình phát triển Liên hiệp quốc) Food and Agriculture Organization (Tổ chức nông nghiệp giới) FAO GSO Food and Agriculture Organization’s Statistics Office (Cơ quan thống kê tổ chức nông nghiệp giới) General Statistics Office (Tổng cục thống kê Việt Nam) IRRI International Rice Research Institute (Viện nghiên cứu lúa gạo quốc tế) FAOSTAT FE The United Nations office for Disaster Risk Reduction (Tổ chức liên hiệp quốc giảm nhẹ rủi ro thiên tai) Vietnam Household Living Standard Survey (Khảo sát mức sống hộ gia đình, Việt Nam) Fixed Effects (Tác động cố định) FGLS Feasible Generalized Least Square (Bình phương tối thiểu tổng quát hóa khả thi) PCSE Panel Corrected Standard Error (Sai số hiệu chỉnh cho liệu mảng) RCP Representative Concentration Pathways (Kịch đường nồng độ khí nhà kính đại diện) Ministry of Natural Resource and Environment (Bộ Tài nguyên môi trường, Việt Nam) Cooperative for American Remittances to Europe (Hợp tác trao đổi Mỹ - Âu) UNDRR VHLSS MONRE CARE SRD MARD Center of Sustainable Rural Development (Trung tâm phát triển nông thôn bền vững) Ministry of Agriculture and Rural Development (Bộ Nông nghiệp Phát triển nông thôn) v DANH MỤC BẢNG Bảng 1.1 Hệ thống tiêu thống kê chủ yếu sử dụng để đo lường khí hậu 19 Bảng 2.1: Tổng hợp biến mơ hình kỳ vọng tác động biến tới suất lúa 63 Bảng 3.1 Lực lượng lao động từ 15 tuổi trở lên khu vực miền Bắc giai đoạn 2010 – 2020 80 Bảng 3.2 Tỷ lệ lao động từ 15 tuổi trở lên làm việc kinh tế qua đào tạo khu vực miền Bắc (%) 81 Bảng 3.3 Phân bố diện tích gieo trồng lúa ĐBSH năm 2020 83 Bảng 3.4 Phân bố diện tích gieo trồng lúa TD&MNPB năm 2020 84 Bảng 4.1 Thống kê mô tả biến xác định cho vụ Đơng Xn tồn miền Bắc (tính từ liệu cấp tỉnh giai đoạn 2001 – 2020) 111 Bảng 4.2 Thống kê mơ tả biến mơ hình vụ Đơng Xn tính cho khu vực (TD-MNPB ĐBSH) .113 Bảng 4.3 Thống kê mô tả biến xác định cho vụ Mùa (toàn miền Bắc) giai đoạn 2001 – 2020 (tính từ liệu cấp tỉnh) .114 Bảng 4.4 Thống kê mô tả biến mơ hình cho vụ Mùa khu vực TD-MNPB ĐBSH giai đoạn 2001 - 2020 .116 Bảng 4.5 Kết ước lượng tác động BĐKH đến suất lúa vụ Đông Xuân Miền Bắc .120 Bảng 4.6 Kết luận tác động yếu tố khí hậu đến suất lúa Đơng Xuân khu vực miền Bắc 123 Bảng 4.7 Kết ước lượng tác động BĐKH đến suất lúa Mùa miền Bắc 124 Bảng 4.8 Kết luận tác động yếu tố khí hậu đến suất lúa Mùa khu vực miền Bắc 127 Bảng 4.9 Dự đoán suất lúa Đông Xuân tỉnh miền Bắc năm 2021 128 Bảng 4.10 Kết dự đoán suất lúa Mùa cho tỉnh miền Bắc với liệu 2021 129 vi DANH MỤC BIỂU ĐỒ Biểu đồ 3.1 Diện tích lúa năm miền Bắc giai đoạn 2001 - 2020 83 Biểu đồ 3.2 Biến động sản lượng lúa hàng năm khu vực Miền Bắc giai đoạn 2001 - 2020 85 Biểu đồ 3.3 Biến động suất lúa chung năm miền Bắc giai đoạn 2001 - 2020 86 Biểu đồ 3.4 Biến động suất lúa Vụ đông xuân với vụ Mùa Khu vực miền Bắc giai đoạn 2001 - 2020 87 Biểu đồ 3.5 Biến động suất lúa vụ Đông Xuân khu vực miền Bắc giai đoạn 2001 - 2020 .90 Biểu đồ 3.6 Biến động diện tích sản lượng lúa Mùa miền Bắc giai đoạn 2001 2020 91 Biểu đồ 3.7 Biến động suất lúa Mùa khu vực miền Bắc giai đoạn 2001 – 2020 91 Biểu đồ 3.8 Biến động nhiệt độ trung bình năm khu vực miền Bắc giai đoạn 2001 2020 92 Biểu đồ 3.9 Biến động nhiệt độ tháng cao khu vực miền Bắc giai đoạn 2001 - 2020 93 Biểu đồ 3.10 Biến động nhiệt độ tháng thấp khu vực miền Bắc giai đoạn 2001 2020 94 Biểu đồ 3.11 Biến động nhiệt độ trung bình vụ Đơng Xn khu vực miền Bắc giai đoạn 2001 - 2020 95 Biểu đồ 3.12 Biến động nhiệt độ trung bình vụ Mùa khu vực miền Bắc giai đoạn 2001 - 2020 .96 Biểu đồ 3.13 Biến động lượng mưa khu vực miền Bắc giai đoạn 2001 - 2020 97 Biểu đồ 3.14 Biến động lượng mưa vụ Đông Xuân khu vực miền Bắc giai đoạn 2001 - 2020 98 Biểu đồ 3.15 Biến động lượng mưa vụ Mùa khu vực miền Bắc giai đoạn 2001 - 2020 99 Biểu đồ 3.16 Biến động số lượng bão (từ cấp 10 trở lên đổ bộ) miền Bắc giai đoạn 2001 - 2020 .100 Biểu đồ 3.17 Số đợt lũ/ lụt trung bình khu vực miền Bắc giai đoạn 2001 – 2020 (tính trung bình hàng năm từ liệu địa phương) 101 Biểu đồ 3.18 Số đợt rét đậm - rét hại kéo dài liên tục ngày tổng số ngày rét đậm rét hại khu vực TD&MNPB giai đoạn 2001 – 2020 102 vii Biểu đồ 3.19 Ngày rét đậm -rét hại số ngày bình quân đợt rét đậm rét hại khu vực ĐBSH giai đoạn 2001 – 2020 102 Biểu đồ 3.20 Số ngày nắng nóng kéo dài số ngày bình quân đợt nóng khu vực TD&MNPB giai đoạn 2001 – 2020 .104 Biểu đồ 3.21 Số đợt nắng nóng kéo dài số ngày bình qn đợt nóng khu vực ĐBSH giai đoạn 2001 – 2020 105 Biểu đồ 4.1 Quan hệ biến nhiệt độ suất lúa vụ Đông Xuân 118 Biểu đồ 4.2 Quan hệ biến lượng mưa suất lúa Đông Xuân 118 Biểu đồ 4.3 Quan hệ biến nhiệt độ suất lúa vụ Mùa .119 Biểu đồ 4.4 Quan hệ biến lượng mưa suất lúa vụ Mùa 119 viii DANH MỤC HÌNH Hình 1.1 Biểu đồ minh họa biến đổi nhiệt độ trung bình tượng thời tiết cực đoan 13 Hình 1.2 Thay đổi nồng độ số thành phần khí khí trái đất khoảng 200 kỷ trước .17 Hình 1.3 Khung lý thuyết nghiên cứu ảnh hưởng BĐKH đến suất miền Bắc Việt Nam 54 Hình 3.1 Bản đồ khu vực miền Bắc – Việt Nam 76 155 PHỤ LỤC Phụ lục 1: Kịch biến đổi nhiệt độ đến cuối kỷ 21 cho khu vực miền Bắc theo kịch RPC4.5 RPC8.5 156 157 Phụ lục 2: Kịch biến đổi lượng mưa đến cuối kỷ 21 cho khu vực miền Bắc theo kịch RPC4.5 RPC8.5 158 159 Phụ lục 3: Kịch biến đổi số ngày nắng nóng gay gắt đến cuối kỷ 21 cho khu vực miền Bắc 160 Phụ lục Hệ số tương quan Pearson biến khí hậu suất lúa Đông Xuân Ln_NS Ln_NS 1.000 T_Mean -0.267 T_Mean T_min T_max R_Total R_min R_max Extr_cold_days 1.000 (0.000) T_min T_max R_Total R_min R_max Extr_cold_days -0.296 0.541 1.000 (0.000) (0.000) 0.130 0.382 0.097 (0.005) (0.000) (0.034) 0.060 -0.253 -0.128 -0.024 (0.192) (0.000) (0.005) (0.603) 0.011 -0.038 -0.047 -0.069 0.343 (0.820) (0.415) (0.309) (0.135) (0.000) 0.092 -0.144 -0.075 -0.033 0.647 0.020 (0.045) (0.002) (0.103) (0.468) (0.000) (0.670) -0.257 -0.119 -0.052 0.118 0.088 -0.052 0.034 (0.001) (0.009) (0.256) (0.010) (0.055) (0.259) (0.458) 1.000 1.000 1.000 1.000 1.000 (Nguồn: tính tốn tác giả) Ghi chú: giá trị (…) giá trị P_value kiểm định quan hệ tương quan hai biến 161 Phụ lục Hệ số tương quan Pearson kiểm định tương quan biến khí hậu suất lúa Mùa Ln_NS Ln_NS T_Mean T_min T_max R_Total R_min R_max Extr_hot_days Storm Flood T_Mean T_min T_max R_Total R_min R_max Extr_hot_days 0.0811 (0.077) 0.0539 (0.241) 0.0333 (0.070) -0.230 (0.000) -0.0118 (0.7976) -0.156 (0.0006) -0.3326 (0.0000) 0.6161 (0.0000) 0.6062 (0.0000) -0.2331 (0.0000) -0.0758 (0.0989) -0.0381 (0.4071) -0.2054 (0.000) 0.1424 (0.0019) 0.0592 (0.1976) 0.1123 (0.0143) 0.1272 (0.0055) -0.1541 (0.0008) -0.336 (0.0000) -0.2595 (0.000) -0.1557 (0.0007) 0.1085 (0.018) 0.3828 (0.000) 0.6562 (0.000) 0.1898 (0.000) 0.0489 (0.288) -0.0756 (0.0998) 0.0542 (0.2386) -0.3638 (0.000) -0.5495 (0.000) 0.0077 (0.8666) -0.1636 (0.0003) 0.0889 (0.0527) -0.0164 (0.7213) -0.1082 (0.0183 -0.2071 (0.000) 0.2125 (0.000) 0.4132 (0.000) 0.1917 (0.000) 0.1845 0.0001 0.1177 (0.0102) 0.2238 (0.000) 0.177 (0.0001) 0.3304 (0.000) Storm Flood 1 1 1 0.2574 (0.000) (Nguồn: tính tốn tác giả) Ghi chú: giá trị (…) giá trị P_value kiểm định quan hệ tương quan hai biến 162 Phụ lục 6: Kết phân tích tượng đa cộng tuyến với biến mơ hình vụ Đơng Xn 163 Phụ lục 7: Kết phân tích tượng đa cộng tuyến với biến mơ hình vụ Mùa 164 Phụ lục 8: Kết ước lượng mơ hình cho vụ Đơng Xn • Mơ hình 1A (MH1A) Ln_NS ln_l ln_s ln_m ln_p ln_f Irr Ln_Eco Extr_cold_days T_Mean R_Total _cons sigma_u sigma_e Rho • Coef 0.111638 0.158321 0.026204 -0.04439 -0.00601 0.001955 0.009362 -0.00367 -0.00401 -0.00308 3.27733 0.077944 0.037695 0.810448 Robust Std Err 0.0245291 0.0126069 0.0087563 0.010070 0.0146967 0.0006132 0.0082077 0.0006958 0.0021438 0.0015842 0.1179998 (fraction t 4.55 12.56 2.99 -4.41 -0.41 3.19 1.14 -5.28 -1.87 -1.94 27.81 P 0.000 0.000 0.003 0.000 0.710 0.002 0.255 0.000 0.062 0.053 0.000 of variance due [95% Conf 0.0634367 0.1335462 0.0089936 -0.0641813 -0.0348115 0.0007495 -0.0067669 -0.0050406 -0.0082256 -0.0061907 3.02404 to Interval] 0.1598399 0.1830934 0.0434072 -0.0246048 0.02279961 0.0031596 0.0254909 -0.002306 0.0002 0.0000355 3.483946 u_i) Mơ hình 2A (MH2A): Ln_NS ln_l ln_s Coef 0.110122 0.156683 Robust Std Err 0.024635 0.0127427 ln_m ln_p ln_f Irr Ln_Eco Extr_cold_days 0.026287 -0.04385 -0.00588 0.002002 0.008294 0.0087997 0.0101341 0.0147762 0.0006161 0.0083085 2.99 -4.33 -0.40 3.25 1.00 -0.00381 T_mean R_Total [95% Conf 0.0617117 0.1316427 Interval] 0.1585321 0.1817242 0.003 0.000 0.718 0.001 0.319 0.0089945 -0.0637691 -0.0348462 0.0007915 -0.0080326 0.0435789 -0.0239402 0.0230765 0.0032128 0.0246215 0.0007128 -5.35 0.000 -0.0052142 -0.0024128 -0.00533 0.0316411 -0.17 0.866 -0.0675115 0.0568445 T_mean_square R_total_square _cons 0.000745 3.77E-05 -0.00044 3.2863463 0.0043855 0.000823 0.0004689 0.3233022 0.17 0.05 -0.93 10.17 -0.0078734 -0.0015797 -0.0013597 2.612641 0.0093624 0.001655 0.0004831 3.909485 sigma_u 0.0778492 sigma_e 0.0377406 Rho 0.8097015 (fraction t P 4.47 0.000 12.30 0.000 0.865 0.964 0.350 0.000 of variance due to u_i) 165 • Mơ hình 3A (MH3A): Robust Ln_NS ln_l ln_s ln_m ln_p ln_f Irr Ln_Eco Extr_cold_days T_min T_max R_min R_max _cons Coef Std Err 0.10319 0.02414 0.15529 0.01235 0.02330 0.00873 -0.03982 0.00989 -0.00519 0.01447 0.00200 0.00060 0.00467 0.00809 t 4.27 12.58 2.67 -4.03 -0.36 3.32 0.58 P 0.000 0.000 0.008 0.000 0.756 0.001 0.564 [95% Conf 0.05576 0.13102 0.00615 -0.05926 -0.03355 0.00082 -0.01123 Interval] 0.15063 0.17955 0.04045 -0.02039 0.02317 0.00318 0.02058 -0.00353 -0.00473 0.00661 -0.03217 0.00068 3.08340 -5.08 -3.79 3.86 -2.41 0.20 26.37 0.000 0.000 0.000 0.016 0.839 0.000 -0.00490 -0.00717 0.00325 -0.05838 -0.00593 2.74244 -0.00217 -0.00228 0.00997 -0.00596 0.00730 3.31855 0.00070 0.00125 0.00171 0.01334 0.00337 0.11692 sigma_u sigma_e 0.073156 0.036858 rho 0.797548 (fraction • of variance due to u_i) Mơ hình 4A (MH4A): Ln_NS ln_l ln_s ln_m ln_p ln_f Irr Ln_Eco T_min T_max R_min R_max T_max_square _cons Coef 0.102639 0.155079 0.024187 -0.0401 -0.00551 0.00201 0.004644 -0.00368 -0.00473 0.052497 -0.03283 0.001076 -0.0010 2.526709 sigma_u sigma_e 0.0767352 0.0367254 rho 0.8136321 Extr_cold_days Robust Std Err 0.0240536 0.0123035 0.0087054 0.0098549 0.0144994 0.0005988 0.008064 0.000697 0.0012407 0.0220502 0.0132943 0.0033585 0.000478 0.293601 t 4.27 12.60 2.78 -4.07 -0.38 3.34 0.58 -5.28 -3.81 2.38 -2.47 0.32 -2.09 8.61 P>t 0.000 0.000 0.006 0.000 0.748 0.001 0.565 0.000 0.000 0.018 0.014 0.749 0.037 0.000 of variance (fraction due [95% Conf 0.0553712 0.1309012 0.0070802 -0.0594708 -0.0339261 0.0008234 -0.0112024 -0.0050518 -0.0071666 0.0091663 -0.0589536 -0.005524 -0.0017694 1.882553 to u_i) Interval] 0.1499071 0.1792568 0.0412944 -0.0207387 0.02291152 0.0031768 0.0204907 -0.0023123 -0.0022905 0.0958285 -0.0067039 0.0076758 -0.0000533 3.155066 166 Phụ lục 9: Kết ước lượng mơ hình cho vụ Mùa • Mơ hình 1B( MH1B) Ln_NS ln_l ln_s ln_m ln_p ln_f Irr Ln_Eco Storm Flood Extr_hot-_days _cons Coef 0.1154969 0.023257 0.0232488 -0.0309658 0.101667 0.0074441 0.0037376 -0.0129618 -0.0239695 -0.0027206 -0.0020439 -0.0002292 2.556033 sigma_u sigma_e 0.12143994 0.06701481 rho 0.76656442 T_mean R_Total • Robust Std Err 0.0804948 0.0267767 0.0296126 0.021002 0.0387051 0.0017753 0.0209436 0.0061008 0.0073731 0.0010702 0.0066375 0.0012868 0.3715132 (fraction of t 1.43 0.87 0.79 -1.47 2.63 4.19 0.18 -2.12 -3.25 -2.54 -0.31 -0.18 6.88 P 0.164 0.394 0.440 0.153 0.015 0.000 0.860 0.044 0.003 0.018 0.761 0.860 0.000 variance due [95% Conf -0.0506362 -0.0320074 -0.0378687 -0.0743119 0.0217835 0.00378 -0.0394879 -0.0255533 -0.0391867 -0.0049293 -0.0157431 -0.002885 1.789267 Interval] 0.28163 0.078522 0.084366 0.01238 0.181551 0.011108 0.046963 -0.00037 -0.00875 -0.00051 0.011655 0.002427 3.322798 to u_i) Mơ hình 2B (MH2B) Robust Ln_NS_VUMUA ln_l ln_s ln_m ln_p ln_f Irr Ln_Eco Storm Flood Extr_hot_days Coef Std Err 0.0668634 0.0251772 0.0275759 0.0199984 0.0365595 0.0014373 0.0196826 0.0060037 0.0075847 0.0011158 0.2291162 0.0043711 0.0044461 0.0001614 3.002003 sigma_u sigma_e 0.052709 0.0282485 0.0289667 -0.0341508 0.0919244 0.0066764 0.0083759 -0.0137498 -0.0242763 -0.0029626 0.5934401 -0.0062104 -0.0110462 0.0001957 -4.849603 0.1108546 0.06482021 rho 0.74520574 (fraction T_mean R_Total T_mean_square R_Total_square _cons t 0.79 1.12 1.05 -1.71 2.51 4.65 0.43 -2.29 -3.20 -2.66 2.59 -1.42 -2.68 1.21 -1.62 Of P 0.438 0.273 0.304 0.10 0.019 0.000 0.674 0.031 0.004 0.014 0.016 0.168 0.013 0.237 0.119 [95% Conf Interval] -0.0852902 0.190708 -0.0237147 0.080212 -0.0279472 0.085881 -0.0754255 0.007124 0.0164693 0.16738 0.0037099 0.009643 -0.032247 0.048999 -0.0261409 -0.00136 -0.0399304 -0.00862 -0.0052655 -0.00066 0.1204745 1.06622 -0.0152319 0.002811 -0.0207783 -0.00243 -0.0001373 0.000529 -11.04543 1.346227 variance due to u_i) 167 • Mơ hình 3B (MH3B): Ln_NS ln_l ln_s ln_m ln_p ln_f Irr Ln_Eco Storm Flood Extr_hot_days T_min T_max R_min R_max _cons Coef 0.1078738 0.0212454 0.0229205 -0.0354483 0.0993048 0.0076158 0.0102851 -0.0132503 -0.0227851 -0.0026703 0.0014272 -0.002617 -0.0004457 -0.0085349 2.608306 sigma_u sigma_e 0.12188726 0.0665192 rho 0.77051327 • Robust Std Err 0.0814021 0.026125 0.0291458 0.019902 0.038258 0.0017624 0.0208956 0.006264 0.0075164 0.0010775 0.0035721 0.0066272 0.0071302 0.0023881 0.3527344 (fraction t 1.33 0.81 0.79 -1.78 2.60 4.32 0.49 -2.12 -3.03 -2.48 0.40 -0.39 -0.06 -3.57 7.39 P 0.198 0.424 0.439 0.088 0.016 0.000 0.627 0.045 0.006 0.021 0.693 0.696 0.951 0.002 0.000 of [95% Conf -0.060132 -0.032674 -0.0372336 -0.0765256 0.0203442 0.0039784 -0.0328412 -0.0261786 -0.0382982 -0.0048942 -0.0059452 -0.0162948 -0.0151617 -0.0134636 1.880298 Interval] 0.27588 0.075165 0.083075 0.005629 0.178265 0.011253 0.053411 -0.00032 -0.00727 -0.00045 0.0088 0.011061 0.01427 -0.00361 3.336313 variance due to u_i) Mơ hình 4B (MH4B) Ln_NS ln_l ln_s ln_m ln_p ln_f Irr Ln_Eco Storm Flood Extr_hot_days T_min T_max R_min R_max T_max_square T_min_square _cons Coef 0.0706888 0.0227615 0.0229749 -0.0355323 0.1025922 0.0073796 0.0121041 -0.0150224 -0.0230882 -0.0030682 -0.0101322 0.3490143 -0.0002304 -0.0064472 -0.0060925 0.0002027 -2.088856 Robust Std Err 0.0691323 0.0249397 0.0274957 0.018766 0.0367897 0.0015915 0.0198471 0.0061068 0.0075515 0.0011623 0.0490983 0.1960157 0.0073708 0.0020166 0.0034071 0.0010981 2.708812 sigma_u sigma_e rho 0.11131107 0.06519096 0.74459995 (fraction t 1.02 0.91 0.84 -1.89 2.79 4.64 0.61 -2.46 -3.06 -2.64 -0.21 1.78 -0.03 -3.20 -1.77 0.18 0.77 P 0.317 0.370 0.412 0.070 0.010 0.000 0.548 0.021 0.005 0.014 0.838 0.091 0.975 0.004 0.086 0.855 0.448 of [95% Conf -0.0719932 -0.0287114 -0.0337735 -0.0742635 0.0266621 0.004095 -0.0288582 -0.0276262 -0.0386738 -0.0054671 -0.1114661 -0.0606422 -0.015443 -0.0106092 -0.0131244 -0.0020637 -7.679569 Interval] 0.213371 0.074235 0.079723 0.003199 0.178522 0.010664 0.053066 -0.00242 -0.0075 -0.00067 0.091202 0.748471 0.014982 -0.00229 0.00094 0.002469 3.501856 variance due to u_i) 168 Phụ lục 10 Phân tích đặc điểm biến động lựa chọn phương pháp dự đoán cho biến độc lập chưa có liệu 2021 (mơ hình dự báo suất lúa Đông Xuân) Đặc điểm biến động Chỉ định phương pháp dự đoán Sai số phương pháp dự đoán BQT trọng số (MSE) Chi hạt giống Xu tăng tất định Hàm xu tuyến tính - - Hàm xu tuyến tính Chi máy móc-sức kéo Xu tăng tất định Hàm xu tuyến tính - - Hàm xu tuyến tính Chi thuốc sâu-diệt cỏ Xu tăng tất định Hàm xu tuyến tính - - Hàm xu tuyến tính Khối lượng phân bón 1ha Xu tăng tất định Hàm xu tuyến tính - - Hàm xu tuyến tính Tăng/giảm khơng ổn định BQT giản đơn/BQT trọng số 0.071 0.108 BQT trọng số Biến động tăng BQT giản đơn/BQT trọng số 5.286 9.643 BQT trọng số Biến Lao động Tỷ lệ diện tích đất NN tưới Sai số phương pháp dự đoán BQT giản đơn (MSE) Lựa chọn phương pháp dự đốn (Nguồn: Phân tích Tổng hợp tác giả) 169 Phụ lục 11 Phân tích đặc điểm biến động lựa chọn phương pháp dự đốn cho biến độc lập chưa có liệu 2021 (mơ hình dự báo suất lúa Mùa) Biến Đặc điểm biến động Chỉ định phương pháp dự đoán Sai số phương pháp dự đoán BQT trọng số (MSE) Chi hạt giống Xu tăng tất định Hàm xu tuyến tính - - Hàm xu tuyến tính Chi máy móc-sức kéo Xu tăng tất định Hàm xu tuyến tính - - Hàm xu tuyến tính Chi thuốc sâu-diệt cỏ Xu tăng tất định Hàm xu tuyến tính - - Hàm xu tuyến tính Khối lượng phân bón 1ha Xu tăng tất định Hàm xu tuyến tính - - Hàm xu tuyến tính Tăng/giảm khơng ổn định BQT giản đơn/BQT trọng số 0.071 0.108 BQT trọng số Biến động tăng, xu ngẫu nhiên BQT giản đơn/BQT trọng số 5.286 9.643 BQT trọng số Lao động Tỷ lệ diện tích đất tưới Sai số phương pháp dự đoán BQT giản đơn (MSE) Lựa chọn phương pháp dự đốn (Nguồn: Phân tích Tổng hợp tác giả)