NH3

r4164

Evaluation of ammonia emissions results calculated by the model (Fertilizer =0)

• Sites and PFT :

PFT	Sites
2	BR-Sa3 ; ID-Pag
4	IL-Yat ; ES-ES1 ; IT-SRo ; US-Blo ; US-Me4 ; US-SP2 ; US-SP3 ; US-SP4 ; DE-Bay ; DE-Tha ; FR-LBr
5	FR-Pue ; IT-Cpz ; IT-Lec ; PT-Esp ; PT-Mi1 ; AU-Tum ; AU-Wac
6	IT-Col ; IT-Non ; IT-PT1 ; IT-Ro2 ; US-MMS ; US-MOz ; DK-Sor ; FR-Fon ; UK-PL3 ; IS-Gun ; IT-Ro1 ; US-WBW ; DE-Hai ; FR-Hes ; UK-Ham ; US-Bar ; US-Ha1 ; US-Oho ; US-UMB ; US-WCr
7	CA-Man ; CA-NS2 ; CA-NS4 ; CA-Qcu ; CA-Qfo ; CA-SF1 ; CA-SF2 ; FI-Hyy ; FI-Sod ; SE-Fla
12	FR-Gri ; DK-Ris ; NL-Lut
13	NL-Lan ; US-Bo1 ; US-Ne1 ; US-Ne2 ; US-Ne3

Sites we have not looked at :

AT-Neu ; CA-NS1 ; CA-NS3 ; CA-NS5 ; CH-Oe1 ; CZ-BK1 ; DE-Gri ; DE-Meh ; DE-Wet ; DK-Lva ; ES-VDA ; FR-Lq1 ; FR-Lqr2 ; HU-Bug ; ID-Dri ; IT-Amp ; IT-Lav ; IT-Lma ; IT-Lmal ; IT-Mbor ; IT-Ren ; NL-Ca1 ; NL-Haar ; NL-Hor ; NL-Loo ; PT-Mi2 ; RU-Fyo ; RU-Ha1 ; RU-Ha2 ; RU-Ha3 ; RU-Zot ; SE-Nor ; SE-SK1 ; SE-Sk2 ; SK-Tat ; UK-EBu ; UK-Gri ; UK-Tad ; US-Aud ; US-Bkg ; US-FPe ; US-Goo ; US-Ho1 ; US-Ho2 ; US-SP1 ; US-Var ; US-Wci0 ; US-Wi1 ; US-Wi2 ; US-Wi4 ; US-Wi5 ; US-Wi8 ; US-Wi9

• Parameters studied

Temperature, Humidity and pH
path : /home/users/taudoux/documents/Results/Parameters/

Nitrogen Pools (Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O)
path : /home/users/taudoux/documents/Results/Pools/

Emissions (NH3) and Depositions (NHX)
path : /home/users/taudoux/documents/Results/Emission-Deposition/

• Files

. The script "script_calcul_parameters.ksh" used in order to create files of Temperature, Humidity and pH of each site in function of time listed in :
/home/users/taudoux/documents/Results/Parameters/T2M/SITES/;
/home/users/taudoux/documents/Results/Parameters/LITTERHUM/SITES/;
/home/users/taudoux/documents/Results/Parameters/pH/SITES/.

. The script "script_calcul_pools.ksh" used in order to create files of soil_NH4+, Soil_NO3-, Soil_NOX and Soil_N2O for each site in function of time listed in :
/home/users/taudoux/documents/Results/Pools/SOIL_NH4/SITES/ ;
/home/users/taudoux/documents/Results/Pools/SOIL_NO3/SITES/ ;
/home/users/taudoux/documents/Results/Pools/SOIL_NOX/SITES/ ;
/home/users/taudoux/documents/Results/Pools/SOIL_N2O/SITES/ .

. The Script "script_calcul_emission-deposition.ksh" used in order to create files of NH3 emission and NHX deposition for each site in function of time in :
/home/users/taudoux/documents/Results/Emission-Deposition/NH3_EMISSION/Sites/ ;
/home/users/taudoux/documents/Results/Emission-Deposition/NHX_DEPOSITION/Sites/.

For each parameter are listed :

. Figure of correlation with NH3 emissions (except pH) obtained using "script_calcul_r2.ksh" to obtain the R2 and then bar_plot.py in order to create the figure ;
. Average emission of ammonia on each site in function of the parameter studied using "scatterplot.py" ;
. Text documents in which are listed the different values.

• Comparison with data used by LMDz-INCA model

path : /home/users/taudoux/documents/Results/Comparison/

Figure, created using "bar_plot_donnees.py", in which for each site are listed the average emission of NH3 for both model.

Average emission of ammonia (g(NH3-N).m⁻².day⁻¹) on each site in function of the parameter considered
Temperature (K)	Humidity (-)	pH (-)	Deposition g(NHx).m-2
SOIL NH4 g(NH4).m-2	SOIL N2O g(N2O).m-2	SOIL NOX g(NOx).m-2	SOIL NO3 g(NO3).m-2

Correlation with NH3 emissions (g(NH3-N).m⁻².day⁻¹)
Temperature (K)	Humidity (-)	Deposition g(NHx).m-2
SOIL NH4 g(NH4).m-2	SOIL N2O g(N2O).m-2	SOIL NOX g(NOx).m-2	SOIL NO3 g(NO3).m-2

Comparison with LMDz-INCA model's data

---

Evaluation of ammonia emissions results calculated by the model (Fertilizer added)

Code was changed : separation of nitrogen input. Initially fertilizer, bnf and atmospheric deposition were in the same loop.

Multiple simulations were run with different amount of fertilizer (0 , 0.015, 0.0225, 0.03, 0.0375, 0.045 , 0.0525, 0.06 and 0.1 g.m-2.day-1). Each simulation were run with two different ratio of NH4 in the fertilizer : 0.5 and 0.8.

We obtain same results when nfert=0 or when nfert is not precised.

• Sites we looked at and their PFT :

PFT	Sites
12	FR-Gri ; DK-Ris ; NL-Lut
13	NL-Lan ; US-Bo1 ; US-Ne1 ; US-Ne2 ; US-Ne3

• Parameters studied

General path : path : /home/users/taudoux/documents/Results/B/Simul_ratio0.5_v1/ or /home/users/taudoux/documents/Results/B/Simul_ratio0.8_v1/
Temperature, Humidity in the folder Parameters/

Nitrogen Pools (Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O) and Leaching (NH4 and NO3) in the folder Pools_fluxes/

Emissions (NH3, N2, N2O, NOx) in the folder NH3_EM/

GPP in the folder GPP/

Nitrogen input summed (fertiliser, NHx and NOx depositions) in the folder DEP_FERT/

We calculated emission factor (EF) of different fluxes of nitrogen (NH3, N2, N2O and NOx emission, but also NH4 and NO3 leaching) with different methods :
. The first (A.) consisted to normalized the fluxes minus emissions when nfert=0 by the amount of fertiliser added ;
. The second (B.) consisted to normalized the fluxes by the amount of fertiliser added ;
. The third (C.) consisted to normalized the fluxes by all nitrogen inputs.

Script python used : barplot.py using files.txt in which values of nitrogen emissions and leaching are listed.

1. EF of each parameters on each site for NH4%=0.5 and 0.8

DK-Ris
0.5	0.8

FR-Gri
0.5	0.8

NL-Lan
0.5	0.8

NL-Lut
0.5	0.8

US-Bo1
0.5	0.8

US-Ne1
0.5	0.8

US-Ne2
0.5	0.8

US-Ne3
0.5	0.8

The sum of those EF is not equal to 1 for each site.