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, 16th November

With Nicolas, Juliette, Didier and Anne.

Regarding the venue of Anne, speaking about the coupled version of models. As done with Maureen, the coupling orchidee/inca is available, remains to be checked. A meeting is scheduled for the following monday between Karine and Anne to see what to do.

Invitation to participate to the "Retraite of CMC" the next two days.

Explanation and summary of the reconstruction using NMVOC as a proxy for N₂O emissions. Two methods were tested: the first one from 1970 to 2019 using the csv data, reconstruction per sector and per region; the second one from 1850 to 2019 using the proxy, reconstruction per region. It exists a slight difference in total emissions between the two.

=> Need to check where it comes from. 2 hypothesis: reconstruction from proxy uses the 2000's netCDF file as a reference and not the 1970's file (because 1970 is the year used for the reconstruction); the average of the eight factors calculated per sector doesn't correspond to the factor computed per region.

TO DO

• See the coupled configuration;
• Check where the difference between reconstructed emissions comes from;
• With Didier, black box model with land, ocean, ant emissions.

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.