wiki:ESM2025-N-cycle

Version 46 (modified by klaurent, 16 months ago) (diff)

--

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:

Last meeting

On Wednesday, 30th November

With Nicolas, Didier and Marion.

Two main points:

  1. Difference in the transfer of information into a netCDF file. Two methods are used, it is necessary to understand why they give different results. A document will be written to explain them. It is important to know where the problem comes from because we want to use good reconstruction emissions in INCA after.
  1. Black box model: from emissions, use a little model in order to reconstruct concentration. This works well from 1950 to the end of the period but we get a plateau before that. This plateau is explained by the constant emissions we have for this period. By making an inverted model, it is possible to reconstruct emissions from concentration (from observations in CMIP6). It shows we need to get two or three small bumps. The absence of these emissions in ours could result in an absence of ENSO for the oceanic part or in imprecise industrial emissions for this period or a modification in land use surfaces...

TO DO

  • Write a document explaining reconstruction methods and the transformation into netCDF file;
  • Investigate the difference between the two methods (either by region or by region+sector), a solution could by to use a weighted factor for the sector part;
  • Look on the new inventory, found by Joeran;
  • Continue development on the coupled model.

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 N2O loss: 10.73 TgN/y
  • Total N2O 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 N2O loss: 11.78 TgN/y
  • Total N2O 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

  1. Proxy_reconstruction_preindustriel.ipynb: With NMVOC emissions from pre industrial, analyses evolution and correlation between NMVOC and N2O. (not very useful except to have dataframe for the next step)
  2. 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.
  3. 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.

Attachments (17)

Download all attachments as: .zip