En considérant qu’environ un quart de l’énergie solaire contribue à l’évaporation de l’eau :
Données : Puissance de l’énergie solaire reçue par la terre : P=1,74 $10^{17}$ W, Rayon de la terre : R=6370 km, Masse d’eau sous forme de vapeur d’eau dans l’atmosphère : $M_{eau}=13~10^{15}$ kg, Chaleur latente d’évaporation : $\Delta H_{vap}=2.5~MJ/kg$
$$Q=\frac{0,25 P}{\Delta H_{vap}}$$
- La masse d'eau évaporée par s (kg/s) -qu'on peut aussi appeler flux ou débit- peut se calculer en divisant la puissance de l'énergie solaire recue par la terre et contribuant à l'évaporation en (J/s) par l'énergie nécessaire pour évaporer un kg d'eau (J/kg)
#Données
P=1.74e17 #W ou J/s
R=6370 #km
M_vap_eau=13e15 #kg
DH_eau=2.5e6 #J/kg
pourc=0.25
rho_eau=1000 #kg m-3
#Changements d'unités
sec_par_an=365*24*60*60
m3_par_km3=1000**3
#Puissance utilisée pour l'évaporation
P_evap=P*pourc #J/s
#Quantité évaporée
Q_evap=P_evap/DH_eau #kg/s
Q_evap1=Q_evap*sec_par_an #kg/an
Q_evap2=Q_evap*sec_par_an/rho_eau/m3_par_km3 #km3/an
print ('Le débit d eau évaporée est {:.2e} kg/s soit {:.2e} kg/an soit {:.2e} km3/an'.format(Q_evap,Q_evap1,Q_evap2))
Le débit d eau évaporée est 1.74e+10 kg/s soit 5.49e+17 kg/an soit 5.49e+05 km3/an
$$P_{par m2}=\frac{P}{4 \pi R^2}$$
- En moyenne, la puissance de l'énergie solaire reçue sur la terre par m2 est la puissance totale reçue (en J/s) divisée par la surface de la terre :
à noter qu'il faut utiliser la surface de toute la terre (et pas seulement la partie éclairée) puisqu'on effectue une moyenne de la puissance dans l'espace (la quantité n'est pas la même aux poles qu'à l'équateur) et dans le temps (la quantité n'est pas la même la nuit que le jour !).
import numpy as np
P_parm2=P/(4*np.pi*R*R*1000*1000)
print ('La puissance solaire reçue est {:.2f} W/m2'.format(P_parm2))
La puissance solaire reçue est 341.24 W/m2
$$j=\frac{0,25 P_{par m2}}{\Delta H_{vap}}$$
- La masse d’eau évaporée en kg.s-1.m-2 qu'on peut aussi appeler densité de flux est :
j_evap=pourc*P_parm2/DH_eau #kg/s m2
j_evap1=j_evap*1000*3600
print ('La densité de flux d eau évaporée est {:.2e} kg s-1 m-2 soit {:.2f} g h-1 m-2'.format(j_evap,j_evap1))
La densité de flux d eau évaporée est 3.41e-05 kg s-1 m-2 soit 122.85 g h-1 m-2
- Le temps de séjour est défini comme le rapport du volume du système sur le débit qui parourt le système :
$$t_s=\frac{V}{Q_{vol}}$$
On pourra considérer (en faisant l'hypothèse d'une masse volumique de la vapeur d'eau constante -elle peut varier selon la pression et donc l'altitude-) que le rapport du volume sur le débit volumique est égal au rapport de la masse sur le débit massique.
t_s=M_vap_eau/Q_evap #s
t_s1=t_s/3600/24 #jours
print ('Le temps de sejour est {:.2e} secondes soit {:.2f} jours'.format(t_s,t_s1))
Le temps de sejour est 7.47e+05 secondes soit 8.65 jours
Ce temps correspond au temps de "vie" moyen d'une goutte de vapeur d eau entre le moment où elle s'évapore et le moment où elle retombe sous forme de précipitation. C'est un temps de vie moyen : certaines gouttes n’auront pas cette chance de voyager dans l’atmosphère autant de temps et se recondenseront juste après s’être évaporée (comme dans les bouteilles d’eau que vous avez devant vous !). D’autres voyageront beaucoup longtemps surtout si elles se font ingurgiter par un goeland ! Cette valeur est donc une moyenne et vous verrrez en master comment on peut avoir une description plus précise de la distribution des temps de séjours.