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

Xâm lấn sinh học và biến đổi khí hậu là hai nhân tố chính thúc đẩy sự thay đổi môi trường toàn cầu, tác động trực tiếp và đa chiều đến tính toàn vẹn của các hệ sinh thái. Các loài thực vật ngoại lai xâm lấn cạnh tranh tài nguyên, làm suy giảm tính đa dạng của thảm thực vật bản địa, biến đổi chu trình dinh dưỡng và phá vỡ cấu trúc cảnh quan tự nhiên. Sự tương tác giữa biến đổi khí hậu và sinh vật xâm lấn tạo nên vòng phản hồi tích cực (positive feedback loop): biến đổi khí hậu mở ra các sinh cảnh mới thích hợp cho loài xâm lấn, trong khi sự phát triển của các loài này làm tăng tính tổn thương của hệ sinh thái trước biến động thời tiết cực đoan.

Chi Keo (Acacia sensu lato) thuộc họ Đậu (Fabaceae) là một trong những chi thực vật lớn nhất ở Úc với gần 1.000 loài bản địa. Nhiều loài keo được đưa vào trồng trọt rộng rãi trên thế giới và tại các bang của Úc nhằm mục đích lâm nghiệp, cải tạo đất, chắn gió và làm cây cảnh. Tuy nhiên, một số loài đã thoát khỏi vùng canh tác, tự nhiên hóa và trở thành loài xâm lấn nguy hiểm. Tại bang Tây Úc (Western Australia - WA), đặc biệt là vùng đô thị mở rộng Perth và đồng bằng ven biển Swan (Swan Coastal Plain), sự xâm lấn của các loài keo có nguồn gốc từ miền Đông và miền Nam nước Úc đang đe dọa nghiêm trọng thảm thực vật bản địa còn sót lại trong các khu bảo tồn.

Đề tài khóa luận tập trung nghiên cứu hai loài keo xâm lấn cụ thể: Acacia longifolia subsp. longifolia (Sydney golden wattle) và Acacia iteaphylla (Flinders Range wattle). Nghiên cứu được thực hiện nhằm giải quyết 3 mục tiêu và các nhiệm vụ cụ thể sau:

  1. Mục tiêu 1: Xác định quy luật và tiềm năng xâm lấn hiện tại của A. longifoliaA. iteaphylla dựa trên các yếu tố khí hậu và môi trường ở quy mô quốc gia và khu vực tại Úc bằng mô hình phân bố loài (Species Distribution Modelling - SDM).
  2. Mục tiêu 2: Dự báo vùng phân bố tiềm năng trong tương lai của hai loài keo này dưới tác động của các kịch bản biến đổi khí hậu toàn cầu đến năm 2045.
  3. Mục tiêu 3: Đánh giá ảnh hưởng của các yếu tố môi trường cục bộ (khoảng cách từ rìa thảm thực vật, độ sâu mực nước ngầm, độ che phủ tầng dưới/mức độ xáo trộn, cháy rừng và sự hiện diện của cây gỗ lớn) đến mật độ và cấu trúc xâm lấn của hai loài keo tại vùng Perth, Tây Úc.

Đối tượng nghiên cứu là hai phân loài/loài keo Acacia longifolia subsp. longifoliaAcacia iteaphylla. Phạm vi không gian bao gồm toàn bộ lãnh thổ nước Úc đối với mô hình hóa phân bố vĩ mô (Tây Úc - WA, Nam Úc - SA, Đông Úc - EA và Tasmania), kết hợp khảo sát thực địa chi tiết tại khu vực Swan Coastal Plain (Perth, Tây Úc). Phạm vi thời gian nghiên cứu thực địa diễn ra trong 3 tháng vào năm 2018, kết hợp chuỗi dữ liệu lịch sử phân bố và dự báo khí hậu đến năm 2045.

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

Nghiên cứu áp dụng khung lý thuyết sinh thái học xâm lấn và lý thuyết ổ sinh thái (ecological niche theory). Các khái niệm cốt lõi bao gồm: đặc điểm thích nghi của loài cỏ dại (Baker, 1965), giả thuyết giải phóng kẻ thù tự nhiên (enemy release hypothesis), hội chứng phát tán hạt nhờ kiến (myrmecochory) và chim, cũng như cơ chế phản ứng của thực vật với chế độ cháy rừng (fire regimes).

Để thực hiện các mục tiêu đề ra, tác giả kết hợp hai phương pháp nghiên cứu chính:

1. Phương pháp mô hình hóa phân bố loài (SDM)

  • Thuật toán và nền tảng: Sử dụng thuật toán Entropy Cực đại (MaxEnt) thông qua hai hệ thống trực tuyến là Phòng thí nghiệm ảo về Đa dạng sinh học và Biến đổi khí hậu (Biodiversity and Climate Change Virtual Laboratory - BCCVL) và Atlas of Living Australia (ALA).
  • Dữ liệu phân bố loài: Trích xuất tọa độ điểm xuất hiện tự nhiên (native range) và điểm đã tự nhiên hóa (naturalised range) của A. longifoliaA. iteaphylla từ cơ sở dữ liệu quốc gia ALA.
  • Biến số môi trường và khí hậu: Bao gồm các biến sinh khí hậu (bioclimatic variables như nhiệt độ tuần ấm nhất B05, lượng mưa quý khô nhất B17, lượng mưa năm B12, tính đẳng nhiệt B03, tính mùa của nhiệt độ B04/B15) và các biến môi trường vật lý (chỉ số ẩm el864/el873/el891, khoảng cách đến bờ biển el819, hàm lượng sét el814, độ sâu tầng đất el816, bề mặt xói mòn el823, chỉ số cách biệt el647).
  • Kịch bản biến đổi khí hậu tương lai: Mô phỏng đến năm 2045 dựa trên mô hình hoàn lưu chung CSIRO Mark 3 kết hợp hai kịch bản nồng độ khí nhà kính đại diện: Kịch bản trung bình (RCP 4.5) và Kịch bản phát thải cao (RCP 8.5).
  • Đánh giá mô hình: Sử dụng diện tích dưới đường cong ROC (Area Under Curve - AUC) để kiểm định độ chính xác dự báo của mô hình.

2. Phương pháp khảo sát thực địa và phân tích thống kê cục bộ

  • Địa điểm khảo sát: Khu vực đồng bằng ven biển Swan, thành phố Perth, bang Tây Úc (bao gồm các khu bảo tồn và tuyến đường như Thomsons Lake Nature Reserve, Anketell, De Haer Road, Piara, Shirley Balla Swamp Reserve).
  • Chỉ tiêu đo đếm thực địa: Thiết lập các ô tiêu chuẩn theo tuyến (transects) từ rìa thảm thực vật vào sâu bên trong để thu thập dữ liệu về:
    • Mật độ cây (Plant Density - số cây/ha).
    • Diện tích gốc (Basal Area - $\text{m}^2/\text{ha}$).
    • Thể tích sinh khối (Biovolume - $\text{m}^3/\text{ha}$).
    • Khoảng cách từ rìa thảm thực vật tự nhiên (Edge effect).
    • Độ sâu mực nước ngầm (Depth to groundwater).
    • Phân loại độ che phủ tầng dưới/mức độ xáo trộn (Understorey cover).
    • Lịch sử cháy rừng (vùng bị cháy - burnt và vùng không bị cháy - unburnt).
    • Sự hiện diện của cây gỗ lớn tầng cao (Occurrence of larger trees).
  • Phương pháp xử lý số liệu: Sử dụng kiểm định t hai mẫu giả định phương sai khác nhau (Two-sample t-Test assuming unequal variances, mức ý nghĩa $\alpha = 0{,}05$) để so sánh sinh trưởng giữa vùng cháy và không cháy; sử dụng kiểm định Chi bình phương ($\chi^2$) để xác định mức độ tương quan giữa vị trí cây keo và bóng tán cây gỗ lớn.

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

Phần I: Giới thiệu (PART I: INTRODUCTION)

Chương mở đầu phân tích tính cấp thiết của đề tài trong bối cảnh toàn cầu hóa và biến đổi khí hậu. Tác giả chỉ ra lỗ hổng nghiên cứu: các công trình trước đây về Acacia xâm lấn tập trung chủ yếu ở Nam Phi hoặc không phân biệt rõ ràng hai phân loài của A. longifolia (subsp. longifoliasubsp. sophorae), đồng thời thiếu các khảo sát định lượng chi tiết ở cấp độ vi mô tại khu vực đô thị Perth. Chương này cũng công bố hệ thống mục tiêu, 3 nhóm câu hỏi nghiên cứu kèm các giả thuyết $H_0/H_1$ tương ứng, xác định giới hạn nghiên cứu (thời gian thực địa 3 tháng, tính chính xác của bản đồ cháy lịch sử và độ che phủ cỏ dại) và định nghĩa các thuật ngữ khoa học then chốt.

Phần II: Tổng quan tài liệu (PART II: LITERATURE REVIEW)

Tổng quan hệ thống hóa các cơ sở lý luận sinh thái:

  • Cơ chế xâm lấn sinh học và tác động: Đánh giá thiệt hại kinh tế hàng tỷ USD và sự suy thoái đa dạng sinh học toàn cầu.
  • Các nhân tố điều khiển xâm lấn: Tác động của biến đổi khí hậu; các đặc tính cỏ dại theo tiêu chí Baker (1965) như nảy mầm không liên tục, hạt tồn tại lâu trong đất, khả năng tự thụ hoặc thụ phấn chéo; các yếu tố môi trường (khoảng cách phát tán qua trung gian kiến, chim, thú có túi; độ sâu nước ngầm; chế độ cháy).
  • Phân loại học và lịch sử chi Acacia: Phân tích lịch sử danh pháp phức tạp của chi Acacia, quyết định tái định danh của Đại hội Thực vật Quốc tế (IBC) và đặc điểm của phân chi đơn ngành Phyllodineae với hơn 950 loài bản địa Úc.
  • Đặc điểm hai loài nghiên cứu:
    • Acacia longifolia subsp. longifolia: Cây bụi đứng cao đến 10m, phiến lá mỏng dài (5–20 cm), chịu sương giá đến $-6^\circ\text{C}$, nhu cầu lượng mưa tối thiểu 550 mm/năm, cố định đạm tốt trên đất nghèo dinh dưỡng.
    • Acacia iteaphylla: Cây bụi rậm nguồn gốc dãy Flinders (Nam Úc), cao 2–5m, chịu hạn rất tốt, chịu lạnh đến $-7^\circ\text{C}$, được phân loại là loài xâm lấn xuyên lục địa (trans-continental invader).
  • Phát tán hạt và phản ứng với lửa: Khái quát vai trò của thể dầu (elaiosome) trong việc thu hút kiến phát tán cự ly ngắn (trung bình 5,4–12,2 m; tối đa 180 m) và chim phát tán cự ly xa (trên 200 m); cơ chế phá vỡ trạng thái ngủ nghỉ của hạt vỏ cứng nhờ nhiệt độ đám cháy.

Phần III: Phương pháp nghiên cứu (PART III: METHODS)

Mô tả quy trình kỹ thuật trích xuất dữ liệu tọa độ từ ALA, xây dựng tập biến sinh khí hậu và môi trường số trên hệ thống BCCVL và ALA. Quy định phương pháp khảo sát thực địa tại các điểm đại diện trên đồng bằng Swan, phương thức đo đường kính gốc, chiều cao tán để quy đổi thành diện tích gốc và thể tích sinh khối cây. Trình bày chi tiết công thức và các kiểm định thống kê sinh học áp dụng trong mô hình.

Phần IV: Kết quả và thảo luận (PART IV: RESULTS AND DISCUSSIONS)

Chương này tích hợp toàn bộ kết quả phân tích mô hình hóa và thực địa:

  • Độ chính xác của mô hình SDM: Giá trị AUC của các mô hình dự báo đều đạt mức tin cậy cao theo Bảng 1.
  • Phân tích phản ứng sinh thái: Các đường cong phản ứng (response curves) chỉ ra ngưỡng chịu đựng nhiệt độ, lượng mưa và các nhân tố thổ nhưỡng của từng loài.
  • Biến động phân bố vĩ mô: So sánh bản đồ phân bố hiện tại và tương lai (năm 2045) theo hai kịch bản RCP 4.5 và RCP 8.5 (Hình 1 đến Hình 16).
  • Dữ liệu thực địa tại Perth: Trình bày đường cong tái sinh (recruitment curves), biểu đồ mật độ theo độ sâu mực nước ngầm, biểu đồ tròn thể hiện tỷ lệ độ che phủ tầng dưới và kết quả kiểm định thống kê so sánh giữa vùng cháy/không cháy và ảnh hưởng của cây gỗ lớn.

Phần V: Kết luận và khuyến nghị (PART V: CONCLUSIONS AND RECOMMENDATIONS)

Tổng hợp các phát hiện cốt lõi từ SDM và thực địa; đưa ra các kiến nghị quản lý thảm thực vật tự nhiên tại Perth, kiểm soát việc buôn bán cây cảnh sân vườn và định hướng các nghiên cứu tiếp theo.

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

Khóa luận đã định lượng hóa sự khác biệt rõ rệt về mặt sinh thái học và phản ứng môi trường giữa hai loài keo xâm lấn:

Tiêu chí so sánh Acacia longifolia subsp. longifolia Acacia iteaphylla
Phân bố tự nhiên (Native range) Đông Nam nước Úc (NSW, Victoria, Queensland) Nam Úc (Flinders Ranges, Gawler Ranges, Eyre Peninsula)
Ảnh hưởng của biến môi trường vào SDM hiện tại Vùng phân bố tiềm năng bị thu hẹp so với chỉ dùng biến khí hậu Vùng phân bố tiềm năng được mở rộng khi tích hợp thêm biến môi trường
Dự báo phân bố tương lai (năm 2045 - RCP 4.5) Tiếp tục mở rộng về phía Nam Đông Úc và đảo Tasmania Ngừng mở rộng và bắt đầu suy giảm mật độ tại WA, SA, EA
Dự báo phân bố tương lai (năm 2045 - RCP 8.5) Mở rộng liên tục và mạnh mẽ Thu hẹp liên tục trên toàn bộ phạm vi phân bố
Môi trường sống ưa thích tại thực địa (Perth) Khu vực đất ẩm ướt, ít bị xáo trộn (undisturbed wet areas) Khu vực đất cát khô ráo, xáo trộn mạnh (disturbed dry areas)
Phản ứng với cháy rừng Mật độ, diện tích gốc và thể tích cây ở vùng bị cháy cao hơn đáng kể ($p < 0{,}05$) Thích nghi trên đất xáo trộn, tái sinh mạnh sau tác động vật lý
Mối quan hệ với cây gỗ lớn tầng cao Ưu tiên xuất hiện dưới bóng tán cây lớn ($\chi^2$ có ý nghĩa thống kê) Ưu tiên xuất hiện dưới bóng tán cây lớn ($\chi^2$ có ý nghĩa thống kê)
     [ Acacia longifolia subsp. longifolia ]              [ Acacia iteaphylla ]
     • Khí hậu tương lai: MỞ RỘNG (RCP 4.5 & 8.5)        • Khí hậu tương lai: THU HẸP (RCP 4.5 & 8.5)
     • Điều kiện lập địa: Khu vực ẩm, ít xáo trộn         • Điều kiện lập địa: Khu vực khô, xáo trộn mạnh
     • Tác động của lửa: Tăng mật độ sau cháy            • Tác động cạnh tranh: Chịu hạn tốt trên đất cát

Đóng góp khoa học của nghiên cứu:

  • Tác giả đã phân lập độc lập phân loài Acacia longifolia subsp. longifolia, khắc phục sự nhầm lẫn học thuật kéo dài khi gộp chung với subsp. sophorae.
  • Chứng minh rằng việc bổ sung các biến môi trường vật lý (độ sâu đất, khoảng cách bờ biển, thành phần sét) làm thay đổi đáng kể độ chính xác của mô hình phân bố loài so với việc chỉ sử dụng dữ liệu sinh khí hậu đơn thuần.
  • Cung cấp bằng chứng thực địa định lượng đầu tiên về mối liên hệ giữa chiều sâu mực nước ngầm, chế độ đốt thực bì kiểm soát và sự phân bố của hai loài keo xâm lấn tại các dải thực vật đô thị ven biển Perth.

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

Khóa luận được tác giả ghi nhận một số hạn chế khách quan trong quá trình thực hiện:

  • Thời gian và mẫu điều tra: Thời gian khảo sát thực địa tại bang Tây Úc bị giới hạn trong 3 tháng, chịu ảnh hưởng bởi thời tiết mùa mưa khiến số lượng địa điểm lấy mẫu chưa thể bao phủ toàn diện.
  • Độ chính xác của dữ liệu nền: Bản đồ dữ liệu lịch sử cháy rừng hàng năm và bản đồ phân bố độ che phủ cỏ dại có sai số nhất định so với hiện trạng quan sát thực tế ngoài thực địa tại một số điểm khảo sát.
  • Biến số trong mô hình SDM: Hệ thống BCCVL và ALA bị hạn chế về số lượng lớp dữ liệu môi trường có sẵn; các mô hình SDM chưa tích hợp được các biến tương tác sinh học (biotic interactions) do thiếu nguồn dữ liệu không gian tương ứng.

Hướng nghiên cứu tiếp theo được gợi mở bao gồm: mở rộng mạng lưới ô tiêu chuẩn thực địa sang các mùa sinh trưởng khác nhau; tích hợp các biến số tương tác sinh học và vi khí hậu vào mô hình SDM; đánh giá chi tiết ngân hàng hạt trong đất (soil seed bank) sau các đợt đốt thực bì có kiểm soát.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo chuyên khảo hữu ích cho:

  • Giảng viên, sinh viên và học viên cao học chuyên ngành Quản lý Môi trường, Khoa học Môi trường, Lâm nghiệp, Sinh thái học và Đa dạng sinh học.
  • Các nhà nghiên cứu quan tâm đến quy trình ứng dụng thuật toán MaxEnt trên nền tảng BCCVL và ALA để dự báo phân bố loài sinh vật dưới các kịch bản biến đổi khí hậu (RCP 4.5, RCP 8.5).
  • Các cán bộ quy hoạch bảo tồn thiên nhiên, kiểm lâm và nhà quản lý tài nguyên tại các khu bảo tồn thực vật đô thị trong việc xây dựng kế hoạch kiểm soát loài ngoại lai xâm lấn và quản lý chế độ cháy rừng phòng hộ.

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

1. Hai loài keo trong nghiên cứu có phản ứng khác nhau như thế nào dưới kịch bản biến đổi khí hậu tương lai (năm 2045)?
Theo kết quả mô hình hóa CSIRO Mark 3, dưới cả hai kịch bản RCP 4.5 và RCP 8.5, Acacia longifolia có xu hướng tiếp tục mở rộng phạm vi phân bố về phía Nam vùng Đông Úc và đảo Tasmania. Ngược lại, Acacia iteaphylla ngừng mở rộng và suy giảm mật độ ở tất cả các vùng phân bố hiện tại dưới kịch bản RCP 4.5, đồng thời liên tục bị thu hẹp diện tích dưới kịch bản phát thải cao RCP 8.5.

2. Khảo sát thực địa tại Perth cho thấy sự khác biệt về sinh cảnh ưa thích của hai loài keo này là gì?
Acacia longifolia có xu hướng phân bố và sinh trưởng nhiều hơn tại các khu vực ẩm ướt, ít bị xáo trộn (undisturbed wet areas). Trong khi đó, Acacia iteaphylla xuất hiện và phát triển ưu thế tại các khu vực khô ráo, có mức độ xáo trộn cao (well disturbed dry areas).

3. Cháy rừng có tác động như thế nào đến sự phát triển của Acacia longifolia?
Kết quả kiểm định t hai mẫu cho thấy các điểm khảo sát bị cháy (burnt sites) có mật độ cây (plant density), diện tích gốc (basal area) và thể tích sinh khối (biovolume) của Acacia longifolia cao hơn đáng kể so với các điểm không bị cháy (unburnt sites), do nhiệt lượng đám cháy kích thích sự nảy mầm hàng loạt từ ngân hàng hạt trong đất.

4. Hai loài keo này có đặc điểm chung nào khi phân bố ở quy mô cục bộ?
Kết quả phân tích kiểm định Chi bình phương ($\chi^2$) cho thấy cả hai loài keo A. longifoliaA. iteaphylla đều thể hiện sự ưu tiên phân bố đáng kể dưới tán của các cây gỗ lớn hơn trong thảm thực vật tự nhiên.

Kết luận

Khóa luận tốt nghiệp của tác giả Aleisa Hannah Tadios Ayson đã giải quyết trọn vẹn mục tiêu đánh giá tác động của biến đổi khí hậu và các nhân tố môi trường đến sự phân bố của hai loài keo xâm lấn Acacia longifolia subsp. longifoliaAcacia iteaphylla tại Úc. Nghiên cứu đã kết hợp thành công mô hình hóa sinh thái quy mô vĩ mô (MaxEnt, BCCVL, ALA) với phương pháp điều tra thực địa định lượng tại vùng Perth. Các kết quả đạt được cung cấp cơ sở khoa học xác thực cho công tác dự báo lan truyền sinh vật ngoại lai và hoạch định chiến lược quản lý thảm thực vật bản địa trong bối cảnh biến đổi khí hậu toàn cầu.