| 35 | | * //Bouwman Inventory// : |
| 36 | | 1. 846 N,,2,,O emission measurements in agricultural fields and 99 measurements for NO emissions |
| 37 | | 2. The data set includes literature reference; location of the Measurement; climate; soil type, texture, organic C content, N content, drainage, and pH; residues left in the field; crop; fertilizer type; N application rate; method and timing of fertilizer application; NH,,4,,+ application rate (for organic fertilizers), N,,2,,O/NO emission/denitrification (expressed as total over the measurement period, as % of N rate, and as % of N rate accounting for control); measurement technique; length of measurement period; frequency of the measurements; and additional information, such as year/season of measurement, information on soil, crop or fertilizer management, specific characteristics of the fertilizer used, and specific weather events important for explaining the measured emissions. |
| 38 | | 3. global gridded (1°x1° resolution) data bases of soil type, soil texture, NDVI (vegetation indices) and climate. |
| 39 | | 4. global emission thus calculated is 6.8 Tg N,,2,,O-N y-1. The tropics (± 30° of the equator) contribute 5.4 Tg N,,2,,O-N y-1 and the emission from extra-tropical regions (poleward of 30°) is 1.4 Tg N,,2,,O-N y-1 . |
| 40 | | |
| 41 | | |
| 42 | | * //Transcom Inventory// : |
| 43 | | 1. "Emissions from natural soils (6–7 TgN yr −1 ) account for 60–70 % of global N2O emissions (Syakila and Kroeze, 2011; Zaehle et al., 2011). The remaining 30–40 % of emissions is from oceans (4.5 TgN yr −1 )" |
| 44 | | 2. Five different inversion frameworks (chemistry transport model) : MOZART4 (2.5° × 1.88°), ACMt42167 (2.8° × 2.8°), TM3 (5.0° × 3.75°), TM5 (6.0° × 4.0°), LMDZ4 (3.75° × 2.5°). |
| 45 | | 3. Data set from Orchidee O-CN, Pisces, edgar-4.1, gfed-2 and from different category (terrestrial biosphere, ocean, waste water, solid waste, solvents, fuel prod, ground transport, industry combustion, residential and other combustion, shipping, biomass burning) |
| 46 | | |
| 47 | | * Values of N,,2,,O emissions per inventory : |
| 48 | | |
| 49 | | || (MtN/yr) || Ocean || Land || Total || (MtN,,2,,O/yr) || Ocean || Land || Total || |
| 50 | | ||**Bouwman** || 3.595 || 7.531 || 11.126 || xxx || 5.650 || 11.835 || 17.485 || |
| 51 | | ||**Transcom** || 4.222 || 10.573 || 14.795 || xxx || 6.634 || 16.614 || 23.248 || |
| 52 | | ||**Piscidee** || 3.962 || 7.107 || 11.069 || xxx ||6.226 || 11.168 || 17.394 || |
| 53 | | |
| 54 | | * Comparison between Incasflx outputs and masks on NetCDF files: |
| 55 | | |
| 56 | | || |||||| - Output INCASFLX - |||||| - Notebook computing (FracMask) - |||||| - Notebook computing (0/1 Mask) - || |
| 57 | | || (MtN2O/yr) || Ocean || Land || Ratio O/L || Ocean || Land || Ratio O/L || Ocean || Land || Ratio O/L || |
| 58 | | ||**Bouwman** || 5.6505 || 11.835 || 0.478 || 6.58 || 10.945 || 0.6012 || 5.973 || 11.552 || 0.517 || |
| 59 | | || || 32.3% || 67.7% || || 37.6% || 62.4% || || 34.1% || 65.9% || || |
| 60 | | ||**Transcom** || 6.634 || 16.614 || 0.399 || 8.879 || 14.421 || 0.616 || 8.118 || 15.184 || 0.535 || |
| 61 | | || || 28.5% || 71.5% || || 38.1% || 61.9% || || 34.8% || 65.2% || || |
| 62 | | ||**Piscidee** || 6.226 || 11.168 || 0.557 || 6.651 || 10.782 || 0.617 || 6.095 || 11.339 || 0.538 || |
| 63 | | || || 35.8% || 64.2% || || 38.2% || 61.8% || || 34.9% || 65.1% || || |
| 64 | | |
| 65 | | |
| 66 | | __General remarks__: Variable for Bouwman = fn2o_oce & fn2o_soil. Land emissions > Ocean emissions. Ocean emissions ~ 6 Mt/yr. Land emissions -> differences because of different period (pre-industrial & nowadays). Per line, Ocean+Land are equal independently of the calculation made. |
| 67 | | |
| 68 | | * CEDS emissions for N,,2,,O (and NH,,3,,): |
| 69 | | [[Image(N2O_CEDS_supplementary_avion.png, 700px)]] |
| 70 | | [[Image(NH3_CEDS_supplementary.png, 700px)]] |
| 71 | | |
| 72 | | == Meeting Reports == |
| 73 | | |
| 74 | | |
| 75 | | === On Wednesday, 24th August === |
| 76 | | |
| 77 | | In person with Didier and Nicolas. PDF support [see attached file[[Image(meeting_220824.pdf)]]]. |
| 78 | | |
| 79 | | Resume of august simulations and work. |
| 80 | | |
| 81 | | * Transformation of CEDS emissions: from 1970 to 1999, use of netCDF file 2000 with a factor per sector to create netcdf emissions. Notebooks **Transform_2000nc_from1970_to1999_via_CEDScsv.ipynb** and **work_on_data_to_complete_sectors_with_unknows_categories.ipynb** |
| 82 | | * 2 simulations: L39.v03 as a base (pre-industrial, constrained to 285ppm), then pi_piscideeL (free, pre industrial, 40 years). |
| 83 | | |
| 84 | | Results seem good regarding vertical pressure profile, total loss N,,2,,O, pattern of emissions. Need to investigate more and continue. |
| 85 | | |
| 86 | | |
| 87 | | __//**TO DO**//__ |
| 88 | | * Make graph of pressure for those two simulations, |
| 89 | | * Continue piscideeL 5 years and compute mean of loss, |
| 90 | | * Resume differences between former simulations, |
| 91 | | * Re-create netCDF files with only manure_management (and not soil_emissions) for agricultural sector. |
| 92 | | |
| 93 | | |
| 94 | | |
| 95 | | === On Tuesday, 12th July === |
| 96 | | |
| 97 | | In person with Didier. |
| 98 | | |
| 99 | | Discussion of simulations and set-up to make. |
| 100 | | |
| 101 | | Check N,,2,,O loss in 1999 for v03 simulation (to be continued) and create pi_piscideeL run from it to turn 40 years (1980-2020, L means libre/free/not rescaled at 285 ppm). |
| 102 | | |
| 103 | | Explanation of how to create ncfiles from 1970 to 1999 from 2000's geographical repartition and ratio computed with total emissions per sector (csv file). |
| 104 | | Once done, adapt ciclad's code to transform them on Incasflx 144x142 grid. |
| 105 | | |
| 106 | | |
| 107 | | === On Monday, 11th July === |
| 108 | | |
| 109 | | In person with Didier, Nicolas and Juliette. |
| 110 | | |
| 111 | | Conversation on Workshop N,,2,,O source in Toulouse (6-8 July) : very interesting to see different communities (observation, inversion, modellers) |
| 112 | | |
| 113 | | First result with a decreasing total emission (from 285 ppm to 280 ppm in 20 years) : seems to be promising. |
| 114 | | |
| 115 | | Reflection on what is needed in order to **move to coupled system**. Remembrance of deadlines and expectations for ESM2025 project. Land emission per year ok, Nicolas will ask to Sarah or Laurent to have same emission for ocean. After that, work with Thibault Lurton and/or Anne Cozic to adapt code. |
| 116 | | |
| 117 | | **How to reconstruct anthropologic N,,2,,O emission between 1850 and 1970 ?** Try to find some proxy (population or chemical species that exists during this period) |
| 118 | | |
| 119 | | |
| 120 | | __//**TO DO**//__ |
| 121 | | * Reconstruct N,,2,,O emission with CEDS inventory with rescaling year by year (height with total budget and spatial with year 2000) |
| 122 | | * When all is fixed, simulation from 1970 to 2015 in order to see if the model is coherent or not (with a rescaling at ~300ppm : CMIP6). |
| 123 | | * Find information on Japanese model (the one of Prabir Patra) |
| 124 | | * Get more information on Pisces simulation (does it use flag or many change in order to have it? Is it possible to have emissions per year from pre-industrial to present time?) |
| 125 | | |
| 126 | | __//**Next**//__ |
| 127 | | * Computation of loss rate |
| 128 | | * Coupled model (with Anne and Thibaut) |
| 129 | | |
| 130 | | |
| 131 | | |
| 132 | | === On Wednesday, 29th June === |
| 133 | | |
| 134 | | In person with Didier and Juliette + in visio with Marion. |
| 135 | | |
| 136 | | Explanation of the remapping problem with pisces' file. The common CDO functions doesn't work (not conservative + white columns or column with outliers). |
| 137 | | |
| 138 | | Solution : See Germain if he has some same files (grid 294x362 and not 292x292) + in correspondence with O. Marti for a better interpolation. //work in progress// |
| 139 | | |
| 140 | | |
| 141 | | __//**TO DO**//__ |
| 142 | | * Understand and find a solution to the remapping problem. |
| 143 | | * Send an e-mail to Sarah Berthet, in charge of different projects on N,,2,,O : could be seen during the N,,2,,O workshop in Toulouse (from 6th to 8th July). |
| 144 | | * Simulations in 39 layers with invariants outputs in two different configurations. |
| 145 | | * On profile picture, try to invert axis + use logarithm scale. |
| 146 | | * Check availabilities for next meetings. |
| 147 | | |
| 148 | | ==== + Meeting with Ddiier ==== |
| 149 | | * See different scripts to compute the N,,2,,O-loss and N,,2,,O-burden. |
| 150 | | * Show what and how to change Incasflx files in order to use a irregular grid as input (avoiding the remapping problem). |
| 151 | | |
| 152 | | |
| 153 | | === On Wednesday, 15th June === |
| 154 | | |
| 155 | | In person with Nicolas and Didier. |
| 156 | | |
| 157 | | Waiting for the first result with the 39-layers configuration. |
| 158 | | |
| 159 | | After the General Assembly of ESM2025, pisces emissions seems to be too high... In fact, there was a problem of units conversion (mostly because of a misunderstanding of input units (mol N rather than mol N,,2,,O). |
| 160 | | |
| 161 | | Nicolas asked if lightning is implemented because it seems to have an total emission between 0.5 and 2 Tg/Nyr. |
| 162 | | |
| 163 | | Some works will be presented during the N,,2,,O workshop in Toulouse from July, 6th to 8th. |
| 164 | | |
| 165 | | __//**To DO**//__ |
| 166 | | * Analyse the vertical profile when the 39-layers configuration of one year is done. |
| 167 | | * + See Didier for the next steps. |
| 168 | | * Re-compute the calculation with the good value of pisces emissions. |
| 169 | | |
| 170 | | === On Wednesday, 1st June === |
| 171 | | |
| 172 | | In person with Nicolas, Juliette, Didier. Special guest : Anne. |
| 173 | | |
| 174 | | Special demands to Anne : |
| 175 | | * Passage from 79 layers to 39 : is it easy and feasible ? Seems to change some configuration's files (oxidants for example). |
| 176 | | * svn for INCA ok, if need to change orchidee's sources, ask to Nicolas (but no need for this moment). |
| 177 | | * Make some time series for long run test (feasible with the training course), |
| 178 | | * Add output variables that would be computing into the code (have to see Anne when it will be time). |
| 179 | | |
| 180 | | Show results on Bouwman (with aqua and soil flux): slight difference again. Using mask 0/1 (and not fractional one) is better. |
| 181 | | |
| 182 | | Show on-going result of CEDS inventory (from McDuffie): problem with units, needs to calculate annual mean to correspond to article's images. |
| 183 | | |
| 184 | | __//**To DO**//__ |
| 185 | | * Finish charts of CEDS inventory. Notebooks used: CEDS_charts_supplementary_data_N2O.ipynb and Charts_N2O_CEDS.ipynb. |
| 186 | | * Schedule a meeting with Didier when 39 layers' configuration is ready. |
| 187 | | * Make simulation with oce and soil flux for Bouwman (not with aqua !!). |
| 188 | | |
| 189 | | === On Wednesday, 18th May === |
| 190 | | |
| 191 | | In person with Juliette, Didier, Nicolas. |
| 192 | | |
| 193 | | * Notebook **Budget_inventories** : total budget of N,,2,,O emissions (in tables), differences found between outputs INCASFLX and calculations in notebook (with mask, area and change units), different charts (by latitude, by land/ocean part, with histograms...) |
| 194 | | |
| 195 | | |
| 196 | | Using the mask with fraction: |
| 197 | | |
| 198 | | || FracMask |||||| - Output INCASFLX - |||||| - Notebook computing - || |
| 199 | | || (Mt/yr) || Ocean || Land || Ratio O/L || Ocean || Land || Ratio O/L || |
| 200 | | ||**Bouwman** || 5.650486 || 11.83547 || 0.478 || 13.967 || 16.726 || 0.835 || |
| 201 | | || || 32.3% || 67.7% || || 45.5% || 54.5% || || |
| 202 | | ||**Transcom** || 6.634 || 16.614 || 0.399 || 8.879 || 14.421 || 0.616 || |
| 203 | | || || 28.5% || 71.5% || || 38.1% || 61.9% || || |
| 204 | | ||**Piscidee** || 6.432 || 7.107 || 0.905 || 6.579 || 6.989 || 0.941 || |
| 205 | | || ||47.5% || 52.5% || || 48.5% || 51.5% || || |
| 206 | | |
| 207 | | Using the mask with 0/1 values: |
| 208 | | |
| 209 | | || 0/1 Mask |||||| - Output INCASFLX - |||||| - Notebook computing - || |
| 210 | | || (Mt/yr) || Ocean || Land || Ratio O/L || Ocean || Land || Ratio O/L || |
| 211 | | ||**Bouwman** || 5.650486 || 11.83547 || 0.478 || 12.797 || 17.899 || 0.714 || |
| 212 | | || || 32.3% || 67.7%|| || 41.7% || 58.3% || || |
| 213 | | ||**Transcom** || 6.634 || 16.614 || 0.399 || 8.118 || 15.184 || 0.535 || |
| 214 | | || || 28.5% || 71.5% || || 34.8% || 65.2% || || |
| 215 | | ||**Piscidee** || 6.432 || 7.107 || 0.905 || 6.174 || 7.395 || 0.834 || |
| 216 | | || || 47.5% || 52.5% || || 45.5% || 54.5% || || |
| 217 | | |
| 218 | | |
| 219 | | Those differences can be explained by the different power that may exist between marine or continental fluxes at coasts notably. With mask, these two fluxes have the same importance so it changes total flux. |
| 220 | | |
| 221 | | Values per latitudes: |
| 222 | | || (Tg/yr) ||90-30S|| 30S-0||0-30N||30-90N||Total|| |
| 223 | | || **Bouwman** ||5.728||7.115||9.948||8.392||31.183|| |
| 224 | | || **Trancsom** ||2.714||7.432||8.762||4.857||23.765|| |
| 225 | | || **Piscidee** ||1.736||4.413||4.392||3.302||13.843|| |
| 226 | | |
| 227 | | |
| 228 | | Calcul file by file: |
| 229 | | From Transcom (orchidee+pisces) : (no mask + not output of sflx) -> correspond to values of table 2 thompson part 1. |
| 230 | | |
| 231 | | || |||||| - Output INCASFLX - |||||| - Notebook computing - || |
| 232 | | || (Mt/yr) || Ocean || Land || Ratio O/L || Ocean || Land || Ratio O/L || |
| 233 | | ||**Transcom** || 6.634174 || 16.61481 || 0.399293 || 4.23125 || 10.59692 || 0.399291 || |
| 234 | | |
| 235 | | |
| 236 | | |
| 237 | | Transcom is not a pre-industrial inventory, it's actual. So, it's better to work with piscidee inventory. |
| 238 | | |
| 239 | | |
| 240 | | __//**To DO**//__ |
| 241 | | * Modify Bouwman run to take only 'oce' and 'soil' variables -> relaunch notebook to compare if it's better, |
| 242 | | * Continue modification in INCA's code. |
| 243 | | |
| 244 | | === On Friday, 13th May === |
| 245 | | |
| 246 | | Informal with Didier. |
| 247 | | |
| 248 | | Preparation to change INCA's code in Irene. |
| 249 | | |
| 250 | | * Modify inca.card to take new NAT files and BBG from ciclad. |
| 251 | | * Modify mksflx_p2p.F90 to put flx_n2o_ant to zero and add aflux terms. |
| 252 | | * Modify sflx_inti.F90 to take into account BBG emissions. |
| 253 | | * Modify set_ub_vals.F90 to rescale N,,2,,O flux (with parallel part and constant N,,2,,O flux at 285 ppb) |
| 254 | | |
| 255 | | Launch different runs in order to check each step. |
| 256 | | |
| 257 | | __//**Next**//__ |
| 258 | | * Change ANT in ciclad, |
| 259 | | * Put flx_n2o_ant to 0, |
| 260 | | * Erase the rescaling, |
| 261 | | * Add N2O loss as an output... |
| 262 | | |
| 263 | | === On Wednesday, 4th May === |
| 264 | | |
| 265 | | In person with Juliette and Marion. |
| 266 | | |
| 267 | | * Presentation of the Wiki page. |
| 268 | | * Presentation of the two notebooks. |
| 269 | | * Discussion on the different conceptions used in ocean and land communities (different grids, different units in the flux (mol/m2/s vs kg/m2/s)...). |
| 270 | | |
| 271 | | On Analyse_3_sflx_NAT_dyn.ipynb, discussions on visualization with scale, calculation of the full emission on earth per year or per month. |
| 272 | | |
| 273 | | __//**To DO**//__ |
| 274 | | * Make an analysis per month, per hemisphere, in order to get the total flux emissions for each inventory. |
| 275 | | * Take information on where the inventories come from (especially Bouwman and Transcom). |
| 276 | | * Draw tables which resume datas used, total emissions, models used in files... |
| 277 | | * Compare native grid and regular grid for pisces. |
| 278 | | * Compare values of each file used with outputs of INCASFLX. |
| 279 | | * Compare total emissions on ocean and land with literature (ciais, ippc) and Bouwman and Trancsom inventories. |
| 280 | | |
| 281 | | === On Wednesday, 20th April === |
| 282 | | |
| 283 | | In person with Nicolas. |
| 284 | | |
| 285 | | Discussion on: |
| 286 | | |
| 287 | | * New page on Wiki/igcmg (this one !) |
| 288 | | * Explanation of the created notebooks (Data_viewer and Analyse_3_sflx_NAT_dyn) |
| 289 | | * Data_viewer : interactive ok but add choices/texts... |
| 290 | | * Analyse_3_sflx_NAT_dyn : change the seasonal analysis made. |
| 291 | | * Test with Nicolas' network to know about access to my notebooks. |
| 292 | | * Speak about dods/orchidee to avoid downloading images on wiki. |
| 293 | | |
| 294 | | Now, focus on N2O with long living time, but after, we can test on other compounds. |
| 295 | | |
| 296 | | __//**To DO**//__ |
| 297 | | * Change colorbar with mean+/-std, add some texts on Notebooks and add units on charts |
| 298 | | * Erase rescaling in Bouwman to have a better comparison |
| 299 | | * Analyse the 3 datasets in a different way : global tables with sum per year, per region, per type (ocean/land) (watch out units) |
| 300 | | * Give "procedure" and access to notebooks to all |
| 301 | | * Compare with Auburn data (and Tian 2020) |
| 302 | | * Look at Mendeley (or Zottero) for bibliography |
| 303 | | |
| 304 | | Nicolas: Paths to have regional masks (from Transcom) |
| 305 | | |
| 306 | | |
| 307 | | |
| 308 | | === On Wednesday, 6th April === |
| 309 | | |
| 310 | | In Visio with Nicolas, Juliette, Didier, Marion. |
| 311 | | |
| 312 | | Mainly, this meeting was to talk a bit more about the project and the work I(=Karine) can make before the training session (April 14th & 15th). |
| 313 | | |
| 314 | | Discussion on: |
| 315 | | * the first result made with INCAFLX and 3 different inventories, |
| 316 | | * the planning we can set up (which kind of simulation...) |
| 317 | | * the impact of N2O on chemistry (if everything is coupled), |
| 318 | | * the certainty (or uncertainty) of inventories used (to determine an appropriate scaling), |
| 319 | | * the work of Pierre for aquatic emissions (importance of river emissions), |
| 320 | | |
| 321 | | //Planning :// Run simulation without chemistry and no rescaling, first for pre-industrial, then until now to see/understand the correction we can make on both period (a different correction for each period). |
| 322 | | |
| 323 | | => This have to be detailed and confirmed |
| 324 | | |
| 325 | | Maybe run only 10 years to know where are sinks, determine some lifetime/seasonal cycle in order to calculate the scaling we may use. |
| 326 | | |
| 327 | | Suggestions to work with the GES branch for the code (branch for coupled model with CO,,2,, , CH,,4,, and N,,2,,O) |
| 328 | | |
| 329 | | **Meeting every 2 weeks same day, same time (ie, Wednesday at 10:00)** |
| 330 | | |
| 331 | | |
| 332 | | __//**To DO**//__ |
| 333 | | * Ask for writing a Wiki page on IGCMG |
| 334 | | * Log on JupyterHub |
| 335 | | * Make an "experiment plan" for future simulations |
| 336 | | * Analyse the results that I've already (seasonal analyses, comparison between inventories with Bouwman as a reference, recap origin/grid resolution...) |
| 337 | | * Show the aquatic variable of Bouwman (and maybe some others) |
| 338 | | |
| 339 | | Didier : emissions file for pre-industrial fires to be send to Karine |
| 340 | | |
| 341 | | == Bibliography == |
| 342 | | |
