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, 25th January

With Juliette, Anne, Nicolas, Marion and Didier.

Simulations with coupled model are good: modifications on the differents files are done. Comparison tests have been done.

Orchidee v3 has been installed and run. It exists a slight difference between Orchidee3 and Orchidee2_2 (Quest) due to new nitrogen computed fluxes.

TO DO

• Transmit nitrogen fluxes from Orchidee to Inca (using Maureen's code);
• With Didier, change names and variables in sflx_noBio: each species need to have a land, ocean and miscellaneous component (on Obelix);
• With Anne, include those changes in the INCA code to read good fields;
• Test the reproducibility of the code by doing 2*2 years or 5* 2years simulation;
• With Nicolas, create input states from Orchidee (with LMDZOR to achieve an equilibrium state): once for pre industrial period and an other for 2000's.

Last intermediate conclusions

PRE INDUSTRIAL RUN

Some miscellaneous characteristics:

• GES configuration,
• 39 vertical layers, resolution 144x142,
• from fictives years 1970 to 2015,
• no anthropogenic emissions,
• BBG emissions stayed at 1850 (as well as aircraft).

Quantities to remember: (mean of 5 years)

• Total N₂O loss: 10.73 TgN/y
• Total N₂O emissions: 10.73 TgN/y
• Burden : 1368 TgN (ref Tian= 1462 TgN in the 1980s)

PRESENT-DAY RUN

Some miscellaneous characteristics:

• GES configuration,
• 39 vertical layers, resolution 144x142,
• 'real' years 2000 to 2007,
• anthropogenic emissions.

Quantities to remember: For 2000:

• Total N₂O loss: 11.78 TgN/y
• Total N₂O emissions: 10.73 TgN/y (ref Tian= 9.7 TgN/yr natural + 0.6 bbg + 0.3 waste water)
• Burden : 1514 TgN (ref Tian= 1555 TgN in 2007-2016)

Last notebooks/working documents

13. Proxy_reconstruction_preindustriel.ipynb: With NMVOC emissions from pre industrial, analyses evolution and correlation between NMVOC and N₂O. (not very useful except to have dataframe for the next step)
14. Proxy_reconstruction_preindustriel-NOx.ipynb: same as the previous one but using NOx emissions rather than NMVOC. Made because emissions converge to zero in 1750, which is not the case for NMVOC.
15. create_proxy_n2o_via_nmvoc.ipynb: using NMVOC data, reconstruction of a csv files with emissions per country from 1850 to 1969. Also, reconstruction without RCO sector (in order to begin more around 0 Tg/Yr? in 1850).

+ on Irene:

• on workdir: scripts_nco/prepare_N2OLOSS.job: will compute N2O total loss, N2O phtolyse loss, N2O o1d loss and N2O burden. This file creates some nc file on scratch directory (in Loss_files) and some txt file in place.

+ 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_country_1970_1999_ncfiles.sh: final step on nc files while transforming excel data into nc files.
• Transform_into_ncfile_preindustrial.sh: script corresponding to the notebook Transform_country_csv_to_nc_1970_2019.ipynb for preindustrial period.