(Luận văn thạc sĩ) tác động kinh tế của biến đổi khí hậu đến thu nhập trồng lúa ở việt nam

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ TP HỒ CHÍ MINH PHẠM THỊ KIM PHỤNG TÁC ĐỘNG KINH TẾ CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN THU NHẬP TRỒNG LÚA Ở VIỆT NAM LUẬN VĂN THẠC SĨ KINH TẾ TP Hồ Chí Minh – Năm 2012 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ TP HỒ CHÍ MINH -CHƯƠNG TRÌNH GIẢNG DẠY KINH TẾ FULBRIGHT PHẠM THỊ KIM PHỤNG TÁC ĐỘNG KINH TẾ CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN THU NHẬP TRỒNG LÚA Ở VIỆT NAM Chun ngành: Chính sách cơng Mã số: 603114 LUẬN VĂN THẠC SĨ KINH TẾ NGƯỜI HƯỚNG DẪN KHOA HỌC TS Đinh Cơng Khải TP Hồ Chí Minh – Năm 2012 i LỜI CAM ĐOAN Tôi cam đoan luận văn hồn tồn tơi thực Các đoạn trích dẫn số liệu sử dụng luận văn dẫn nguồn có độ xác cao phạm vi hiểu biết Luận văn không thiết phản ánh quan điểm Trường Đại học Kinh tế thành phố Hồ Chí Minh hay Chương trình Giảng dạy Kinh tế Fulbright TP Hồ Chí Minh, ngày 30 tháng năm 2012 Tác giả luận văn Phạm Thị Kim Phụng ii LỜI CẢM ƠN Hồn thành luận văn này, tơi xin cảm ơn đến bậc sinh thành gia đình ln ủng hộ, khích lệ tơi chun tâm học tập Tơi biết ơn Chương trình Giảng dạy Kinh tế Fulbright tạo hội truyền dạy kiến thức kinh tế cho vốn học viên thuộc khối ngành Kinh tế theo học chương trình sau đại học Tôi xin cảm ơn cá nhân tổ chức giúp đỡ tơi nhiệt trình q trình thực luận văn Đặc biệt, xin bày tỏ lời cảm ơn chân thành tới cá nhân giúp tơi lớn q trình thực luận văn: Thầy Đinh Công Khải – Viện Trưởng Viện Chính sách cơng, Đại học Kinh tế thành phố Hồ Chí Minh nhiệt tình hướng dẫn có góp ý quan trọng cho tơi q trình thực hồn thiện luận văn Ơng Nguyễn Châu Thoại – Khoa Kinh tế Tài nguyên, Đại học Tài ngun Mơi trường thành phố Hồ Chí Minh, Việt Nam khơi nguồn ý tưởng cho thực luận văn Ông Salvatore Virdis – Khoa Sinh thái thực vật học Địa chất, Đại học Sassari, Ý Cô Lê Thị Hạnh – chuyên viên GIS dự án IMOLA, Huế, Việt Nam giúp nhiều chun mơn GIS Ơng Ronald Vargas Rojas – chuyên viên quản lý đất FAO người bạn ông cung cấp tài liệu liệu số đồ phân loại đất giới Ông Richard Williams – Khoa Xã hội học, Đại học Notre Dame, Mỹ giúp lớn mặt kỹ thuật tính tốn tác động biên dự báo mơ hình hồi quy thơng thường hồi quy phân vị phần mềm xử lý phân tích số liệu Stata Đồng hành tơi, ln có bạn bè đồng nghiệp, người động viên, chia sẻ giúp đỡ nhiều thời gian vừa qua Chân thành cảm ơn tất người iii TĨM TẮT ĐỀ TÀI Có câu nói: “Nắng mưa chuyện trời, lúa gạo chuyện đời đời phải lo” Quả vậy, việc trồng lúa gạo vốn phụ thuộc nhiều vào yếu tố khách quan khí hậu điều kiện tự nhiên vùng Trong bối cảnh biến đổi khí hậu (BĐKH) nay, thật BĐKH nỗi lo đặc biệt liên quan đến việc làm thu nhập nông dân Việt Nam chiếm khoảng 80% dân số Do vậy, nghiên cứu áp dụng tiếp cận Ricardian để đánh giá mối quan hệ chủ yếu yếu tố khí hậu (nhiệt độ lượng mưa) tác động đến thu nhập trồng lúa (TNTL) nông hộ Việt Nam Sau đó, dự báo mức thiệt hại BĐKH gây cho hộ trồng lúa Việt Nam tương lai GIS xem công cụ hỗ trợ lớn việc kết nối liệu đưa vào mơ hình kinh tế lượng liệu VHLSS 2008, khí hậu, loại đất, độ cao Kết nghiên cứu 4279 hộ trồng lúa từ liệu VHLSS 2008 cho thấy yếu tố khí hậu có tác động tuyến tính phi tuyến tính đến TNTL Nếu nhiệt độ tăng 1oC/tháng TNTL hộ trung bình giảm 425 nghìn đồng/ha/oC lượng mưa tăng 1mm/tháng TNTL hộ tăng nghìn đồng/ha/mm Tuy nhiên, có tăng lượng mưa vào mùa khơ có tác động tích cực TNTL Giả định yếu tố khác không đổi ngoại trừ yếu tố khí hậu thay đổi; rõ ràng, hộ trồng lúa Việt Nam chịu thiệt hại BĐKH nhiệt độ có xu hướng tăng mùa lượng mưa tăng lên vào mùa mưa giảm vào mùa khô Dự báo vào năm 2050 2100, mức độ thiệt hại hộ trồng lúa dao động trung bình từ 15 nghìn đồng/ha đến 1,6 triệu đồng/ha hay dao động từ 0,1% đến 14% so với thu nhập hộ trồng lúa vào năm 2008 Trong đó, hai vùng ảnh hưởng BĐKH vùng Duyên hải Nam Trung Bộ Nam Bộ Ngoài ra, khu vực Trung Du miền núi phía Bắc trở thành vùng cứu cánh cho nông nghiệp lúa gạo Nhà nước chuyển hướng phát triển cho vùng mức cao TNTL phản ánh loại đất trồng, đặc điểm kinh tế hành vi hộ Hầu hết, loại đất trồng lúa có tác động tích cực đến TNTL, Luật bảo vệ đất trồng lúa trở nên khả thi tình Hộ chủ động tưới tiêu, hộ bán lẻ cho tiêu dùng bán sỉ cho tư thương, hộ sống xã có trạm khuyến nơng có tác động tích cực lớn đến TNTL Bên cạnh, cần cân nhắc vấn đề hỗ trợ tín dụng cho nơng hộ, quan tâm TNTL nhóm hộ dân tộc thiểu số, nghiên cứu lại số vụ năm để vừa đạt hiệu kinh tế đảm bảo tính mặt nơng học yếu tố tác động tiêu cực đến TNTL Đề tài hạn chế liệu khí hậu dài hạn thiếu vài liệu để phân tích sâu cho ngành trồng lúa Việt Nam Việc định vị vị trí trồng lúa nằm ngồi phạm vi bao phủ trạm khí tượng lệch lạc mặt khơng gian, ảnh hưởng kết đề tài iv MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii TÓM TẮT ĐỀ TÀI .iii MỤC LỤC iv DANH MỤC CÁC CHỮ VIẾT TẮT VÀ KÝ HIỆU vi DANH MỤC CÁC BẢNG vii DANH MỤC CÁC HÌNH viii DANH MỤC CÁC PHỤ LỤC ix Chương 1: GIỚI THIỆU ĐỀ TÀI 1.1 Bối cảnh vấn đề nghiên cứu 1.2 Mục tiêu nghiên cứu 1.3 Câu hỏi nghiên cứu 1.4 Đối tượng phạm vi nghiên cứu 1.5 Phương pháp nghiên cứu 1.6 Bố cục nghiên cứu Chương 2: CƠ SỞ LÝ THUYẾT 2.1 Tiếp cận Ricardian đánh giá tác động biến đổi khí hậu 2.1.1 2.1.2 Vận dụng lý thuyết lợi so sánh sử dụng đất nơng nghiệp Mơ hình Ricardian 2.2 Tổng quan kết nghiên cứu trước 2.2.1 2.2.2 2.2.3 Các yếu tố khí hậu ảnh hưởng đến thu nhập ròng Các nhóm đất ảnh hưởng đến thu nhập ròng 11 Các đặc điểm kinh tế hộ ảnh hưởng đến thu nhập ròng 12 Chương 3: CÁC YẾU TỐ ẢNH HƯỞNG ĐẾN THU NHẬP TRỒNG LÚA Ở VIỆT NAM VÀ MƠ HÌNH ƯỚC LƯỢNG 14 3.1 Các yếu tố ảnh hưởng đến thu nhập trồng lúa Việt Nam 14 3.1.1 3.1.2 3.1.3 Các yếu tố khí hậu 14 Các yếu tố loại đất 16 Các yếu tố liên quan đặc điểm kinh tế - xã hội hộ 17 3.2 Mơ hình ước lượng, chiến lược ước lượng liệu nghiên cứu 19 3.2.1 3.2.2 3.2.3 Mơ hình ước lượng 19 Chiến lược ước lượng mơ hình 20 Dữ liệu nghiên cứu 20 v Chương 4: KẾT QUẢ NGHIÊN CỨU 22 4.1 Thống kê mô tả 22 4.2 Ma trận tương quan 23 4.3 Kết hồi quy mơ hình 23 4.3.1 4.3.2 4.3.3 Các biến khí hậu tác động đến thu nhập trồng lúa 25 Các nhóm đất tác động đến thu nhập trồng lúa 25 Các biến đặc điểm kinh tế - xã hội tác động đến thu nhập trồng lúa 26 4.4 Tác động biên biến khí hậu tác động đến thu nhập trồng lúa 27 4.5 Dự báo ảnh hưởng biến đổi khí hậu đến thu nhập trồng lúa 29 Chương 5: KẾT LUẬN VÀ GỢI Ý CHÍNH SÁCH 32 5.1 Kết luận đề tài 32 5.2 Gợi ý sách 33 5.2.1 5.2.2 Giảm thiểu BĐKH 33 Thích ứng với BĐKH 34 5.3 Đóng góp nghiên cứu 37 5.4 Hạn chế nghiên cứu 38 5.5 Đề xuất hướng nghiên cứu 39 TÀI LIỆU THAM KHẢO 40 PHỤ LỤC 44 vi DANH MỤC CÁC CHỮ VIẾT TẮT VÀ KÝ HIỆU ANLT : An ninh lương thực BĐKH : Biến đổi khí hậu : Tổ chức Lương thực Nông nghiệp Liên FAO Food and Agriculture Organization of the United Nations ĐBSCL Hiệp Quốc : Đồng sông Cửu Long GIS Geographic Information System : Hệ thống thông tin địa lý GSO General Statistics Office : Tổng cục Thống kê IMHEN Institute Of Meteorology, : Viện Khoa học Khí tượng Thủy văn Môi trường Hydrology and Environment MARD Ministry of Agriculture and Rural : Bộ Nông nghiệp Phát triển Nông thôn : Bộ Tài nguyên Môi trường TNR : Thu nhập ròng TNTL : Thu nhập ròng từ trồng lúa UBND : Ủy ban nhân dân : Điều tra mức sống hộ gia đình Việt Nam : Ngân hàng giới Development MONRE Ministry of Natural Resources and Environment VHLSS Vietnam Household Living Standard Survey WB World Bank vii DANH MỤC CÁC BẢNG Bảng 3.1 Thời gian mùa theo phân vùng khí hậu 15 Bảng 3.2 Số lượng mẫu nghiên cứu theo nhóm đất 16 Bảng 3.3 Các biến giải thích sử dụng mơ hình 19 Bảng 4.1 Kết hồi quy mơ hình hiệu chỉnh 24 Bảng 4.2 Kết phân tích tác động biên mơ hình hiệu chỉnh 27 Bảng 4.3 Mức biến đổi TNTL theo khả kịch BĐKH 29 viii DANH MỤC CÁC HÌNH Hình 2.1 Giá trị đất thay đổi nhiệt độ Hình 3.1 Tổng lượng mưa tháng theo mùa vùng khí hậu 16 Hình 3.2 Nhiệt độ trung bình tháng theo mùa vùng khí hậu 16 Hình 3.3 Các yếu tố ảnh hưởng TNTL 18 Hình 4.1 Tác động biên nhiệt độ lượng mưa đến TNTL theo vùng khí hậu 28 Hình 4.2 Dự báo mức độ thiệt hại hộ TNTL BĐKH theo vùng khí hậu 30 45 Phụ lục Nhiệt độ trung bình giai đoạn 2001 – 2010 Nguồn: Tác giả tính tốn vẽ theo nguồn liệu MARD Phụ lục Tổng lượng mưa giai đoạn 2001 – 2010 Nguồn: Tác giả tính tốn vẽ theo nguồn liệu MARD 46 Phụ lục Các mốc nhiệt độ giai đoan sinh trưởng phát triển lúa Nhiệt độ (oC) Các giai đoạn phát triển lúa ngưỡng thấp 16-19 ngưỡng cao 45 Khoảng tối ưu 18-40 Nhú mầm hạt giống 12 35 25-30 Ra rễ 16 35 25-28 9-16 33 25-31 15 - - Sự phân hóa hoa 15-20 30 - Lúa làm đòng 22 35-36 30-33 12-18 > 30 20-29 Nảy mầm Đẻ nhánh Bắt đầu phân hóa hoa Trổ bơng Nguồn: Yoshida, 1978 trích lại từ Nguu Nguyen Van, 2004 47 Phụ lục Bản đồ phân bố điểm trồng lúa nhóm đất 48 Phụ lục Các mùa vụ trồng lúa vùng năm Vụ vùng Vụ Đông xuân - Đồng sông Hồng - Đồng sông Cửu long - Miền núi - Duyên hải miền Trung - Tây Nguyên, Đông Nam Vụ Hè thu - Đồng sông Hồng - Đồng sông Cửu long - Miền núi - Duyên hải miền Trung - Tây Nguyên, Đông Nam Vụ Mùa/Thu đông - Đồng sông Hồng - Đồng sông Cửu long - Miền núi - Duyên hải miền Trung - Tây Nguyên, Đông Nam Thời gian gieo cấy Thời gian thu hoạch Tháng 2, Tháng 12, Tháng 3, Tháng 1, Tháng 2, Tháng 5, Tháng 3, Tháng 6,7 Tháng 4, Tháng 5, Tháng 4, Tháng 5, Tháng 6, Tháng 7, Tháng 7, Tháng 9, 10 Tháng 7, Tháng Tháng 8, Tháng 8, Tháng 8, Tháng 10, 11 Tháng 12, Tháng 12, Tháng 11, 12 Tháng 10, 11 Nguồn: GSO (2008) 49 Phụ lục Kịch phát thải trung bình B2 Mức tăng nhiệt độ trung bình (oC) thay đổi lượng mưa (%) so với thời kỳ 1980-1999 vùng khí hậu Việt Nam theo kịch phát thải trung bình B2 Vùng Trung bình thời kỳ tháng 12, 1, tháng 3, 4, Tây Bắc tháng 6, 7, tháng 9, 10, 11 tháng 12, 1, tháng 3, 4, Ðông Bắc tháng 6, 7, tháng 9, 10, 11 tháng 12, 1, Đồng tháng 3, 4, Bằng Trung Du tháng 6, 7, Bắc tháng 9, 10, 11 tháng 12, 1, Bắc Trung tháng 3, 4, Bộ tháng 6, 7, tháng 9, 10, 11 tháng 12, 1, Duyên Hải tháng 3, 4, Nam tháng 6, 7, Trung Bộ tháng 9, 10, 11 tháng 12, 1, tháng 3, 4, Tây Nguyên tháng 6, 7, tháng 9, 10, 11 tháng 12, 1, tháng 3, 4, Nam Bộ tháng 6, 7, tháng 9, 10, 11 tháng đến tháng 11 (*) tháng 12 đến tháng (*) 2050 2100 Nhiệt độ Lượng mưa Nhiệt độ Lượng mưa (oC) (%) (oC) (%) 1,5 2,9 3,1 3,7 1,5 -2,8 -3,7 0,8 5,9 1,7 7,8 1,20 1,1 2,5 1,4 1,4 1,9 3,1 2,5 1,4 -2,2 2,8 -2,9 0,8 6,3 1,6 8,3 1,3 1,5 2,6 1,3 2,1 2,8 2,8 1,7 -3,4 3,1 -4,5 0,8 7,5 1,7 9,9 1,1 2,4 2,2 3,1 1,5 2,9 1,4 1,8 -4,9 3,2 -6,5 1,3 7,2 2,6 9,5 1,4 4,2 2,7 5,6 -5,1 -6,7 1 -7,1 2,2 -9,3 0,7 1,9 1,4 2,6 2,1 7,9 0,8 -7,4 1,8 -9,7 0,8 -8,7 1,8 -11,4 0,7 0,1 1,4 0,2 0,7 6,1 1,5 8,1 0,8 -7,7 1,7 -10,1 0,9 -7,2 1,9 -9,4 1,2 0,8 2,1 1,1 1,2 6,3 2,3 8,5 1.01 1.24 4,09 -3,44 2,03 2,53 (*): tác giả tính giá trị trung bình cho nước theo trung bình thời kỳ Nguồn: MONRE (2009), trích từ Phụ lục Phụ lục 10 5,43 -4,51 50 Phụ lục Các bước xây dựng mẫu nghiên cứu cho đề tài GIS hệ thống bao gồm phần mềm, phần cứng, liệu giúp người việc tính tốn, phân tích, thể thơng tin gắn kết với vị trí không gian bề mặt đất (Đỗ Nguyên Hải Hoàng Văn Mùa, 2007) Dưới bước thực tóm tắt nhằm kết nối liệu mà đề tài sử dụng cấp hộ nghiên cứu: Bước 1: Số hóa liệu đặc điểm kinh tế hộ vào GIS Công việc bao gồm bước nhỏ sau: - Xác định mẫu nghiên cứu hộ trồng lúa, bao gồm 4691 mẫu nghiên cứu - Định vị tạo code cho mẫu nghiên cứu dựa theo địa hộ kê khai theo cấp tỉnh, huyện, xã, địa bàn, hộ số Trong đó, code tỉnh, huyện mẫu tương ứng code tỉnh, huyện danh mục hành 2001 cịn code xã mẫu tương ứng với code xã VHLSS 2006 (GSO, 2008) - Kết nối liệu đặc điểm kinh tế hộ tạo code tương ứng với code đồ số cấp xã dạng điểm vào phần mềm ArcGis 10 (cấp thấp công tác số hóa đồ hành tồn quốc biết nay) Lúc này, vị trí mẫu đồ tương ứng điểm mà mẫu thuộc địa phận xã Kết cho thấy 4691 mẫu tương ứng với 1976 điểm đồ, gọi điểm nghiên cứu Bước 2: Xác định loại đất cho điểm nghiên cứu Công việc bao gồm bước nhỏ sau: - Chuyển liệu đất dạng raster thành dạng vector - Dạng vector dạng rời rạc loại đất tính đất Qua đó, xác định nhóm đất cho điểm nghiên cứu dựa theo vị trí đồ (Phụ lục 5) Bước 3: Số hóa liệu khí hậu vào GIS xác định thơng tin khí hậu cho điểm nghiên cứu Công việc bao gồm bước nhỏ sau: - Tạo code phân loại trạm khí tượng theo bảy vùng khí hậu - Tính tốn giá trị trung bình tháng phân theo mùa khơ mùa mưa cho vùng theo đề cập Bảng 3.1 - Chuyển liệu khí hậu vào phần mềm ArcGis 10 dựa theo tọa độ địa lý trạm - Định vị điểm nghiên cứu thuộc phạm vi trạm khí tượng bao phủ cách gán code trạm cho điểm nghiên cứu theo mức ưu tiên sau: 51 (i) Gán code trạm cho điểm trạm nằm địa phận hành xã24 Kết có 26 mẫu gán code 14 trạm (ii) Gán code trạm cho điểm nằm vùng bán kính bao phủ trạm gần độ cao tương ứng (iii) Nếu điểm nhận code trạm thỏa mãn tiêu chí (ii) ưu tiên lấy code trạm nằm địa phận huyện, tỉnh với điểm, theo hướng từ Bắc xuống Nam Đông sang Tây25 Bước 4: Chọn lại mẫu nghiên cứu nhằm loại bỏ quan sát sai lệch khơng gian Quan sát vị trí điểm nghiên cứu so với bán kính trạm khí tượng, cho thấy điểm phân bố cách trạm có bán kính tối đa 120 km, có khoảng 90% điểm nghiên cứu nằm phạm vi bán kính cách trạm khí tượng 35 km (Bảng cho Phụ lục 8) Nhiều nghiên cứu nước sử dụng liệu thu từ vệ tinh, ưu điểm liệu liên tục có phạm vi bao phổ rộng đến vùng trái đất (Mendelsohn đ.t.g, 2004) Nguyễn Hữu Dũng Nguyễn Châu Thoại (2011) chọn mẫu nơng hộ trồng trọt nằm bán kính 50km cách trạm khí tượng gần Tuy nhiên, địa hình khơng đồng vùng miền nên bán kính bao phủ trạm nên khác Trong đề tài này, tác giả áp dụng bán kính bao phủ trạm cho điểm miền Bắc, Trung Nam 30, 35 40km (Phụ lục 9) Việc đưa giá trị đồng thời có xét đến số lượng trạm số lượng mẫu miền hạn chế làm giảm mẫu cho mơ hình Kết cịn lại 1798 điểm nghiên cứu tương ứng 4279 mẫu nghiên cứu Như vậy, kết trình xử lý số liệu nhằm gắn kết liệu cấp độ khác (VHLSS 2008 cấp hộ, liệu đất cấp xã, liệu khí tượng cấp huyện) thành liệu cấp hộ Ví dụ điểm nghiên cứu nằm phạm vi nhóm đất xám thuộc trạm khí tượng Huế hộ trồng lúa thuộc điểm nghiên cứu mang nhóm đất xám mang thơng tin khí hậu trạm khí tượng Huế 24 Khí hậu đồng xã Việc lấy số liệu đồng huyện hẳn không xác dù khơng thấy có trạm nằm huyện Nếu trường hợp mẫu nằm địa bàn huyện có đặt trạm lấy thơng tin trạm nằm ngồi huyện đạt tiêu chí đề 25 Dựa theo hướng gió thổi từ hướng Bắc xuống từ biển vào 52 Bảng cho Phụ lục Khoảng cách trạm khí tượng đến điểm mẫu nghiên cứu Bán kính trạm (km) xã có mẫu (0) 10 20 25 30 35 40 45 50 55 60 65 70 75 120 Số lượng điểm tỷ lệ (%) so với điểm Tương ứng với số lượng mẫu tỷ lệ (%) so với mẫu 14 0,71% 26 0,55% 423 1194 1489 1671 1783 1867 1919 1944 1956 1962 1973 1974 1975 1976 21,41% 60,43% 75,35% 84,56% 90,23% 94,48% 97,12% 98,38% 98,99% 99,29% 99,85% 99,90% 99,95% 100,00% 982 2893 3618 4016 4268 4463 4583 4630 4651 4663 4685 4687 4690 4691 20,93% 61,67% 77,13% 85,61% 90,98% 95,14% 97,70% 98,70% 99,15% 99,40% 99,87% 99,91% 99,98% 100,00% Nguồn: Tác giả tính tốn theo nguồn MARD, GSO (VHLSS, 2008), Cục đồ (2009) 53 Phụ lục Bản đồ chọn điểm nghiên cứu theo vùng khí hậu 54 Phụ lục 10 Bảng thống kê mô tả biến mơ hình hồi quy Biến PLE Đơn vị tính Triệu đồng/ha TeRa o TeRaSq o C/tháng ( C/tháng) 16 28 Độ lệch chuẩn 4,71 1,44 Giá trị nhỏ -18 19 Giá trị lớn 69 29 759 776 72,78 350 814 22 20 20 3,34 12 28 Trung bình 11,75 27 Trung vị 11,70 28 731 Yếu vị TeDr o TeDrSq o ( C/tháng) 493 410 396 152,44 146 787 RaRa mm/tháng 249 230 202 122 708 (mm/tháng) 67391 52931 40866 14976 501512 mm/tháng 47 42 45 74,82 54024,2 20,86 136 2621 1782 1996 2770,49 85 18558 Sac (1: nhóm đất xám) 0,38 0 0,49 SGl (1: nhóm đất glây) SFL (1: nhóm đất phù sa) SAR (1: nhóm đất cát ven biển) 0,36 0 0,48 0,17 0 0,38 0,02 0 0,15 Sex (1: nam) 0,83 1 0,38 48,57 47 45 12,61 20 91 0,21 0 0,41 C/tháng 2 RaRaSq RaDr RaDrSq (mm/tháng) Age Ethn (1: người dân tộc thiểu số) năm Educ lớp 6,97 3,35 12 HoSi người 4,34 4 1,61 14 0,18 0 0,39 0,89 1 0,32 ha/năm 0,77 0,40 0,40 1,92 0.01 64 ha/năm -0,37 -0,40 -0,40 0,40 -2 0,46 0 0,50 0,15 0,47 0 0 0,35 0,50 0 1 0,30 0 0,46 0,05 0 0,22 Rice (1: trồng gạo) Irr (1: tưới tiêu chủ động) Area LgAr (logarit số 10 Area) MiMa (1: bán cho tư thương) Reta (1: bán lẻ) Cred (1: cịn nợ) NoFa (1: có việc làm phi nơng nghiệp) Exte (1: xã có trạm khuyến nơng) 55 Phụ lục 11 Bảng ma trận tương quan biến hồi quy Biến Ple TeRa Ple 1.00 TeRa 0.15 1.00 TeRaSq TeRaSq TeDr TeDrSq RaRa RaRaSq RaDr RaDrSq SAc SGl SFL SAR 0.16 1.00 TeDr (0.33) 0.05 0.04 1.00 TeDrSq (0.34) 0.01 (0.00) 1.00 1.00 RaRa (0.15) (0.37) (0.37) 0.20 0.21 1.00 RaRaSq (0.12) (0.33) (0.33) 0.17 0.18 0.97 1.00 RaDr (0.02) (0.19) (0.18) (0.06) (0.06) 0.71 0.70 1.00 RaDrSq (0.07) (0.19) (0.19) 0.01 0.01 0.74 0.76 0.96 1.00 SAc (0.02) (0.30) (0.30) (0.25) (0.24) 0.11 0.12 0.10 0.08 1.00 SGl 0.01 0.22 0.22 0.11 0.10 (0.03) (0.04) 0.04 0.04 (0.59) 1.00 SFL 0.06 0.24 0.24 0.17 0.17 (0.16) (0.14) (0.24) (0.19) (0.36) (0.34) 1.00 SAR 0.05 (0.01) (0.01) 0.01 0.00 0.08 0.07 0.14 0.13 (0.12) (0.12) (0.07) 1.00 Age Sex Educ Ethn HoSi Rice Irr lgArea MiMa Reta NoFa Cred 1.00 Age (0.00) 0.13 0.13 0.15 0.14 (0.00) (0.01) 0.00 0.00 (0.15) 0.12 0.07 0.02 1.00 Sex 0.00 (0.05) (0.05) (0.06) (0.06) (0.01) 0.00 (0.01) 0.00 0.04 (0.05) (0.02) 0.01 (0.21) Educ 0.18 0.24 0.24 (0.17) (0.19) (0.12) (0.12) 0.08 0.02 (0.04) 0.08 (0.00) 0.01 (0.20) 0.19 1.00 Ethn (0.09) (0.39) (0.39) (0.33) (0.31) 0.08 0.10 (0.05) (0.02) 0.45 (0.35) (0.20) (0.08) (0.19) 0.11 (0.24) 1.00 HoSi (0.08) (0.22) (0.22) 0.02 0.03 0.05 0.05 (0.06) (0.03) 0.10 (0.13) (0.02) (0.04) (0.11) 0.21 (0.11) 0.24 Rice (0.06) 0.11 0.12 0.33 0.33 (0.02) (0.01) (0.10) (0.07) (0.24) 0.09 0.21 0.02 0.03 (0.03) (0.05) (0.18) (0.03) 1.00 0.15 0.25 0.26 0.08 0.07 (0.08) (0.07) (0.01) (0.04) (0.19) 0.18 0.10 0.01 0.09 (0.06) 0.17 (0.33) (0.17) 0.09 1.00 lgArea (0.18) 0.01 0.01 0.38 0.39 (0.04) (0.03) (0.22) (0.14) (0.16) 0.02 0.21 (0.05) 0.05 0.11 (0.10) (0.02) 0.24 0.13 0.02 1.00 MiMa (0.03) 0.14 0.14 0.34 0.34 (0.08) (0.07) (0.15) (0.11) (0.34) 0.21 0.19 0.03 0.10 (0.02) (0.03) (0.24) (0.09) 0.16 0.17 0.47 1.00 Reta 0.09 (0.08) (0.08) (0.18) (0.18) 0.03 0.03 0.07 0.05 0.13 (0.05) (0.07) (0.02) (0.02) 0.03 0.04 0.11 (0.00) (0.07) (0.00) 0.04 (0.38) 1.00 NoFa 0.03 0.00 0.01 (0.04) (0.04) 0.01 0.00 0.01 (0.00) (0.07) 0.03 0.08 (0.01) (0.08) 0.04 0.08 (0.07) 0.06 0.11 0.02 (0.03) (0.04) (0.03) Cred (0.05) (0.03) (0.03) (0.02) (0.02) (0.03) (0.03) (0.03) (0.02) 0.06 (0.02) (0.05) (0.03) (0.12) 0.03 0.01 0.09 0.09 (0.02) 0.01 0.04 0.01 (0.01) 0.04 1.00 Exte (0.00) (0.04) (0.04) 0.03 0.03 0.00 (0.00) (0.01) (0.01) 0.03 0.00 (0.04) (0.04) (0.01) (0.01) (0.03) 0.05 0.01 (0.03) (0.03) (0.01) (0.01) 0.01 (0.02) (0.01) Irr Exte 1.00 1.00 1.00 1.00 56 Phụ lục 12 Bảng kết hồi quy mơ hình ước lượng Median regression Raw sum of deviations 15599.39 (about 11.703042) Min sum of deviations 13446.8 Ple Coef TeRa TeRaSq TeDr TeDrSq RaRa RaRaSq RaDr RaDrSq SAc SGl SFL SAR Age Sex Educ Ethn HoSi Rice Irr lgArea MiMa Reta NoFa Cred Exte _cons -.8894989 0148655 3.374692 -.0844052 -.0251702 0000414 0833387 -.0008085 824266 8712108 2.204541 2.626193 4.47e-06 1925245 0277345 -1.08241 0890716 1808162 1.494209 -.353774 1.113753 1.249175 026207 -.3779701 6308474 -9.899883 Std Err 1.182527 0231449 4792874 0103837 0049301 6.84e-06 0158117 0001214 3167793 339414 3747758 5619853 0062408 2048535 0254323 2613165 0506582 20897 2495866 1099708 2020676 2371534 1630386 1479522 3356435 12.93524 t -0.75 0.64 7.04 -8.13 -5.11 6.05 5.27 -6.66 2.60 2.57 5.88 4.67 0.00 0.94 1.09 -4.14 1.76 0.87 5.99 -3.22 5.51 5.27 0.16 -2.55 1.88 -0.77 P>|t| 0.452 0.521 0.000 0.000 0.000 0.000 0.000 0.000 0.009 0.010 0.000 0.000 0.999 0.347 0.276 0.000 0.079 0.387 0.000 0.001 0.000 0.000 0.872 0.011 0.060 0.444 Number of obs = Pseudo R2 = 4279 0.1380 [95% Conf Interval] -3.207868 -.0305107 2.435038 -.1047626 -.0348357 000028 0523396 -.0010466 2032132 2057823 1.469784 1.524408 -.0122307 -.2090952 -.0221262 -1.594727 -.0102448 -.2288741 1.004889 -.5693741 7175949 7842302 -.2934338 -.6680336 -.027189 -35.2597 1.428871 0602416 4.314345 -.0640479 -.0155047 0000548 1143378 -.0005705 1.445319 1.536639 2.939297 3.727977 0122396 5941443 0775951 -.5700934 1883881 5905065 1.98353 -.1381739 1.509911 1.714119 3458477 -.0879066 1.288884 15.45993 Phụ lục 13 Bảng kết hồi quy mơ hình hiệu chỉnh Median regression Raw sum of deviations 15599.39 (about 11.703042) Min sum of deviations 13446.87 Ple Coef TeRa TeRaSq TeDr TeDrSq RaRa RaRaSq RaDr RaDrSq SAc SGl SFL SAR Sex Educ Ethn HoSi Rice Irr lgArea MiMa Reta Cred Exte _cons -.8638214 0143138 3.392905 -.0848489 -.0253489 0000417 0836931 -.0008119 8517163 9196837 2.251346 2.673046 1938217 0288392 -1.090182 0934391 1838715 1.465654 -.3498604 1.103479 1.228078 -.3801926 6455524 -10.37657 Std Err 1.09462 0214242 4428999 0095967 0045724 6.36e-06 0146996 0001135 2934806 3136572 345931 5198943 1878727 0228903 2405472 0467614 1920222 2312032 1016843 18653 2190833 1360373 3103646 11.96219 t -0.79 0.67 7.66 -8.84 -5.54 6.55 5.69 -7.15 2.90 2.93 6.51 5.14 1.03 1.26 -4.53 2.00 0.96 6.34 -3.44 5.92 5.61 -2.79 2.08 -0.87 P>|t| 0.430 0.504 0.000 0.000 0.000 0.000 0.000 0.000 0.004 0.003 0.000 0.000 0.302 0.208 0.000 0.046 0.338 0.000 0.001 0.000 0.000 0.005 0.038 0.386 Number of obs = Pseudo R2 = 4279 0.1380 [95% Conf Interval] -3.009847 -.0276888 2.52459 -.1036635 -.0343131 0000292 0548741 -.0010345 2763412 304752 1.57314 1.653782 -.1745068 -.0160377 -1.56178 0017624 -.1925923 1.012375 -.5492147 7377828 7985605 -.6468966 0370759 -33.82871 1.282204 0563164 4.26122 -.0660343 -.0163847 0000542 112512 -.0005893 1.427091 1.534615 2.929551 3.692311 5621501 0737162 -.6185838 1851158 5603353 1.918933 -.1505061 1.469175 1.657596 -.1134886 1.254029 13.07556 57 Phụ lục 14 Xu hướng tác động yếu tố khí hậu đến TNTL Nhiệt độ mùa mưa (TeRa) Thu nhap rong tu lua (trieu dong/ha) 11.5 12 12.5 13 Thu nhap rong tu lua (trieu dong/ha) 12 14 16 18 13.5 Lượng mưa mùa mưa (RaRa) 10 x 200 400 600 Tong Luong mua thang trung binh cua mua mua mm/thang x 18 22 24 Nhiet mua mua C/thang 26 28 Nhiệt độ mùa khô (TeDr) 14 Lượng mưa mùa khô (RaDr) 12 20 x 6 thu nhap rong lua (trieu/ha) 10 12 Thu nhap rong tu lua (trieu dong/ha) 10 x 50 100 Tong Luong mua thang trung binh cua mua kho mm/thang x: giá trị 10 15 20 Nhiet mua kho C/thang 25 Giá trị tới ngưỡng Bảng cho Phụ lục 14 Biến số khí hậu TeRa TeDr RaRa RaDr Hệ số hồi quy bậc -0.86382 3.392905 -0.02535 0.083693 Hệ số hồi quy bậc 0.0143138 -0.0848489 0.0000417 -0.0008119 Giá trị trung bình tới ngưỡng (*) 30.17oC/tháng 19.99oC/tháng 303.94 mm/tháng 51.54 mm/tháng Giá trị trung bình 27.00oC/tháng 21.96oC/tháng 248.59mm/tháng 46.76 mm/tháng Giá trị yếu tố khí hậu tới ngưỡng đạo hàm biến phụ thuộc theo biến giải thích cho yếu tố khí hậu (vế trái công thức (2.4) 0) 30 58 Phụ lục 15 Dự báo TNTL hộ theo khả kịch BĐKH (triệu đồng/ha/hộ) sum p0 p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11 p12 p13 p14 Variable Obs Mean p0 p1 p2 p3 p4 4279 4279 4279 4279 4279 11.51188 11.43439 11.46918 10.97028 11.4973 p5 p6 p7 p8 p9 4279 4279 4279 4279 4279 p10 p11 p12 p13 p14 4279 4279 4279 4279 4279 Std Dev Min Max 2.22565 2.240419 2.207548 2.738258 2.223672 4.905142 4.849225 4.908757 3.222991 4.854744 15.94703 15.89492 15.845 15.83951 15.91256 11.39169 10.9557 10.85999 11.38635 11.45711 2.222001 2.737514 2.741644 2.255865 2.202004 4.852839 3.172593 3.117091 4.644962 4.911805 15.79289 15.80505 15.71086 15.87226 15.81357 10.12626 11.49186 11.33158 10.10624 9.925933 3.347981 2.223325 2.231385 3.348335 3.343829 1.185678 4.838058 4.761222 1.118593 1.042872 15.44964 15.90088 15.7388 15.40821 15.20249 Giả sử yếu tố khác không đổi, ngoại trừ yếu tố khí hậu thay đổi Các kết dự báo dựa theo mơ hình hồi quy hiệu chỉnh Phụ lục 13 với khả BĐKH xảy sau: p0: TNTL dự báo theo yếu tố khí hậu (2001-2010) Khả p1: TNTL dự báo nhiệt độ mùa mưa tăng lên 1.01oC/tháng vào năm 2050 Khả p2: TNTL dự báo lượng mưa mùa mưa tăng lên 4.1% (10.2mm/tháng) vào năm 2050 Khả p3: TNTL dự báo nhiệt độ mùa khô tăng lên 1.24oC/tháng vào năm 2050 Khả p4: TNTL dự báo lượng mưa mùa khô giảm 3.4% (1.6mm/tháng) vào năm 2050 Khả p5: TNTL dự báo nhiệt độ mùa mưa tăng lên 1.01oC/tháng lượng mưa mùa mưa tăng lên 4.1% (10.2mm/tháng) vào năm 2050 Khả p6: TNTL dự báo nhiệt độ mùa khô tăng lên 1.24oC/tháng lượng mưa mùa khô giảm 3.4% (1.6mm/tháng) vào năm 2050 Khả p7: TNTL dự báo thay đổi nhiệt độ lượng mưa đồng thời vào mùa theo khả vào năm 2050 Khả p8: TNTL dự báo nhiệt độ mùa mưa tăng lên 2.03oC/tháng vào năm 2100 Khả p9: TNTL dự báo lượng mưa mùa mưa tăng lên 5.43% (13.5mm/tháng) vào năm 2100 Khả p10: TNTL dự báo nhiệt độ mùa khô tăng lên 2.53oC/tháng vào năm 2100 Khả p11: TNTL dự báo lượng mưa mùa khô giảm 4.51% (2.1mm/tháng) vào năm 2100 Khả p12: TNTL dự báo nhiệt độ mùa mưa tăng lên 2.03oC/tháng lượng mưa mùa mưa tăng lên 5.43% (13.5mm/tháng) vào năm 2100 Khả p13: TNTL dự báo nhiệt độ mùa khô tăng lên 2.53oC/tháng lượng mưa mùa khô giảm 4.51% (2.1mm/tháng) vào năm 2100 Khả p14: TNTL dự báo thay đổi nhiệt độ lượng mưa đồng thời vào mùa theo khả 12 13 vào năm 2100 59 Phụ lục 16 Bản đồ dự báo tác động BĐKH vào năm 2100 đến TNTL hộ ... biến khí hậu tác động đến thu nhập trồng lúa 25 Các nhóm đất tác động đến thu nhập trồng lúa 25 Các biến đặc điểm kinh tế - xã hội tác động đến thu nhập trồng lúa 26 4.4 Tác động. .. TRƯỜNG ĐẠI HỌC KINH TẾ TP HỒ CHÍ MINH -CHƯƠNG TRÌNH GIẢNG DẠY KINH TẾ FULBRIGHT PHẠM THỊ KIM PHỤNG TÁC ĐỘNG KINH TẾ CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN THU NHẬP TRỒNG LÚA Ở VIỆT NAM Chun ngành:... Chương CÁC YẾU TỐ ẢNH HƯỞNG ĐẾN THU NHẬP TRỒNG LÚA Ở VIỆT NAM VÀ MƠ HÌNH ƯỚC LƯỢNG 3.1 Các yếu tố ảnh hưởng đến thu nhập trồng lúa Việt Nam Ở Việt Nam, TNTL bị ảnh hưởng yếu tố đề cập chương