Transport Working Group: Difference between revisions

From Geos-chem
Jump to navigation Jump to search
Line 107: Line 107:
<span style="color:darkorange">'''''This update will be included in [[GEOS-Chem 12#12.2.0|GEOS-Chem 12.2.0]].'''''</span>
<span style="color:darkorange">'''''This update will be included in [[GEOS-Chem 12#12.2.0|GEOS-Chem 12.2.0]].'''''</span>


The transport Working Group has proposed implementing additional tracers in GEOS-Chem for benchmarking purposes and diagnosing transport in GEOS-Chem compared to GEOS-5.
The Transport Working Group has proposed implementing additional tracers in GEOS-Chem for benchmarking purposes and diagnosing transport in GEOS-Chem compared to GEOS-5. The [[#List of tracers in GEOS|list of tracers in GEOS]] is provided below and is implemented via the tracer gridded component "TR_GridComp" in GEOS. In GEOS-Chem, the tracer suite has been added as an extension of the [[Rn-Pb-Be simulation]] (simulation type "1") by default, but may be added to any GEOS-Chem simulations by following the [[#Implementation in GEOS-Chem|instructions below]].
 
'''''Hongyu Liu wrote:'''''
 
:Andrea and I would like to seek your comment on the effort needed to implement additional tracers in GCC and GCHP. See below for the [[#List of tracers in GEOS|list of tracers in GEOS]] that Andrea sent earlier. 
 
:It would be straightforward to implement all these tracers in GCHP. Per Andrea, all we need to do is to download some additional code (called the tracer gridded component "Tracer_GridComp" in GEOS).  
 
:If we want the same tracer suites in GCC (e.g., an extension of current [[Rn-Pb-Be_simulation|radionuclide tracers]], simulation type "1") as in GCHP, how much effort would it be? Is there any practical reason that a reduced number of tracers is preferred in GCC?
 
'''''[[User:Bmy|Bob Yantosca]] replied:'''''
 
:It looks as if most of these are passive species with a fixed atmospheric lifetime, right?  If so, then it is very straightforward to add these to any GEOS-Chem simulation.  See our [[HEMCO examples#Passive_species|instructions for adding passive species to GEOS-Chem]].
 
:So basically in the passive species menu of input.geos you would list the names of each passive species, its default concentrations, and atmospheric lifetime is.  If the species has an emissions flux, you can specify that with a corresponding entry in the HEMCO_Config.rc file, as described on the wiki.
 
:I think the only modification we might need to make is that you are currently limited to 20 passive species but that is easy to extend by just changing a parameter value.
 
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 18:42, 23 February 2018 (UTC)


=== List of tracers in GEOS ===
=== List of tracers in GEOS ===
Line 159: Line 141:
--[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 15:41, 15 June 2018 (UTC)
--[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 15:41, 15 June 2018 (UTC)


=== Tracer suggestions ===
=== Origin tracers ===


'''''Chris Holmes wrote:'''''
'''''Chris Holmes wrote:'''''
Line 171: Line 153:


--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 18:42, 23 February 2018 (UTC)
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 18:42, 23 February 2018 (UTC)
=== Implementation in GEOS-Chem ===
Several of the GEOS tracers listed above have been added as [[HEMCO_examples#Passive_species|passive species in GEOS-Chem]]. The new tracers are listed in the advected species menu of input.geos as:
%%% ADVECTED SPECIES MENU %%%:
Type of simulation      : 1
Species Entries ------->: Name
Species name            : Rn222
Species name            : Pb210
Species name            : Pb210Strat
Species name            : Be7
Species name            : Be7Strat
Species name            : Be10
Species name            : Be10Strat
Species name            : Passive
<span style="color:green">Species name            : TR_SF6
Species name            : TR_CH3I
Species name            : TR_CO_25
Species name            : TR_CO_50
Species name            : TR_e90
Species name            : TR_e90_N
Species name            : TR_e90_S
Species name            : TR_Uniform
Species name            : TR_Clock</span>
where <tt>TR_</tt> denotes that these species are from TR_GridComp of GEOS. In the passive species menu of input.geos, each species is defined along with the molecular weight, lifetime, default background concentration, and long name.
%%% PASSIVE SPECIES MENU %%%:
Passive species name    : Passive    1.0    -1      1.0e-7  Passive_tracer_for_mass_conservation_eval
<span style="color:green">Passive species name    : TR_SF6    146.0  1.0e11  1.0e-20  Sulfur_hexafluoride
Passive species name    : TR_CH3I    142.0  4.32e5  1.0e-20  Methyl_iodide
Passive species name    : TR_CO_25  28.0  2.16e6  1.0e-20  Anthro_CO_25day_tracer
Passive species name    : TR_CO_50  28.0  4.32e6  1.0e-20  Anthro_CO_50day_tracer
Passive species name    : TR_e90    1.0    7.776e6 1.0e-20  Constant_emission_tracer
Passive species name    : TR_e90_N  1.0    7.776e6 1.0e-20  Constant_emission_tracer_for_NH
Passive species name    : TR_e90_S  1.0    7.776e6 1.0e-20  Constant_emission_tracer_for_SH
Passive species name    : TR_Uniform 28.0  2.16e6  1.0e-20  CO-like_tracer_with_uniform_emission
Passive species name    : TR_Clock  1.0    -1      0.0      Age_of_air_in_seconds</span>
Species with an emissions flux also have a corresponding entry in the HEMCO_Config.rc file.
#==============================================================================
# SF6 emissions from EDGAR v4.2
#==============================================================================
0 SRC_2D_SF6  /n/regal/jacob_lab/msulprizio/v42_SF6_2008-1970_IPCC_2_nc/v2019-01/EDGAR_v42_SF6_IPCC_2.generic.01x01.nc emi_sf6 1970-2008/1-12/1/0 C xy kg/m2/s TR_SF6 - 1 1
#==============================================================================
# CH3I emitted over the oceans at rate of 1 molec/cm2/s
#==============================================================================
0 SRC_2D_CH3I 1.0 - - - xy molec/cm2/s TR_CH3I 1000 1 1
#==============================================================================
# Anthropogenic CO 25-day and 50-day tracers
#==============================================================================
0 CEDS_CO_25_ENE  $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_ene  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_ENE  -                                                          -      -                  - -  -      TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_IND  $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_ind  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_IND  -                                                          -      -                  - -  -      TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_TRA  $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_tra  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_TRA  -                                                          -      -                  - -  -      TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_RCO  $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_rco  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_RCO  -                                                          -      -                  - -  -      TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_SLV  $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_slv  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_SLV  -                                                          -      -                  - -  -      TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_WST  $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_wst  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_WST  -                                                          -      -                  - -  -      TR_CO_50  26 1/2/12 5
#==============================================================================
# e90 tracer is emitted with a uniform flux of 100 ppbv
#==============================================================================
0 SRC_2D_e90  100 - - - xy ppbv TR_e90  -    1 1
0 SRC_2D_e90n  100 - - - xy ppbv TR_e90_N 1001 1 1
0 SRC_2D_e90s  100 - - - xy ppbv TR_e90_S 1002 1 1
#==============================================================================
# Uniform surface emission of 2400 Tg/yr
#==============================================================================
0 SRC_2D_Uniform 1.5e-7 - - - xy kg/m2/s TR_Uniform - 1 1
...
1000 OCEAN_MASK $METDIR/$CNYR/01/$MET.$CNYR0101.CN.$RES.$NC FROCEAN 2000/1/1/0 C xy 1 1 -180/-90/180/90
1001 NH_MASK    -180/0/180/90  - 2000/1/1/0 C xy 1 1 -180/0/180/90
1002 SH_MASK    -180/-90/180/0 - 2000/1/1/0 C xy 1 1 -180/-90/180/0
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 13:12, 4 January 2019 (UTC)

Revision as of 13:12, 4 January 2019

All users interested in the GEOS-Chem transport scheme are encouraged to subscribe to the transport email list (click on the link in the contact information section section below).

Contact information

Transport Working Group Co-Chairs
Transport Working Group email list geos-chem-transport [at] g.harvard.edu
To subscribe to email list Either
  • Send an email to geos-chem-transport+subscribe [at] g.harvard.edu

Or

To unsubscribe from email list Either
  • Send an email to geos-chem-transport+unsubscribe [at] g.harvard.edu

Or

--Bob Y. (talk) 18:38, 21 August 2015 (UTC)

Current GEOS-Chem Transport Projects (please add yours!)

User Group Description Contact Person
U. Toronto Assessing impact of transport biases on methane simulation and applying weak constraint 4D-Var to correct them Ilya Stanevich
U. Toronto Characterizing model transport errors using weak constraint 4D-Var Martin Keller
U. Toronto Quantifying the impact of errors in PBL mixing on CO2 Wei Lu
Harvard Examining the effect of grid resolution on vertical transport and convective mass fluxes Karen Yu
Harvard Quantifying the impact of numerical diffusion on plume transport and comparing to observed atmospheric diffusion Sebastian D. Eastham
Harvard Assessing the co-benefits of increasing vertical and horizontal grid resolutions for tracer transport Jiawei Zhuang
Dalhousie Implementation of a new PBL mixing scheme Aaron van Donkelaar
NIA / NASA Langley Constraints on Precipitation Scavenging and Transport of Soluble Trace Species From Airborne and Surface Radionuclide (222Rn, 210Pb, 7Be, 10Be) Observations Bo Zhang Hongyu Liu
Harvard Vertical resolution needs for GEOS-Chem global transport Ada Shaw
Harvard
U. Toronto 		Re-diagnosed convection in off-line GEOS-Chem Bob Yantosca
Tailong He
Add yours here!    

Current GEOS-Chem Transport Issues

Text to be added

Common set of tracers in GEOS-Chem and GEOS

This update will be included in GEOS-Chem 12.2.0.

The Transport Working Group has proposed implementing additional tracers in GEOS-Chem for benchmarking purposes and diagnosing transport in GEOS-Chem compared to GEOS-5. The list of tracers in GEOS is provided below and is implemented via the tracer gridded component "TR_GridComp" in GEOS. In GEOS-Chem, the tracer suite has been added as an extension of the Rn-Pb-Be simulation (simulation type "1") by default, but may be added to any GEOS-Chem simulations by following the instructions below.

List of tracers in GEOS

17 GMI tracers: https://gmi.gsfc.nasa.gov/uploads/files/gmi_tracersuite.pdf

Additional Tracers in GEOS - used in part to get an age spectrum 

# Name     Units        Long Name
# -----    ------       --------------------------------
nh_5      'mol mol-1'   Northern Hemisphere 5 day tracer
nh_50     'mol mol-1'   Northern Hemisphere 50 day tracer
sh_5      'mol mol-1'   Southern Hemisphere 5 day tracer
sh_50     'mol mol-1'   Southern Hemisphere 50 day tracer
ntr_5     'mol mol-1'   Northern Hemisphere Subtropical 5 day tracer
ntr_50    'mol mol-1'   Northern Hemisphere Subtropical 50 day tracer
str_5     'mol mol-1'   Southern Hemisphere Subtropical 5 day tracer
str_50    'mol mol-1'   Southern Hemisphere Subtropical 50 day tracer
tr_5      'mol mol-1'   Tropical 5 day tracer
tr_50     'mol mol-1'   Tropical 50 day tracer
e90       'mol mol-1'   Constant emission 90 day tracer
e90_n     'mol mol-1'   Constant emission Northern Hemisphere 90 day tracer
e90_s     'mol mol-1'   Constant emission Southern Hemisphere 90 day tracer
st80_25   'mol mol-1'   Stratospheric source 25 day tracer
aoa_nh    'days'        Age of air tracer (Northern Hemisphere surface source)
::

and a few others 

CH3I      'mol mol-1'   Methyl iodide
CO_50_na  'mol mol-1'   Anthro CO North America 50 day tracer

--Lizzie Lundgren (talk) 15:41, 15 June 2018 (UTC)

Origin tracers

Chris Holmes wrote:

Another tracer suggestion: I like and use the "origin" tracers defined by Clara Orbe et al. (2013). These tracers give the fraction of air that that originates in specified source regions.
In addition to diagnosing transport, I think the tracers will be useful as templates for users to modify for their own purposes (e.g. customizing source regions). For that reason, a comprehensive set of tracers would be useful, if they can be optional.
Reference
Orbe, C., M. Holzer, L. M. Polvani, and D. Waugh (2013), Air-mass origin as a diagnostic of tropospheric transport, J. Geophys. Res. Atmos., 118, 1459–1470, doi:10.1002/jgrd.50133.

--Melissa Sulprizio (talk) 18:42, 23 February 2018 (UTC)

Implementation in GEOS-Chem

Several of the GEOS tracers listed above have been added as passive species in GEOS-Chem. The new tracers are listed in the advected species menu of input.geos as: 

%%% ADVECTED SPECIES MENU %%%:
Type of simulation      : 1
Species Entries ------->: Name
Species name            : Rn222
Species name            : Pb210
Species name            : Pb210Strat
Species name            : Be7
Species name            : Be7Strat
Species name            : Be10
Species name            : Be10Strat
Species name            : Passive
Species name            : TR_SF6
Species name            : TR_CH3I
Species name            : TR_CO_25
Species name            : TR_CO_50
Species name            : TR_e90
Species name            : TR_e90_N
Species name            : TR_e90_S
Species name            : TR_Uniform
Species name            : TR_Clock

where TR_ denotes that these species are from TR_GridComp of GEOS. In the passive species menu of input.geos, each species is defined along with the molecular weight, lifetime, default background concentration, and long name. 

%%% PASSIVE SPECIES MENU %%%:
Passive species name    : Passive    1.0    -1      1.0e-7   Passive_tracer_for_mass_conservation_eval
Passive species name    : TR_SF6     146.0  1.0e11  1.0e-20  Sulfur_hexafluoride
Passive species name    : TR_CH3I    142.0  4.32e5  1.0e-20  Methyl_iodide
Passive species name    : TR_CO_25   28.0   2.16e6  1.0e-20  Anthro_CO_25day_tracer
Passive species name    : TR_CO_50   28.0   4.32e6  1.0e-20  Anthro_CO_50day_tracer
Passive species name    : TR_e90     1.0    7.776e6 1.0e-20  Constant_emission_tracer
Passive species name    : TR_e90_N   1.0    7.776e6 1.0e-20  Constant_emission_tracer_for_NH
Passive species name    : TR_e90_S   1.0    7.776e6 1.0e-20  Constant_emission_tracer_for_SH
Passive species name    : TR_Uniform 28.0   2.16e6  1.0e-20  CO-like_tracer_with_uniform_emission
Passive species name    : TR_Clock   1.0    -1      0.0      Age_of_air_in_seconds

Species with an emissions flux also have a corresponding entry in the HEMCO_Config.rc file. 

#==============================================================================
# SF6 emissions from EDGAR v4.2
#==============================================================================
0 SRC_2D_SF6  /n/regal/jacob_lab/msulprizio/v42_SF6_2008-1970_IPCC_2_nc/v2019-01/EDGAR_v42_SF6_IPCC_2.generic.01x01.nc emi_sf6 1970-2008/1-12/1/0 C xy kg/m2/s TR_SF6 - 1 1

#==============================================================================
# CH3I emitted over the oceans at rate of 1 molec/cm2/s
#==============================================================================
0 SRC_2D_CH3I 1.0 - - - xy molec/cm2/s TR_CH3I 1000 1 1

#==============================================================================
# Anthropogenic CO 25-day and 50-day tracers
#==============================================================================
0 CEDS_CO_25_ENE   $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_ene  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_ENE   -                                                          -       -                  - -  -       TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_IND   $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_ind  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_IND   -                                                          -       -                  - -  -       TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_TRA   $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_tra  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_TRA   -                                                          -       -                  - -  -       TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_RCO   $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_rco  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_RCO   -                                                          -       -                  - -  -       TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_SLV   $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_slv  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_SLV   -                                                          -       -                  - -  -       TR_CO_50  26 1/2/12 5
0 CEDS_CO_25_WST   $ROOT/CEDS/v2018-08/$YYYY/CO-em-anthro_CMIP_CEDS_$YYYY.nc  CO_wst  1950-2014/1-12/1/0 C xy kg/m2/s TR_CO_25  26 1/2/12 5
0 CEDS_CO_50_WST   -                                                          -       -                  - -  -       TR_CO_50  26 1/2/12 5

#==============================================================================
# e90 tracer is emitted with a uniform flux of 100 ppbv
#==============================================================================
0 SRC_2D_e90   100 - - - xy ppbv TR_e90   -    1 1
0 SRC_2D_e90n  100 - - - xy ppbv TR_e90_N 1001 1 1
0 SRC_2D_e90s  100 - - - xy ppbv TR_e90_S 1002 1 1

#==============================================================================
# Uniform surface emission of 2400 Tg/yr
#==============================================================================
0 SRC_2D_Uniform 1.5e-7 - - - xy kg/m2/s TR_Uniform - 1 1

...

1000 OCEAN_MASK $METDIR/$CNYR/01/$MET.$CNYR0101.CN.$RES.$NC FROCEAN 2000/1/1/0 C xy 1 1 -180/-90/180/90
1001 NH_MASK    -180/0/180/90  - 2000/1/1/0 C xy 1 1 -180/0/180/90
1002 SH_MASK    -180/-90/180/0 - 2000/1/1/0 C xy 1 1 -180/-90/180/0

--Melissa Sulprizio (talk) 13:12, 4 January 2019 (UTC)

Retrieved from "https://wiki.seas.harvard.edu/geos-chem/index.php?title=Transport_Working_Group&oldid=40154"