Shapiro-Wilk-Tests zeigten, dass die Verteilungen aller gemessenen Variablen in beiden Anbausaisons nicht signifikant von der Normalverteilung abwichen (p > 0,05). Die Varianzanalyse ergab, dass die Wechselwirkung zwischen Komposttee-Konzentration (Ct) und Blattdüngung mit Biostimulanzien (F) die Wachstumseigenschaften, die Zahl der Dolden sowie ölbasierte Merkmale (Ölanteil, Ölproduktion pro Pflanze und Ölproduktion pro Hektar) in beiden Anbausaisons signifikant beeinflusste.
Ertragsreaktion im Vergleich zur Kontrollbehandlung
In der ersten Saison (Abbildung 2) stieg der Ölertrag pro Hektar (Oy ha⁻1) stetig mit zunehmender Kompostteemenge an. Alle Ct3-Behandlungen (600 L ha⁻1) erzielten höhere Ölerträge als die Kontrollbehandlung (T1: Ct1 × F1, 360 L ha⁻1 ohne Blattsprühung). Den höchsten Ölertrag (55,58 kg ha⁻1) erhielt man bei Kombination von Ct3 mit Spirulina in einer Konzentration von 2 cm3 L⁻1 (T13), der signifikant über der Kontrolle und allen Ct1-Behandlungen lag und den LSD-Wert (2,93 kg ha⁻1) übertraf. Im Vergleich zur Kontrolle zeigte diese Behandlung eine deutliche Steigerung des Ölertrags pro Hektar (Abbildung 3).

Abbildung 3: Biplot-Hauptkomponentenanalyse für die fünfzehn kombinierten Behandlungen, basierend auf vier vegetativen Merkmalen, der Anzahl der Dolden pro Pflanze und drei ölbasierten Merkmalen, unter Verwendung der ersten und zweiten Hauptkomponenten; jede Behandlung wird durch einen Punkt dargestellt. Die Position innerhalb des Diagramms zeigt das Ausmaß der Ähnlichkeit oder Unterschiedlichkeit zwischen den Behandlungen gemäß den Eingabefaktoren (8 Wachstums- und Ölmerkmale); je näher beieinander, desto ähnlicher, und je weiter auseinander, desto unterschiedlicher. PH: Pflanzenhöhe (cm), NB: Anzahl der Seitentriebe/Pflanze, HFW: Gesamtfrischgewicht des Krauts (g), HDW: Gesamttrockengewicht des Krauts (g), O%: Gehalt an ätherischem Öl, Oy/p: Ölproduktion pro Pflanze (cm3), Oy/ha: Ölproduktion pro Hektar (L) und NU: Anzahl der Dolden/Pflanze T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4) und T15 (Ct3 × F5). Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Ähnliche Trends zeigten sich bezüglich des Ölertrags pro Pflanze und der Anzahl der Dolden, was darauf hinweist, dass Behandlungen mit höherem Ölertrag pro Hektar auch mit einer stärkeren Fortpflanzungsleistung und größerer Biomasseakkumulation verbunden waren. Diese Zusammenhänge werden als statistische Assoziationen zwischen gemessenen Merkmalen dargestellt, ohne dass kausale Wirkmechanismen abgeleitet werden.
In der zweiten Saison erzielte Ct3 × F3 (T13) erneut den höchsten Ölgehalt pro Hektar (58,85 kg ha⁻1), der signifikant über der Kontrolle und den meisten Ct2-Behandlungen lag (LSD = 3,17 kg ha⁻1). Die Rangfolge der Behandlungen im Vergleich zur Kontrolle war über beide Saisons weitgehend konsistent, was auf reproduzierbare Behandlungseffekte unter den getesteten Umweltbedingungen hinweist.
Hauptkomponentenanalyse (PCA)
Die Hauptkomponentenanalyse wurde als exploratives Werkzeug angewandt, um multivariate Beziehungen zwischen den acht gemessenen Merkmalen über die 15 Behandlungen hinweg zusammenzufassen. Die ersten beiden Hauptkomponenten erklärten 98,4 % der gesamten Varianz, wobei PC₁ 96,5 % und PC₂ 1,9 % ausmachte (Tabelle 2;Abbildung 4). Die Dominanz von PC₁ spiegelt eine starke Kovarianz zwischen wachstums- und ertragsbezogenen Merkmalen im Datensatz wider, anstatt klar voneinander getrennte zugrundeliegende biologische Prozesse.
| Parameter | PC1 | PC2 | PC3 | PC4 |
| NU | 0.3516 | 0.3064 | -0.3732 | 0.772 |
| O% | 0.349 | -0.6068 | -0.1873 | -0.0634 |
| Oy/p | 0.3591 | -0.11 | -0.0028 | -0.0974 |
| Oy/ha | 0.3583 | -0.1579 | -0.1046 | -0.2021 |
| PH | 0.3476 | 0.5757 | 0.3834 | -0.3296 |
| NB | 0.3495 | -0.26 | 0.7504 | 0.3693 |
| HFW | 0.3551 | 0.3138 | -0.1821 | -0.1535 |
| HDW | 0.358 | -0.0575 | -0.2675 | -0.2847 |
| Eigenwert | 7.7201 | 0.1518 | 0.0721 | 0.0346 |
| % der Varianz | 96,5 | 1,9 | 0,9 | 0,43 |
| Kumulativ (%) | 96,5 | 98,4 | 99,3 | 99,73 |
| PH: Pflanzenhöhe (cm), NB: Anzahl der Äste/Pflanze, HFW: gesamte frische Krautmasse (g), HDW: gesamte trockene Krautmasse (g), O%: Gehalt an ätherischem Öl, Oy/p: Öl-Ausbeute pro Pflanze (cm³), Oy/ha: Öl-Ausbeute pro Hektar (L) und NU: Anzahl der Dolden/Pflanze |
Tabelle 2: Eigenwerte, Eigenvektoren, prozentuale Varianz und kumulierte Varianz der ersten vier Hauptkomponenten basierend auf acht Wachstums- und ölbasierten Merkmalen über 15 Behandlungskombinationen hinweg. Die Abkürzungen für die Merkmale lauten wie folgt: PH (Pflanzenhöhe, cm), NB (Anzahl der Seitenzweige pro Pflanze), HFW (Frischmasse des Krauts, g), HDW (Trockenmasse des Krauts, g), O% (Gehalt an ätherischem Öl), Oy/p (Ölausbeute pro Pflanze, cm3), Oy/ha (Ölausbeute pro Hektar, L) und NU (Anzahl der Dolden pro Pflanze).

Abbildung 4: Ladeplot des Eigenvektors für 8 untersuchte Merkmale, behandelt mit 15 wechselwirkenden Behandlungen: Die Balken für PC₁ verwenden die linke Achse; PC₂, PC₃, PC₄ verwenden unterschiedliche Symbolformen und Farben auf der rechten Achse. PH: Pflanzenhöhe (cm), NB: Anzahl der Seitentriebe/Pflanze, HFW: Gesamtfrischmasse der Pflanze (g), HDW: Gesamttrockenmasse der Pflanze (g), O%: Gehalt an ätherischem Öl, Oy/p: Öl-Ausbeute pro Pflanze (cm3), Oy/ha: Öl-Ausbeute pro Hektar (L) und NU: Anzahl der Dolden/Pflanze. Fehlerbalken zeigen die Standardfehler des Mittelwerts (SEM). Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Alle Merkmale wiesen starke und ähnliche positive Ladungen auf der PC₁ auf, einschließlich der Anzahl der Schirme, des Ölanteils, des Ölertrags pro Pflanze, des Ölertrags pro Hektar, der Pflanzenhöhe, der Anzahl der Seitenzweige, der frischen und der trockenen Krautmasse (Tabelle 2; Abbildung 5). Demnach kann die PC₁ als allgemeiner Produktivitätsgradient interpretiert werden, der die Variation in Biomasse- und Ölertragsmerkmalen über die Behandlungen hinweg zusammenfasst. Die PC₂ erklärte nur einen geringen Anteil der Varianz und trug nur begrenzt zur zusätzlichen Differenzierung zwischen den Behandlungen bei.

Abbildung 5: PCA-Score-Plot zur Veranschaulichung der Verteilung von 15 Behandlungskombinationen basierend auf den gemessenen Merkmalen. Die Behandlungen sind wie folgt definiert: T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4) und T15 (Ct3 × F5). Die Fehlerbalken zeigen die Standardfehler des Mittelwerts (SEM). Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
PCA-Scores und Vergleich mit der Kontrolle
Behandlungen mit hohen positiven PC₁-Werten (insbesondere T13, T12, T15 und T14) wiesen ebenfalls die größten prozentualen Zunahmen im Vergleich zur Kontrolle (T1) bei den meisten Merkmalen auf (Abbildung 5 und Abbildung 6). Beispielsweise zeigte T13 eine Zunahme von 94,44 % bei der Anzahl der Dolden und von 98,18 % beim Ölertrag pro Hektar im Vergleich zur Kontrolle.

Abbildung 6: Relative prozentuale Zunahme der gemessenen Merkmale im Vergleich zur Kontrollbehandlung (T1: Ct1 × F1). Die Behandlungen sind wie folgt definiert: T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4) und T15 (Ct3 × F5). Bitte klicken Sie hier, um eine vergrößerte Ansicht dieser Abbildung einzusehen.
Behandlungen mit intermediären PC₁-Werten (z. B. T8 und T10) zeigten moderate Verbesserungen gegenüber der Kontrolle, während Behandlungen mit negativen PC₁-Werten, einschließlich der Kontrolle und anderer Ct1-Kombinationen, nur minimale relative Veränderungen aufwiesen. Diese Muster deuten darauf hin, dass die PC₁-Werte eng mit den baselinerichtigten Behandlungsantworten korrelierten, was die Ergebnisse der univariaten ANOVA untermauert. Die PCA wurde jedoch ausschließlich zur Datenreduktion und Visualisierung multivariater Muster verwendet und ersetzt keine inferenzstatistischen Tests.
Merkmalsassoziationen im PCA-Biplot
Die PCA-Biplot-Darstellung (Abbildung 4) verdeutlicht enge Zusammenhänge zwischen Öl-Ausbeute pro Pflanze, Öl-Ausbeute pro Hektar und Trockengewicht des Krauts, was auf starke Korrelationen zwischen diesen Variablen hinweist. Pflanzenhöhe, Frischgewicht des Krauts und Anzahl der Dolden wiesen ebenfalls ähnliche Vektororientierungen auf, was auf eine koordinierte Variation dieser Merkmale innerhalb des Datensatzes schließen lässt. Die Anzahl der Seitentriebe zeigte eine teilweise Trennung von anderen Ertragskomponenten, was auf schwächere Assoziationen unter den vorliegenden Versuchsbedingungen hindeutet. Diese Beziehungen beschreiben lediglich beobachtete Korrelationen und implizieren keine mechanistischen Verknüpfungen.
Behandlungen (T12–T15), die Komposttee in einer Menge von 600 L ha⁻1 in Kombination mit Blattdüngern erhielten, zeigten eine überlegene multivariate Leistung hinsichtlich Wachstums- und Ölmerkmalen. Im Gegensatz dazu wiesen die Behandlungen (T7–T11) vergleichsweise niedrigere Werte bei den meisten untersuchten Parametern auf. Behandlungen, die in entgegengesetzter Richtung zu den meisten Vektoren lagen, wie beispielsweise Ct1 (360 L ha⁻1) ohne Blattdüngung oder in Kombination mit Nostoc, waren mit niedrigeren Werten bei den meisten gemessenen Merkmalen assoziiert.
Insgesamt lieferte die PCA eine prägnante multivariate Zusammenfassung, die die ANOVA ergänzt, indem sie verdeutlicht, wie sich die Behandlungen hinsichtlich mehrerer Merkmale gleichzeitig vom Kontrollansatz unterschieden. Die Übereinstimmung zwischen den Mustern der PCA und den univariaten Ergebnissen nach Basislinienanpassung untermauert die Robustheit der Behandlungseffekte, ohne die Interpretation über die gemessenen Variablen hinaus auszuweiten.
VERFÜGBARKEIT VON DATEN:
Die Rohdaten für diese Arbeit sind in den ergänzenden Dateien verfügbar.
Zusatz-Tabelle 1: Rohdaten der Experimente zu den Ergebnissen. Bitte klicken Sie hier, um diese Datei herunterzuladen.