Changes between Version 33 and Version 34 of DevelopmentActivities/ORCHIDEE-CN/NH3
- Timestamp:
- 2017-07-11T15:13:11+02:00 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
DevelopmentActivities/ORCHIDEE-CN/NH3
v33 v34 1 = Work of Thomas Audoux carried out during an internship on NH3 =1 = Work of Thomas Audoux carried out during an internship on NH3 fluxes in ORCHIDEE model = 2 2 This has been done with rev 4164. My work is divided in 3 parts, the first part corresponds to the work done for the internship report. Then, the second part is a brief bibliography about emission factors. Finally, the last part corresponds to the second part of my internship: analyzing ammonia emission factor. 3 3 == I. Evaluation of ammonia emissions results calculated by the model ( Fertilizer =0 ) == … … 25 25 * Humidity 26 26 * pH 27 * Nitrogen Pools (Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O) [[BR]]27 * Nitrogen Pools (Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O) 28 28 * Emissions (NH3) and Depositions (NHX) [[BR]] 29 29 path : /home/users/taudoux/documents/Results/Parameters/ [[BR]] … … 52 52 * Text documents in which are listed the different values. 53 53 54 === Comparison with data used by LMDz-INCA model ===55 path : /home/users/taudoux/documents/Results/Comparison/ [[BR]]56 57 Figure, created using "bar_plot_donnees.py", in which for each site are listed the average emission of NH3 for both model.58 54 59 55 … … 64 60 || [[Image(nh3_em_nh4_soil.png, 50%)]] || [[Image(nh3_em_n2o_soil.png, 50%)]] || [[Image(nh3_em_nox_soil.png, 50%)]] || [[Image(nh3_em_no3_soil.png, 50%)]] || 65 61 66 [[BR]] 67 62 Figure, created using "bar_plot_donnees.py", in which for each site are listed the average emission of NH3 for both model. 68 63 69 64 ||||||||= Correlation with NH3 emissions (g(NH3-N).m⁻².day⁻¹) =|| … … 79 74 || [[Image(Comparison.png, 50%)]] || 80 75 81 76 path : /home/users/taudoux/documents/Results/Comparison/ 82 77 83 78 [[BR]] … … 125 120 126 121 127 === Code was changed ===122 === Code was changed (v1) === 128 123 '''Separation of nitrogen input. Initially fertilizer, bnf and atmospheric deposition were in the same loop ''' 129 124 … … 145 140 * Temperature 146 141 * Humidity 147 * Nitrogen Pools ( Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O) and Leaching (NH4 and NO3)148 * Emissions (NH3, N2, N2O, NO x)142 * Nitrogen Pools (soil_nh4, soil_no3, Soil_nox and soil_n2o) and Leaching (NH4 and NO3) 143 * Emissions (NH3, N2, N2O, NOX) 149 144 * GPP 150 * Nitrogen input summed (fertilizer, NHx and NOx depositions) 145 * Nitrogen input summed (fertilizer, NHX and NOX depositions) 146 * SOIL_ACTIVE_N , SOIL_SURF_N , SOIL_SLOW_N , SOIL_PASSIVE_N , TOTAL_M_N LITTER_STR_AB_N , LITTER_MET_AB_N , LITTER_STR_BE_N , LITTER_MET_BE_N , LITTER_WOD_AB_N , LITTER_WOD_BE_N 151 147 152 148 Respectively in the folders Parameters (T° and Hum), Pools_fluxes, NH3_EM, DEP_FERT and GPP. … … 270 266 271 267 === C. EF of each parameters on each site for NH4%=0.5 and 0.8 === 268 272 269 path : /home/users/taudoux/documents/Results/B/EFc/ [[BR]] 273 270 Python : barplot_c.py (arg1= Name of the site) using files.txt in which values of nitrogen emissions and leaching are listed in order to plot EF in function of nitrogen inputs (fertilizer, NHx and NOx depositions). … … 355 352 356 353 357 In order to know if the loops are closed : [[BR]] 358 Loop at inorganic nitrogen scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; NH4 and NO3 Plant Uptake ; difference between nitrogen pools shown previously) divided by all nitrogen inputs (NOX and NHX deposition; Mineralisation but also Fertilizer) [[BR]] 354 The values we found seem to be a little odd... So, in order to know if they are exploitable, we looked at the closure of different loops : 355 356 * Loop at inorganic nitrogen scale : 357 358 We stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; NH4 and NO3 Plant Uptake ; difference between nitrogen pools shown previously) divided by all nitrogen inputs (NOX and NHX deposition; Mineralisation but also Fertilizer) [[BR]] 359 359 360 path : /home/users/taudoux/documents/Results/B/EFe/ [[BR]] 360 361 python : barplot_e.py arg1=name of the site [[BR]] 361 362 362 Values in abscissa are all thenitrogen inputs but the fertilization is for each point : 0 , 54.8 , 82.2 , 109.6 , 140 , 164.4 , 191.8 , 219.2 , 365.25 kg/ha/y.363 Values in abscissa correspond to all nitrogen inputs but the fertilization is for each point : 0 , 54.8 , 82.2 , 109.6 , 140 , 164.4 , 191.8 , 219.2 , 365.25 kg/ha/y. 363 364 364 365 ||||= DK-Ris =|| … … 397 398 398 399 399 Loop at ORCHIDEE scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; ; difference between nitrogen pools shown previously + difference between SOIL_ACTIVE_N , SOIL_SURF_N , SOIL_SLOW_N , SOIL_PASSIVE_N , TOTAL_M_N LITTER_STR_AB_N , LITTER_MET_AB_N , LITTER_STR_BE_N , LITTER_MET_BE_N , LITTER_WOD_AB_N , LITTER_WOD_BE_N); divided by all nitrogen inputs (NOX and NHX deposition and Fertilizer) [[BR]] 400 * Loop at ORCHIDEE scale : 401 We stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; ; difference between nitrogen pools shown previously + difference between SOIL_ACTIVE_N , SOIL_SURF_N , SOIL_SLOW_N , SOIL_PASSIVE_N , TOTAL_M_N LITTER_STR_AB_N , LITTER_MET_AB_N , LITTER_STR_BE_N , LITTER_MET_BE_N , LITTER_WOD_AB_N , LITTER_WOD_BE_N); divided by all nitrogen inputs (NOX and NHX deposition and fertilizer) [[BR]] 402 400 403 path : /home/users/taudoux/documents/Results/B/EFebis/ [[BR]] 401 404 barplot_ebis.py arg1=name of the site … … 436 439 437 440 438 We can see that there is a problem with NO3 stocks that are continually increasing or leached. It must be a problem with the denitrification. So we looked at the parameterization and saw that a variable called "anvf" takes an important role about denitrifiers. In order to know its effect, we ran a simulation in which this variable was deleted. [[BR]]441 We can see that there is a '''problem with NO3 stocks that are continually increasing or leached'''. It must be a problem with the denitrification. So we looked at the parameterization and saw that a variable called "anvf" takes an important role about denitrifiers. In order to know its effect, we ran a simulation in which this variable was deleted. [[BR]] 439 442 Here are the results : [[BR]] 443 440 444 path : /home/users/taudoux/documents/Results/B/EFi/ [[BR]] 441 445 python : barplot_init.py arg1=name of the site … … 460 464 461 465 462 Code was whanged : run_off, the variable "harvest_above_n" were added and ''the denitrifiers activity was changed'' (UPPER BOUND). [[BR]] 463 464 path of files.txt : /home/users/taudoux/documents/Results/B/Simul 466 === Code was changed (v2) === 467 Run_off, the variable "harvest_above_n" were added and ''the denitrifiers activity was changed'' (UPPER BOUND). [[BR]] 468 469 path of files.txt : /home/users/taudoux/documents/Results/B/Simul_ratio0.5_v2[[BR]] 465 470 We looked again at the differents loops. [[BR]] 466 471 467 Loop at inorganic nitrogen scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; NH4 and NO3 Plant Uptake ; difference between nitrogen pools shown previously) divided by all nitrogen inputs (NOX and NHX deposition; Mineralisation but also Fertilizer) [[BR]] 472 * Loop at inorganic nitrogen scale 468 473 469 474 path : /home/users/taudoux/documents/Results/B/EFf/ 470 475 python : barplot_f.py arg1=name of the site 471 476 472 Values in abscissa are all the nitrogen inputs but the fertilization is for each point : 0 , 54.8 , 82.2 , 109.6 , 140 , 164.4 , 191.8 , 219.2 , 365.25 kg/ha/y.473 474 477 ||= DK-Ris =|| 475 478 ||= 0.5 =|| … … 506 509 ||= [[Image(US-Ne3_0.5f.png, 50%)]] =|| 507 510 508 Loop at ORCHIDEE scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; ; difference between nitrogen pools shown previously + difference between SOIL_ACTIVE_N , SOIL_SURF_N , SOIL_SLOW_N , SOIL_PASSIVE_N , TOTAL_M_N LITTER_STR_AB_N , LITTER_MET_AB_N , LITTER_STR_BE_N , LITTER_MET_BE_N , LITTER_WOD_AB_N , LITTER_WOD_BE_N) and the variable harvest_above_n; divided by all nitrogen inputs (NOX and NHX deposition and Fertilizer) [[BR]] 511 * Loop at ORCHIDEE scale 512 509 513 path : /home/users/taudoux/documents/Results/B/EFg/[[BR]] 510 514 python : barplot_g.py arg1=name of the site 515 511 516 ||= DK-Ris =|| 512 517 ||= 0.5 =|| … … 544 549 545 550 546 Code was changed : ''the denitrifiers activity was changed'' (LOWER BOUND) 551 === Code was changed === 552 553 ''the denitrifiers activity was changed'' (LOWER BOUND) 547 554 path of files.txt : /home/users/taudoux/documents/Results/B/RundBis [[BR]] 548 555 … … 550 557 We looked again at the differents loops. [[BR]] 551 558 552 Loop at inorganic nitrogen scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; NH4 and NO3 Plant Uptake ; difference between nitrogen pools shown previously) divided by all nitrogen inputs (NOX and NHX deposition; Mineralisation but also Fertilizer) [[BR]] 553 559 * Loop at inorganic nitrogen scale 554 560 path : /home/users/taudoux/documents/Results/B/EFfbis/ [[BR]] 555 561 python : barplot_fbis.py arg1=name of the site 556 562 557 Values in abscissa are all the nitrogen inputs but the fertilization is for each point : 0 , 54.8 , 82.2 , 109.6 , 140 , 164.4 , 191.8 , 219.2 , 365.25 kg/ha/y.558 559 563 ||= DK-Ris =|| 560 564 ||= 0.5 =|| … … 591 595 ||= [[Image(US-Ne3_0.5fbis.png, 50%)]] =|| 592 596 593 Loop at ORCHIDEE scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; ; difference between nitrogen pools shown previously + difference between SOIL_ACTIVE_N , SOIL_SURF_N , SOIL_SLOW_N , SOIL_PASSIVE_N , TOTAL_M_N LITTER_STR_AB_N , LITTER_MET_AB_N , LITTER_STR_BE_N , LITTER_MET_BE_N , LITTER_WOD_AB_N , LITTER_WOD_BE_N) and the variable harvest_above_n; divided by all nitrogen inputs (NOX and NHX deposition and Fertilizer) [[BR]] 597 * Loop at ORCHIDEE scale 594 598 595 599 path : /home/users/taudoux/documents/Results/B/EFgbis/ [[BR]] 596 600 python : barplot_gbis.py arg1=name of the site 601 597 602 ||= DK-Ris =|| 598 603 ||= 0.5 =|| 599 ||= [[Image(DK-Ris_0.5gbis.png, 50%)]] =|604 ||= [[Image(DK-Ris_0.5gbis.png, 50%)]] 600 605 601 606 ||= FR-Gri =|| … … 629 634 ||= [[Image(US-Ne3_0.5gbis.png, 50%)]] =|| 630 635 631 Those results seem to suit good compared to others. So, in order to compare values of the emission with data we can find in the literature we calculate the emission factors of different fluxes. Plot EF in function of nitrogen inputs (fertilizer, NHx and NOx depositions)636 Those results seem to suit better compared to others. So, in order to compare values of the emission with data we can find in the literature we calculate the emission factors of different fluxes. Plot EF in function of nitrogen inputs (fertilizer, NHx and NOx depositions) [[BR]] 632 637 633 638 path : /home/users/taudoux/documents/Results/B/EFz/ [[BR]] … … 669 674 670 675 Globally, by comparing with the first simulation, we can see that NH3 and NOx Emission Factor decreased in favor of NO3 and NH4 leaching. Moreover, the nitrogen seems to no longer be blocked at the nitrate state making it possible to be denitrified and emitted by N2, N2O, NOx emissions or to be leached as NOx. 676 677 ---- 678 679 === Brief look at the global simulation done by Palmira MESSINA === 680 This has been done with rev ?. [[BR]] 681 682 I looked how ammonia emissions respond in this global simulation in order to see if the behavior previously seen site to site is the same on a larger scale. 683 684 ||||||||= Average emission of ammonia (g(NH3-N).m⁻².day⁻¹) in function of the parameter considered =|| 685 ||= Temperature (K) =||= Humidity (-) =||= pH (-) =||= Deposition g(NHx).m^-2 =|| 686 ||= [[Image(nh3_em_t2m_g.gif, 50%)]] =||= [[Image(nh3_em_hum_g.gif, 50%)]] =||= [[Image(nh3_em_pH_g.gif, 50%)]] =||= [[Image(nh3_em_nhx_dep_g.gif, 50%)]] =|| 687 ||= SOIL NH4 g(NH4).m^-2 =||= SOIL N2O g(N2O).m^-2 =||= SOIL NOX g(NOx).m^-2 =||= SOIL NO3 g(NO3).m^-2 =|| 688 ||= [[Image(nh3_em_soil_nh4_g.gif, 50%)]] =||= [[Image(nh3_em_soil_n2o_g.gif, 50%)]] =||= [[Image(nh3_em_soil_nox_g.gif, 50%)]] =||= [[Image(nh3_em_soil_no3_g.gif, 50%)]] =||
