Differences between Morton (1983) and the Excel sheet in the Oudin's PET formula
- De : Alexandre Lhosmot [univ-fcomte.fr]
- Envoyé : mercredi 10 juin 2020 10:01
- À : @charles.perrin ; @vazken.andreassian
- Cc : Guillaume Bertrand [univ-fcomte.fr]
Je me permets de vous contacter à propos du calcul de l'ETP journalier que vous avez effectué ("Potential evapotranspiration (PE) calculation using the formulation proposed by Oudin et al., (2005)"). Mon directeur de thèse (Guillaume Bertrand, Chrono-Environnement) m'a transmis la feuille Excel de calcul que vous avez réalisé. Je souhaite calculer l'ETP journalier pour mon site d'étude (Tourbière de Frasne) où je possède déjà une chronique de température et de précipitation journalière sur 10 ans.
En regardant "Appendix C of the article by Morton, F,I, 1983, Operational estimates of areal evapotranspiration and their significance to the science and practice of hydrology, Journal of Hydrology 66 (1/4), 1-76," où les calculs de la radiation solaire sont explicités, je m’aperçois que certaines formules diffèrent entre la feuille Excel et l'article que vous mentionnez.
Je vous contacte donc dans le but d'éclaircir ces différences et de mieux comprendre les calculs. 'J' correspond au Julian day.
Pour le calcul de la déclinaison solaire :
-
Teta = 23.2 sin (29.5 * J - 94)
en degrés dans Morton (1983) -
Teta = 0.4093 * SIN(J / 58.1 - 1.405)
(Radian) dans la feuille Excel En divisant Teta de l'article par180/Pi
, je tombe quasiment sur la même valeur que celle de la feuille Excel, mais je ne comprends pas d'où vient la formule que vous utilisez.
Pour le calcul de Eta :
-
Eta = 1 + (1/60)*sin(29.5*J - 106)
dans Morton (1983) -
Eta = 1 + COS(J / 58.1) / 30
dans la feuille Excel PourJ = 1
, Eta est respectivement égal à0.985
et1.033
. Je ne comprends pas d'où vient cette différence.
Pour le calcul de la radiation solaire :
-
G = (1354 / Eta²) * (OM / 180) * cos(Pz)
dans l'article,341 W/m²
-
G = 446 * Eta * OM * cos(Pz)
dans la feuille Excel,116 W/m²/day
Je ne vois pas d'où vient le446
dans votre formule.
Enfin, pour le calcul de l'ETP :
- dans Oudin,
PE = (G / Lambda * Rho) * (T + 5) / 100
- dans la feuille,
PE = (G * (T + 5)) / 28.5 / 100
, j'imagine que vous avez directement mis les valeurs de Lambda et Rho dans l'équation. Cependant, en appliquant la formule de Oudin je ne tombe pas sur le même résultat, pourLambda = 2.45 MJ/kg
,Rho = 1000 Kg/m3
et en convertissant G en MJ/m²/day.
Mes incompréhensions viennent peut être de problèmes d'unité ou de conversion.