Bối cảnh và vấn đề nghiên cứu

Đồng bằng sông Cửu Long (ĐBSCL) là vùng trọng điểm nông nghiệp của Việt Nam nhưng đang đối mặt với những tác động nghiêm trọng từ biến đổi khí hậu. Trong các hình thái thiên tai, hạn hán và xâm nhập mặn diễn biến phức tạp, đe dọa trực tiếp đến nguồn nước ngọt phục vụ sản xuất và sinh kế dân cư. Bến Tre là tỉnh ven biển với chiều dài bờ biển khoảng 65 km, địa hình trũng thấp và chịu ảnh hưởng trực tiếp của xâm nhập mặn qua ba nhánh sông lớn là sông Cửa Đại, sông Hàm Luông và sông Cổ Chiên. Mùa khô hạn mặn tại đây thường kéo dài từ 6 đến 7 tháng mỗi năm; các kịch bản biến đổi khí hậu dự báo khi mực nước biển dâng thêm 1 m, có tới 50% diện tích các huyện của tỉnh Bến Tre đứng trước nguy cơ ngập mặn.

Cây bưởi da xanh là loại cây ăn trái đặc sản, có nguồn gốc từ huyện Mỏ Cày Bắc, được xác định là cây trồng chủ lực của tỉnh Bến Tre với diện tích canh tác toàn tỉnh đạt 8.824 ha (chiếm khoảng 30% tổng diện tích cây ăn trái). Cây bưởi mang lại giá trị kinh tế cao và phục vụ xuất khẩu, tuy nhiên đây là loại cây mẫn cảm với độ mặn. Thực trạng xâm nhập mặn ngày càng gay gắt tại xã Bình Phú (thành phố Bến Tre) – khu vực cửa ngõ chịu ảnh hưởng mặn từ sông Hàm Luông – đặt ra yêu cầu phải đánh giá định lượng thiệt hại năng suất và năng lực thích ứng của nông hộ.

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

  • Mục tiêu chung: Đánh giá tác động của xâm nhập mặn đến năng suất các nông hộ trồng bưởi tại xã Bình Phú, thành phố Bến Tre, tỉnh Bến Tre.
  • Mục tiêu cụ thể:
    1. Mô tả thực trạng tình hình xâm nhập mặn tại xã Bình Phú, thành phố Bến Tre, tỉnh Bến Tre.
    2. Phân tích mức độ ảnh hưởng của xâm nhập mặn đến các hộ trồng bưởi tại xã Bình Phú, thành phố Bến Tre, tỉnh Bến Tre.
    3. Đánh giá tác động xâm nhập mặn đến hiệu quả sản xuất của nông hộ trồng bưởi tại xã Bình Phú, thành phố Bến Tre, tỉnh Bến Tre.
    4. Đề xuất một số giải pháp nhằm nâng cao khả năng thích ứng với xâm nhập mặn của nông hộ trồng bưởi tại xã Bình Phú, thành phố Bến Tre, tỉnh Bến Tre.

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

  • Đối tượng nghiên cứu: Tác động của xâm nhập mặn đến năng suất của các hộ trồng bưởi da xanh.
  • Đối tượng khảo sát: Các nông hộ trực tiếp canh tác bưởi da xanh tại địa bàn xã Bình Phú.
  • Phạm vi không gian: Xã Bình Phú, thành phố Bến Tre, tỉnh Bến Tre.
  • Phạm vi thời gian: Thu thập số liệu sơ cấp và thứ cấp trong năm 2022 (đánh giá chu kỳ sản xuất năm 2021). Thời gian triển khai nghiên cứu từ tháng 09/2022 đến tháng 02/2023.

Cơ sở lý thuyết và phương pháp

Khóa luận vận dụng lý thuyết kinh tế nông hộ, các khái niệm về xâm nhập mặn và biến đổi khí hậu (theo định nghĩa của Công ước khung Liên Hợp Quốc về Biến đổi khí hậu - UNFCCC). Các chỉ tiêu kinh tế cơ bản được sử dụng gồm: tổng chi phí sản xuất ($TC = CPVC + CPLĐ$), doanh thu ($TR = P \times Q$), lợi nhuận ($\pi = TR - TC$), thu nhập hỗn hợp ($I = \pi + \text{chi phí lao động gia đình}$); cùng hệ thống tỷ suất tài chính: tỷ suất doanh thu/chi phí ($TR/TC$), tỷ suất lợi nhuận/chi phí ($\pi/TC$), và tỷ suất thu nhập/chi phí ($I/TC$).

Phương pháp thu thập dữ liệu

  • Dữ liệu thứ cấp: Thu thập từ các văn bản, báo cáo tổng kết của Ủy ban nhân dân xã Bình Phú, Sở Nông nghiệp và Phát triển nông thôn tỉnh Bến Tre, Đài Khí tượng Thủy văn tỉnh Bến Tre và các tài liệu chuyên ngành đã công bố.
  • Dữ liệu sơ cấp: Thu thập thông qua khảo sát trực tiếp 60 nông hộ trồng bưởi da xanh tại xã Bình Phú bằng phiếu điều tra cấu trúc sẵn.
  • Quy mô chọn mẫu: Áp dụng công thức chọn mẫu hồi quy của Harri (1985): $N = 50 + m$ (trong đó $N$ là cỡ mẫu tối thiểu, $m$ là số biến độc lập). Với $m = 9$ biến độc lập, mẫu nghiên cứu đạt $N = 60$ hộ, thỏa mãn yêu cầu phân tích thực nghiệm.

Phương pháp phân tích và mô hình kinh tế lượng

Đề tài sử dụng phương pháp thống kê mô tả, phương pháp so sánh (tuyệt đối và tương đối) và thang đo Likert 5 mức độ (từ 1: Rất ít nghiêm trọng đến 5: Rất nghiêm trọng) để đo lường mức độ ảnh hưởng của xâm nhập mặn.

Để lượng hóa các nhân tố ảnh hưởng đến năng suất bưởi da xanh, tác giả thiết lập hàm sản xuất Cobb-Douglas và ước lượng bằng phương pháp bình phương nhỏ nhất (OLS) ở dạng log-log qua phần mềm EViews 8:

$$\ln(Y) = \beta_0 + \beta_1 \ln(X_1) + \beta_2 \ln(X_2) + \beta_3 \ln(X_3) + \beta_4 \ln(X_4) + \beta_5 \ln(X_5) + \beta_6 \ln(X_6) + \beta_7 \ln(X_7) + \beta_8 X_8 + \beta_9 D_1$$

Tên biến Ký hiệu Đơn vị tính Kỳ vọng dấu Cơ sở kế thừa
Năng suất bưởi (Biến phụ thuộc) $Y$ kg/1.000 m² Khung hàm sản xuất
Trình độ học vấn của chủ hộ $X_1$ Năm $(+)$ Nguyễn Thị Tuyết Nhung (2016)
Kinh nghiệm canh tác bưởi $X_2$ Năm $(+)$ Nguyễn Thị Tuyết Nhung (2016), Lê Anh Khoa (2018)
Lượng nước tưới $X_3$ m³/1.000 m²/năm $(-)$ Mai Văn Xuân & Nguyễn Văn Hóa (2011), Trương Quỳnh Hoa và cs (2022)
Lượng phân bón hữu cơ $X_4$ kg/1.000 m² $(+)$ Đặng Thị Tuyết Ngân (2015), Nguyễn Thị Tuyết Nhung (2016)
Lượng phân bón vô cơ $X_5$ kg/1.000 m² $(+)$ Đặng Thị Tuyết Ngân (2015), Trương Quỳnh Hoa và cs (2022)
Chi phí thuốc bảo vệ thực vật $X_6$ 1.000 đồng/1.000 m² $(+)$ Đặng Thị Tuyết Ngân (2015), Trương Quỳnh Hoa và cs (2022)
Công lao động (nhà và thuê) $X_7$ Giờ/năm $(+)$ Quãng Trọng Chinh (2013), Nguyễn Thị Tươi & Nguyễn Phú Sơn (2021)
Mức độ nhiễm mặn của vườn $X_8$ Thang đo 1 - 5 $(-)$ Quãng Trọng Chinh (2013), Lê Anh Khoa (2018)
Khuyến nông (Biến giả) $D_1$ 1: Có tham gia; 0: Không $(+)$ Mai Văn Xuân & Nguyễn Văn Hóa (2011), Đặng Thị Tuyết Ngân (2015)

Mô hình được tiến hành kiểm định các khuyết tật: hiện tượng đa cộng tuyến (thông qua hồi quy phụ và so sánh $R^2_{aux}$ với $R^2_{gốc}$), hiện tượng phương sai sai số thay đổi và tự tương quan.


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

Chương 1: Mở đầu

Chương 1 thiết lập tính cấp thiết của đề tài xuất phát từ tác động tiêu cực của biến đổi khí hậu và hạn mặn tại ĐBSCL, đặc biệt tại tỉnh Bến Tre – nơi có mạng lưới sông rạch chịu ảnh hưởng mặn sâu vào đất liền. Tác giả trình bày vị trí kinh tế của cây bưởi da xanh, xác định rõ mục tiêu chung, 4 mục tiêu cụ thể, đối tượng khảo sát là 60 hộ trồng bưởi tại xã Bình Phú, phạm vi không gian - thời gian nghiên cứu và kết cấu 5 chương của khóa luận.

Chương 2: Tổng quan

Chương 2 tổng hợp các công trình nghiên cứu trong và ngoài nước liên quan đến xâm nhập mặn và năng suất cây trồng:

  • Nghiên cứu quốc tế: Paul M. Barlow và Eric G. Reichard (2010) về xâm nhập mặn tầng nước ngầm tại Bắc Mỹ; Anna Bergqvist và cộng sự (2012) về tác động của hạn mặn tại Bến Tre; Susmita Dasgupta (2018) về độ mặn của đất làm suy giảm 15,6% sản lượng lúa tại Bangladesh.
  • Nghiên cứu trong nước: Các nghiên cứu hàm sản xuất Cobb-Douglas và hồi quy OLS của Lầu Tắc Phu (2012), Quãng Trọng Chinh (2013), Lê Anh Khoa (2018); phân tích chuỗi giá trị và thị trường của Đoàn Thị Hương (2013); các khảo sát thích ứng hạn mặn của Lê Thị Phương Mai (2017), Ngô Quốc Hậu và cộng sự (2017), Dương Hoàng Long (2018) và Võ Văn Bình và cộng sự (2021).
  • Tổng quan địa bàn: Xã Bình Phú có tổng diện tích tự nhiên 965,55 ha (đất nông nghiệp chiếm 621,8 ha), gồm 6 ấp. Tỷ lệ lao động có việc làm đạt 90,79%, thu nhập bình quân 54 triệu đồng/người/năm. Khí hậu chia thành 2 mùa, mùa mưa chiếm 94 - 98% tổng lượng mưa cả năm.
  • Tổng quan về xâm nhập mặn và cây bưởi da xanh: Cung cấp cơ sở sinh học của cây bưởi da xanh (tuổi thọ kinh tế, giai đoạn kiến thiết cơ bản 4 năm đầu, chu kỳ cho trái 8 - 12 năm) và diễn biến mặn tại Bến Tre qua các đợt hạn mặn lịch sử 2015-2016, 2019-2020.

Chương 3: Nội dung và phương pháp nghiên cứu

Chương này hệ thống hóa cơ sở lý luận về kinh tế nông hộ, cơ chế dịch chuyển của khối nước mặn vào hệ thống tầng trữ nước ngọt và sông rạch vùng cửa sông. Tác giả chi tiết hóa các công thức tính toán chỉ tiêu tài chính ($TC, TR, \pi, I$), các hệ số tỷ suất hiệu quả tài chính, quy trình chọn mẫu điều tra 60 hộ theo công thức Harri (1985), thang đo Likert 5 điểm, thiết kế hàm hồi quy OLS Cobb-Douglas dạng log-log với 9 biến độc lập và các tiêu chuẩn kiểm định thống kê liên quan.

Chương 4: Kết quả nghiên cứu và thảo luận

Chương 4 phản ánh các phát hiện từ khảo sát thực địa:

  • Thực trạng mặn: Đợt hạn mặn 2019-2020 ghi nhận độ mặn 4‰ vào sâu tới 86 km trên các sông chính, xâm nhập sâu vào xã Bình Phú từ tháng 12/2019 đến tháng 6/2020; đợt hạn mặn mùa khô 2021-2022 xuất hiện ranh mặn 4‰ cách cửa sông 48 - 58 km.
  • Đặc điểm nông hộ: Người ra quyết định sản xuất chủ yếu là nam giới (88,33%), nữ giới chiếm 11,67% (53 nam và 7 nữ trên tổng 60 quan sát).
  • Mức độ ảnh hưởng: Đánh giá qua thang đo Likert cho thấy mức độ nhiễm mặn bình quân của các vườn bưởi khảo sát đạt 2,75 điểm.
  • Hiệu quả kinh tế năm 2021: Canh tác bưởi da xanh mang lại lợi nhuận bình quân 25,3 triệu đồng/1.000 m² (1 công đất).
  • Kết quả hồi quy mô hình: Trong 9 biến độc lập tham gia hồi quy, có 7 biến tác động có ý nghĩa thống kê đến năng suất bưởi da xanh. Trong đó, 6 yếu tố mang dấu dương (tác động thuận chiều) gồm: kinh nghiệm, lượng nước tưới, lượng phân hữu cơ, lượng phân vô cơ, công lao động và tham gia khuyến nông. Yếu tố mức độ nhiễm mặn mang dấu âm (tác động nghịch chiều), phản ánh mức độ nhiễm mặn càng tăng thì năng suất bưởi càng suy giảm rõ rệt.

Chương 5: Kết luận và kiến nghị

Chương 5 tóm lược các kết quả chính về tình hình sản xuất, mức độ tổn thất năng suất do mặn và hiệu quả kinh tế của nông hộ. Từ các phát hiện nghiên cứu, tác giả đề xuất 3 nhóm giải pháp trọng tâm: nâng cao nhận thức của nông hộ về xâm nhập mặn; chủ động nguồn nước ngọt tưới tiêu (nạo vét ao mương, sử dụng túi chứa nước, hồ bồn trữ nước); địa phương tăng cường hệ thống cảnh báo, quan trắc độ mặn định kỳ và điều tiết công trình thủy lợi ngăn mặn.


Kết quả và đóng góp

Kết quả nghiên cứu chính

  • Hiệu quả kinh tế: Năm 2021, mô hình trồng bưởi da xanh tại xã Bình Phú đạt mức lợi nhuận trung bình là 25,3 triệu đồng/1.000 m² (tương đương 25,3 triệu đồng/công đất).
  • Mức độ nhiễm mặn: Điểm đánh giá mức độ nhiễm mặn trung bình của 60 hộ khảo sát theo thang đo Likert 5 mức độ là 2,75, chứng minh xâm nhập mặn diễn biến ở mức độ đáng kể và gây tác động trực tiếp đến vườn cây.
  • Yếu tố tác động đến năng suất: Ước lượng mô hình OLS xác nhận 7/9 biến có tác động có ý nghĩa thống kê:
    • 6 biến tác động tích cực: Kinh nghiệm chủ hộ ($X_2$), Lượng nước tưới hợp lý ($X_3$), Lượng phân bón hữu cơ ($X_4$), Lượng phân bón vô cơ ($X_5$), Thời gian công lao động ($X_7$), và Việc tham gia các lớp tập huấn khuyến nông ($D_1$).
    • 1 biến tác động tiêu cực: Mức độ nhiễm mặn ($X_8$). Vườn bưởi chịu mức độ nhiễm mặn càng cao thì năng suất thu hoạch càng giảm sút nghiêm trọng.
Nhóm kết quả Chỉ tiêu / Biến số cụ thể Dữ liệu / Tác động ghi nhận
Kinh tế nông hộ (2021) Lợi nhuận trung bình 25,3 triệu đồng/1.000 m²
Khảo sát mặn Mức độ nhiễm mặn vườn bưởi 2,75 / 5 (Thang đo Likert)
Nhân khẩu học mẫu Giới tính người ra quyết định Nam: 88,33% (53 hộ); Nữ: 11,67% (7 hộ)
Hồi quy OLS (Năng suất) Yếu tố tác động thuận chiều $(+)$ Kinh nghiệm, Lượng nước tưới, Phân hữu cơ, Phân vô cơ, Công lao động, Khuyến nông
Yếu tố tác động ngược chiều $(-)$ Mức độ nhiễm mặn ($X_8$)
Yếu tố không có ý nghĩa thống kê Trình độ học vấn ($X_1$), Chi phí thuốc BVTV ($X_6$)

Giải pháp và khuyến nghị của tác giả

  • Về phía nông hộ: Cần chủ động trang bị thiết bị trữ nước ngọt như bồn chứa, túi nhựa, nạo vét mương vườn; theo dõi chặt chẽ nồng độ mặn trước khi bơm tưới; gia tăng tỷ lệ sử dụng phân bón hữu cơ và vôi để cải tạo, rửa mặn cho đất; tuân thủ kỹ thuật canh tác do cán bộ nông nghiệp hướng dẫn.
  • Về phía chính quyền và cơ quan quản lý: Đài Khí tượng Thủy văn và Phòng Nông nghiệp cần tăng cường tần suất dự báo độ mặn trên các sông chính (nhất là sông Hàm Luông); vận hành kịp thời hệ thống cống đập ngăn mặn; mở rộng các lớp tập huấn kỹ thuật chăm sóc cây ăn trái trong điều kiện hạn mặn.

Đóng góp của đề tài

Công trình cung cấp bằng chứng định lượng cụ thể tại địa bàn cấp xã (Bình Phú, thành phố Bến Tre) về mức độ suy giảm năng suất bưởi da xanh dưới tác động của mặn. Nghiên cứu kết hợp giữa phân tích hiệu quả tài chính nông hộ và lượng hóa kinh tế lượng, chỉ rõ vai trò của hoạt động khuyến nông và bón phân hữu cơ trong việc giảm thiểu tổn thất năng suất do thiên tai mặn xâm nhập.


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

  • Hạn chế của nghiên cứu: Dữ liệu sơ cấp giới hạn ở quy mô 60 nông hộ trên địa bàn một xã duy nhất (xã Bình Phú), chưa bao quát toàn bộ các vùng chuyên canh bưởi da xanh khác của tỉnh Bến Tre (như huyện Mỏ Cày Bắc, Châu Thành hay Chợ Lách). Số liệu tài chính thu thập tập trung vào chu kỳ sản xuất năm 2021, chưa phản ánh được biến động đa niên của các chu kỳ xâm nhập mặn có cường độ khác nhau.
  • Hướng nghiên cứu tiếp: Mở rộng phạm vi khảo sát liên huyện nhằm so sánh năng lực thích ứng giữa các vùng sinh thái mặn khác nhau; áp dụng dữ liệu chuỗi thời gian (time-series) hoặc dữ liệu bảng (panel data) để theo dõi sự phục hồi năng suất của vườn bưởi qua nhiều năm; tích hợp các phân tích chi phí - lợi ích (CBA) đối với từng mô hình trữ nước và công nghệ tưới tiết kiệm.

Giá trị tham khảo

  • Đối tượng tham khảo:
    • Sinh viên và học viên cao học: Thuộc các chuyên ngành Kinh tế nông nghiệp, Kinh tế phát triển, Quản lý tài nguyên và Môi trường; người thực hiện khóa luận tốt nghiệp về phân tích sản xuất nông hộ và tác động của biến đổi khí hậu.
    • Cán bộ khuyến nông và nhà quản lý: Cán bộ quản lý nông nghiệp tại các huyện, thị vùng ĐBSCL trong việc xây dựng kế hoạch ứng phó hạn mặn cho vùng cây ăn trái.
  • Nội dung giá trị nhất:
    • Quy trình xây dựng bảng câu hỏi và thiết kế hàm hồi quy OLS Cobb-Douglas log-log ứng dụng cho cây ăn trái lâu năm chịu tác động ngoại cảnh.
    • Khung hệ thống chỉ tiêu tính toán tài chính nông hộ phân tách giữa chi phí vật chất và chi phí lao động (lao động gia đình và lao động thuê ngoài).
    • Các dẫn chứng tổng quan tài liệu về hạn mặn tại ĐBSCL được trích dẫn chi tiết từ các tác giả đi trước.

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

1. Cỡ mẫu 60 nông hộ trong khóa luận được xác định dựa trên công thức khoa học nào?

Cỡ mẫu được xác định dựa trên nguyên tắc chọn mẫu hồi quy của Harri (1985): $N = 50 + m$. Trong đó $N$ là số mẫu tối thiểu cần tìm và $m$ là số biến độc lập trong mô hình. Vì mô hình hồi quy của đề tài sử dụng 9 biến độc lập ($m = 9$), cỡ mẫu lý thuyết tối thiểu là 59 quan sát. Do đó, việc tác giả khảo sát 60 nông hộ trồng bưởi tại xã Bình Phú đảm bảo độ tin cậy để ước lượng mô hình hồi quy tuyến tính.

2. Mô hình kinh tế lượng nào được tác giả sử dụng để đánh giá tác động đến năng suất bưởi?

Tác giả sử dụng hàm sản xuất Cobb-Douglas được chuyển đổi sang dạng tuyến tính log-log và ước lượng bằng phương pháp bình phương nhỏ nhất (OLS) thông qua phần mềm EViews 8. Biến phụ thuộc là logarit tự nhiên của năng suất bưởi ($\ln Y$), chịu sự tác động của 8 biến độc lập liên tục (được logarit hóa hoặc dùng thang đo) và 1 biến giả thể hiện việc tham gia khuyến nông.

3. Kết quả mô hình OLS cho thấy xâm nhập mặn tác động đến năng suất bưởi như thế nào?

Biến mức độ nhiễm mặn của vườn bưởi ($X_8$, đo bằng thang đo Likert 1 - 5) có hệ số hồi quy mang dấu âm và đạt mức ý nghĩa thống kê trong mô hình OLS. Kết quả này khẳng định mức độ nhiễm mặn tác động tiêu cực (nghịch chiều) đến năng suất bưởi da xanh: các vườn bưởi ghi nhận mức độ nhiễm mặn càng nghiêm trọng thì năng suất thu hoạch càng bị sụt giảm.

4. Hiệu quả kinh tế của các nông hộ trồng bưởi tại xã Bình Phú trong năm 2021 đạt mức bao nhiêu?

Theo kết quả điều tra thực tế được tổng kết trong khóa luận, trong năm sản xuất 2021, trung bình mỗi hộ trồng bưởi da xanh tại địa bàn xã Bình Phú thu được mức lợi nhuận thuần là 25,3 triệu đồng trên diện tích 1.000 m² (tương đương 1 công đất canh tác).

5. Những yếu tố nào tác động cùng chiều làm tăng năng suất bưởi da xanh theo kết quả nghiên cứu?

Mô hình chỉ ra 6 yếu tố tác động tích cực có ý nghĩa thống kê đến năng suất bưởi gồm: kinh nghiệm canh tác của nông hộ ($X_2$), lượng nước tưới hợp lý ($X_3$), khối lượng phân bón hữu cơ ($X_4$), khối lượng phân bón vô cơ ($X_5$), thời gian công lao động chăm sóc ($X_7$), và việc tham gia các hoạt động tập huấn khuyến nông của địa phương ($D_1$).


Kết luận

Khóa luận tốt nghiệp của tác giả Trần Thị Yến Nhi đã hoàn thành việc mô tả thực trạng xâm nhập mặn và đánh giá định lượng tác động của hiện tượng này đến các nông hộ trồng bưởi da xanh tại xã Bình Phú, thành phố Bến Tre. Thông qua phương pháp ước lượng OLS mô hình Cobb-Douglas từ dữ liệu khảo sát 60 hộ, công trình chứng minh mức độ nhiễm mặn làm sụt giảm rõ rệt năng suất bưởi, đồng thời làm nổi bật vai trò giảm thiểu rủi ro của phân bón hữu cơ, công chăm sóc và công tác khuyến nông. Các kiến nghị đề xuất tập trung vào việc nâng cao nhận thức, chủ động dự trữ nước ngọt và cải thiện hạ tầng thủy lợi, mang lại giá trị tham khảo thiết thực cho công tác quản lý và thích ứng biến đổi khí hậu tại địa phương.