2  Données d’observation

2.1 Structure des données brutes

Les fichiers IVD contiennent:

• le rendement
• la date épiaison

2.2 Hypothèses sous jacente à la paramétrisation

Pour paramétrer la variété Apache, plusieurs hypothèses sont formulées vis à vis de la dépendance des stades phénologiques entre eux chez le blé tendre. En particulier, on pose des hypothèses pour estimer la date d’autres stades phénologiques, utiles à la calibration, à partir de la date d’épiaison.

L’ensemble des hypothèses est détaillé dans data_raw/obs/obs_hypotheses.yml.

Les dates estimées concernent :

• le stade épi 1 cm,
• stade feuille étendard (= LAI maximal)
• l’anthèse

Pour réaliser ces estimations, on s’appuie sur la notion de phyllochrone, à l’instar de la démarche de modélisation retenue dans Sirius. On s’appuie sur ces hypothèses pour lier les deux phénologies : celle du cycle végétatif et celle du cycle reproductif (qui, sous STICS, peuvent être parfaitement indépendants).

2.2.1 Variable directement observée

Le rendement grain est observé directement dans le dataset IVD sous la variable yield.obs, puis exporté vers data_YIELD_MAFRUIT_OBS.csv sous le nom STICS mafruit. L’unité d’origine est le q/ha. Au chargement de l’étape 4, un facteur de conversion × 0.1 est appliqué via obs_scale_factors.mafruit: 0.1 afin d’assurer la cohérence avec l’unité de simulation utilisée pour mafruit (t/ha).

2.2.2 Point d’ancrage observé pour la phénologie

La variable de départ pour la reconstruction des stades phénologiques est heading.obs, décrite dans le fichier d’hypothèses comme la date d’épiaison, codée en DOY de l’année year. Cette variable sert d’ancrage commun aux dérivations thermiques réalisées par scripts/00_prepare_obs_apache.R.

2.2.3 Estimation du stade “épi 1 cm” (iamfs)

iamfs n’est pas observé directement. Il est calculé en remontant dans le temps depuis heading.obs jusqu’à atteindre 4 phyllochrones, soit 4 × 112 = 448 °C.j.

Cette hypothèse correspond à une simplification de travail : le stade épi 1 cm est supposé survenir environ 4 phyllochrones avant l’épiaison pour Apache. Comme pour ilaxs, le cumul thermique repose sur une température journalière moyenne (tmin + tmax) / 2.

2.2.4 Estimation du stade “feuille étendard” (ilaxs)

ilaxs n’est pas observé directement non plus. Il est calculé en remontant dans le temps depuis heading.obs jusqu’à atteindre un cumul thermique de 2 phyllochrones, soit 2 × 112 = 224 °C.j.

Les hypothèses associées sont :

• phyllochrone Apache fixé à 112 °C.j (He et al. 2012) - phyllochrone moyen/médian chez le blé tendre
• la feuille-drapeau apparaît 2 phyllochrones avant l’épiaison
• la température journalière utilisée pour le cumul est (tmin + tmax) / 2

2.2.5 Estimation de l’anthèse

Dans STICS, il est suggéré que iflos = idrps pour une céréale. On propose donc une estimation de ces stades en avançant dans le temps depuis l’épiaison (heading.obs) jusqu’à atteindre 1.3 phyllochrone, soit 145.6 °C.j, avec le même phyllochrone Apache de 112 °C.j.

Dans l’implémentation actuelle, cette estimation est exportée dans data_PHENO_IDRPS_OBS.csv sous le nom de variable idrps pour alimenter l’étape 3 de calibration. Le document de calibration doit donc être lu comme une approximation opérationnelle, et non comme une observation directe de idrps.

2.3 Synthèse des fichiers

Fichier source	Variable	Étape
data_PHENO_IAMFS_OBS.csv	iamfs — stade épi 1 cm	Étape 1
data_PHENO_ILAXS_OBS.csv	ilaxs — stade feuille étendard	Étape 2
data_PHENO_IDRPS_OBS.csv	idrps — début remplissage grain	Étape 3
data_YIELD_MAFRUIT_OBS.csv	mafruit — rendement grain	Étape 4

2.3.1 Météo et couverture temporelle

Les dérivations thermiques utilisent la feuille Weather du classeur inputs/base/inputs_stics.xlsx. Ce choix garantit la cohérence entre la météo utilisée pour préparer les observations et celle utilisée ensuite dans les simulations STICS.