Tổng quan luận án

Luận án tiến sĩ "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" được nghiên cứu sinh Nguyễn Đăng Khoa thực hiện dưới sự hướng dẫn khoa học của PGS. Trần Thị Kim Thu, bảo vệ tại Trường Đại học Kinh tế Quốc dân năm 2023, thuộc chuyên ngành Kinh tế học (Thống kê kinh tế), mã số 9310101.

+-----------------------------------------------------------------------------------+
| THÔNG TIN TỔNG QUAN LUẬN ÁN                                                       |
+----------------------+------------------------------------------------------------+
| Đề tài               | 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                      |
| Nghiên cứu sinh      | Nguyễn Đăng Khoa                                           |
| Người hướng dẫn      | PGS. Trần Thị Kim Thu                                      |
| Cơ sở đào tạo        | Trường Đại học Kinh tế Quốc dân                            |
| Chuyên ngành         | Kinh tế học (Thống kê kinh tế)                             |
| Mã số chuyên ngành   | 9310101                                                    |
| Năm hoàn thành       | 2023                                                       |
+----------------------+------------------------------------------------------------+

Tính cấp thiết và khoảng trống nghiên cứu

Biến đổi khí hậu (BĐKH) đang là thách thức toàn cầu đe dọa an ninh lương thực và hoạt động sinh kế, trong đó nông nghiệp là ngành dễ bị tổn thương nhất do phụ thuộc trực tiếp vào điều kiện tự nhiên (Adams và cộng sự, 1998; IPCC, 2014). Tại Việt Nam, lúa gạo đóng vai trò trụ cột trong an ninh lương thực quốc gia, chiếm khoảng 30% tổng giá trị sản xuất nông nghiệp và 86% sản lượng lương thực có hạt (FAO, 2018; MARD, 2020). Khu vực miền Bắc có trên 70% diện tích đất cây hàng năm dành cho trồng lúa với hơn 50% dân cư nông thôn tham gia canh tác (GSO, 2018).

Tuy nhiên, tốc độ tăng năng suất lúa tại miền Bắc có xu hướng chững lại hoặc suy giảm trong giai đoạn gần đây. Miền Bắc sở hữu địa hình phức tạp, phân hóa mùa vụ sâu sắc và chịu tác động của các hiện tượng thời tiết cực đoan. Theo kịch bản RCP 4.5 của Bộ Tài nguyên và Môi trường (MONRE, 2020), nhiệt độ trung bình tại miền Bắc dự báo tăng khoảng 1,3°C vào giữa thế kỷ 21, đi kèm sự gia tăng của các cực trị thời tiết.

Khoảng trống nghiên cứu được chỉ ra: Các công trình thực nghiệm về tác động của BĐKH đến sản xuất lúa tại Việt Nam phần lớn tập trung vào khu vực phía Nam, đặc biệt là Đồng bằng sông Cửu Long (ĐBSCL) với nguy cơ hạn hán và xâm nhập mặn. Các nghiên cứu hiện có chủ yếu tiếp cận qua các biến số trung bình (nhiệt độ trung bình, tổng lượng mưa), để trống mảng nghiên cứu định lượng phân tích tác động của các giá trị cực trị (nhiệt độ tối đa, nhiệt độ tối thiểu), các hiện tượng thời tiết cực đoan (rét đậm, rét hại, nắng nóng gay gắt kéo dài) và thiên tai (bão, lũ lụt) theo từng mùa vụ sản xuất cụ thể tại miền Bắc.

Mục tiêu nghiên cứu

Mục tiêu chung của luận án là đánh giá ảnh hưởng của BĐKH đến năng suất lúa ở miền Bắc Việt Nam xét trên nhiều bình diện của khí hậu (trung bình, cực trị, tính cực đoan và thiên tai), từ đó cung cấp bằng chứng định lượng phục vụ công tác dự báo và hoạch định chính sách thích ứng cho sản xuất lúa.

Hệ thống câu hỏi nghiên cứu được xác định gồm 4 câu hỏi:

  1. Thực trạng BĐKH ở miền Bắc Việt Nam diễn ra như thế nào? Những nguy cơ của nó đối với sản xuất lúa trong từng mùa vụ sản xuất lúa?
  2. Biến đổi về nhiệt độ ở các khía cạnh trung bình, cực đại, cực tiểu có ảnh hưởng như thế nào lên năng suất lúa trong mỗi vụ sản xuất lúa ở miền Bắc Việt Nam?
  3. Biến đổi về lượng mưa ở các khía cạnh trung bình, cực đại, cực tiểu có ảnh hưởng như thế nào lên năng suất lúa trong mỗi vụ sản xuất ở miền Bắc Việt Nam?
  4. Biến đổi về tần số các hiện tượng thời tiết cực đoan (nắng nóng gay gắt, rét đậm - rét hại) và thiên tai có ảnh hưởng như thế nào đến năng suất lúa trong từng vụ sản xuất lúa ở miền Bắc Việt Nam?

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Tác động của BĐKH đến năng suất lúa ở khu vực miền Bắc Việt Nam.
  • Phạm vi nội dung: Tập trung vào chỉ tiêu năng suất lúa hiện vật tính theo diện tích gieo trồng (đơn vị: tạ/ha). Các yếu tố khí hậu khảo sát gồm nhiệt độ (trung bình, tối đa, tối thiểu), lượng mưa (tổng lượng mưa, tối đa, tối thiểu), thời tiết cực đoan (số đợt, số ngày rét đậm - rét hại, nắng nóng gay gắt) và thiên tai (bão từ cấp 10 trở lên, lũ lụt). Luận án không đi sâu vào xâm nhập mặn và nước biển dâng do đây không phải nguy cơ phổ biến trên diện rộng tại miền Bắc.
  • Phạm vi không gian: 25 tỉnh/thành phố thuộc 3 tiểu vùng khí hậu miền Bắc gồm Đồng bằng sông Hồng (ĐBSH), Miền núi Tây Bắc và Trung du - miền núi Đông Bắc (gọi chung là Trung du và miền núi phía Bắc - TD-MNPB).
  • Phạm vi thời gian: Giai đoạn quan sát chính từ năm 2001 đến năm 2020 (chuỗi số liệu 20 năm), kết hợp dữ liệu năm 2021 để dự đoán thống kê.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trên thế giới và tại Việt Nam

Tổng quan tài liệu của luận án hệ thống hóa các công trình theo các nhóm chủ đề:

Hướng nghiên cứu Tác giả và năm công bố Phát hiện chính được trích dẫn
Tác động nông nghiệp phạm vi toàn cầu Rosenzweig và Parry (1994); Darwin và cộng sự (1995); Mendelsohn và Dinar (2009); Nelson và cộng sự (2014); Burke và Emerick (2016); Altieri và Nicholl (2017); Ding (2019) BĐKH tác động không đồng đều giữa các vùng; khu vực nhiệt đới và các nước đang phát triển chịu tổn thương nặng hơn; tồn tại tính phức tạp theo vùng miền và loại cây trồng.
Kênh tác động của nhiệt độ tới lúa Yoshida (1978, 1981); Matsui và Horie (1992); Mendelsohn và cộng sự (1994, 2001); Welch và cộng sự (2010); IRRI (2012); Sanchez và cộng sự (2013); Ahmad và cộng sự (2014); Fridaus và cộng sự (2014, 2020); Haque và cộng sự (2016) Ngưỡng nhiệt thích hợp của lúa từ 20°C - 32°C; trên 30°C vào ban đêm làm giảm 10% năng suất mỗi độ tăng thêm (IRRI, 2012); tồn tại quan hệ phi tuyến giữa nhiệt độ và năng suất lúa.
Kênh tác động của lượng mưa và thiên tai Moomaw và Vergara (1965); Basneyake và cộng sự (2006); Auffhammer và cộng sự (2012); Koc và Ceylan (2013) Lượng mưa tác động tùy thời điểm sinh trưởng; mưa lớn giai đoạn trổ bông làm giảm thụ phấn; bão, hạn hán, lũ lụt làm giảm trực tiếp sản lượng thu hoạch.
Nghiên cứu thực nghiệm tại Việt Nam The (2009); Yu và cộng sự (2010); Nguyen (2011); The và cộng sự (2013); Đặng Thị Thanh Hoa (2013); Balica và cộng sự (2014); Lê Anh Tuấn và cộng sự (2014); Nguyen và Hoang (2015); Le (2016); Trinh (2018); Hương và cộng sự (2020) Phần lớn nghiên cứu thực hiện tại ĐBSCL hoặc tính toán chung toàn quốc; Le (2016) cho thấy lượng mưa có quan hệ thuận với năng suất lúa vụ Đông Xuân và Hè Thu, trái ngược với phát hiện giảm sản lượng của Nguyen (2011).

Khoảng trống luận án giải quyết

Luận án định vị trọng tâm giải quyết ba khoảng trống cụ thể:

  1. Thiếu các công trình nghiên cứu định lượng chuyên sâu cho khu vực miền Bắc - nơi có chế độ nhiệt và phân hóa mùa vụ khác biệt so với phía Nam.
  2. Thiếu sự kết hợp đồng thời giữa các đại lượng trung bình với các cực trị thời tiết (nhiệt độ cao nhất/thấp nhất, lượng mưa cực đại) và các chỉ số thời tiết cực đoan (rét đậm - rét hại, nắng nóng kéo dài, bão cấp cao).
  3. Thiếu mô hình hóa tách biệt cho từng mùa vụ canh tác (vụ Đông Xuân và vụ Mùa) trên nền tảng dữ liệu mảng cấp tỉnh dài hạn.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý luận và các kênh tác động

Cơ sở sinh học của hoạt động canh tác lúa trải qua ba giai đoạn chính:

  • Giai đoạn sinh dưỡng: Từ khi nảy mầm đến khi phân hóa đòng;
  • Giai đoạn sinh thực: Từ khi hình thành đòng đến khi trổ bông và thụ phấn;
  • Giai đoạn chín: Từ khi thụ phấn thành công đến khi thu hoạch (kéo dài khoảng 30 ngày ở vùng nhiệt đới).

Năng suất lúa chịu sự chi phối của yếu tố giống (di truyền, chiếm khoảng 40% tiềm năng năng suất), đất đai, phân bón, kỹ thuật canh tác và các yếu tố khí hậu. Luận án phân định hai công thức đo lường năng suất dạng hiện vật:

$$\text{Năng suất gieo trồng} = \frac{\text{Sản lượng thu hoạch (tạ)}}{\text{Tổng diện tích gieo trồng (ha)}}$$

$$\text{Năng suất thu hoạch} = \frac{\text{Sản lượng thu hoạch (tạ)}}{\text{Tổng diện tích thu hoạch thực tế (ha)}}$$

Luận án lựa chọn biến phụ thuộc là năng suất lúa tính trên diện tích gieo trồng (tạ/ha) nhằm phản ánh chân thực hiệu quả sử dụng đất trước tác động khí hậu, loại trừ biến động giá cả thị trường.

+-----------------------------------------------------------------------------------+
| HỆ THỐNG CHỈ TIÊU THỐNG KÊ ĐO LƯỜNG KHÍ HẬU TRONG NGHIÊN CỨU                      |
+---------------------+-------------------------------+-----------------------------+
| Nhóm yếu tố         | Chỉ tiêu thống kê             | Phạm vi thời gian xác định  |
+---------------------+-------------------------------+-----------------------------+
| Nhiệt độ            | - Nhiệt độ trung bình         | Tháng, vụ, năm              |
|                     | - Nhiệt độ cao nhất           | Tháng, vụ, năm              |
|                     | - Nhiệt độ thấp nhất          | Tháng, vụ, năm              |
+---------------------+-------------------------------+-----------------------------+
| Lượng mưa           | - Tổng lượng mưa              | Tháng, vụ, năm              |
|                     | - Lượng mưa trung bình        | Tháng, vụ, năm              |
|                     | - Lượng mưa cao nhất          | Tháng, vụ                   |
|                     | - Lượng mưa thấp nhất         | Tháng, vụ                   |
+---------------------+-------------------------------+-----------------------------+
| Thời tiết cực đoan  | - Số đợt rét đậm, rét hại     | Vụ Đông Xuân, năm           |
| và thiên tai        | - Tổng số ngày rét đậm, hại   | Vụ Đông Xuân, năm           |
|                     | - Số đợt nắng nóng kéo dài    | Vụ Mùa, năm                 |
|                     | - Tổng số ngày nắng nóng      | Vụ Mùa, năm                 |
|                     | - Số cơn bão (từ cấp 10+)     | Vụ Mùa, năm                 |
|                     | - Số đợt lũ, lụt              | Vụ Mùa, năm                 |
+---------------------+-------------------------------+-----------------------------+

Quy chuẩn kỹ thuật của các hiện tượng cực đoan:

  • Rét đậm: Nhiệt độ trung bình ngày $T_{tb} < 15^\circ\text{C}$.
  • Rét hại: Nhiệt độ trung bình ngày $T_{tb} \le 13^\circ\text{C}$. Đợt rét đậm - rét hại xác định khi kéo dài liên tục từ 3 ngày trở lên.
  • Nắng nóng: Nhiệt độ cao nhất ngày $T_{max} \ge 35^\circ\text{C}$ ($35^\circ\text{C} \le T_{max} < 37^\circ\text{C}$ là nắng nóng; $37^\circ\text{C} \le T_{max} < 39^\circ\text{C}$ là nắng nóng gay gắt; $T_{max} \ge 39^\circ\text{C}$ là đặc biệt gay gắt). Đợt diện rộng khi có ít nhất 2/3 số trạm đạt $T_{max} \ge 35^\circ\text{C}$ và 1/2 số trạm đạt $T_{max} \ge 37^\circ\text{C}$, kéo dài $\ge 3$ ngày.

Phương pháp nghiên cứu định lượng

  • Phương pháp thu thập dữ liệu: Thu thập dữ liệu thứ cấp từ trạm quan trắc của Tổng cục Khí tượng Thủy văn (Bộ TN&MT), số liệu thiên tai từ UNDRR (https://www.preventionweb.net/), số liệu diện tích, sản lượng, năng suất và lao động từ Tổng cục Thống kê (GSO), Niên giám Thống kê tỉnh, Khảo sát Mức sống Dân cư (VHLSS) và Tổng điều tra nông thôn, nông nghiệp và thủy sản (2000 - 2020).
  • Cấu trúc bộ dữ liệu: Dữ liệu mảng cân bằng cấp tỉnh gồm 25 tỉnh/thành phố miền Bắc trong 20 năm (2001 - 2020), tổng cộng 500 quan sát.
  • Phương pháp phân tích:
    1. Thống kê mô tả: Phân tích giá trị trung tâm, độ biến thiên, tần suất phân phối bằng bảng biểu thống kê.
    2. Phân tích dãy số thời gian: Xác định đường xu thế tuyến tính $Y = b_0 + b_1 t$ của các biến khí hậu và năng suất lúa qua 20 năm.
    3. Phân tích hồi quy dữ liệu mảng: Sử dụng mô hình Tác động cố định (Fixed Effects - FE), kiểm định và khắc phục khuyết tật mô hình (phương sai sai số thay đổi, tự tương quan) bằng sai số chuẩn mạnh (Robust Standard Errors), ước lượng PCSE (Panel Corrected Standard Errors) và FGLS (Feasible Generalized Least Squares).
    4. Dự đoán thống kê: Sử dụng mô hình ước lượng để dự báo năng suất lúa vụ Đông Xuân và vụ Mùa cho 25 tỉnh năm 2021.

Nội dung chính theo từng chương

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

Chương 1 hệ thống hóa khái niệm về khí quyển, thời tiết, khí hậu và BĐKH. Tác giả phân tích nguyên nhân tự nhiên (biến đổi quỹ đạo Trái Đất: độ lệch tâm, độ nghiêng trục, tiến động; hoạt động núi lửa; phát xạ mặt trời) và nguyên nhân nhân tạo (phát thải khí nhà kính $\text{CO}_2$, $\text{CH}_4$, $\text{N}_2\text{O}$ từ năng lượng hóa thạch và nông nghiệp).

Chương này luận giải đặc điểm sinh trưởng của cây lúa qua 3 giai đoạn và xác lập các kênh tác động vật lý - sinh học:

  • Nhiệt độ tác động trực tiếp đến hô hấp, quang hợp, chuyển hóa dinh dưỡng và tốc độ sinh trưởng. Ngưỡng nhiệt vượt quá mức tối ưu ($>32^\circ\text{C}$) gây thoái hóa hạt đòng, giảm tỷ lệ thụ phấn.
  • Lượng mưa quyết định độ ẩm đất và cung cấp nước tưới, nhưng mưa cực đoan vào giai đoạn thụ phấn hoặc thu hoạch gây thất thoát lớn.
  • Khí hậu cực đoan và thiên tai làm gián đoạn sinh trưởng, gây chết cây con (rét hại) hoặc đổ rạp, úng ngập (bão, lũ).

Chương 1 khép lại bằng việc tổng quan các công trình nghiên cứu trong nước và quốc tế, chỉ ra khoảng trống về dữ liệu và phương pháp, từ đó xây dựng khung phân tích cho luận án.

Chương 2: Phương pháp nghiên cứu

Chương 2 trình bày quy trình thiết kế nghiên cứu kinh tế lượng thực nghiệm:

  • Xây dựng dạng hàm sản xuất nông nghiệp mở rộng có chứa các yếu tố khí hậu:

$$\text{Yield}{it} = \alpha_i + \beta \mathbf{Climate}{it} + \gamma \mathbf{Control}{it} + \varepsilon{it}$$

  • Trong đó, nhóm biến khí hậu ($\mathbf{Climate}$) gồm nhiệt độ trung bình, nhiệt độ cực đại, nhiệt độ cực tiểu, lượng mưa, cực trị mưa, số ngày rét đậm - rét hại (vụ Đông Xuân), số ngày nắng nóng kéo dài và số cơn bão/đợt lũ (vụ Mùa).
  • Luận án đưa thêm các số hạng bậc hai (bình phương của nhiệt độ cực đại, nhiệt độ trung bình) nhằm kiểm định giả thuyết về mối quan hệ phi tuyến (hình chữ U ngược).
  • Nhóm biến kiểm soát ($\mathbf{Control}$) bao gồm các biến đại diện cho trình độ canh tác, lao động nông nghiệp đã qua đào tạo, diện tích gieo trồng và đặc điểm vùng.
  • Quy trình ước lượng giải trình chi tiết việc thu thập 500 quan sát cấp tỉnh giai đoạn 2001 - 2020 và kiểm định lựa chọn mô hình phù hợp (FE so với RE, kiểm định Wooldridge về tự tương quan, kiểm định Modified Wald về phương sai thay đổi).

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

Chương 3 tổng quan về điều kiện tự nhiên, kinh tế - xã hội và canh tác lúa tại 2 tiểu vùng: ĐBSH và TD-MNPB.

+-----------------------------------------------------------------------------------+
| TỔNG HỢP DIỄN BIẾN SẢN XUẤT LÚA VÀ KHÍ HẬU MIỀN BẮC (2001 - 2020)                |
+----------------------+------------------------------------------------------------+
| Chỉ tiêu             | Thực trạng và xu thế ghi nhận giai đoạn 2001 - 2020        |
+----------------------+------------------------------------------------------------+
| Diện tích gieo trồng | Xu hướng giảm dần do chuyển đổi mục đích sử dụng đất       |
|                      | sang công nghiệp, đô thị hóa và cơ sở hạ tầng.             |
+----------------------+------------------------------------------------------------+
| Sản lượng và         | Năng suất trung bình tăng nhưng tốc độ tăng chậm dần;     |
| Năng suất            | Năng suất vụ Đông Xuân cao hơn đáng kể so với vụ Mùa.      |
|                      | ĐBSH có năng suất bình quân cao hơn TD-MNPB.               |
+----------------------+------------------------------------------------------------+
| Nhiệt độ             | Nhiệt độ trung bình năm và các vụ có xu thế tăng nhẹ;      |
|                      | Gia tăng biên độ dao động nhiệt giữa các năm.              |
+----------------------+------------------------------------------------------------+
| Lượng mưa            | Phân bố không đều theo không gian và thời gian;            |
|                      | Lượng mưa tập trung cao điểm vào các tháng giữa vụ Mùa.    |
+----------------------+------------------------------------------------------------+
| Thiên tai &          | Xuất hiện các đợt rét đậm - rét hại kéo dài ở vụ Đông Xuân |
| Hiện tượng cực đoan  | (đặc biệt tại TD-MNPB); Nắng nóng gay gắt và bão lũ        |
|                      | tập trung nhiều ở vụ Mùa, gây thiệt hại cục bộ.            |
+----------------------+------------------------------------------------------------+

Chương này cũng rà soát thực trạng chính sách thích ứng BĐKH trong sản xuất nông nghiệp tại các địa phương, chỉ ra hạn chế về hạ tầng thủy lợi và chuyển đổi cơ cấu giống.

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

Chương 4 là chương trọng tâm chứa kết quả ước lượng thực nghiệm:

  • Thống kê mô tả: Xác lập bảng tham số mô tả (trung bình, độ lệch chuẩn, giá trị lớn nhất, giá trị nhỏ nhất) cho các biến trong mô hình vụ Đông Xuân và vụ Mùa cho toàn miền Bắc và tách riêng 2 tiểu vùng (ĐBSH và TD-MNPB).
+-----------------------------------------------------------------------------------+
| SO SÁNH TÁC ĐỘNG CỦA CÁC YẾU TỐ KHÍ HẬU GIỮA HAI MÙA VỤ                           |
+--------------------+--------------------------------+-----------------------------+
| Yếu tố khí hậu     | Vụ Đông Xuân                   | Vụ Mùa                      |
+--------------------+--------------------------------+-----------------------------+
| Nhiệt độ trung bình| Tác động tích cực ở mức độ vừa | Tăng nhiệt độ trung bình có |
|                    | phải (giúp giảm lạnh đầu vụ).  | xu hướng tác động bất lợi.  |
+--------------------+--------------------------------+-----------------------------+
| Nhiệt độ cực đại   | Quan hệ phi tuyến chữ U ngược; | Quan hệ phi tuyến chữ U     |
| ($T_{max}$)        | Vượt ngưỡng tối ưu gây giảm    | ngược rõ nét; Nắng nóng     |
|                    | năng suất cuối vụ.             | gay gắt làm giảm đỗ hạt.    |
+--------------------+--------------------------------+-----------------------------+
| Lượng mưa          | Thiếu mưa đầu vụ gây hạn; Mưa  | Lượng mưa quá lớn gây ngập  |
|                    | vừa phải có lợi cho đẻ nhánh.  | úng; Mưa cực đại giảm mùa.  |
+--------------------+--------------------------------+-----------------------------+
| Thiên tai &        | Rét đậm - rét hại kéo dài làm  | Bão từ cấp 10 trở lên và lũ |
| Hiện tượng cực đoan| chết mạ, kéo dài thời gian vụ. | lụt làm giảm mạnh năng suất.|
+--------------------+--------------------------------+-----------------------------+
  • Kết quả ước lượng:
    • Mô hình FE hiệu chỉnh sai số chuẩn khẳng định tác động phi tuyến của nhiệt độ cực đại lên năng suất lúa ở cả hai vụ. Hệ số của $T_{max}$ mang dấu dương trong khi hệ số của $(T_{max})^2$ mang dấu âm có ý nghĩa thống kê.
    • Vụ Đông Xuân: Hiện tượng rét đậm - rét hại kéo dài $\ge 3$ ngày có hệ số hồi quy âm, phản ánh tác động tiêu cực rõ nét đến giai đoạn sinh dưỡng đầu vụ.
    • Vụ Mùa: Tần suất các đợt nắng nóng gay gắt và số lượng cơn bão đổ bộ (cấp 10 trở lên) có mối liên hệ nghịch chiều với năng suất thu hoạch.
  • Dự đoán thống kê năm 2021: Ứng dụng mô hình ước lượng để tính toán năng suất lúa dự báo cho 25 tỉnh thành miền Bắc trong năm 2021, so sánh với số liệu thống kê thực tế để kiểm tra độ tin cậy của phương trình hồi quy.

Kết quả và những đóng góp mới

Đóng góp mới về mặt học thuật và lý luận

  1. Luận án mở rộng khung phân tích ảnh hưởng của BĐKH bằng cách tích hợp đồng thời các chỉ tiêu trung bình, cực trị (nhiệt độ tối đa, nhiệt độ tối thiểu, lượng mưa cực đại) và các chỉ số phản ánh tính cực đoan của thời tiết (đợt rét đậm - rét hại kéo dài, đợt nắng nóng gay gắt) cùng thiên tai (bão, lũ).
  2. Thiết lập cơ sở lý luận và kiểm định thực nghiệm về sự khác biệt trong cơ chế tác động của khí hậu theo từng mùa vụ canh tác (vụ Đông Xuân và vụ Mùa), chứng minh tính không đồng nhất về phản ứng của cây lúa trước BĐKH giữa các thời điểm trong năm.
  3. Đề xuất quy trình xử lý, tổng hợp và chuẩn hóa dữ liệu khí tượng thủy văn kết hợp dữ liệu kinh tế - xã hội thành bộ dữ liệu mảng cấp tỉnh qua 20 năm, khắc phục hạn chế thiếu hụt dữ liệu dài hạn trong nghiên cứu thống kê nông nghiệp tại miền Bắc.

Những phát hiện và kết luận rút ra từ nghiên cứu

  1. Tác động của BĐKH đến năng suất lúa biểu hiện đa chiều qua các giá trị cực trị và hiện tượng thời tiết cực đoan chứ không thuần túy qua nhiệt độ trung bình hay tổng lượng mưa cả năm.
  2. Tồn tại mối quan hệ phi tuyến dạng hình chữ U ngược giữa nhiệt độ (đặc biệt là nhiệt độ tối đa) và năng suất lúa ở cả hai vụ sản xuất, chứng minh rằng sự gia tăng nhiệt độ vượt ngưỡng sinh học trong dài hạn sẽ gây tổn hại đến năng suất lúa miền Bắc.
  3. Mức độ tác động của các biến khí hậu có sự phân hóa mùa vụ sâu sắc: Vụ Đông Xuân chịu rủi ro lớn từ rét đậm - rét hại đầu vụ và biến động nhiệt độ tối thiểu; vụ Mùa chịu tổn thương nghiêm trọng từ các đợt nắng nóng gay gắt, mưa lớn tập trung và bão lũ.
  4. Thiên tai và hiện tượng thời tiết cực đoan có tác động tiêu cực tức thời và đáng kể đến năng suất gieo trồng, đặc biệt ở các địa phương ven biển và miền núi có cơ sở hạ tầng chống chịu kém.

Kiến nghị và giải pháp đề xuất

  • Quy hoạch mùa vụ và cơ cấu giống: Chuyển dịch lịch gieo cấy vụ Đông Xuân để tránh các đợt rét đậm - rét hại cực đoan; đưa vào canh tác các giống lúa thuần ngắn ngày, giống chịu nhiệt và chịu úng cho vụ Mùa.
  • Nâng cấp cơ sở hạ tầng: Tăng cường năng lực tiêu thoát nước cho vùng ĐBSH trong vụ Mùa và hệ thống trữ nước tưới chống hạn cho vùng TD-MNPB trong vụ Đông Xuân.
  • Hoàn thiện công tác dự báo và cảnh báo: Xây dựng hệ thống cảnh báo sớm rủi ro thời tiết nông vụ cấp địa phương dựa trên dữ liệu khí tượng quan trắc.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

  • Dữ liệu sử dụng dừng lại ở cấp tỉnh (macro-level panel data), chưa phân rã được đến cấp huyện, xã hoặc số liệu vi mô cấp hộ nông dân, do đó chưa kiểm soát hoàn toàn các sai biệt về tập quán canh tác cục bộ.
  • Chuỗi dữ liệu quan sát dừng ở mốc năm 2020 (với dự đoán năm 2021), chưa cập nhật các biến động khí hậu bất thường giai đoạn sau năm 2021.
  • Mô hình tập trung vào chỉ tiêu năng suất hiện vật (tạ/ha) mà chưa lượng hóa đầy đủ các khía cạnh về tổn thất giá trị kinh tế, chi phí phòng chống thiên tai và chất lượng hạt lúa.

Hướng nghiên cứu tiếp

  • Mở rộng phân tích sử dụng dữ liệu vi mô cấp hộ (kết hợp các kỳ điều tra VHLSS chuyên sâu) để đánh giá năng lực thích ứng và hành vi kinh tế của nông dân.
  • Bổ sung các biến số vi khí hậu như độ ẩm tương đối, bức xạ mặt trời và tốc độ gió khi hệ thống trạm quan trắc tự động cung cấp dữ liệu liên tục.
  • Ứng dụng các kịch bản BĐKH mới nhất để mô phỏng định lượng biến động năng suất lúa đến năm 2030 và 2050 theo từng tiểu vùng sinh thái nông nghiệp.

Giá trị tham khảo

Luận án cung cấp nguồn tư liệu học thuật và thực tiễn cho nhiều nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học chuyên ngành Kinh tế học, Thống kê kinh tế, Kinh tế nông nghiệp: Tham khảo quy trình thiết kế mô hình dữ liệu mảng (FE, Robust SE, PCSE) trong xử lý biến khí hậu và phương pháp xây dựng bộ dữ liệu thứ cấp liên ngành.
  • Giảng viên và nhà nghiên cứu kinh tế lượng ứng dụng: Sử dụng cấu trúc mô hình phi tuyến và phương pháp kiểm định tác động của BĐKH làm tài liệu giảng dạy và nghiên cứu thực nghiệm.
  • Cơ quan quản lý nhà nước và hoạch định chính sách (Bộ Nông nghiệp và Phát triển Nông thôn, Bộ Tài nguyên và Môi trường, Sở NN&PTNT 25 tỉnh miền Bắc): Sử dụng các bằng chứng định lượng về tác động mùa vụ và ngưỡng nhiệt để xây dựng khung lịch thời vụ, quy hoạch cơ cấu giống lúa và phân bổ ngân sách đầu tư thủy lợi thích ứng BĐKH.

Câu hỏi thường gặp

1. Tại sao luận án sử dụng chỉ tiêu năng suất hiện vật (tạ/ha) thay vì chỉ tiêu giá trị (doanh thu/ha)?

Luận án giải thích rằng chỉ tiêu giá trị (doanh thu, thu nhập trên 1 ha) chịu sự chi phối mạnh mẽ bởi biến động cung cầu và giá bán lúa gạo trên thị trường. Trong những năm khí hậu khắc nghiệt làm giảm sản lượng, nguồn cung khan hiếm có thể đẩy giá bán lên cao, khiến thu nhập trên 1 ha không đổi hoặc tăng giả tạo. Do đó, chỉ tiêu năng suất hiện vật (tạ/ha) phản ánh trực tiếp và chính xác hơn phản ứng sinh học của cây lúa và hiệu quả sử dụng đất trước tác động của thời tiết.

2. Dữ liệu mảng của luận án được xây dựng từ những nguồn nào và quy mô quan sát là bao nhiêu?

Dữ liệu mảng được tổng hợp từ hệ thống trạm quan trắc của Tổng cục Khí tượng Thủy văn (Bộ TN&MT), cổng dữ liệu thiên tai UNDRR, Niên giám Thống kê tỉnh của Tổng cục Thống kê (GSO), điều tra VHLSS và Tổng điều tra nông thôn, nông nghiệp, thủy sản. Quy mô bộ dữ liệu gồm 25 tỉnh/thành phố miền Bắc trong 20 năm (2001 - 2020), đạt 500 quan sát.

3. Mối quan hệ giữa nhiệt độ cực đại ($T_{max}$) và năng suất lúa ở miền Bắc có đặc điểm gì?

Kết quả ước lượng thực nghiệm chỉ ra mối quan hệ phi tuyến dạng hình chữ U ngược giữa nhiệt độ tối đa và năng suất lúa ở cả hai vụ canh tác. Khi nhiệt độ tối đa tăng trong giới hạn sinh học, năng suất có thể được thúc đẩy; tuy nhiên khi vượt qua ngưỡng tối ưu, nhiệt độ cao gay gắt làm gián đoạn quá trình quang hợp, thụ phấn và làm hạt lép, dẫn đến suy giảm năng suất nghiêm trọng.

4. Luận án định nghĩa thế nào là đợt rét đậm - rét hại và nắng nóng gay gắt?

  • Rét đậm - rét hại: Rét đậm xảy ra khi $T_{tb} < 15^\circ\text{C}$, rét hại khi $T_{tb} \le 13^\circ\text{C}$. Đợt rét đậm - rét hại được tính khi hiện tượng này kéo dài liên tục từ 3 ngày trở lên.
  • Nắng nóng gay gắt: Xảy ra khi nhiệt độ cao nhất ngày đạt $37^\circ\text{C} \le T_{max} < 39^\circ\text{C}$ (và đặc biệt gay gắt khi $T_{max} \ge 39^\circ\text{C}$). Đợt diện rộng được xác định khi có ít nhất 2/3 số trạm trong khu vực đạt $T_{max} \ge 35^\circ\text{C}$ và 1/2 số trạm đạt $T_{max} \ge 37^\circ\text{C}$, kéo dài từ 3 ngày liên tục trở lên.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Đăng Khoa đã giải quyết bài toán định lượng về ảnh hưởng của BĐKH đến năng suất lúa tại 25 tỉnh miền Bắc giai đoạn 2001 - 2020 bằng mô hình kinh tế lượng dữ liệu mảng. Nghiên cứu đã chứng minh rằng các yếu tố cực trị thời tiết (nhiệt độ tối đa, tối thiểu), các hiện tượng cực đoan (rét đậm - rét hại, nắng nóng kéo dài) và thiên tai (bão lũ) có tác động phân hóa rõ rệt theo từng mùa vụ sản xuất. Kết quả nghiên cứu cung cấp cơ sở khoa học và bằng chứng thực nghiệm phục vụ việc điều chỉnh lịch thời vụ, cơ cấu giống lúa và hoạch định chính sách nông nghiệp thích ứng với BĐKH tại miền Bắc Việt Nam.