OBJECTIVES: Coupling with LMDZ

The objective of this page is to present and discuss the existing coupling scheme between ORC and LMDZ and to propose new solution to improve physical/numerical consistency in the coupling.

The main points concern i) the treatment of the Evaporation and how the "beta" factor (soil water stress) is treated as well as how the potential evaporation is calculated and ii) the treatment of the radiation in ORC (longwave and shortwave) to ensure energy conservation when coupled with LMDz.

Several documents have been produced and are usefull to understand such coupling and the major points:

1. General Coupling with GCM (Polcher et al., 1998): ​Article on Coupling
2. Updated-note on the coupling by Polcher: ​Updated note (Polcher)
3. Note by Dufresne and Ghattas on coupling with LMDZ: ​Dufresne & Ghattas
4. New note by Fuxing on recent implementations: ​Fuxing note
5. Short note by Agnes on the "beta" coefficient showing how the energy budget formulation in the current ORC-TRUNK is not consistent with the implicit scheme of Polcher et al., 1998:​Agnes note.
6. Milly, 1992, Potential evaporation and Soil Moisture in GCMs:​Milly's paper on Milly's correction
7. Derivation of Penman equation for the evaporation of a wet surface, based on the bulk aerodynamic approach and the linearization of qsat(Ts) around Ta: ​Penman equation

EVAPO-TRANSPIRATION

Evaporation during coupling (beta problem)

1. Problem:

The derivation of the surface temperature as it is done in ORCHIDEE-TRUNK is not consistent with surface energy budget equation proposed by Polcher (Polcher et al. 1998; Document 1) for implicit coupling with the vertial diffusion in the boundary layer. Document 3 by Dufresne and Ghattas gives the corresponding equations in LMDZ.

In the trunk of ORCHIDEE, the stress factor beta is missing in the denominator in the sensitivity of the latent heat flux to the surface temperature at the old time step. This is coded in enerbil as:

larsub_old = chalsu0 * vbeta1 * (un - vbeta5) * (peqBcoef -  qsol_sat) / (zikq - peqAcoef)

This is not consistent with a fully implicit surface energy budget equation (cf Documents 4 and 5) and could lead to energy conservation problem.

2. Action decided

It has been decided to re-introduce the beta and take the following formulation for the ORCHIDEE-TRUNK:

larsub_old = chalsu0 * vbeta1 * (un - vbeta5) * (peqBcoef -  qsol_sat) / (zikq - vbeta1 * (un - vbeta5) * peqAcoef)

Note that the correction does not change significantly the results of ORC-LMDZ simulations.

3. Complementary notes

As mentioned in the Updated note by Polcher (Document 2), the current formulation can be seen as a mix between an implicit scheme for net radiation, H, potential ET, and an explicit scheme for beta. This has been introduced by Jan to support the implementation of a new energy budget based on the temperature of a wet surface (the hypothetically wet temperature of Milly, 1992, cf Document 6), further noted T*. The above change can easily be reverted when the new energy budget will be implemented (which is not yet the case).

Overestimation of Potential Evaporation

The implementation of the potential evaporation formula based on a bulk formulation that considers qsat(Ts) is not satisfactory if one refers to Budyko's framework, which defines E=beta*Ep(T*), with T* defined as above, and the potential evaporation Ep = rau/ra*[qsat(T*)-qa].

Currently the code considers qsat(Ts) and not qsat(T*), and thus over-estimates the evaporation, since Ts is warmer on average than T*.

In the present version, this is tackled by correcting the bare soil evaporation afterwards following Milly's approach.

Jan proposed improved approaches based on the Penman-Monteith method to calculate the ETP (formulation that provide smaller ETP than the bulk formulation that is used). They have been developed by Anais Barella during her PhD thesis and tested in forced mode. The correct surface energy budget equation for the implicit coupling is not yet available and work has to be done before implementing this approach in Orchidee.trunk. ​http://www.lmd.polytechnique.fr/~intro/Files/2014_These_Barella.pdf).

It is stressed by all the participants that improvements on the stress factor are important and that Anais work is a first step in this direction. However, more work is needed and a clear formulation for the new approach must be well documented.

Note that other GCM groups also face the same issues and have used various approaches not to overestimate ETP.

RADIATION

The problem of energy conservation

In order to save computing time, the time step of the radiation code is longer than that of the other components of the model (rest of the atmospheric physics and call to LSM). If the radiation budget of each component is updated at a shorter time step, the energy is not conserved. This is the case in the ORCHIDEE.trunk version where the LWnet (through lwup) is calculated separately in LMDZ and ORCHIDEE.

Note that the same is true for the SWnet. To ensure energy conservation, SWnet should be constant between two calls of the RT model of LMDZ. This implies that the surface albedo does not vary between these two calls.

An option which favors the energy conservation has been introduced by Fuxing and Jean-Louis: Practically the LWup is updated only when the radiative code is called. Note that in this case, fixing the Lwup is equivalent to fix the surface temperature between two call of the RT model.

=> It is recognized that none of the option is satisfactory and that improvement is needed. Accurate simulations of the diurnal cycle of the surface temperature and fluxes requires to be able to change the surface temperature each time ORC is called and also to change the albedo (maybe less critical) However, in the mean time it is required that both option are available in Orchidee.trunk to be able to check for energy conservation problems (the option favoring energy conservation has to be committed in the trunk).

Anticipated Improvements

To go further some work is necessary on the LMDZ and ORCHIDEE sides as well as on the interface.

• LMDZ is currently implementing a new radiation code (RRTM) for the SW radiation, which will allow to increase the number of spectral bands for the albedo (expecting that the spectral bands are consistent with the one adopted for the albedo calculation in Orchidee). A post-doc student is working at LMD on the implementation of the new radiative code.

• A proposition is to follow the approach described in Polcher et al. 1998 (see above link): ORCHIDEE would provide to LMDZ an averaged (over the time-steps at which the radiation is not called) radiative temperature (averaged to the power 4), emissivity and albedo (simple averages). ORCHIDEE would thus be able to have varying albedo and temperature between calls to RT model but it will still provide to LMDz the mean

• The interface variables have to be redefined and they should include : net LW and net SW radiations, in several bands ?; emissivity and albedo (several bands ?) and the mean radiative temperature.

But first these changes have to be documented very clearly. It is stressed that these developments would suppress the backward compatibility (should we keep the backward compatibility with a flag ?)