Explicit representation of impurities

In order to use the model configuration including light-absorbing impurities in snow presented in Tuzet et al. (2017) you have to install the branch cen_dev from the SURFEX_Git2 repository.

To activate the explicit representation of impurities as well as TARTES radiative transfer scheme, you only have to make some changes in the default namelist.
These changes are detailed hereafter for each namelist section and a default namelist activating the explicit representation of impurities is attached to the wiki page.
It is important to note that the s2m command to launch is the same as a classic Crocus run, only the namelist must be modified.

Necessary changes in namelist


NIMPUR=0/1/2              ; default=0

NIMPUR is the number of impurity you want to use in your simulation.
For now you can set NIMPUR=1 if you want to run simulations with BC only and NIMPUR=2 if you want to run, simulations with BC and Dust.
If you want to run a simulation with dust only you can set NIMPUR=2 and prescribe no BC deposition.


LSPECSNOW = .TRUE.        ; default=.FALSE. ,   #Enable spectral computation inside SURFEX/Crocus, necessary to run TARTES radiative transfer scheme    
NIMPUROF=0/1/2            ; default=0       #Initialize the number of impurities in the OFFLINE module, you have to set NIMPUROF to the same value as NIMPUR. 


LPROBANDS= .TRUE.         ; default=.FALSE.  #Enable the spectral resolution of Crocus diagnostics, necessary if you want to get spectral albedo and spectral direct/diffuse ratio diagnostics


CSNOWRAD='T17'            ; default=B92   #Set the radiative transfer scheme to TARTES+impurity as in Tuzet et al. (2017)     
#The following option is currently unavailable, setting it to .TRUE. will cause some errors 
#LATMORAD=.FALSE.          ; default=.FALSE. # This option is not stable yet, but it is supposed to compute the direct/diffuse ratio directly from atmospheric informations (AOD, Ozone column, Water column..)


CSELECT: You have to add 'SNOWIMP1','SNOWIMP2','SPEC_ALB', 'DIFF_RATIO', SPEC_TOT in this field.

CSELECT makes a selection of the diagnostics you want to compute in your output NetCDF file (PRO file).
SNOWIMP1 is the BC concentration in each snow layer (g/g_of_snow)
SNOWIMP2 is the Dust concentration in each snow layer (g/gg_of_snow)
SPEC_ALB is the spectral albedo for the 186 Crocus spectral bands (from 300 to 4000nm included by step of 20nm; i.e [300,320...,3980,4000])
DIFF_RATIO is the spectral direct to diffuse ratio for the 186 Crocus spectral bands (from 300 to 4000nm included by step of 20nm; i.e [300,320...,3980,4000])
SPEC_TOT is the total incoming irradiance after spectral repartition i.e the spectral irradiance used by TARTES for radiative transfer computations

Optional changes (to activate impurity deposition)

The following changes set the way you want to prescribe your deposition fluxes in input of SURFEX/Crocus.
You can either feed the model with prescribed and constant deposition fluxes or introduce a wet and dry deposition field directly in the forcing file.

  • If you want to prescribe constant deposition fluxes over time:


XIMPUR_WET(1)=5.e-11     ; default=0.
XIMPUR_WET(2)=5.e-9      ; default=0.
XIMPUR_DRY(1)=1.e-11     ; default=0.
XIMPUR_DRY(2)=1.e-9     ; default=0.

These variable set the different depositions fluxes and are expressed in g m^(-2) s^(-1).

XIMPUR_WET corresponds to the initial amount of impurity present in the falling snow (wet deposition), activated in case of rain or snow.
XIMPUR_DRY corresponds to the dry deposition coefficient always activated.
You can precise the type of impurity into parenthesis as done above. e.g XIMPUR_WET(1)=1.e-9 set the wet depostion coefficient of BC to 1.e-9 g m^(-2) s^(-1)

  • If you want to prescribe directly the deposition fluxes at each model time step from a forcing file:


LFORCIMP = .TRUE.         ; default=.FALSE.
#The following option is currently unavailable, setting it to .TRUE. will cause some errors
#LFORCATMOTARTES =.FALSE.  ; default=.FALSE. # Activate the prescription of aerosol optical depth and ozone column from forcing file if LATMORAD=.TRUE.

When you activate LFORCIMP by setting it to .TRUE., you have to add new variables in your forcing file.
You have two variables to add for each type of impurity: IMPWET (wet deposition coefficient) and IMPDRY (dry deposition coefficient). Both these deposition fluxes have to be in g m^(-2) s^(-1).
If you have one type of impurity you will need IMPWET1 and IMPDRY1.
If you have two type of impurity you will need IMPWET1 and IMPDRY1 (BC) and IMPWET2 and IMPDRY2 (Dust). etc...
These new variables must be defined with the same dimension as the snowfall rate for exemple.
If those variables are not defined in your forcing files you will get an error when running the simulation.
Note that if you activate LFORCIMP the deposition values contained in NAM_SURF_SNOW_CSTS are ignored.

##Nominal behaviour of LFORCATMOTARTES when it will be repared##
When you activate LFORCATMOTARTES by setting it to .TRUE., you have to add new variables in your forcing file : AODTOT (total aerosol optical depth) and OZONE (total ozone column).
These new variables must be defined with the same dimension as the snowfall rate for exemple. If those variables are not defined in your forcing files you will get an error when running the simulation. ###

P.S: In the attached namelist the other parameters have default values, you might need to tune this parameters for your study or to use the multiphysics ensemble of Crocus, ESCROC described in Lafaysse et al. (2017).

If you need to add another type of impurity in the model please open a new ticket on the support page and someone will get back to you.

Initialization of an user prescribed snowpack with impurities

In the case you want to initialize an user-prescribed snowpack by defining manually it's layers, it's possible to initialize the concentration of impurities in each layer by adding the values in NAM_PREP_ISBA_SNOW by adapting the following exemple:

&NAM_PREP_ISBA_SNOW XIMPURSNOW= 0,0,5e-9,0,0 #BC concentration (g/g_of_snow) in each layer (5 layers here)
XIMPURSNOW= 0,0,5000e-9,0,0 #Dust concentration (g/g_of_snow) in each layer(5 layers here)

And so on if you have more impurity types. By default the concentration in all layers will be 0.

Information on default values used in TARTES radiative transfer scheme

The refractive of ice is taken from Warren and Brandt. 2008.

The values of the shape parameters B and g used in Kokhanovsky and Zege(2004) theory are set to 1.6 and 0.845 respectively (according to Libois et al. 2013 and Dumont et al. 2017)

The Mass Absoprtion Coefficient (MAC) of black carbon is based on refractive index advised by Bond and Bergstrom 2006, with a value of 11.25 kg m-1 at 550nm as in Tuzet et al. 2019.
The Mass Absoprtion Coefficient (MAC) of mineral dust is based on refractive index found in Caponi et al. 2017 for Lybian dust of PM 2.5

Important information about new diagnostics

By activating these new developments, you will have new diagnostics in your PRO file (if you correctly modified LSELCT in &NAM_WRITE_DIAG_SURFn).
The impurity concentration diagnostics (SNOWIMP) have classic dimensions (time,number of point,number of snow layers).
The spectral diagnostics (SPEC_ALB,DIFF_TOT and SPEC_TOT) have special dimensions (time,number of spectral bands,number of snow layers).
For now the spectral bands run from 300 to 4000 (included) by step of 20, meaning 186 spectral bands. (300,320,340,......,3980,4000)


Bond, T. C. and Bergstrom, R. W.: Light absorption by carbonaceous particles: An investigative review, Aerosol science and technology, 40,
27–67, 2006

Caponi, L., Formenti, P., Massabo, D., Biagio, C. D., Cazaunau, M., Pangui, E., Chevaillier, S., Landrot, G., Andreae, M. O., Kandler, K.,
et al.: Spectral-and size-resolved mass absorption efficiency of mineral dust aerosols in the shortwave spectrum: a simulation chamber
study, Atmospheric Chemistry and Physics, 17, 7175–7191, 2017.

Dumont, M., Arnaud, L., Picard, G., Libois, Q., Lejeune, Y., Nabat, P., Voisin, D., and Morin, S.: In situ continuous visible and near-infrared
spectroscopy of an alpine snowpack, The Cryosphere, 11, 1091–1110,, http://www.the-cryosphere.
net/11/1091/2017/, 2017.

Kokhanovsky, A. and Zege, E.: Scattering optics of snow, Applied Optics, 43(7), 1589–1602,,

Libois, Q., Picard, G., France, J. L., Arnaud, L., Dumont, D., Carmagnola, C. M., and King, M. D.: Influence of grain shape on light penetration
in snow, The Cryosphere, 7, 1803–1818,, 2013.

Tuzet, F., Dumont, M., Lafaysse, M., Picard, G., Arnaud, L., Voisin, D., Lejeune, Y., Charrois, L., Nabat, P., and Morin, S.: A multilayer
physically based snowpack model simulating direct and indirect radiative impacts of light-absorbing impurities in snow, The Cryosphere,
11, 2633–2653, 2017

Tuzet, F. et al. (In discussion) :Influence of light absorbing particles on snow spectral irradiance profiles, TCD 2019

Warren S, Brandt R. : Optical constants of ice from the ultraviolet to the microwave: A revised compilation. Journal of Geophysical Research: Atmospheres. 2008

Tartes_Imp.nam - Running Crocus with TARTES and explicit forcing of dust and black carbon deposition (12.7 KB) francois tuzet, 06/11/2019 05:14 PM

Tartes_NoImp.nam - Default namelist including impurity scheme (12.7 KB) francois tuzet, 06/11/2019 05:15 PM