Changeset 9995


Ignore:
Timestamp:
2018-07-24T15:08:22+02:00 (2 years ago)
Author:
vancop
Message:

Radiative transfer SI3 documentation

Location:
NEMO/trunk/doc/si3_doc
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/doc/si3_doc/tex_main/SI3_manual.tex

    r9983 r9995  
    8484 
    8585\subfile{../tex_sub/chap_dynamics} 
    86 % 
     86 
    8787\subfile{../tex_sub/chap_transport} 
    88 % 
     88 
    8989\subfile{../tex_sub/chap_ridging_rafting} 
    90 % 
     90 
    9191\subfile{../tex_sub/chap_radiative_transfer} 
    92 % 
     92 
    9393\subfile{../tex_sub/chap_thermo} 
    94 % 
     94 
    9595\subfile{../tex_sub/chap_interfaces} 
    96 % 
     96 
    9797\subfile{../tex_sub/chap_output_diagnostics} 
    98 % 
     98 
    9999\subfile{../tex_sub/chap_single_category_use} 
    100 % 
     100 
    101101\subfile{../tex_sub/chap_bdy_agrif} 
    102 % 
     102 
    103103\subfile{../tex_sub/chap_miscellaneous} 
    104104 
  • NEMO/trunk/doc/si3_doc/tex_sub/chap_radiative_transfer.tex

    r9974 r9995  
    33 
    44\begin{document} 
     5 
     6 
     7 
    58 
    69% ================================================================ 
     
    1518$\ $\newline    % force a new line 
    1619 
    17 \section{Surface albedo} 
    18  
    19 \section{Transmission} 
    20  
    21 \section{Attenuation} 
     20Radiative transfer in SI$^3$ currently reduces to the parameterization of solar radiation partitionning through the snow/ice/open water system, treated using a single wavelength band. This will likely be improved in future versions of the code. In this chapter, we first explain how solar radiation is partionned in the snow-ice system, then describe how, solar radiation-wise, the snow-ice system is framed in the context of the atmosphere-ice-ocean boundary. 
     21 
     22\section{Solar radiation partitionning in the snow-ice system} 
     23 
     24%-------------------------------------------------------------------------------------------------------------------- 
     25% 
     26% FIG x : Radiation balance % 
     27\begin{figure}[!ht] 
     28\begin{center} 
     29\vspace{0cm} 
     30\includegraphics[height=6cm,angle=-00]{../Figures/radiative_transfer.png} 
     31\caption{Partitionning of solar radiation in the snow-ice system, as represented in SI$^3$.} 
     32\label{fig_radiative_transfer} 
     33\end{center} 
     34\end{figure} 
     35% 
     36%-------------------------------------------------------------------------------------------------------------------- 
     37 
     38Solar radiation in the snow-ice system is represented following the principles of \cite{MaykutUntersteiner71}, see Fig.\ref{fig_radiative_transfer}, using a unique band of solar radiation. Incident solar radiation (W/m$^2$, counted per unit ice area - not per grid cell area) is specified in the SBC routines and is a priori category dependent, because multiple atmosphere-surface reflexions are frequent in polar regions imply that incident radiation depends on the surface albedo and therefore surface state.  
     39 
     40Net solar radiation qsr\_ice(i,j,l) is obtained by substracting the reflected part of the incident radiation using the surface albedo $\alpha(i,j,l)$, parameterized as a function of environmental conditions.  
     41 
     42The subsequent attenuation of solar radiation through the snow-ice system is represented assuming the presence of a highly diffusive surface scattering layer, absorbing a fraction $i_o$ of net solar radiation, which is transformed into sensible heat, contributing to the surface energy balance.  
     43 
     44The remainder of solar radiation, qtr\_ice\_top(i,j,l), is transmitted below the surface and attenuates following Beer-Lambert law. The part of solar radiation that is absorbed on its path to the base of the ice is given as sensible heat to the snow/ice system, via a source term in the heat diffusion equation. The rest of solar radiation that reaches the ice base, qtr\_ice\_bot(i,j,l), is transmitted to the ocean.  
     45 
     46In the rest of this section, we describe how the albedo, the surface transmission parameter $i_o$ and the attenuation of solar radiation are parameterized. 
     47 
     48\subsection{Surface albedo} 
     49 
     50The surface albedo determines the amount of solar radiation that is reflected by the ice surface, hence also net solar radiation. The philosophy of the parameterization of surface albedo is the following: each ice category has its own albedo value $\alpha(i,j,l)$, determined as a function of cloud fraction, ice thickness, snow depth, melt pond fraction and depth, using observation-based empirical fits.  
     51 
     52The original \cite{ShineHenderson85} parameterization had a few inconsistencies and flaws that the revisited parameterization described hereafter fixes. In particular, the dependencies of albedos on ice thickness, snow depth and cloud fraction have been revised in the light of recent observational constraints \citep{Brandtetal05,GrenfellPerovich04}. In addition, the asymptotic properties of albedo are better specified and now fully consistent with oceanic values. Finally, the effect of melt ponds has been included \citep{Lecomteetal15}. 
     53 
     54The user has control on 5 reference namelist values, which describe the asymptotic values of albedo of snow and ice for dry and wet conditions, as well as the deep ponded-ice albedo. Observational surveys, in particular during SHEBA in the Arctic \citep{Perovichetal02alb} and further additional experiments \citep{GrenfellPerovich04}, as well as by \cite{Brandtetal05} in the Antarctic, have provided relatively strong constraints on the surface albedo. In this context, the albedo can hardly be used as the main model tuning parameter, at least outside of these observation-based bounds (see namalb for reference values). 
     55 
     56\forfile{../namelists/namalb} 
     57 
     58%-------------------------------------------------------------------------------------------------------------------- 
     59% 
     60% FIG x : cloud correction 
     61% 
     62\begin{figure}[ht] 
     63\begin{center} 
     64\vspace{0cm} 
     65\includegraphics[height=10cm,angle=-00]{../Figures/albedo_cloud_correction.png} 
     66\caption{Albedo correction $\Delta \alpha$ as a function of overcast sky (diffuse light) albedo $\alpha_os$, from field observations \cite[][their Table 3]{GrenfellPerovich04} (squares) and 2nd-order fit (Eq. \ref{eq_albedo_cloud_correction}). Red squares represent the irrelevant data points excluded from the fit. For indication, the amplitude of the correction used in the ocean component is also depicted (blue circle).} 
     67% ocean uses 0.06 for overcast sky (Payne 74) and Briegleb and Ramanathan parameterization 
     68\label{fig_albedo_cloud_correction} 
     69\end{center} 
     70\end{figure} 
     71% 
     72%-------------------------------------------------------------------------------------------------------------------- 
     73 
     74Because the albedo is not an intrinsic optical property, it depends on the type of light (diffuse of direct), which is practically handled by weighting the clear (cs) and overcast (os) skies values by cloud fraction $c(i,j)$ \citep{FichefetMaqueda97}: 
     75\begin{equation} 
     76\alpha(i,j,l) = [ 1 - c(i,j) ] \cdot \alpha_{cs} (i,j,l) + c (i,j) \cdot \alpha_{os}(i,j,l). 
     77\end{equation} 
     78For concision, we drop the spatial and category indices hereafter. \cite{GrenfellPerovich04} observations at Point Barrow, on the Alaskan Coast, suggest that clear and overcast sky albedos are directly related through 
     79\begin{equation} 
     80\alpha_{cs} = \alpha_{os} - \Delta \alpha(\alpha_{os}). 
     81\end{equation} 
     82The relation between $\Delta \alpha$ and $ \alpha_{os}$ can well be handled using a 2$^{nd}$-order polynomial fit (Fig. \ref{fig_albedo_cloud_correction}): 
     83\begin{equation} 
     84\Delta \alpha =  ( -0.1010 \cdot \alpha_{os}^2 + 0.1933 \cdot \alpha_{os} - 0.0148 ). 
     85\label{eq_albedo_cloud_correction} 
     86\end{equation} 
     87Overcast sky surface albedo is used as a reference, from which the clear-sky value is derived. 
     88 
     89%-------------------------------------------------------------------------------------------------------------------- 
     90% 
     91% FIG x : thickness, snow depth, pond depth dependencies 
     92% 
     93\begin{figure}[ht] 
     94\begin{center} 
     95\vspace{0cm} 
     96\includegraphics[height=4cm,angle=-00]{../Figures/albedo_dependencies.png} 
     97\caption{Example albedo dependencies on ice thickness, snow depth and pond depth, as parameterized in SI$^3$.} 
     98\label{fig_albedo_dependencies} 
     99\end{center} 
     100\end{figure} 
     101% 
     102%-------------------------------------------------------------------------------------------------------------------- 
     103 
     104The second important parameter that controls surface albedo is surface type. In each category, we assume that three types of surfaces can coexist (bare, snow-covered and ponded ice), with respective fractions $f_{ice}$, $f_{snw}$ and $f_{pnd}$ summing to 1. Then the overcast albedo is expressed as 
     105\begin{equation} 
     106\alpha_{os}(i,j,l) = f_{ice} \cdot \alpha_{ice} + f_{snw} \cdot \alpha_{snw} + f_{pnd}\cdot \alpha_{pnd} 
     107\end{equation} 
     108with a specific albedo value for each surface type.  
     109 
     110The surface fractions $f_{ice}$, $f_{snw}$ and $f_{pnd}$ are currently crudely parameterized: if snow is present ($h_s>0$), then $f_{snw}=1$ and $f_{ice}=f_{pnd}=0$. In the absence of snow, $f_{pnd}$ is either specified or calculated (depending on melt pond options in nampnd), and $f_{ice}=1.-f_{pnd}$. Admittedly, more refined parameterizations of $f_{snw}$ could improve the realism of the model. Note finally that the dependence of surface albedo on the presence of melt ponds can be included or not (namelist parameter ln\_pnd\_alb). If the latter is set to false, $f_{pnd}$ is always assumed zero in the albedo computations. 
     111 
     112Works by \cite{Brandtetal05} and references therein, indicate that the dependence of the albedo of bare ice on ice thickness depends is linear/logarithmic/constant from thin to thick ice. Hence, the following expressions capture  the essence of their works: 
     113\begin{eqnarray} 
     114\alpha_{ice} =  
     115\begin{cases} 
     116\alpha_{ice}^{\infty} & \text{ if } h_i > 1.5, \\ 
     117\alpha_{ice}^{\infty} + ( 0.18 - \alpha_{ice}^{\infty} ) \cdot \frac{ln(1.5) - ln(h_i)}{ln(1.5) - ln(0.05)} & \text{ if } 0.05 < h_i ,<= 1.5 \\ 
     118\alpha_{oce}  + ( 0.18 - \alpha_{oce} ) h_i /0.05 & \text{ if } h_i < 0.05. 
     119\end{cases} 
     120\end{eqnarray} 
     121The thick-ice constant albedo value depends on whether the surface is dry or melting: 
     122\begin{eqnarray} 
     123\alpha_{ice}^{\infty} = 
     124\begin{cases} 
     125\alpha_{i,dry} & \text{ if } T_{su} < T_{fr} \\ 
     126\alpha_{i,mlt} & \text{ if } T_{su} = T_{fr}, 
     127\end{cases} 
     128\end{eqnarray} 
     129values that are to be specified from the namelist. 
     130 
     131\cite{GrenfellPerovich04} suggest that the dependence of surface albedo on snow depth is exponential, 
     132\begin{eqnarray} 
     133\alpha_{snw} = \alpha_{snw}^{\infty} - ( \alpha_{snw}^{\infty} - \alpha_{ice} ) * exp( -h_s / h_s^{ref} ), 
     134\end{eqnarray} 
     135where $h_s^{ref} = 0.02$ $(0.03)$ m for dry (wet) snow. As for bare ice, the deep-snow asymptotic albedo also depends on whether the surface is dry or melting: 
     136\begin{eqnarray} 
     137\alpha_{snw}^{\infty} = 
     138\begin{cases} 
     139\alpha_{s,dry} & \text{ if } T_{su} < T_{fr} \\ 
     140\alpha_{s,mlt} & \text{ if } T_{su} = T_{fr}, 
     141\end{cases} 
     142\end{eqnarray} 
     143values that are to be specified from the namelist. 
     144 
     145Based on ideas developed from melt ponds on continental ice \citep{ZuoOerlemans96}, the albedo of ponded ice was proposed to follow \citep{Lecomteetal11}: 
     146\begin{eqnarray} 
     147\alpha_{pnd} = \alpha_{dpnd} - ( \alpha_{dpnd} -\alpha_{ice} ) \cdot exp( -h_{pnd} / 0.05 ) 
     148\end{eqnarray} 
     149$\alpha_{dpnd}$ is a namelist parameter. \cite{EbertCurry93} also use such dependency for their multi-spectral albedo. 
     150 
     151The dependencies of surface albedo on ice thickness, snow depth and pond depth are illustrated in Fig. \ref{fig_albedo_dependencies}. 
     152 
     153\subsection{Transmission below the snow/ice surface} 
     154 
     155The transmitted solar radiation below the surface is represented following \cite{FichefetMaqueda97} and \cite{MaykutUntersteiner71}: 
     156\begin{eqnarray} 
     157qtr\_ice\_top(i,j,l) = i_o(i,j) qsr\_ice(i,j,l), 
     158\end{eqnarray} 
     159where $i_o=0$ in presence of snow, and depends on cloud fraction otherwise, based on works of \cite{GrenfellMaykut77}. This parameterization needs to be re-evaluated and likely updated. 
     160 
     161\subsection{Attenuation and transmission below the ice/ocean interface} 
     162 
     163Attenuation of solar radiation through the ice follows Beer-Lambert law. In practise, we assume that irradiance below layer $k$ is given by 
     164 
     165\begin{eqnarray} 
     166radtr\_i(i,j,k,l) = qtr\_ice\_top(i,j,l) \cdot exp(-\kappa_i z), 
     167\end{eqnarray} 
     168where $\kappa_i = 1$ m$^{-1}$ is the exponential attenuation coefficient (namelist parameter rn\_kappa\_i). Hence, at the ice base, remains below the $l^{th}$ category a transmitted flux:  
     169\begin{eqnarray} 
     170qtr\_ice\_bot(i,j,l) = qtr\_ice\_top(i,j,l) \cdot exp(-\kappa_i h_i). 
     171\end{eqnarray} 
     172 
     173\section{Solar radiation: framing sea ice at the ocean-atmosphere boundary} 
     174 
     175How solar radiation transfer through sea ice is framed into the atmosphere-ice-ocean is nearly identical but not exactly the same in forced and coupled mode (see Fig. \ref{fig_radiative_transfer}. 
     176 
     177The basic principle of the computation is that the irradiant flux given to the ocean model (qsr) is computed as the average flux per grid cell area (qsr\_tot) minus what is given to the sea ice ($\sum a(l) qsr\_ice(l)$), plus what is transmitted below sea ice $\sum$ qtr\_ice\_bot(jl) (see at the base of Fig. \ref{fig_radiative_transfer}). Such formulation ensures heat conservation by construction. 
     178 
     179%-------------------------------------------------------------------------------------------------------------------- 
     180% 
     181% FIG x : Radiation balance % 
     182\begin{figure}[!ht] 
     183\begin{center} 
     184\vspace{0cm} 
     185\includegraphics[height=8cm,angle=-00]{../Figures/radiation_atm_ice_oce.png} 
     186\caption{Framing solar radiation transfer through sea ice into the atmosphere-ice-ocean context.} 
     187% ocean uses 0.06 for overcast sky (Payne 74) and Briegleb and Ramanathan parameterization 
     188\label{fig_radiative_transfer} 
     189\end{center} 
     190\end{figure} 
     191% 
     192%-------------------------------------------------------------------------------------------------------------------- 
     193 
     194\subsection{Forced mode} 
     195 
     196In forced-atmosphere mode, it is the incoming solar irradiance fluxes above the ocean and sea ice (categories) that are specified (from files) or computed (from bulk formulae), and constitute the basis of solar radiation transfer computations. Then the net solar fluxes above open water (qsr\_oce) and ice categories (qsr\_ice) are obtained by multiplication by $1-\alpha$. qsr\_tot is then diagnosed as a weighted sum of qsr\_oce and the qsr\_ice(jl)'s. 
     197 
     198\subsection{Coupled mode} 
     199 
     200In coupled-atmosphere mode, qsr\_tot and qsr\_ice have to be provided by the atmospheric model, whereas qsr\_oce is diagnosed from qsr\_ice and qsr\_tot. 
     201 
     202Some atmospheric models enable \textit{tiling} and can provide solar fluxes over individual ice categories. For such atmospheric models, net solar radiation fluxes are directly useable by SI$^3$ (nn\_flxdist =  -1). Other models cannot do tiling, being only able to provide a net solar flux above all ice categories, seen as a single surface type. For such models a first option is to give the net solar flux above sea ice identically to all sea ice categories (nn\_flxdist = 0). Yet a better option is to redistribute the mean solar flux above sea ice $< qsr_ice >$ above categories (nn\_flxdist = 2) using the following scaling, conserving heat by construction: 
     203 
     204\begin{eqnarray} 
     205qsr\_ice(jl) = < qsr\_ice > \frac{1 - \alpha(jl)}{1 - < \alpha >} 
     206\end{eqnarray} 
     207 
     208where $ < \alpha >$ is the albedo averaged over the ice categories. Note that for testing, the flux redistributor can be emulated in forced mode (nn\_flxdist = 1). 
    22209 
    23210\end{document} 
  • NEMO/trunk/doc/si3_doc/tex_sub/todolist.tex

    r9983 r9995  
    3232\item 5. Ridging and rafting: V1 from previous doc. Should be refreshed. A bit too long. All namelist parameters are not explained (namdyn\_rdgrft). 
    3333 
    34 \item 6. Radiative transfer: Full rewriting required (Martin \& Clem, 3-6h). Available elements from previous documentation, but mostly obsolete. 
     34\item \textcolor{gray}{6. Radiative transfer: Full rewriting required (Martin \& Clem, 3-6h). Available elements from previous documentation, but mostly obsolete.} 
    3535 
    3636\item 7. Thermodynamics: Lots of rewriting required. Available elements from previous documentation, but mostly obsolete (Martin, 6h) 
Note: See TracChangeset for help on using the changeset viewer.