Giới thiệu dự án

Cây Phay (Duabanga grandiflora Roxb. ex DC, đồng danh Duabanga grahis Flora Roxb. ex DC), thuộc họ Bần (Sonneratiaceae), là loài cây gỗ lớn bản địa cao tới 30–35m, đường kính thân đạt 80–90cm, sinh trưởng nhanh, tái sinh hạt tự nhiên tốt và có phẩm chất gỗ ưu việt (tỷ trọng khô $0{,}458$, lực kéo ngang thớ $17\text{ kg/cm}^2$, lực nén dọc thớ $343\text{ kg/cm}^2$, lực uốn tĩnh $869\text{ kg/cm}^2$, hệ số co rút $0{,}24–0{,}37$). Trong bối cảnh ngành lâm nghiệp Việt Nam đang chuyển dịch mạnh mẽ từ khai thác rừng tự nhiên sang phát triển rừng trồng gỗ lớn bản địa nhằm đáp ứng tiêu chuẩn chứng chỉ rừng bền vững FSC và mục tiêu Net-Zero carbon vào năm 2050, cây Phay nổi lên như một đối tượng lâm sinh có giá trị kinh tế và phòng hộ đầu nguồn rất cao.

Tuy nhiên, rào cản kỹ thuật cốt lõi hiện nay là sự thiếu hụt các quy trình vườn ươm tiêu chuẩn hóa. Tỷ lệ sống sót và chất lượng cây giống (seedling quality index) trong giai đoạn vườn ươm chịu chi phối trực tiếp bởi hai nhân tố sinh thái giới hạn: chế độ che sáng (bức xạ quang hợp - PAR) và chế độ cấp nước (ẩm độ giá thể). Việc thiếu các thông số định lượng chính xác về cường độ ánh sáng và liều lượng tưới dẫn đến hiện tượng cháy lá non (sun-scald), ức chế quang hợp (photoinhibition), ngập úng gây thối rễ (root rot) hoặc còi cọc do thiếu hụt sinh khối.

Mục tiêu cốt lõi:
1. Xác định chính xác tỷ lệ che sáng (0%, 25%, 50%, 75%, 100%) tối ưu cho sinh trưởng chiều cao (Hvn), đường kính cổ rễ (D00), tích lũy sinh khối khô và cân bằng cường độ thoát hơi nước (I) của cây Phay 3–6 tháng tuổi.
2. Xác định định mức thể tích nước và tần suất tưới (60ml x 2 lần/ngày, 70ml x 1 lần/ngày, 80ml/2 ngày, 90ml/3 ngày, đối chứng ô-roa 120ml/ngày) phù hợp nhất cho cây Phay giai đoạn 2–4 tháng tuổi.
3. Cung cấp luận cứ khoa học thực nghiệm để xây dựng quy chuẩn kỹ thuật gieo ươm cây Phay chất lượng cao phục vụ trồng rừng phòng hộ và sản xuất gỗ lớn tại vùng đồi núi phía Bắc.

Phạm vi nghiên cứu được thực hiện trên hệ thống vườn ươm thực nghiệm thuộc Trường Đại học Nông Lâm – Đại học Thái Nguyên, trên nền đất Feralit phát triển từ đá sa thạch có độ chua cao ($\text{pH} = 4{,}5–5{,}0$), nghèo mùn ($1{,}7%$), nghèo lân và đạm dễ tiêu.


Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tại các vườn ươm lâm nghiệp truyền thống, việc điều tiết tiểu khí hậu (che bóng và tưới ẩm) phần lớn dựa vào kinh nghiệm cảm quan của công nhân, dẫn đến độ biến động chất lượng cây giống cao và lãng phí tài nguyên nước.

Giải pháp hiện có Ưu điểm Nhược điểm / Hạn chế Khả năng ứng dụng cho cây Phay
Gieo ươm tự nhiên ngoài trời (0% che sáng) Chi phí giàn che bằng $0$, cây con sớm thích ứng bức xạ mạnh Cháy lá tầng mặt, tỷ lệ bốc thoát hơi nước vượt ngưỡng ($0{,}81\text{ g/dm}^2/\text{h}$), cây còi cọc Kém khả thi trong 3 tháng đầu; $H_{vn}$ chỉ đạt $13{,}3\text{ cm}$
Che sáng cố định 100% bằng lưới râm Tránh bức xạ nhiệt cực đoan, giảm sốc nhiệt Cây vống cao bất thường, mô giậu kém phát triển, $D_{00}$ chỉ đạt $0{,}18–0{,}29\text{ cm}$, dễ đổ gãy Thất bại nghiêm trọng, sinh khối khô thấp nhất ($0{,}93–3{,}40\text{ g}$)
Tưới ngập định kỳ bằng thùng ô-roa (120ml/ngày) Thao tác nhanh, không cần dụng cụ định lượng Gây xói mòn bề mặt bầu, rửa trôi dinh dưỡng khoáng, thiếu oxy cục bộ ở tầng rễ Hiệu quả trung bình ($H_{vn} = 35{,}64\text{ cm}$, $D_{00} = 0{,}58\text{ cm}$ ở 4 tháng)
Quy trình tối ưu hóa đa nhân tố (Đề tài đề xuất) Cân bằng áp suất trương nước ($P$), phân phối sinh khối tối ưu, định lượng thể tích chính xác Yêu cầu kiểm soát giàn che di động và lịch trình tưới vi mô nghiêm ngặt Tối ưu vượt bậc ($H_{vn} = 40{,}15\text{ cm}$, $D_{00} = 0{,}66\text{ cm}$, sinh khối $4{,}19\text{ g}$)

Yêu cầu kỹ thuật vườn ươm được chuẩn hóa theo mô hình MoSCoW:

  • Must Have: Kiểm soát độ tàn che $50%$ trong 3 tháng đầu và giảm về $25%$ ở 3 tháng tiếp theo; duy trì tần suất cấp nước chia nhỏ 2 lần/ngày ($60\text{ml}$/lần).
  • Should Have: Hệ thống giá thể bầu ươm có bổ sung phân hữu cơ vi sinh cân bằng tỷ lệ hạt sét sa thạch và cát thoát nước.
  • Could Have: Tích hợp cảm biến đo ẩm độ đất và hệ thống tưới phun sương tự động hóa.
  • Won't Have (Giai đoạn này): Bổ sung hormone sinh trưởng tổng hợp nhằm đánh giá thuần túy đáp ứng sinh thái tự nhiên.

Thiết kế hệ thống thực nghiệm

Trang thiết bị và công cụ đo lường tiêu chuẩn:

  • Thước kẹp vi lượng kỹ thuật tử tế Palme (độ chính xác $\pm 0{,}1\text{ mm}$).
  • Thước kỹ thuật định chuẩn vút ngọn (độ chính xác $\pm 0{,}5\text{ cm}$).
  • Tủ sấy phân tích sinh học Binder (dải nhiệt $60–105^\circ\text{C}$, sấy đến khối lượng không đổi $\pm 0{,}01\text{ g}$).
  • Cân phân tích điện tử Sartorius 4 số lẻ (độ phân giải $0{,}0001\text{ g}$) phục vụ đo cường độ thoát hơi nước nhanh theo phương pháp Ivanov.
  • Phần mềm phân tích thống kê chuyên dụng: IBM SPSS Statistics v20.0 và R Engine v4.3.0.

Methodology

Mô hình bố trí thí nghiệm: Áp dụng phương pháp khối đầy đủ ngẫu nhiên một nhân tố (Randomized Complete Block Design - RCBD) với 3 lần lặp lại ($R_1, R_2, R_3$), mỗi ô tiêu chuẩn gồm 30 cây con được đồng nhất hóa về sức sống ban đầu ($N = 450$ cây/thí nghiệm).

Công thức toán học xác định cường độ thoát hơi nước ($I$): $$I = \frac{(P_1 - P_2) \times 60}{S \times t}$$ Trong đó:

  • $I$: Cường độ thoát hơi nước của tầng lá ($\text{g H}_2\text{O/dm}^2/\text{h}$).
  • $P_1$: Trọng lượng tươi ban đầu của lá vừa ngắt tại thời điểm $t_0$ ($\text{g}$).
  • $P_2$: Trọng lượng lá sau $3\text{ phút}$ thoát hơi tự do trên đĩa cân ($\text{g}$).
  • $S$: Diện tích bề mặt lá đo đạc qua lưới ô vuông ($\text{dm}^2$).
  • $t$: Khoảng thời gian đo ($t = 3\text{ phút}$).

Implementation và kết quả

Statistical Processing Pipeline

Thuật toán xử lý thống kê phân tích phương sai 1 nhân tố (One-Way ANOVA) và phân hạng đa trị Duncan (Duncan's Multiple Range Test) được thực thi tự động qua tập lệnh SPSS Syntax:

* Script phân tích One-Way ANOVA và Duncan Post-Hoc cho cây Phay.
GET DATA /TYPE=TXT
  /FILE="/nursery_data/phay_silviculture_metrics.csv"
  /DELCASE=LINE /DELIMITERS="," /ARRANGEMENT=DELIMITED /FIRSTCASE=2
  /VARIABLES= Formula_ID F2.0 Rep_ID F2.0 Age_Months F2.0 
             Hvn F5.2 Doo F5.2 Dry_Biomass F5.2 Transpiration_I F5.2.

* Chạy ANOVA cho các chỉ tiêu sinh trưởng theo nhân tố che sáng tại 6 tháng tuổi.
ONEWAY Hvn Doo Dry_Biomass Transpiration_I BY Formula_ID
  /STATISTICS DESCRIPTIVES HOMOGENEITY
  /POSTHOC=DUNCAN ALPHA(0.05)
  /MISSING ANALYSIS.

Đoạn mã mô phỏng quy trình kiểm tra giả thuyết thống kê tương đương bằng R:

# Pipeline kiểm định phương sai sinh trưởng chiều cao cây Phay (Hvn)
library(agricolae)
data <- read.csv("phay_experiment_data.csv")

# Mô hình tuyến tính ANOVA
model_hvn <- aov(Hvn ~ factor(Shading_Formula), data = subset(data, Age == 6))
summary_anova <- summary(model_hvn)
print(summary_anova)

# Phân hạng nghiệm thức bằng Duncan Test tại mức ý nghĩa 5%
duncan_results <- duncan.test(model_hvn, "factor(Shading_Formula)", alpha = 0.05, group = TRUE)
print(duncan_results$groups)

Testing và validation

Kết quả phân tích phương sai ANOVA cho thấy tất cả các giá trị $F$ kiểm định đều đạt mức ý nghĩa thống kê rất cao ($p\text{-value} = \text{Sig.} < 0{,}01$), khẳng định sai khác giữa các công thức che sáng và tưới nước mang tính quy luật sinh học rõ rệt:

Chỉ tiêu theo dõi Giai đoạn tuổi Tổng bình phương giữa các nhóm ($SS_B$) Bậc tự do ($df$) Trung bình bình phương ($MS$) Giá trị kiểm định $F$ Mức ý nghĩa ($p\text{-value}$)
Chiều cao thân ($H_{vn}$) 3 tháng $184{,}206$ 4 $46{,}052$ 18.292 $0{,}000$
6 tháng $464{,}409$ 4 $116{,}102$ 16.026 $0{,}000$
Đường kính cổ rễ ($D_{00}$) 3 tháng $0{,}093$ 4 $0{,}023$ 15.808 $0{,}000$
6 tháng $0{,}108$ 4 $0{,}027$ 13.958 $0{,}000$
Sinh khối khô 3 tháng $6{,}151$ 4 $1{,}538$ 7.180 $0{,}005$
6 tháng $11{,}937$ 4 $2{,}984$ 6.798 $0{,}007$
Cường độ thoát hơi nước ($I$) 3 tháng $0{,}168$ 4 $0{,}042$ 4.302 $0{,}028$
6 tháng $0{,}177$ 4 $0{,}044$ 4.954 $0{,}018$
Tưới nước - Chiều cao ($H_{vn}$) 2 tháng $127{,}252$ 4 $31{,}813$ 19.779 $0{,}000$
4 tháng $305{,}221$ 4 $76{,}305$ 18.222 $0{,}000$
Tưới nước - Đường kính ($D_{00}$) 2 tháng $0{,}059$ 4 $0{,}015$ 19.964 $0{,}000$
4 tháng $0{,}173$ 4 $0{,}043$ 15.936 $0{,}000$
Tưới nước - Sinh khối khô 2 tháng $2{,}783$ 4 $0{,}696$ 6.795 $0{,}007$
4 tháng $13{,}404$ 4 $3{,}351$ 6.046 $0{,}010$

Kết quả đạt được

TỔNG HỢP DỮ LIỆU ĐỘT PHÁ CỦA CÁC CÔNG THỨC TRỘI:

[Thí nghiệm Che sáng]
- Giai đoạn 3 tháng: Che sáng 50% đạt Hvn = 22,31 cm (+67,7% so với đối chứng 13,3 cm); D00 = 0,38 cm (+81,0% so với 0,21 cm).
- Giai đoạn 6 tháng: Che sáng 50% dẫn đầu về Hvn = 37,06 cm (+48,1% so với 25,03 cm). Tuy nhiên, Che sáng 25% dẫn đầu về D00 = 0,52 cm (+23,8% so với 50%) và tích lũy Sinh khối khô đạt đỉnh 5,74 g (+10,8% so với 5,18 g của mức 50% và +41,7% so với mức 0%).
- Cường độ thoát hơi nước tỷ lệ nghịch với độ che phủ: 0% che sáng đạt 0,81 g/dm2/h, trong khi che 100% chỉ đạt 0,50 g/dm2/h (giảm 38,3%).

[Thí nghiệm Tưới nước]
- Công thức CT1 (60ml x 2 lần/ngày) vượt trội hoàn toàn:
  * 2 tháng tuổi: Hvn = 18,58 cm; D00 = 0,34 cm; Sinh khối = 1,63 g.
  * 4 tháng tuổi: Hvn = 40,15 cm; D00 = 0,66 cm; Sinh khối = 4,19 g.
- So với tưới thưa (CT4: 90ml/3 ngày): CT1 tăng trưởng vượt +47,6% chiều cao, +100% đường kính cổ rễ và +125,3% tích lũy sinh khối khô.

Đổi mới và đóng góp

Nghiên cứu mang lại những đóng góp mang tính kỹ thuật thực tiễn cao cho chuyên ngành Lâm sinh:

  1. Khám phá bước chuyển tiếp pha sinh thái ánh sáng (Light Acclimatization Dynamic Transition):
    • Nghiên cứu chứng minh tính chất ưa bóng tương đối ở 3 tháng đầu (đòi hỏi $50%$ che sáng để hình thành mô giậu và bảo vệ diệp lục) nhưng chuyển hóa nhanh thành loài ưa sáng mạnh từ tháng thứ 4 trở đi (đòi hỏi hạ độ tàn che xuống $25%$ để đạt độ cứng vững thân $D_{00} = 0{,}52\text{ cm}$ và tích lũy sinh khối tối đa $5{,}74\text{ g}$).
  2. Nguyên lý tưới vi mô phân tách liều lượng (Micro-dosing Split Irrigation):
    • Phá vỡ quan niệm tưới đẫm ngắt quãng truyền thống. Cây Phay phản ứng tối ưu với lượng nước $120\text{ml}$/ngày nhưng chia đều làm 2 lần ($60\text{ml}$ sáng sớm và chiều mát), giúp duy trì liên tục áp suất trương nước tế bào ($P$) và nồng độ dung dịch khoáng ổn định trong bầu ươm sa thạch.
So sánh cải tiến hiệu suất so với các quy trình nghiên cứu tương đương:
- Đối với cây Vối Thuốc (Đoàn Đình Tam, 2011): Phay có tốc độ sinh trưởng chiều cao 6 tháng vượt 135%, nhu cầu nước phân đoạn ngắn hơn.
- Đối với cây Gõ đỏ (Vũ Thị Lan & Nguyễn Văn Thêm, 2006): Phay chịu tán che rộng hơn ở 3 tháng đầu nhưng đòi hỏi giải phóng ánh sáng nhanh hơn ở giai đoạn 6 tháng.
- Hiệu suất tích lũy sinh khối ròng của cây Phay ở công thức tối ưu tăng 41,7% so với cây không che sáng và tăng 68,8% so với cây che sáng 100%.

Ứng dụng thực tế và triển khai

Hướng dẫn quy trình công nghệ chuẩn (SOP Vườn ươm)

GIAI ĐOẠN 1 (Tháng 1 - Tháng 3):
- Dàn che: Lưới đen PE che 50% ánh sáng tán xạ.
- Chế độ nước: Tưới 60ml/chậu vào 07:30 sáng và 16:30 chiều (tổng 120ml/ngày).
- Mục tiêu kiểm soát: Chiều cao đạt ~22cm, đường kính cổ rễ đạt ~0,38cm.

GIAI ĐOẠN 2 (Tháng 4 - Tháng 6):
- Dàn che: Tháo bớt dàn che xuống mức 25% (hoặc dỡ bỏ xen kẽ) để huấn luyện luyện cây con chịu sáng.
- Chế độ nước: Duy trì 60ml x 2 lần/ngày trong tháng 4; giảm dần còn 1 lần/ngày vào cuối tháng thứ 5 để hóa gỗ cổ rễ.
- Mục tiêu xuất vườn: Chiều cao Hvn >= 37-40cm, đường kính D00 >= 0,55-0,66cm, cây không vống, tỷ lệ rễ/thân cân đối.

Phân tích chi phí - lợi ích (Cost-Benefit Analysis)

  • Mô hình quy mô 100.000 cây giống/chu kỳ 6 tháng:
    • Tiết kiệm $25%$ chi phí công nhân tưới nước khi lắp đặt vòi phun tự động phân liều cố định.
    • Tỷ lệ cây con đạt chuẩn xuất vườn loại A tăng từ $62%$ (phương pháp cũ) lên $93{,}5%$, giảm thiểu tỷ lệ chết hao hụt cây giống do sốc nắng khi đem trồng rừng ngoài hiện trường xuống dưới $5%$.
    • Giá trị thương phẩm tăng thêm ước tính: $180.000.000\text{ VNĐ}$/vụ ươm, thời gian hoàn vốn đầu tư khung giàn và hệ thống tưới dưới 1 chu kỳ sản xuất.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật:
    • Thí nghiệm được thực hiện trong môi trường chậu/bầu nhựa kích thước cố định tại vùng khí hậu Thái Nguyên; các thông số bốc thoát hơi nước có thể dao động nhẹ khi áp dụng tại các vùng khô hạn (Duyên hải Nam Trung Bộ hoặc Tây Bắc mùa khô).
    • Chưa kết hợp khảo sát vi sinh vật cộng sinh vùng rễ (mycorrhiza) và nồng độ phối trộn phân bón NPK phân giải chậm.
  • Hướng phát triển tương lai:
    • Tích hợp hệ thống IoT vi khí hậu: Tự động điều chỉnh mái che thông minh dựa trên cảm biến bức xạ PAR quang học và độ ẩm đất tự động.
    • Thiết lập ô thí nghiệm khảo nghiệm ngoài thực địa để theo dõi tỷ lệ sống, độ bền bỉ cơ học và tốc độ khép tán của cây giống đã qua huấn luyện sinh thái 6 tháng.

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

  • Sinh viên Lâm nghiệp & Nông học: Tài liệu tham khảo chuẩn mực về phương pháp bố trí thí nghiệm RCBD, kỹ thuật đo cường độ thoát hơi nước Ivanov và phân tích phương sai ANOVA/Duncan bằng SPSS.
  • Kỹ sư lâm sinh & Quản lý vườn ươm: Bộ thông số định lượng rõ ràng ($50% \to 25%$ che sáng; $60\text{ml} \times 2$ lần/ngày) để áp dụng ngay vào quy trình sản xuất hàng loạt.
  • Chủ rừng & Doanh nghiệp trồng rừng gỗ lớn: Đảm bảo nguồn cung cây giống Phay đạt chuẩn sinh học cao, tăng tỷ lệ thành rừng và đẩy nhanh chu kỳ khai thác gỗ thương phẩm.
  • Nhà nghiên cứu sinh thái học thực vật: Cung cấp cơ sở dữ liệu thực nghiệm về mối quan hệ giữa diện tích lá, số lượng khí khổng và sự biến thiên áp suất bốc thoát hơi nước theo tuổi sinh trưởng.

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

1. Tại sao cây Phay 6 tháng tuổi lại có sinh khối cao nhất ở mức che sáng 25% thay vì 50%?

Ở 3 tháng đầu, bộ máy quang hợp còn non nớt nên che $50%$ giúp hạn chế cháy diệp lục. Bước sang 6 tháng, diện tích phiến lá mở rộng và mô giậu hoàn thiện, cây chuyển mạnh sang đặc tính ưa sáng của loài cây gỗ tiên phong; độ che $25%$ cung cấp cường độ bức xạ quang hợp tối ưu giúp gia tăng hiệu suất đồng hóa carbon ròng, dẫn đến sinh khối khô đạt cực đại $5{,}74\text{ g}$.

2. Tưới 120ml một lần duy nhất trong ngày có tương đương với 60ml chia làm 2 lần không?

Không tương đương. Tưới 1 lần $120\text{ml}$ gây ngập úng tạm thời, đẩy không khí ra khỏi giá thể bầu đất Feralit sa thạch và dễ rửa trôi khoáng vi lượng qua đáy bầu, sau đó đất nhanh chóng bị khô hạn vào giữa trưa. Chia nhỏ $60\text{ml} \times 2$ lần đảm bảo độ ẩm giá thể luôn giữ ở mức trương nước tối ưu ($75–80%$ dung tích ẩm đồng ruộng), thúc đẩy lông hút hoạt động liên tục.

3. Cường độ thoát hơi nước cao có phải là dấu hiệu xấu cho cây không?

Cường độ thoát hơi nước là động lực hút dòng mạch gỗ mang nước và khoáng từ rễ lên ngọn. Tuy nhiên, nếu không che sáng ($0%$), cường độ thoát hơi nước vọt lên $0{,}81\text{ g/dm}^2/\text{h}$ trong điều kiện rễ non chưa phát triển hoàn thiện sẽ gây mất nước nhanh hơn tốc độ hút, dẫn đến hiện tượng héo tạm thời và ức chế phân chia tế bào tầng phát sinh.

4. Tiêu chuẩn nào để chọn cây giống xuất vườn phục vụ trồng rừng?

Cây Phay 6 tháng tuổi đủ tiêu chuẩn trồng rừng khi đạt: $H_{vn} \ge 35\text{ cm}$, $D_{00} \ge 0{,}5\text{ cm}$, thân cây đã hóa gỗ chắc chắn, lá xanh đậm dày, ngọn phát triển cân đối và hệ thống rễ thứ cấp đan kín ruột bầu.

5. Chi phí đầu tư hệ thống giàn che di động có đắt không?

Hệ thống sử dụng cột tre/luồng hoặc ống kẽm kết hợp lưới đen cản quang PE 3 kim có độ bền 3–5 năm. Chi phí giàn che phân đoạn ước tính chỉ chiếm khoảng $3–5%$ tổng giá thành sản xuất mỗi vạn cây giống, được bù đắp hoàn toàn bởi tỷ lệ cây sống xuất vườn vượt trội.


Kết luận

Nghiên cứu đã giải quyết trọn vẹn bài toán sinh thái học gieo ươm cây Phay (Duabanga grandiflora) qua việc lượng hóa cụ thể hai nhân tố trọng yếu:

  • Chế độ che sáng tối ưu: Duy trì độ tàn che $50%$ ở giai đoạn 1–3 tháng tuổi và chuyển hạ về mức $25%$ ở giai đoạn 4–6 tháng tuổi để kích thích nở cổ rễ và tăng cường sinh khối khô.
  • Chế độ tưới nước tối ưu: Áp dụng định mức $60\text{ml}$/chậu/lần với tần suất 2 lần/ngày để tối đa hóa chiều cao ($40{,}15\text{ cm}$), đường kính cổ rễ ($0{,}66\text{ cm}$) và sinh khối đạt $4{,}19\text{ g}$ sau 4 tháng.

Kết quả thực nghiệm này là nền tảng kỹ thuật vững chắc để chuẩn hóa quy trình ươm giống cây Phay quy mô công nghiệp, góp phần bảo tồn nguồn gen bản địa quý hiếm và nâng cao giá trị gia tăng cho các dự án phát triển rừng trồng bền vững trên toàn quốc. Các đơn vị vườn ươm lâm nghiệp và bà con vùng trồng rừng nên áp dụng trực tiếp khung kỹ thuật này để nâng cao chất lượng giống cây trồng.