ESM2025-N-cycle

This page is dedicated to the ESM2025 project (Earth System Models - European Project), especially for the N-cycle modelling. On this topic are working Marion Gehlen, Juliette Lathière, Didier Hauglustaine, Nicolas Vuichard and Karine Laurent.

You will find information on the evolution of the project on these different pages:

• Meetings: report of meetings every 15 days;
• Intermediate Conclusions: some partial conclusions;
• Notebooks list: description of working documents;
• References: some interesting references.

Last meeting

On Wednesday, 10th May

With Didier and Nicolas. Invited Thibaut and Maureen.

Discussion about simulations for ESM2025/WP1. Two different simulations: one with DMS and an other with nitrogen (NH₃ + N₂O + NO (?)). We have to know if versions of code (especially LMDZ) are similar. Thibaut is in charge of scheduling a meeting with ocean group.

We speak a little on the deliverable for WP6 (some minor correction to be made). The second GA of ESM2025 will be the next week and I have to prepare a poster: correction and finalization to do.

Box model: new results with fluxes from NMIP2: more emissions are emitted then it increases the N₂O concentration. Discussion about a pulse model in order to find the steady-state concentration. (how do they do for carbon dioxide or methane?)

TO DO

• Test the new pulse model and compare the two versions;
• Run pre industrial simulation (when Josefine files ok). First, offline with total emissions, then offline with sum of emissions (land/ocean/misc), then coupled version.
• Continue the poster.

Last intermediate conclusions

Reconstruction of anthropogenic N₂O emissions in pre industrial

1. from 1850 to 1969: proxy with NMVOC and linear regression, only total emissions per region (no information about sectors);
2. from 1850 to 1969: reconstruction per region using 1970's emissions as reference;
3. from 1970 to 1999 (extended to 2019): reconstruction per region and sector, using 2000's emissions as reference. (+modification of agriculture sector by factor 0.40 in order to drop soil emissions.)

Remark: A slight difference exists between data from csv file and data in netCDF file; it can be shown with the comparison of 'Notebooks/N2O_sector_country_via_ncfile_maks.csv' (data from netCDF file) and 'Notebooks/total_em_CEDS_sector_country_1970_2019_repartEDGAR2.csv'. This might be due to the maks used...

Last notebooks/working documents

16. compare_ant_emissions.ipynb: Different methods are used to reconstruct nc files for N₂O emissions. This notebook served to visualisation of global emissions.
17. check_coeff_reconstruction_proxy.ipynb: with different method of reconstruction, result are different, this notebook explore coefficient used to reduce emissions for pre industrial period and after (coeff only per region and coeff per region and sector).
18. black_box_emissions_model.ipynb: little box model to compute concentration with given emissions and inverted model.

+ on Irene:

• on workdir: scripts_nco/prepare_N2OLOSS.job: will compute N2O total loss, N₂O photolyse loss, N₂O o1d loss and N₂O burden. This file creates some nc file on scratch directory (in Loss_files) and some txt file in place.
• on workdir: scripts_ferret: some scripts to draw more or less beautiful pictures.
• path to simulation directory (on work-gen2201-cont003):
  • with Orchidee_Quest: ESM2025/modipsl/config/LMDZORINCA_v6
  • with Orchidee_3: ESM2025/modipsl_v6.3/config/LMDZORINCA_v6
  • with Maureen's Code: extract_Maureen/modipsl/config/LMDZORINCA_v6

+ on Ciclad/Spirit?:

• INCAFLXANT_1970_1999.sh: bash script to create sflx_lmdz_ANT for irene (source code in ANT/McDuffie). The creation is made on datadir.

+ on Obelix:

• CEDS_transform/nco_transform_1970_1999_ncfiles.sh and nco_transform_country_1970_1999_ncfiles.sh: final step on nc files while transforming excel data into nc files.
• Interpol_With_rmp.ipynb: Notebook of O. Marti that transform eORCA1 grid into regular one (144x142 here).
• Ciclad_N2O_values.ipynb: Notebook to quickly see evolution of BBG and ANT emissions, needs a grep command on Ciclad and file transfer to Obelix.
• Transform_into_ncfile_preindustrial.sh: script corresponding to the notebook Transform_country_csv_to_nc_1970_2019.ipynb for preindustrial period.