Giới thiệu dự án

Rau Tàu bay (Crassocephalum crepidioides (Benth.) S. Moore), thuộc họ Cúc (Asteraceae), là loài cây bản địa giàu giá trị dinh dưỡng và hoạt tính sinh học. Theo Tổ chức Y tế Thế giới (WHO), khoảng 80% dân số tại các quốc gia đang phát triển phụ thuộc vào các nguồn dược liệu tự nhiên trong chăm sóc sức khỏe ban đầu. Tại Việt Nam, hệ thống thực vật ghi nhận 3.948 loài có công dụng làm thuốc và hơn 150 loài rau rừng ăn được. Phân tích hóa thực vật học cho thấy sinh khối rau Tàu bay chứa 91,1% nước, 2,5% protein thực vật, 0,2% lipid, 1,6% cellulose, 3,7% dẫn xuất không protein, 0,9% khoáng tổng số, cùng hàm lượng vi chất chống oxy hóa cao gồm 3,4 mg% Caroten (tiền vitamin A) và 10 mg% Vitamin C.

Thực trạng khai thác tự nhiên mang tính tận diệt kết hợp với tập quán hái lượm phân tán đã khiến trữ lượng tự nhiên suy giảm nghiêm trọng. Việc thương mại hóa và thuần dưỡng cây rau Tàu bay gặp rào cản lớn do thiếu quy trình kỹ thuật chuẩn hóa về dinh dưỡng khoáng, dẫn đến năng suất không ổn định, chất lượng cây giống phân hóa mạnh và tỷ lệ sống thấp khi nhân rộng.

+-------------------------------------------------------------------------+
|                  VẤN ĐỀ VÀ MỤC TIÊU NGHIÊN CỨU                          |
+-------------------------------------------------------------------------+
| [Thực trạng]                                                            |
|  - Khai thác tự nhiên cạn kiệt nguồn gen quý.                          |
|  - Thiếu quy trình bón phân chuẩn hóa cho cây con giai đoạn vườn ươm.   |
|                                                                         |
| [Mục tiêu giải pháp]                                                    |
|  1. Xác định động thái sinh trưởng theo 4 công thức phân bón.           |
|  2. Tối ưu hóa tỷ lệ sống (>95%), chiều cao vút ngọn (Hvn), số lá, mầm.  |
|  3. Thiết lập thông số kỹ thuật chuẩn hóa để nhân rộng quy mô trang trại.|
+-------------------------------------------------------------------------+

Mục tiêu cụ thể của dự án nghiên cứu:

  1. Đánh giá ảnh hưởng định lượng của 3 nhóm phân bón (Phân chuồng hoai mục, Phân khoáng NPK 5:10:3, Phân hữu cơ vi sinh Sông Gianh) so với đối chứng đến tỷ lệ sống của cây giống.
  2. Xác định công thức phân bón tối ưu cho các chỉ tiêu sinh học: Chiều cao vút ngọn ($H_{vn}$), số lượng lá hữu hiệu ($SL$), số mầm nhánh sau chu kỳ 30 ngày theo dõi.
  3. Ứng dụng mô hình phân tích phương sai một nhân tố (One-way ANOVA) để kiểm định ý nghĩa thống kê ($p < 0.01$), thiết lập cơ sở dữ liệu cho quy trình canh tác lâm nghiệp và nông nghiệp công nghệ cao.

Phạm vi nghiên cứu được thực hiện tại Vườn thực nghiệm Khoa Lâm nghiệp – Trường Đại học Nông Lâm Thái Nguyên trên đất Feralit phát triển trên đá sa thạch (tầng canh tác $>40\text{ cm}$, thành phần cơ giới thịt nhẹ pha cát). Giới hạn đề tài tập trung vào giai đoạn vườn ươm cây giống từ hạt trong khung thời gian 30 ngày với chu kỳ theo dõi 10 ngày/lần.


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

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

Phương thức thu hái tự nhiên và canh tác truyền thống bộc lộ nhiều điểm hạn chế về năng suất sinh khối và độ ổn định chất lượng. Bảng đối chiếu dưới đây làm rõ ưu và nhược điểm của các giải pháp hiện hành:

Phương pháp Ưu điểm Nhược điểm Chi phí / Hiệu quả
Khai thác tự nhiên Chi phí đầu tư ban đầu bằng 0, không tốn công chăm sóc. Phụ thuộc mùa vụ, phá vỡ sinh thái, hàm lượng dược tính không đồng đều. Chi phí thu gom cao, rủi ro cạn kiệt nguồn gen.
Canh tác hữu cơ truyền thống (Phân chuồng) Cải tạo cấu trúc đất tốt, cung cấp mùn hữu cơ lâu dài. Tác dụng chậm, dinh dưỡng phân giải không đồng đều, tốn công vận chuyển (100–200 kg/ha). Hiệu quả trung bình ở giai đoạn vườn ươm ngắn ngày.
Canh tác NPK cân đối (Đề xuất) Dinh dưỡng dễ tiêu, hiệu lực tức thì, thúc đẩy phân hóa chồi mầm và kéo dài rễ nhanh. Yêu cầu kiểm soát liều lượng nghiêm ngặt để tránh rửa trôi và thoái hóa đất. Chi phí vật tư thấp ($0,2\text{ kg}/2\text{ m}^2$), hiệu quả kinh tế và sinh khối vượt trội.

Hệ thống yêu cầu kỹ thuật đối với mô hình vườn ươm được phân loại theo mô hình MoSCoW:

  • Must have (Bắt buộc): Khả năng thoát nước luống $>95%$, tỷ lệ sống cây con $\ge 90%$, kiểm soát độ ẩm đất đạt $75 - 80%$.
  • Should have (Nên có): Ứng dụng phân tích phương sai ANOVA để đánh giá sai số thí nghiệm, phân tầng chất lượng cây giống (Tốt, Trung bình, Xấu).
  • Could have (Có thể có): Hệ thống tưới rãnh bán tự động (tưới ngấm qua dõng luống đạt $1/3 - 1/2$ chiều cao gờ).
  • Won't have (Chưa áp dụng): Hệ thống nhà kính kiểm soát quang chu kỳ tự động ở giai đoạn thử nghiệm cơ sở.

Thiết kế hệ thống

Thí nghiệm được thiết kế theo khối ngẫu nhiên đầy đủ (Randomized Complete Block Design - RCBD) với 4 công thức thí nghiệm (CTTN), mỗi công thức lặp lại 3 lần trên ô tiêu chuẩn $2\text{ m}^2$, mật độ 25 cây/ô ($100\text{ cây/công thức}$):

  • CT1: Phân chuồng hoai mục (bón lót $0,4\text{ kg}/\text{ô}$).
  • CT2: Phân khoáng NPK 5:10:3*KS ($5%\text{ N}$, $10%\text{ P}_2\text{O}_5$, $3%\text{ K}_2\text{O}$, $18-20%\text{ CaO}$, $8-10%\text{ S}$, $4-5%\text{ SiO}_2$, $2-2.5%\text{ MgO}$, $20-30\text{ ppm Cu}$, $40-50\text{ ppm Zn}$ - bón lót $0,2\text{ kg}/\text{ô}$).
  • CT3: Phân hữu cơ vi sinh Sông Gianh (Chất hữu cơ $15%$, Axit Humic $1,5%$, $\text{P}_2\text{O}_5\text{ hh } 3%$, mật độ vi sinh vật Azotobacter sp., Bacillus sp., Streptomyces sp. $\ge 1\times 10^6\text{ CFU/g}$ - bón lót $0,15\text{ kg}/\text{ô}$).
  • CT4: Đối chứng (Không bón phân).

Cấu trúc cơ sở dữ liệu quan trắc sinh trắc học được chuẩn hóa theo lược đồ SQL phục vụ quản lý dữ liệu thí nghiệm:

CREATE TABLE nursery_experiment_logs (
    log_id INT PRIMARY KEY AUTO_INCREMENT,
    formula_code VARCHAR(10) NOT NULL, -- CT1, CT2, CT3, CT4
    replication_index TINYINT NOT NULL, -- 1, 2, 3
    day_interval TINYINT NOT NULL,      -- 10, 20, 30
    plant_height_cm DECIMAL(5,2) NOT NULL,
    leaf_count INT NOT NULL,
    shoot_count INT NOT NULL,
    survival_status BOOLEAN DEFAULT TRUE,
    quality_grade ENUM('GOOD', 'MEDIUM', 'POOR') NOT NULL,
    recorded_at DATE NOT NULL
);

Methodology

Quy trình nghiên cứu áp dụng phương pháp nghiên cứu thực nghiệm nông học kết hợp mô hình phân tích sinh trắc học hiện đại:

  • Chuẩn bị nền đất: Làm đất tơi xốp, dọn sạch tàn dư cỏ dại, lên luống cao $25 - 30\text{ cm}$, bề mặt luống phẳng, rãnh thoát nước sâu $30\text{ cm}$.
  • Khoảng cách gieo cấy: Hàng cách hàng $30\text{ cm}$, cây cách cây $20\text{ cm}$.
  • Đo đếm chỉ tiêu: Chiều cao vút ngọn ($H_{vn}$) đo bằng thước khắc vạch độ chính xác $\pm 0,1\text{ cm}$ từ cổ rễ đến đỉnh sinh trưởng; đếm toàn bộ số lá thật và số mầm nhánh cấp 1.
  • Xử lý thống kê: Sử dụng mô hình ANOVA một nhân tố để kiểm định sự sai khác giữa các giá trị trung bình theo công thức:

$$\text{SS}{\text{total}} = \text{SS}{\text{between}} + \text{SS}_{\text{within}}$$

$$F_{\text{calc}} = \frac{\text{MS}{\text{between}}}{\text{MS}{\text{within}}} = \frac{\text{SS}{\text{between}} / (k - 1)}{\text{SS}{\text{within}} / (N - k)}$$

Nếu $F_{\text{calc}} > F_{\text{crit}}(\alpha, k-1, N-k)$ với $\alpha = 0.05$, bác bỏ giả thuyết $H_0$ và kết luận có sự khác biệt có ý nghĩa thống kê giữa các công thức phân bón.


Implementation và kết quả

Development process

Quá trình triển khai thực địa được phân chia theo 4 giai đoạn chính:

  • Pha 1 (Ngày 0): Xử lý hạt giống, phân tích hóa tính đất Feralit tầng $0-10\text{ cm}$, bón lót phân bón theo đúng định mức thiết kế của từng công thức.
  • Pha 2 (Ngày 1 - 10): Theo dõi quá trình hồi xanh, đo đếm chỉ tiêu chu kỳ 1.
  • Pha 3 (Ngày 11 - 20): Kiểm soát cỏ dại, xới xáo phá váng tầng mặt, đo đếm chỉ tiêu chu kỳ 2.
  • Pha 4 (Ngày 21 - 30): Đánh giá động thái tăng trưởng cực đại, phân loại phẩm chất và chạy mô hình ANOVA.

Script phân tích tự động dữ liệu thí nghiệm bằng Python:

import numpy as np
import scipy.stats as stats

# Dữ liệu chiều cao (Hvn - cm) tại ngày thứ 30 của 4 công thức (3 lần lặp)
data_ct1 = [32.10, 36.20, 34.72]  # Trung bình: 34.34 cm
data_ct2 = [47.10, 49.50, 48.12]  # Trung bình: 48.24 cm
data_ct3 = [24.50, 26.80, 25.53]  # Trung bình: 25.61 cm
data_ct4 = [21.80, 24.10, 23.34]  # Trung bình: 23.08 cm

# Chạy One-way ANOVA
f_val, p_val = stats.f_oneway(data_ct1, data_ct2, data_ct3, data_ct4)

print(f"--- KẾT QUẢ KIỂM ĐỊNH ANOVA CHIỀU CAO (HVN) ---")
print(f"Giá trị F tính toán: {f_val:.4f}")
print(f"Giá trị P-value: {p_val:.4e}")

if p_val < 0.05:
    print("Kết luận: Chấp nhận đối thuyết H1 (Khác biệt có ý nghĩa thống kê cực kỳ cao)")

Testing và validation

Số liệu quan sát thực nghiệm sau 30 ngày được tổng hợp qua bảng ma trận chỉ tiêu sinh học:

Công thức thí nghiệm Tỷ lệ sống 30 ngày (%) Chiều cao $H_{vn}$ TB (cm) Số lá TB (lá/cây) Số mầm TB (mầm/cây) Tỷ lệ cây Tốt (%) Tỷ lệ cây Xấu (%)
CT1 (Phân chuồng) 93,33 34,34 16,93 5,60 85,71 4,29
CT2 (Phân NPK) 97,33 48,24 23,26 7,13 89,04 4,68
CT3 (Vi sinh) 90,67 25,61 13,22 4,25 83,82 4,41
CT4 (Đối chứng) 89,33 23,08 11,62 3,24 68,65 13,44

Kết quả phân tích phương sai ANOVA một nhân tố cho thấy sự vượt trội tuyệt đối của công thức NPK:

  • Tỷ lệ sống: $F_{\text{calc}} = 10,05 > F_{\text{crit}} = 4,07$ ($p = 0,0043 < 0,01$).
  • Chiều cao vút ngọn ($H_{vn}$): $F_{\text{calc}} = 149,38 > F_{\text{crit}} = 4,07$ ($p = 2,31\times 10^{-7}$).
  • Số lượng lá hữu hiệu: $F_{\text{calc}} = 36,01 > F_{\text{crit}} = 4,07$ ($p = 5,41\times 10^{-5}$).
  • Số mầm nhánh: $F_{\text{calc}} = 53,26 > F_{\text{crit}} = 4,07$ ($p = 1,25\times 10^{-5}$).
[Biểu đồ phân bố động thái tăng trưởng chiều cao Hvn qua các chu kỳ (cm)]
CT2 (NPK)        : |==== 14.65cm ====> 34.70cm ========> 48.24cm
CT1 (Chuồng hoai): |=== 12.56cm =====> 22.12cm ======> 34.34cm
CT3 (Vi sinh)    : |== 10.28cm ====> 13.72cm ====> 25.61cm
CT4 (Đối chứng)  : |== 10.51cm ====> 13.54cm ===> 23.08cm
                   +-------------------+-------------------+
                   Ngày 10             Ngày 20             Ngày 30

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

  1. Hiệu suất sinh trưởng vượt bậc: Công thức CT2 (NPK 5:10:3) đạt chiều cao trung bình 3 chu kỳ là $32,53\text{ cm}$, cao gấp 2,07 lần so với đối chứng CT4 ($15,71\text{ cm}$). Tại thời điểm kết thúc 30 ngày, cây ở công thức CT2 đạt $48,24\text{ cm}$, bứt phá hoàn toàn so với CT1 ($34,34\text{ cm}$) và CT3 ($25,61\text{ cm}$).
  2. Khả năng phân hóa cơ quan sinh dưỡng: Số lá hữu hiệu ở CT2 đạt $23,26\text{ lá/cây}$ (vượt đối chứng $+100,2%$), số mầm nhánh đạt $7,13\text{ mầm/cây}$ (vượt đối chứng $+120,1%$).
  3. Chất lượng thương phẩm cây giống: CT2 đạt tỷ lệ cây loại Tốt $89,04%$, giảm thiểu tỷ lệ cây còi cọc và sâu bệnh xuống mức tối thiểu ($4,68%$), đảm bảo tiêu chuẩn xuất vườn phục vụ trồng đại trà.

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

Nghiên cứu mang lại những đóng góp khoa học và kỹ thuật ứng dụng rõ nét:

  • Thiết lập công thức bón lót tối ưu: Lần đầu tiên định lượng chính xác hiệu quả của dòng phân NPK giàu Lân và Canxi (5:10:3*KS) đối với loài rau dại bản địa Crassocephalum crepidioides trên nền đất Feralit chua nghèo dinh dưỡng. Lượng lân dễ tiêu ($10%\text{ P}_2\text{O}_5$) kết hợp với hàm lượng Canxi khử chua ($18-20%\text{ CaO}$) đã giải phóng bộ rễ, kích thích lông hút phát triển mạnh mẽ ngay từ giai đoạn 10 ngày đầu.
  • Cải thiện hiệu suất định lượng so với các giải pháp hiện hành:
    • Tăng trưởng chiều cao vượt $+108,9%$ so với đối chứng tự nhiên và $+40,5%$ so với phân chuồng truyền thống.
    • Tăng trưởng số mầm nhánh vượt $+120,1%$ so với đối chứng, tạo tiền đề nâng cao năng suất cắt ngọn thương phẩm cho các chu kỳ thu hoạch sau này.
  • Đóng góp học thuật: Cung cấp bộ dữ liệu sinh trắc học có giá trị kiểm định ANOVA chuẩn mực cho nguồn gen cây dược liệu - rau rừng Việt Nam, hỗ trợ công tác bảo tồn gen ex-situ tại các vườn thực vật chuyên ngành.

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

Kịch bản ứng dụng thực địa

  • Mô hình nông lâm kết hợp: Trồng xen cây rau Tàu bay dưới tán rừng trồng (Keo, Bạch đàn, Quế giai đoạn 1–3 năm tuổi) hoặc vườn cây ăn quả để tận dụng bóng râm nhẹ và tăng thu nhập ngắn hạn.
  • Vùng chuyên canh rau dược liệu sạch: Triển khai theo tiêu chuẩn VietGAP phục vụ chuỗi cung ứng thực phẩm đặc sản cho đô thị và làm nguyên liệu chiết xuất tinh dầu, flavonoid trong công nghiệp dược.
+-----------------------------------------------------------------------------------+
|               HƯỚNG DẪN KỸ THUẬT CANH TÁC RAU TÀU BAY (QUY MÔ 1 HA)              |
+-----------------------------------------------------------------------------------+
| 1. Chuẩn bị đất:                                                                  |
|    - Cày bừa sâu 25-30cm, khử trùng vôi bột 300-400 kg/ha.                        |
|    - Lên luống cao 25-30cm, rộng 1.0-1.2m, rãnh rộng 30cm.                        |
| 2. Bón lót (CT2 chuẩn hóa):                                                       |
|    - Phân NPK 5:10:3*KS: 80 - 100 kg/ha (bón rải theo rãnh hoặc hố, lấp đất mỏng). |
| 3. Mật độ cấy:                                                                    |
|    - Khoảng cách: 30cm x 20cm (mật độ tương đương 140.000 - 160.000 cây/ha).      |
| 4. Chế độ nước:                                                                   |
|    - Duy trì độ ẩm đất 75 - 80%. Ưu tiên tưới rãnh dẫn nước ngập 1/3 luống.       |
| 5. Thu hoạch & Cách ly:                                                           |
|    - Cắt tỉa ngọn non khi cây đạt 35-40cm. Ngừng bón phân trước thu hoạch 7 ngày.|
+-----------------------------------------------------------------------------------+

Phân tích hiệu quả kinh tế (Dự toán trên 1 ha canh tác)

Hạng mục chi phí / Doanh thu Định mức kỹ thuật Đơn giá dự kiến (VNĐ) Thành tiền (VNĐ)
Giống cây con từ hạt 150.000 cây 300 đ/cây 45.000.000
Phân bón NPK 5:10:3 100 kg (Bón lót) + 150 kg (Thúc) 8.000 đ/kg 2.000.000
Công làm đất, lên luống, tưới tiêu 60 công 250.000 đ/công 15.000.000
Vật tư phụ & Bảo vệ thực vật sinh học Trọn gói 1 chu kỳ vụ - 8.000.000
TỔNG CHI PHÍ ĐẦU TƯ - - 70.000.000
Sản lượng ngọn thương phẩm 8.000 kg/vụ (3 đợt cắt) 20.000 đ/kg 160.000.000
LỢI NHUẬN THUẦN DỰ KIẾN Chu kỳ 60 - 75 ngày - 90.000.000

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

Hạn chế kỹ thuật

  • Thí nghiệm mới khảo sát giai đoạn vườn ươm 30 ngày đầu, chưa theo dõi kéo dài qua các lứa thu hoạch sinh khối liên tục (lứa 2, lứa 3).
  • Chưa tích hợp khảo nghiệm phân bón qua lá (Foliar feeding) và các tỷ lệ phối trộn hữu cơ - vô cơ phức hợp.
  • Chưa phân tích sâu biến động hàm lượng hoạt chất thứ cấp (Flavonoid, hợp chất phenolic, tinh dầu) dưới tác động của từng loại phân bón.

Hướng phát triển tiếp theo

  • Thiết lập thí nghiệm bón phân đa nhân tố ($N \times P \times K$ theo ma trận trực giao Box-Behnken) để tìm điểm tối ưu biên.
  • Ứng dụng kỹ thuật IoT trong kiểm soát độ ẩm đất và châm phân tự động (Fertigation).
  • Đánh giá khả năng ức chế tế bào ung thư và kháng viêm của dịch chiết từ nguồn rau trồng thuần hóa so với nguồn rau tự nhiên.

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

  • Sinh viên & Học viên Nông Lâm nghiệp: Tiếp cận quy trình bố trí thí nghiệm khối ngẫu nhiên (RCBD) chuẩn mực và phương pháp xử lý dữ liệu sinh học bằng ANOVA một nhân tố.
  • Kỹ sư nông học & Chủ trang trại: Nắm bắt quy trình kỹ thuật bón lót phân NPK 5:10:3 để rút ngắn chu kỳ vườn ươm cây giống xuống còn 20–25 ngày, tỷ lệ cây đạt chuẩn xuất vườn $>89%$.
  • Hợp tác xã nông nghiệp miền núi: Có mô hình khả thi để chuyển đổi đất nương rẫy bạc màu, đa dạng hóa cơ cấu cây trồng đặc sản có giá trị kinh tế cao ($90.000.000\text{ VNĐ/ha/vụ}$).
  • Nhà nghiên cứu Dược liệu: Cơ sở dữ liệu hình thái và sinh thái học làm tiền đề cho các nghiên cứu chuẩn hóa nguồn nguyên liệu đầu vào theo tiêu chuẩn GACP-WHO.

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

1. Yêu cầu kỹ thuật cốt lõi để triển khai thành công vườn ươm rau Tàu bay là gì?

Cần lựa chọn đất Feralit hoặc đất thịt nhẹ tơi xốp, thoát nước tuyệt đối để tránh bệnh lở cổ rễ. Lên luống cao $25-30\text{ cm}$, duy trì ẩm độ đất $75-80%$. Sử dụng phân NPK 5:10:3 dạng hạt bón lót đều với liều lượng $0,2\text{ kg}/2\text{ m}^2$ luống.

2. Phân bón vi sinh (CT3) có thành phần vi sinh vật hữu ích nhưng tại sao sinh trưởng lại thấp hơn NPK (CT2) ở giai đoạn vườn ươm?

Phân hữu cơ vi sinh cần thời gian để hệ vi sinh vật (Bacillus, Azotobacter) thiết lập mật độ và phân giải chất hữu cơ phức tạp thành dạng ion dễ tiêu. Ở chu kỳ 30 ngày ngắn ngủi của cây con, phân NPK cung cấp ngay lập tức các ion khoáng $\text{NO}_3^-$, $\text{NH}_4^+$, $\text{H}_2\text{PO}_4^-$ và $\text{K}^+$, giúp rễ non hấp thu tức thì.

3. Có thể kết hợp trồng xen rau Tàu bay vào các hệ thống canh tác hiện có không?

Rau Tàu bay có biên độ sinh thái rộng, chịu bóng râm nhẹ rất tốt. Cây hoàn toàn thích hợp trồng xen dưới tán rừng keo 1-2 năm tuổi, vườn cam, bưởi, hoặc xen canh gối vụ trong vườn đồi để chống xói mòn đất dốc.

4. Quy trình quản lý sâu bệnh và cách ly phân bón trước thu hoạch như thế nào?

Cần nhổ bỏ và tiêu hủy ngay các cá thể có biểu hiện thối nhũn hoặc sâu ăn lá. Sử dụng các chế phẩm thảo mộc, dầu khoáng sinh học nếu có dịch hại. Tuyệt đối dừng bón phân và các hoạt chất bảo vệ thực vật trước khi thu hái đọt non từ 5 đến 7 ngày.

5. Khả năng mở rộng quy mô (Scalability) và thời gian hoàn vốn đầu tư?

Mô hình có khả năng mở rộng nhanh do nhân giống bằng hạt với hệ số nhân giống cao. Chu kỳ canh tác ngắn ($30\text{ ngày}$ vườn ươm, $45\text{ ngày}$ sau trồng bắt đầu thu hoạch), giúp thu hồi vốn đầu tư ngay trong vụ đầu tiên sau $2,5 - 3\text{ tháng}$.


Kết luận

Đề tài nghiên cứu đã chứng minh một cách khoa học và định lượng vai trò quyết định của chế độ dinh dưỡng khoáng đối với sự sinh trưởng của cây rau Tàu bay (Crassocephalum crepidioides) tại vườn ươm. Việc áp dụng công thức phân bón NPK 5:10:3 mang lại hiệu quả vượt trội: tỷ lệ sống đạt $97,33%$, chiều cao vút ngọn đạt $48,24\text{ cm}$, số lá hữu hiệu đạt $23,26\text{ lá/cây}$, số mầm nhánh đạt $7,13\text{ mầm/cây}$ và tỷ lệ phẩm chất tốt đạt $89,04%$ với độ tin cậy thống kê cao ($p < 0,01$).

Kết quả này mở ra tiềm năng to lớn trong việc thuần hóa, nhân rộng và thương mại hóa các loài rau rừng bản địa giàu dược tính, góp phần phát triển kinh tế sinh thái vùng cao, nâng cao thu nhập cho đồng bào nông lâm dân và bảo tồn nguồn tài nguyên đa dạng sinh học quý giá của Việt Nam.