Tổng quan về luận án

Biến đổi khí hậu (BĐKH) toàn cầu cùng sự gia tăng về tần suất lẫn cường độ của các loại hình thiên tai cực đoan đang đặt ra những thách thức nghiêm trọng đối với an ninh lương thực và sự phát triển bền vững của ngành nông nghiệp. Nằm ở trung tâm khu vực Bắc Trung Bộ với địa hình phân hóa phức tạp từ miền núi, trung du, đồng bằng đến dải ven biển, tỉnh Nghệ An được ví như một "Việt Nam thu nhỏ" nhưng đồng thời cũng là một trong những tâm điểm hứng chịu tổn thương nặng nề nhất do thiên tai. Thống kê khí tượng - thủy văn cho thấy trong hai thập kỷ qua, khu vực Bắc Trung Bộ đã chịu tác động của 38 cơn bão, trong đó có 10 cơn bão đổ bộ trực tiếp vào Nghệ An. Riêng năm 2020, "ước tính thiệt hại về kinh tế khoảng 1.327,019 tỷ đồng, trong đó (giông, lốc, sét, mưa đá 250,12 tỷ đồng, bão số 5: 141,053 tỷ đồng, đợt mưa từ ngày 15 đến ngày 20/10: 149,349 tỷ đồng; hoàn lưu bão số 9: 786,497 tỷ đồng)" [1]. Thực trạng này đe dọa trực tiếp tới sinh kế của hơn 70% số hộ dân làm nông nghiệp trên địa bàn tỉnh.

Khoảng trống nghiên cứu (Research gap) then chốt mà luận án của NCS. Đặng Văn Thanh (Chuyên ngành Kỹ thuật Tài nguyên nước, Mã số: 9580212, Trường Đại học Thủy lợi, năm 2022) tập trung giải quyết xuất phát từ hạn chế của y văn hiện hành: hầu hết các công trình nghiên cứu trước đây tại Việt Nam chỉ dừng lại ở phương pháp thống kê mô tả thiệt hại quá khứ hoặc đánh giá định tính (Hạ Thị Thiều Dao, 2016; Lương Ngọc Thúy, 2015; Phan Sĩ Mẫn, 2013). Một số nghiên cứu kinh tế lượng tiệm cận mô hình Ricardian (Nguyen Thi Lan Huong, 2018; UNDP, 2012) hoặc hồi quy đa biến (Trần Đại Nghĩa, 2015) thì lại chỉ xem xét các biến khí hậu dài hạn (nhiệt độ, lượng mưa) một cách độc lập mà bỏ qua tổ hợp các cú sốc thiên tai cực đoan (bão, hạn hán, xâm nhập mặn) diễn ra đồng thời. Luận án đã thiết lập các câu hỏi nghiên cứu và giả thuyết cụ thể:

  • Câu hỏi nghiên cứu 1 (RQ1): Mức độ tác động và cường độ thiệt hại do tổ hợp thiên tai (bão, hạn, mặn) tương tác với biến đổi khí hậu (nhiệt độ, lượng mưa) làm suy giảm doanh thu ròng của hộ nông dân trồng trọt tại Nghệ An như thế nào?
  • Câu hỏi nghiên cứu 2 (RQ2): Theo các kịch bản BĐKH của Bộ Tài nguyên và Môi trường (kịch bản trung bình RCP 4.5 và kịch bản cao RCP 8.5), quy mô thiệt hại về doanh thu của hộ trồng trọt sẽ biến động ra sao đến các mốc thời gian trong tương lai?
  • Câu hỏi nghiên cứu 3 (RQ3): Các biện pháp thích ứng, ứng phó của nông hộ tạo ra hiệu ứng đệm kinh tế có ý nghĩa thống kê như thế nào trong việc bảo toàn doanh thu?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 (H1): Các cú sốc thiên tai cấp tính (chỉ số khô hạn Sazonov $Sa$, độ mặn, cường độ bão) và biến thiên nhiệt - ẩm cực trị có mối quan hệ phi tuyến, tác động tiêu cực đáng kể đến doanh thu trồng trọt.
  • Giả thuyết 2 (H2): Dưới kịch bản phát thải cao RCP 8.5, mức suy giảm doanh thu của các loại cây trồng hàng năm và lâu năm sẽ nghiêm trọng hơn đáng kể so với kịch bản RCP 4.5 nếu thiếu vắng các can thiệp thích ứng chủ động.
  • Giả thuyết 3 (H3): Việc chủ động áp dụng các giải pháp ứng phó (thủy lợi nội đồng, chuyển dịch thời vụ, thay đổi cơ cấu giống) có hệ số tác động dương, giảm thiểu tổn thất kinh tế ròng của nông hộ.

Khung lý thuyết của luận án tích hợp phương pháp tiếp cận kinh tế lượng Ricardian cắt ngang (Ricardian cross-sectional approach) kết hợp mô hình hồi quy tuyến tính đa biến OLS vi mô. Phạm vi nghiên cứu bao quát chuỗi dữ liệu thứ cấp giai đoạn 2001-2020 và bộ dữ liệu khảo sát sơ cấp trực tiếp hộ nông dân năm 2020 trên 3 loại cây trồng chủ lực đại diện cho 3 phân khúc sản xuất: cây lúa (chiếm 51,2% diện tích gieo trồng cây hàng năm), cây cam (chiếm 57,6% diện tích và 51,5% sản lượng cây ăn quả), và cây chè (chiếm 76,3% diện tích và 95,5% sản lượng cây công nghiệp lâu năm toàn tỉnh).

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về lượng hóa tác động của BĐKH và thiên tai đến nông nghiệp được định hình qua ba luồng tiếp cận chính:

  1. Luồng mô hình sinh lý cây trồng và mô phỏng kỹ thuật: Các nghiên cứu kinh điển của Eitzinger và cộng sự (2003), Torriani và cộng sự (2007), hay tại Việt Nam là Huỳnh Thị Lan Hương (2013) sử dụng mô hình DSSAT để đánh giá nhu cầu nước và năng suất lúa tại Quảng Ngãi. Ưu điểm của luồng này là phản ánh chính xác các phản ứng sinh học - nông học của cây trồng trước thay đổi nhiệt độ, bốc thoát hơi nước và bức xạ. Tuy nhiên, nhược điểm chí mạng là bỏ qua hoàn toàn hành vi thích ứng kinh tế (economic adaptation) của con người, dẫn đến việc ước tính sai lệch hoặc phóng đại quá mức thiệt hại thực tế.
  2. Luồng hàm sản xuất kinh tế vi mô (Production Function Approach): Điển hình qua các công trình của Isik & Devadoss (2006), Poudel & Kotani (2013). Cách tiếp cận này mô hình hóa mối quan hệ trực tiếp giữa các yếu tố đầu vào vật chất, điều kiện thời tiết với sản lượng đầu ra. Dẫu vậy, phương pháp hàm sản xuất vi mô vẫn vướng phải sự chỉ trích do giả định cố định công nghệ và hành vi canh tác, không tính đến chi phí điều chỉnh biên khi thời tiết thay đổi.
  3. Luồng mô hình kinh tế Hedonic Ricardian: Do Mendelsohn, Nordhaus và Shaw khởi xướng năm 1994, sau đó được phát triển hoàn thiện bởi Di Falco & Veronesi (2013), Dinar & Mendelsohn (2011). Mô hình Ricardian phân tích dữ liệu cắt ngang chéo giữa các vùng sinh thái khí hậu khác nhau để đo lường sự thay đổi trong giá trị đất đai hoặc doanh thu ròng nông nghiệp, từ đó nội suy được toàn bộ năng lực tự thích ứng nội sinh của hộ nông dân.

Y văn quốc tế ghi nhận nhiều tranh luận và bằng chứng thực nghiệm phân hóa sâu sắc. Nghiên cứu của Aung Tun Oo, Guido Van Huylenbroeck và Stijn Speelman (2020) tại các khu vực khô hạn của Myanmar chỉ ra rằng nhiệt độ tăng từ 2-4°C trong giai đoạn 2060-2100 làm sụt giảm 18-35% thu nhập trồng trọt theo mô hình PCM, và giảm 15-42% theo mô hình CGCM3 [29]. Ngược lại, nghiên cứu của Eid Helmy M., El-Marsafawy và Ouda (2007 - World Bank Policy Research Working Paper 4239) tại Ai Cập cho thấy nhiệt độ tăng 1°C làm giảm $968,94/ha đất trồng không có vật nuôi, nhưng khi có hệ thống tưới tiêu thích ứng, hiệu ứng tổn thất được chặn đứng và chuyển thành mức tăng thu nhập $26,17/ha [30]. Trong khi đó, Niggol Seo, Mendelsohn và Munasinghe (2005) tại Sri Lanka chứng minh biên độ tác động của BĐKH dao động rất rộng từ thiệt hại 20% đến gia tăng 72% tùy thuộc vào sự tương tác giữa lượng mưa và nhiệt độ trong từng kịch bản cụ thể [31].

Tại Việt Nam, nghiên cứu của Trần Đại Nghĩa (2015) tại Đồng bằng sông Cửu Long phát hiện khi mức độ nghiêm trọng của hạn hán tăng 1 bậc thì năng suất lúa giảm 1,6%, và lụt tăng 1 điểm làm năng suất giảm 1,2% [16]. Nghiên cứu của Nguyen Thi Lan Huong (2018) ứng dụng mô hình Ricardian tại vùng Tây Bắc dựa trên bộ dữ liệu VHLSS 2012 đã phân tích doanh thu vụ mùa theo nhiệt độ và lượng mưa mùa khô/mưa nhưng chưa lồng ghép biến thiên tai bão, hạn, mặn [35].

Luận án của Đặng Văn Thanh đã định vị chính xác khoảng trống khoa học này: Tiên phong tích hợp đồng thời các biến thiên tai cực đoan cấp tính (bão, hạn, mặn) vào khung kinh tế lượng Ricardian cắt ngang cùng với các biến khí hậu dài hạn (nhiệt độ, lượng mưa) và biến hành vi thích ứng vi mô, tạo nên một bước tiến quan trọng trong phương pháp lượng hóa tổn thất kinh tế nông nghiệp tại vùng nhiệt đới gió mùa.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án tạo ra bước phát triển học thuật có ý nghĩa đối với ba lý thuyết nền tảng:

  • Mở rộng Lý thuyết Ricardian trong Kinh tế tài nguyên và Nông nghiệp (Mendelsohn et al., 1994; Dinar, 2011): Bằng cách đưa vector thiên tai cực đoan (Cường độ bão, Chỉ số hạn hán Sazonov $Sa$, Nồng độ xâm nhập mặn) vào hàm doanh thu ròng bên cạnh các biến khí hậu truyền thống (Nhiệt độ tối đa/tối thiểu mùa vụ, Lượng mưa vụ Đông Xuân và Hè Thu), luận án đã khắc phục giới hạn của mô hình Ricardian nguyên bản vốn thường chỉ tập trung vào các biến đổi khí hậu mang tính bình quân, trơn láng.
  • Bổ sung Lý thuyết Hành vi Thích ứng Nông hộ (Farm-level Adaptation Theory - Di Falco & Veronesi, 2013): Khẳng định vai trò của các biến cấu trúc hộ gia đình (tuổi, trình độ học vấn, quy mô diện tích, tiếp cận tín dụng, khuyến nông) và các quyết định kỹ thuật (thay đổi giống, thời vụ, thủy lợi chủ động) như những cơ chế đệm kinh tế (economic shock absorbers) làm biến đổi độ dốc của hàm phản ứng thiệt hại.
  • Xác lập hệ thống mệnh đề lý thuyết thực nghiệm:
    • Mệnh đề P1: Tác động của khí hậu và thiên tai lên doanh thu ròng nông nghiệp mang tính chất bậc hai (phi tuyến); tồn tại ngưỡng nhiệt - ẩm tối ưu cho từng tiểu vùng khí thái nông nghiệp.
    • Mệnh đề P2: Tác động tổng hợp của thiên tai liên hoàn (bão kèm mưa lớn cực đoan, hạn hán kéo theo xâm nhập mặn) tạo ra hiệu ứng suy giảm doanh thu theo cấp số cộng hoặc nhân, vượt trội so với tác động của sự gia tăng nhiệt độ bình quân đơn lẻ.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đa tầng giữa khoa học Thủy lợi - Khí tượng thủy văn và Kinh tế lượng vi mô:

[Biến số Khí hậu Dài hạn & BĐKH]
(Nhiệt độ max/min, Lượng mưa các vụ, Kịch bản RCP 4.5 & RCP 8.5)
                 +                                ===> [MÔ HÌNH RICARDIAN CẮT NGANG] ===> Lượng hóa thiệt hại doanh thu
[Biến số Thiên tai Cực đoan]                            (Phân tích theo Lúa, Cam, Chè)
(Cường độ bão, Chỉ số hạn Sa, Nồng độ xâm nhập mặn)
                 +                                ===> [MÔ HÌNH HỒI QUY OLS ĐA BIẾN] ===> Đánh giá hiệu quả thích ứng
[Đặc tính Nông hộ & Biện pháp Thích ứng Vi mô]
(Hạ tầng tưới, Chuyển dịch mùa vụ, Giống chịu hạn/mặn, Nguồn lực)

Điều kiện biên (Boundary conditions) của khung phân tích được xác lập chặt chẽ: Áp dụng cho các nông hộ quy mô vừa và nhỏ trong điều kiện thị trường nông sản cạnh tranh tại các vùng sinh thái chuyển tiếp nhiệt đới gió mùa, nơi tài nguyên nước và rủi ro thời tiết cực đoan đóng vai trò là yếu tố hạn chế sản xuất hàng đầu.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết lý thực chứng (Positivism), sử dụng phương pháp định lượng thực nghiệm nghiêm ngặt kết hợp dữ liệu sơ cấp và thứ cấp:

  • Triết lý nghiên cứu: Thu thập dữ liệu thực nghiệm khách quan để kiểm định các giả thuyết kinh tế lượng, lượng hóa chính xác mối quan hệ nhân quả giữa các yếu tố thời tiết cực đoan và hiệu quả kinh tế nông hộ.
  • Thiết kế đa cấp độ (Multi-level design):
    • Cấp độ vĩ mô/trung mô: Chuỗi dữ liệu khí tượng thủy văn trạm quan trắc (Trạm khí tượng Vinh, Quỳnh Lưu,...) và số liệu thống kê ngành nông nghiệp Nghệ An giai đoạn 2001-2020.
    • Cấp độ vi mô: Dữ liệu mặt cắt chéo (cross-sectional microdata) thu thập thông qua khảo sát trực tiếp các hộ nông dân trồng lúa, chè, cam tại các huyện trọng điểm năm 2020.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu bao gồm 4 bước liên hoàn:

  1. Thu thập dữ liệu thứ cấp và xử lý chỉ số khí hậu: Khai thác dữ liệu khí tượng (nhiệt độ, lượng mưa, bão) và thủy văn (độ mặn, dòng chảy). Tính toán chỉ số hạn hán Sazonov ($Sa$) dựa trên công thức khí tượng chuẩn để phân loại các cấp độ hạn hán thực nghiệm.
  2. Thiết kế mẫu và khảo sát sơ cấp: Áp dụng phương pháp chọn mẫu phân tầng ngẫu nhiên (Stratified Random Sampling) theo tiểu vùng sinh thái (đồng bằng ven biển, trung du, miền núi) và theo 3 loại cây trồng đại diện. Bảng hỏi khảo sát được chuẩn hóa qua các bước phỏng vấn thử nghiệm (pilot study) và tham vấn chuyên gia theo phương pháp Dinar (2011).
  3. Tam giác giác hóa dữ liệu (Triangulation): Đối chiếu chéo dữ liệu tự báo cáo của nông hộ (doanh thu, chi phí, thiệt hại) với dữ liệu quan trắc khí tượng trạm đo và số liệu thống kê chính thức của Sở Nông nghiệp & PTNT, Chi cục Thủy lợi Nghệ An nhằm loại bỏ sai số hồi ức (recall bias).
  4. Kiểm định độ tin cậy và giá trị: Đánh giá tính đại diện của mẫu khảo sát, kiểm tra hệ số phân bố và chuẩn hóa thang đo các biến kinh tế - xã hội.

Data và phân tích

Mô hình kinh tế lượng được thiết lập và ước lượng thông qua phần mềm chuyên dụng Stata và Microsoft Excel:

1. Mô hình Ricardian lượng hóa thiệt hại doanh thu:

$$\ln(Y_i) = \beta_0 + \sum_{j} \beta_1 X_{ji} + \sum_{k} \beta_2 C_{ki} + \sum_{k} \beta_3 C_{ki}^2 + \sum_{m} \beta_4 D_{mi} + \sum_{n} \beta_5 A_{ni} + \varepsilon_i$$

Trong đó:

  • $Y_i$: Doanh thu ròng từ hoạt động trồng trọt của hộ $i$ (triệu đồng/ha).
  • $X_{ji}$: Vector đặc tính kinh tế - xã hội nông hộ (tuổi chủ hộ, trình độ học vấn, quy mô diện tích đất, lao động gia đình, tiếp cận tín dụng, tập huấn kỹ thuật).
  • $C_{ki}, C_{ki}^2$: Biến khí hậu tuyến tính và bậc hai (Nhiệt độ trung bình/cực đại/cực tiểu vụ Đông Xuân và Hè Thu; Lượng mưa vụ Đông Xuân và Hè Thu).
  • $D_{mi}$: Biến thiên tai cực đoan (Mức độ hạn hán theo chỉ số $Sa$, Nồng độ xâm nhập mặn, Cường độ và tần suất bão).
  • $A_{ni}$: Các biến đại diện cho biện pháp ứng phó của hộ nông dân.
  • $\varepsilon_i$: Sai số ngẫu nhiên của mô hình.

2. Mô hình hồi quy tuyến tính OLS phân tích hiệu quả biện pháp ứng phó:

Đánh giá mức độ đóng góp của từng quyết định thích ứng vi mô (đầu tư hệ thống tưới tiết kiệm, kiên cố hóa kênh mương, nạo vét thủy lợi nội đồng, chuyển đổi cơ cấu giống ngắn ngày/chịu mặn, điều chỉnh khung lịch thời vụ gieo cấy) đến việc phục hồi doanh thu ròng của hộ nông dân.

3. Quy trình kiểm định Robustness và sai số dự báo:

Mô hình được kiểm định nghiêm ngặt hiện tượng đa cộng tuyến thông qua hệ số phóng đại phương sai (VIF < 10), kiểm định phương sai sai số thay đổi bằng kiểm định Breusch-Pagan / White và khắc phục triệt để bằng sai số chuẩn vững (Robust Standard Errors). Độ tin cậy dự báo được đánh giá qua việc so sánh sai số dự báo của mô hình so với thực tế (RMSE, MAPE). Dự báo tương lai được ngoại suy dựa trên kịch bản biến đổi khí hậu chính thức của Bộ Tài nguyên và Môi trường (Kịch bản RCP 4.5 và RCP 8.5) đến giữa thế kỷ và cuối thế kỷ 21.

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả ước lượng thực nghiệm từ bộ dữ liệu của luận án đã đem lại 5 phát hiện cốt lõi:

  1. Hiệu ứng phi tuyến và tính nhạy cảm nhiệt độ cực đoan: Nhiệt độ tối đa vụ Hè Thu và Đông Xuân có tác động tiêu cực mạnh mẽ và có ý nghĩa thống kê ($p < 0,01$) đến doanh thu trồng trọt. Hiện tượng nắng nóng kéo dài làm gia tăng bốc thoát hơi nước và gia tăng nhu cầu tưới, dẫn đến sụt giảm nghiêm trọng năng suất lúa và cam. Ngược lại, lượng mưa vụ Đông Xuân có mối quan hệ dạng parabol ngược với doanh thu: mưa gia tăng ở mức vừa phải mang lại tác động dương, nhưng khi vượt ngưỡng tối ưu sẽ gây ngập úng cục bộ, làm đảo chiều suy giảm doanh thu ròng.

  2. Cú sốc kép từ hạn hán và xâm nhập mặn: Mức độ hạn hán (được lượng hóa qua chỉ số Sazonov $Sa$) và độ mặn tại các vùng ven sông/ven biển (Cửa Lò, Nghi Lộc, Hưng Nguyên) là hai nhân tố hủy hoại doanh thu lớn nhất. Cứ mỗi cấp độ gia tăng của hạn hán và mặn hóa khiến doanh thu trung bình của hộ trồng lúa suy giảm từ 8,5% đến 14,2% nếu không có nguồn nước ngọt bổ sung kịp thời.

  3. Tổn thất đột biến do bão lũ cực đoan: Cường độ bão tương tác với lượng mưa ngày cực đại gây thiệt hại tức thời và nặng nề nhất đối với các cây ăn quả lâu năm (cam) và cây công nghiệp (chè) do gãy đổ, xói mòn đất tầng mặt và ngập úng rễ, củng cố bằng chứng thực tế khi hoàn lưu bão số 9 năm 2020 đã tàn phá 786,497 tỷ đồng sản xuất nông nghiệp tại Nghệ An.

  4. Giá trị bù đắp kinh tế của các biện pháp thích ứng vi mô: Kết quả mô hình hồi quy OLS chỉ ra rằng các hộ chủ động áp dụng các giải pháp ứng phó kỹ thuật (đặc biệt là đầu tư hạ tầng tưới tiêu chủ động kết hợp giống kháng hạn/mặn) có mức sụt giảm doanh thu thấp hơn từ 18% đến 32% so với các hộ thụ động không áp dụng biện pháp phòng vệ trước cùng một cường độ thiên tai.

  5. Dự báo phân hóa sâu sắc theo các kịch bản BĐKH: Dưới kịch bản biến đổi khí hậu trung bình (RCP 4.5), mức tổn thất doanh thu trồng trọt toàn tỉnh dự báo dao động từ 7,2% đến 15,6% vào giữa thế kỷ. Tuy nhiên, dưới kịch bản phát thải cao (RCP 8.5), nếu các yếu tố công nghệ và hạ tầng không thay đổi, mức thiệt hại doanh thu ròng có thể leo thang vượt mức 28,4% vào cuối thế kỷ 21 do sự cộng hưởng khốc liệt giữa gia tăng nhiệt độ cực trị và tần suất hạn - mặn gay gắt.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Luận án mở rộng thành công mô hình kinh tế lượng Ricardian vào lĩnh vực quản lý tài nguyên nước và phòng chống thiên tai tại vùng duyên hải miền Trung, cung cấp bằng chứng thực nghiệm khẳng định tầm quan trọng của việc tích hợp các biến thiên tai cấp tính vào các mô hình kinh tế lượng khí hậu.
  • Đổi mới phương pháp luận: Xây dựng quy trình tích hợp giữa các chỉ số khí tượng thủy văn chuyên sâu (chỉ số $Sa$, độ mặn cửa sông) với vi mô kinh tế nông hộ, có thể chuyển giao và nhân rộng cho các lưu vực sông khác tại Việt Nam và các quốc gia Đông Nam Á.
  • Hàm ý thực tiễn cho hộ nông dân: Khuyến nghị mô hình chuyển dịch cơ cấu cây trồng linh hoạt: ưu tiên mở rộng diện tích các giống lúa ngắn ngày chịu mặn, áp dụng công nghệ tưới nhỏ giọt tiết kiệm nước cho vùng cam Phủ Quỳ, và phát triển thảm phủ hữu cơ chống hạn cho vùng chè Thanh Chương, Anh Sơn.
  • Hàm ý chính sách quản lý nhà nước (Implementation pathway):
    1. Giải pháp hạ tầng thủy lợi: Nâng cấp, kiên cố hóa hệ thống đê ngăn mặn và các công trình tiêu thoát lũ động lực khu vực hạ lưu sông Lam; hiện đại hóa các hồ chứa thủy lợi xung yếu nhằm chủ động tích trữ nước ngọt cho mùa khô.
    2. Giải pháp phi công trình: Nâng cao chất lượng dự báo, cảnh báo sớm thiên tai cục bộ; hoàn thiện thể chế quản lý rủi ro thiên tai dựa vào cộng đồng (CBDRM); thiết lập cơ chế tài chính quỹ dự phòng và xúc tiến các sản phẩm bảo hiểm nông nghiệp chỉ số thời tiết.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  • Giới hạn về dạng dữ liệu: Bộ dữ liệu sơ cấp được khảo sát tại một thời điểm cắt ngang (năm 2020), chưa xây dựng được bộ dữ liệu bảng theo dõi dài hạn (Panel Data) để kiểm soát triệt để các yếu tố bất biến theo thời gian của từng hộ nông dân.
  • Phạm vi loại hình thiên tai: Luận án tập trung sâu vào 3 loại hình thiên tai chủ yếu (bão, hạn, mặn), chưa lượng hóa chi tiết các dạng thiên tai thứ cấp khác như lũ quét, sạt lở đất cục bộ ở miền núi cao hoặc sương muối, rét hại.
  • Giới hạn kịch bản khí hậu: Nghiên cứu dựa trên kịch bản BĐKH của Bộ Tài nguyên và Môi trường với hai kịch bản RCP 4.5 và RCP 8.5 (không có kịch bản RCP 2.6 và chưa cập nhật các kịch bản SSP mới nhất của IPCC AR6).

Chương trình nghị sự cho các nghiên cứu tiếp theo (Future Research Agenda):

  • Phát triển mô hình Ricardian dạng dữ liệu bảng động (Dynamic Panel Ricardian Model) theo dõi mẫu nông hộ qua 5-10 năm.
  • Mở rộng phân tích sang mô hình Cân bằng Tổng thể Khả tính (CGE) nhằm đánh giá tác động lan tỏa liên ngành giữa nông nghiệp, công nghiệp chế biến và thị trường lao động nông thôn.
  • Nghiên cứu sâu về phân tích chi phí - lợi ích biên (Marginal Cost-Benefit Analysis) đối với từng gói công nghệ tưới thông minh thích ứng BĐKH.

Tác động và ảnh hưởng

  • Tác động học thuật: Cung cấp cơ sở dữ liệu và khung lý thuyết chuẩn tắc cho các đề tài nghiên cứu chuyên sâu về kinh tế tài nguyên nước, quản lý rủi ro thiên tai và thích ứng BĐKH tại Việt Nam.
  • Chuyển đổi thực tiễn ngành nông nghiệp: Định hướng cơ cấu lại ngành trồng trọt tỉnh Nghệ An theo hướng nâng cao giá trị gia tăng và khả năng chống chịu khí hậu, bảo vệ chuỗi giá trị các nông sản chủ lực (lúa gạo, cam Vinh, chè Nghệ An).
  • Đóng góp chính sách cấp tỉnh và quốc gia: Là luận cứ khoa học trực tiếp phục vụ Sở Nông nghiệp & PTNT, Chi cục Thủy lợi và UBND tỉnh Nghệ An trong việc rà soát, lập "Quy hoạch thủy lợi phòng chống thiên tai giai đoạn 2021-2030, tầm nhìn 2050" và "Kế hoạch hành động ứng phó với BĐKH ngành nông nghiệp".
  • Giá trị xã hội và quốc tế: Góp phần bảo vệ sinh kế, giảm thiểu rủi ro tái nghèo do thiên tai cho hàng trăm nghìn hộ nông dân; đồng thời là bài học kinh nghiệm điển hình cho các quốc gia đang phát triển trong việc lượng hóa tổn thất kinh tế do biến đổi khí hậu vi mô.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới học giả: Kế thừa khung lý thuyết Ricardian mở rộng tích hợp chỉ số thiên tai cực đoan và phương pháp ước lượng kinh tế lượng vi mô.
  • Cơ quan Hoạch định Chính sách (Bộ NN&PTNT, UBND tỉnh Nghệ An, Chi cục Thủy lợi): Tiếp cận bộ công cụ lượng hóa thiệt hại và ma trận dự báo tổn thất theo kịch bản để phân bổ nguồn vốn đầu tư công trình - phi công trình tối ưu.
  • Doanh nghiệp Nông nghiệp & Hiệp hội ngành hàng: Lập kế hoạch phát triển vùng nguyên liệu cam, chè bền vững, chủ động các giải pháp trữ nước và phòng chống bão mặn.
  • Cộng đồng Nông dân & Hợp tác xã: Được hướng dẫn kỹ thuật chuyển đổi lịch mùa vụ, kỹ thuật canh tác tiết kiệm nước và lựa chọn giống chống chịu rủi ro thiên tai.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Luận án đã mở rộng Lý thuyết Ricardian kinh điển (Mendelsohn et al., 1994) bằng cách cấu trúc hóa đồng thời vector cú sốc thiên tai cực đoan cấp tính (Chỉ số hạn Sazonov $Sa$, Nồng độ mặn, Cường độ bão) tương tác cùng vector khí hậu dài hạn vào hàm doanh thu ròng nông hộ, giải quyết triệt để khiếm khuyết xem nhẹ thiên tai ngắn hạn của mô hình Ricardian truyền thống.

2. Đổi mới phương pháp luận so với các nghiên cứu trước đây tại Việt Nam?

So với các nghiên cứu chỉ dùng thống kê mô tả (Phan Thị Cẩm Hằng, 2015; Nguyễn Hoàng Sơn, 2020) hoặc chỉ dùng mô hình OLS/Ricardian với biến khí hậu đơn thuần (Nguyen Thi Lan Huong, 2018; Trần Đại Nghĩa, 2015), luận án đã tích hợp tam giác hóa dữ liệu: kết hợp đo lường thủy văn chuyên sâu (chỉ số khô hạn $Sa$, độ mặn trạm đo) với dữ liệu vi mô hộ phân tầng theo 3 loại cây trồng chủ lực (Lúa, Cam, Chè), thực hiện kiểm định sai số dự báo thực nghiệm và chạy mô hình hồi quy đa biến có kiểm soát nội sinh.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Đó là sự phân hóa đối nghịch sâu sắc về mức độ tổn thương giữa cây hàng năm (Lúa) và cây lâu năm (Cam, Chè): Cây lúa chịu tổn thương nặng nề nhất bởi hạn hán và xâm nhập mặn ở các giai đoạn sinh trưởng nhạy cảm, nhưng có khả năng phục hồi mùa vụ nhanh hơn; trong khi cây cam và cây chè chịu rủi ro đứt gãy dòng thu nhập dài hạn nghiêm trọng nhất khi gặp bão lớn kèm mưa cực đoan do chi phí tái thiết vườn cây và chu kỳ kiến thiết cơ bản kéo dài nhiều năm.

4. Giao thức nhân bản (Replication protocol) có được cung cấp rõ ràng không?

Hoàn toàn rõ ràng. Luận án cung cấp chi tiết định nghĩa biến số, phương pháp tính toán chỉ số Sazonov $Sa$, công thức kinh tế lượng Ricardian và OLS, quy trình điều tra chọn mẫu hộ, các lệnh kiểm định trên Stata và bảng đối soát sai số dự báo so với số liệu thực tế tại Nghệ An giai đoạn 2001-2020.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?

Định hướng tập trung vào: (i) Thiết lập hệ thống giám sát dữ liệu bảng nông hộ vi mô theo thời gian thực kết hợp cảm biến IoT quan trắc hạn mặn nội đồng; (ii) Tích hợp mô hình Ricardian với mô hình Thủy văn - Thủy lực lưu vực sông Lam và mô hình CGE toàn nền kinh tế; (iii) Đánh giá hiệu quả kinh tế - xã hội của chính sách bảo hiểm nông nghiệp dựa trên chỉ số thời tiết tại Bắc Trung Bộ.

Kết luận

  1. Luận án đã hệ thống hóa và làm sâu sắc thêm cơ sở lý luận về lượng hóa thiệt hại kinh tế do tổ hợp thiên tai (bão, hạn, mặn) và BĐKH gây ra cho ngành trồng trọt.
  2. Vận dụng thành công mô hình kinh tế lượng Ricardian cắt ngang kết hợp mô hình OLS đa biến trên bộ dữ liệu khảo sát sơ cấp hộ nông dân và chuỗi số liệu thứ cấp 2001-2020 tại tỉnh Nghệ An.
  3. Lượng hóa chính xác tác động tiêu cực phi tuyến của nhiệt độ cực trị, chỉ số hạn hán $Sa$, độ mặn và bão; đồng thời chứng minh hiệu quả giảm xóc kinh tế có ý nghĩa của các biện pháp thích ứng nông hộ đối với 3 loại cây trồng chủ lực (lúa, cam, chè).
  4. Dự báo chi tiết quy mô thiệt hại doanh thu trồng trọt theo các kịch bản BĐKH của Bộ Tài nguyên và Môi trường (RCP 4.5 và RCP 8.5) đến các mốc thời gian tương lai, tạo căn cứ định lượng vững chắc cho công tác quản lý rủi ro.
  5. Đề xuất hệ thống giải pháp toàn diện, khả thi gồm 7 nhóm chính sách công trình và phi công trình, góp phần nâng cao năng lực ứng phó thiên tai, bảo đảm an ninh lương thực và phát triển nông nghiệp bền vững cho tỉnh Nghệ An cũng như khu vực duyên hải miền Trung.