Changes between Version 8 and Version 9 of user/clem/Heat_Flux_air-ice-ocean_interface
- Timestamp:
- 2016-05-19T13:10:42+02:00 (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
user/clem/Heat_Flux_air-ice-ocean_interface
v8 v9 4 4 It is time to approach the matter from a different angle. I describe below what we are doing in NEMO and what we think the solution should be, and you can comment on that.[[BR]] 5 5 6 So far, the total non-solar flux (Qns) is imposed by the atmosphere but it excludes (from what I understand) the heat contents of snow precipitation (Qsnow) and rain (Qrain). It also does not take into account the heat content (if any) of evaporated ice and water (Qevap, already in vapor form) that the atmosphere should receive. The total non-solar (conservative) flux should thus be written :[[BR]]6 So far, the total non-solar flux (Qns) is imposed by the atmosphere but it excludes (from what I understand) the heat contents of snow precipitation (Qsnow) and rain (Qrain). It also does not take into account the heat content (if any) of evaporated ice and water (Qevap, already in vapor form) that the atmosphere should receive. The total non-solar (conservative) flux should thus be written (including the snow melting when it reaches the ocean):[[BR]] 7 7 8 {{{Qns_tot = Qns + Qsnow + Qrain + Qevap }}}8 {{{Qns_tot = Qns + Qsnow + Qrain + Qevap + "Snow-melt-in-ocean"}}} 9 9 10 10 … … 13 13 >>''C Rousset'': I agree but I do not think we have access to ice surface temperature at this stage (not sure though)? 14 14 15 {{{Qsnow = sprecip * ice_frac * ( Cp * SST - Lfus) [in W/m2]}}}15 {{{Qsnow = sprecip * ( Cp * SST ) [in W/m2]}}} 16 16 17 17 {{{Qrain = lprecip * Cp * SST}}} … … 20 20 21 21 {{{Qevap (ice) = 0}}} 22 23 {{{"Snow-melt-in-ocean" = sprecip * ice_frac * ( - Lfus )}}} 22 24 23 25 with sprecip = snow precip, lprecip = liquid precip, ice_frac = ice fraction, Cp = specific heat of sea water, Lfus = Latent heat of fusion of pure ice at 0deg, SST = sea surface temperature.[[BR]] … … 35 37 1) The atmosphere considers all these fluxes to be 0, that’s ok but 36 38 NEMO needs to know 37 and we will have to trick the calculation of snow temperature in LIM3 since it is recalculated from Qsnow (see nota bene). You can see that as NEMO business but it is important you know it has consequences on sea ice thermodynamics.[[BR]]39 and we will have to trick the calculation of snow temperature in LIM3 since it is recalculated from Qsnow and Lfus (see nota bene). You can see that as NEMO business but it is important you know it has consequences on sea ice thermodynamics.[[BR]] 38 40 >''O Marti'': Not exactly. The flux from atmosphere is not 0. But atmosphere doesn't need to compute it. You could diagnose the actual flux by computing the total ATM enthalpy before and after the rains leaves ATM. 39 41 40 2) The atmosphere (you, from your atmospheric models) provides non-zero fluxes and then it is better if NEMO knows how Qsnow is calculated in the atmosphere, which Cp, temperature and Lfus is used (see why below).42 2) The atmosphere (you, from your atmospheric models) provides non-zero fluxes and then it is better if NEMO knows how Qsnow is calculated in the atmosphere, which Cp, temperature and also Lfus is used (see why below). 41 43 42 44 43 45 **Nota bene**[[BR]] 44 The most important flux here (by far) is Qsnow. When snow falls on ice, its heat content adds up with the pre-existing snow to give a total heat contentequal to q and snow temperature is recalculated with this expression:[[BR]]46 The most important flux here (by far) is Qsnow. When snow falls on ice, its enthalpy adds up with the pre-existing snow to give a total enthalpy equal to q and snow temperature is recalculated with this expression:[[BR]] 45 47 {{{q = rho(snow) * ( Cp(snow) * T(snow) - Lfus ) [in J/m3]}}}[[BR]] 46 Therefore a large difference between the expressions of Qsnow and q would lead to erroneous snow temperatures. One consequence could be a snow that always melts because its temperature is always positive. The largest difference between Qsnow and qis certainly Lfus which may differ greatly between the atmospheric model and LIM3. [[BR]]48 Therefore a large difference between the expressions of (Qsnow - sprecip*Lfus) and q would lead to erroneous snow temperatures. One consequence could be a snow that always melts because its temperature is always positive. The largest difference is certainly Lfus which may differ greatly between the atmospheric model and LIM3. [[BR]] 47 49 48 At this point I expect to receive from you the Qsnow, Qevap and Qrain (whether you think it should be zero or a value defined by your atmospheric model calculations). If you do not think this is needed please send me comments based on the above equations so that we can move on in a constructive way.50 At this point I expect to receive from you the Qsnow, Qevap(ocean & ice) and Qrain (whether you think it should be zero or a value defined by your atmospheric model calculations). If you do not think this is needed please send me comments based on the above equations so that we can move on in a constructive way. 49 51 50 52 == LIM3 physical constants == 51 53 Heat capacity (ice/snow) = 2067 52 54 53 Heat capacity (sea water) = 3991.867957 55 Heat capacity (sea water) = 3991.86795711963 54 56 55 57 Latent heat of fusion at 0deg (Lfus) = 334,000 J/kg