Geoscientific Model Development www.geosci-model-dev.net 1991-959X 1991-9603 3 2 2010 10.5194/gmd-3-337-2010 http://www.geosci-model-dev.net/3/337/2010/ http://www.geosci-model-dev.net/3/337/2010/gmd-3-337-2010.html http://www.geosci-model-dev.net/3/337/2010/gmd-3-337-2010.pdf 337 364 2010-08-10 A kinetic chemistry tagging technique and its application to modelling the stable isotopic composition of atmospheric trace gases S. Gromov sergey.gromov@mpic.de P. Jöckel R. Sander C. A. M. Brenninkmeijer Atmospheric Chemistry Department, Max Planck Institute for Chemistry, P.O. Box 3060, 55020 Mainz, Germany now at: Deutsches Zentrum für Luft-und Raumfahrt (DLR), Institut für Physik der Atmosphäre, Oberpfaffenhofen, 82234 Weßling, Germany Isotope composition, in many cases, holds unique information on the sources, chemical modification and sinks of atmospheric trace gases. Vital to the interpretation and use of an increasing number of isotope analyses is appropriate modelling. However, the exact implementation of isotopic information in chemistry-climate models is a challenge, and often studies use simplifications which limit their applicability. Here we implement a thorough isotopic extension in MECCA, a comprehensive kinetic chemistry sub-model. To this end, we devise a generic tagging technique for the kinetic chemistry mechanisms implemented as the sub-submodel MECCA-TAG. The technique is diagnostic and numerically efficient and supports the investigation of various aspects of kinetic chemistry schemes. We focus specifically on the application to the modelling of stable isotopic composition. The results of MECCA-TAG are evaluated against the reference sub-submodel MECCA-DBL, which is implicitly full-detailed, but computationally expensive and thus sub-optimal in practical applications. Furthermore, we evaluate the elaborate carbon and oxygen isotopic mechanism by simulating the multi-isotope composition of CO and other trace gases in the CAABA/MECCA box-model. The mechanism realistically simulates the oxygen isotope composition of key species, as well as the carbon isotope signature transfer. The model adequately reproduces the isotope chemistry features for CO, taking into account the limits of the modelling domain. In particular, the mass-independently fractionated (MIF) composition of CO due to reactions of ozone with unsaturated hydrocarbons (a source effect) versus its intrinsic MIF enrichment induced in the removal reaction via oxidation by OH is assessed. The simulated ozone source effect was up to +1‰ in Δ¹⁷O(CO). The versatile modelling framework we employ (the Modular Earth Submodel System, MESSy) opens the way for implementation of the novel detailed isotopic chemistry treatment in the three-dimensional atmospheric-chemistry general circulation model EMAC. We therefore also present estimates of the computational gain obtained by the developed optimisations. Allan, W., Manning, M. R., Lassey, K. R., Lowe, D. C., and Gomez, A. J.: Modeling the variation of delta C-13 in atmospheric methane: Phase ellipses and the kinetic isotope effect, Global Biogeochem. Cy., 15, 467–481, 2001. Barkan, E. and Luz, B.: High precision measurements of 17O/16O and 18O/16O ratios in H₂O, Rapid Commun. Mass Sp., 19, 3737–3742, 2005. Bergamaschi, P., Hein, R., Brenninkmeijer, C. A. M., and Crutzen, P. J.: Inverse modeling of the global CO cycle 2. Inversion of C-13/C-12 and O-18/O-16 isotope ratios, J. Geophys. Res.-Atmos., 105, 1929–1945, 2000. Bloch, R.: A Practical Method to Integrate Some Stiff Systems, Comput. Biomed. Res., 24, 420–428, 1991. Boering, K. A., Jackson, T., Hoag, K. J., Cole, A. S., Perri, M. J., Thiemens, M., and Atlas, E.: Observations of the anomalous oxygen isotopic composition of carbon dioxide in the lower stratosphere and the flux of the anomaly to the troposphere, Geophys. Res. Lett., 31, L03109, doi:10.1029/2003gl018451, 2004. Brenninkmeijer, C. A. M. and Röckmann, T.: Principal factors determining the O-18/O-16 ratio of atmospheric CO as derived from observations in the southern hemispheric troposphere and lowermost stratosphere, J. Geophys. Res.-Atmos., 102, 25477–25485, 1997. Brenninkmeijer, C. A. M., Röckmann, T., Bräunlich, M., Jöckel, P., and Bergamaschi, P.: Review of progress in isotope studies of atmospheric carbon monoxide, Chemosphere - Global Change Science, 1, 33–52, 10.1016/S1465-9972(99)00018-5, 1999. Brenninkmeijer, C. A. M., Janssen, C., Kaiser, J., Röckmann, T., Rhee, T. S., and Assonov, S. S.: Isotope effects in the chemistry of atmospheric trace compounds, Chem. Rev., 103, 5125–5161, doi:10.1021/Cr020644k, 2003. Coplen, T. B., Bohlke, J. K., De Bievre, P., Ding, T., Holden, N. E., Hopple, J. A., Krouse, H. R., Lamberty, A., Peiser, H. S., Revesz, K., Rieder, S. E., Rosman, K. J. R., Roth, E., Taylor, P. D. P., Vocke, R. D., and Xiao, Y. K.: Isotope-abundance variations of selected elements – (IUPAC Technical Report), Pure Appl. Chem., 74, 1987–2017, 2002. Cuntz, M., Ciais, P., Hoffmann, G., and Knorr, W.: A comprehensive global three-dimensional model of delta O-18 in atmospheric CO₂: 1. Validation of surface processes, J. Geophys. Res.-Atmos., 108, 4527, doi:10.1029/2002jd003153, 2003. Dubey, M. K., Mohrschladt, R., Donahue, N. M., and Anderson, J. G.: Isotope specific kinetics of hydroxyl radical (OH) with water (H₂O): Testing models of reactivity and atmospheric fractionation, J. Phys. Chem. A, 101, 1494–1500, 1997. Duncan, B. N., Logan, J. A., Bey, I., Megretskaia, I. A., Yantosca, R. M., Novelli, P. C., Jones, N. B., and Rinsland, C. P.: Global budget of CO, 1988–1997: Source estimates and validation with a global model, J. Geophys. Res.-Atmos., 112, D22301, doi:10.1029/2007jd008459, 2007. Eriksson, K., Johnson, C., and Logg, A.: Explicit time-stepping for stiff ODES, SIAM J. Sci. Stat. Comp., 25, 1142–1157, 2003. Gros, V., Röckmann, T., Jöckel, P., and Brenninkmeijer, C. A. M.: Monitoring of atmospheric carbon monoxide and its isotopic composition, International Symposium on Isotopomers, Yokohama, Japan, 2001. Gros, V., Jöckel, P., Brenninkmeijer, C. A. M., Röckmann, T., Meinhardt, F., and Graul, R.: Characterization of pollution events observed at Schauinsland, Germany, using CO and its stable isotopes, Atmos. Environ., 36, 2831–2840, Pii S1352-2310(02)00144-9, 2002. Jöckel, P., Sander, R., Kerkweg, A., Tost, H., and Lelieveld, J.: Technical Note: The Modular Earth Submodel System (MESSy) – a new approach towards Earth System Modeling, Atmos. Chem. Phys., 5, 433–444, doi:10.5194/acp-5-433-2005, 2005. Jöckel, P., Tost, H., Pozzer, A., Brühl, C., Buchholz, J., Ganzeveld, L., Hoor, P., Kerkweg, A., Lawrence, M. G., Sander, R., Steil, B., Stiller, G., Tanarhte, M., Taraborrelli, D., Van Aardenne, J., and Lelieveld, J.: The atmospheric chemistry general circulation model ECHAM5/MESSy1: consistent simulation of ozone from the surface to the mesosphere, Atmos. Chem. Phys., 6, 5067–5104, doi:10.5194/acp-6-5067-2006, 2006. Jöckel, P., Kerkweg, A., Buchholz-Dietsch, J., Tost, H., Sander, R., and Pozzer, A.: Technical Note: Coupling of chemical processes with the Modular Earth Submodel System (MESSy) submodel TRACER, Atmos. Chem. Phys., 8, 1677–1687, doi:10.5194/acp-8-1677-2008, 2008. Johnson, M. S., Feilberg, K. L., von Hessberg, P., and Nielsen, O. J.: Isotopic processes in atmospheric chemistry, Chem. Soc. Rev., 31, 313–323, doi:10.1039/B108011n, 2002. Johnston, J. C., Röckmann, T., and Brenninkmeijer, C. A. M.: CO2+O(1D) isotopic exchange: Laboratory and modeling studies, J. Geophys. Res.-Atmos., 105, 15213–15229, 2000. Kato, S., Akimoto, H., Röckmann, T., Bräunlich, M., and Brenninkmeijer, C. A. M.: Stable isotopic compositions of carbon monoxide from biomass burning experiments, Atmos. Environ., 33, 4357–4362, 1999. Kaye, J. A.: Mechanisms and Observations for Isotope Fractionation of Molecular-Species in Planetary-Atmospheres, Rev. Geophys., 25, 1609–1658, 1987. Kohen, A. and Limbach, H.-H.: Isotope Effects in Chemistry and Biology, CRC Press/Taylor and Francis Group, Boca Raton, FL 33487-2742, 1074~pp., 2006. Lassey, K. R., Lowe, D. C., Brenninkmeijer, C. A. M., and Gomez, A. J.: Atmospheric Methane and Its Carbon Isotopes in the Southern-Hemisphere – Their Time-Series and an Instructive Model, Chemosphere, 26, 95–109, 1993. Lehmann, R.: An algorithm for the determination of all significant pathways in chemical reaction systems, J. Atmos. Chem., 47, 45–78, 2004. Lelieveld, J. and Crutzen, P. J.: The Role of Clouds in Tropospheric Photochemistry, J. Atmos. Chem., 12, 229–267, 1991. Liang, M. C. and Yung, Y. L.: Sources of the oxygen isotopic anomaly in atmospheric N₂O, J. Geophys. Res.-Atmos., 112, D13307, doi:10.1029/2006jd007876, 2007. Liang, M. C., Blake, G. A., and Yung, Y. L.: Seasonal cycle of (COO)-O-16-O-16, (COO)-O-16-O-17, and (COO)-O-16-O-18 in the middle atmosphere: Implications for mesospheric dynamics and biogeochemical sources and sinks of CO₂, J. Geophys. Res.-Atmos., 113, D12305, doi:10.1029/2007jd008392, 2008. Lyons, J. R.: Transfer of mass-independent fractionation in ozone to other oxygen-containing radicals in the atmosphere, Geophys. Res. Lett., 28, 3231–3234, 2001. Manning, M. R., Brenninkmeijer, C. A. M., and Allan, W.: Atmospheric carbon monoxide budget of the southern hemisphere: Implications of C-13/C-12 measurements, J. Geophys. Res.-Atmos., 102, 10673–10682, 1997. McCarthy, M. C., Boering, K. A., Rice, A. L., Tyler, S. C., Connell, P., and Atlas, E.: Carbon and hydrogen isotopic compositions of stratospheric methane: 2. Two-dimensional model results and implications for kinetic isotope effects, J. Geophys. Res.-Atmos., 108(18), 4461 doi:10.1029/2002jd003183, 2003. McNaught, A. D. and Wilkinson, A.: Compendium of Chemical Terminology, 2nd edn., Blackwell Science, 1997. Pieterse, G., Krol, M. C., and Röckmann, T.: A consistent molecular hydrogen isotope chemistry scheme based on an independent bond approximation, Atmos. Chem. Phys., 9, 8503–8529, doi:10.5194/acp-9-8503-2009, 2009. Pozzer, A., Jöckel, P., and Van Aardenne, J.: The influence of the vertical distribution of emissions on tropospheric chemistry, Atmos. Chem. Phys., 9, 9417–9432, doi:10.5194/acp-9-9417-2009, 2009. Press, W. H., Teukolsky, S. A., Flannery, B. P., and Vetterling, W. T.: Numerical Recipes in FORTRAN: The Art of Scientific Computing, Cambridge University Press, 1992. Quay, P. D., King, S. L., Stutsman, J., Wilbur, D. O., Steele, L. P., Fung, I., Gammon, R. H., Brown, T. A., Farwell, G. W., Grootes, P. M., and Schmidt, F. H.: Carbon isotopic composition of atmospheric CH₄: Fossil and biomass burning source strengths, Global Biogeochem. Cy., 5, 25–47, 1991. Röckmann, T., Brenninkmeijer, C. A. M., Neeb, P., and Crutzen, P. J.: Ozonolysis of nonmethane hydrocarbons as a source of the observed mass independent oxygen isotope enrichment in tropospheric CO, J. Geophys. Res.-Atmos., 103, 1463–1470, 1998a. Röckmann, T., Brenninkmeijer, C. A. M., Saueressig, G., Bergamaschi, P., Crowley, J. N., Fischer, H., and Crutzen, P. J.: Mass-independent oxygen isotope fractionation in atmospheric CO as a result of the reaction CO+OH, Science, 281, 544–546, 1998b. Röckmann, T., Jöckel, P., Gros, V., Bräunlich, M., Possnert, G., and Brenninkmeijer, C. A. M.: Using $^14$C, $^13$C, $^18$O and $^17$O isotopic variations to provide insights into the high northern latitude surface CO inventory, Atmos. Chem. Phys., 2, 147–159, doi:10.5194/acp-2-147-2002, 2002. Sander, R., Kerkweg, A., Jockel, P., and Lelieveld, J.: Technical note: The new comprehensive atmospheric chemistry module MECCA, Atmos. Chem. Phys., 5, 445–450, doi:10.5194/acp-5-445-2005, 2005. Sander, R., Baumgaertner, A. J. G., Gromov, S., Harder, H., Jöckel, P., Kerkweg, A., Kubistin, D., Riede, H., Taraborrelli, D., and Tost, H.: The atmospheric chemistry box model CAABA/MECCA-3.0, Geosci. Model Dev., in preparation, 2010. Sandu, A. and Sander, R.: Technical note: Simulating chemical systems in Fortran90 and Matlab with the Kinetic PreProcessor KPP-2.1, Atmos. Chem. Phys., 6, 187–195, doi:10.5194/acp-6-187-2006, 2006. Saueressig, G., Crowley, J. N., Bergamaschi, P., Bruhl, C., Brenninkmeijer, C. A. M., and Fischer, H.: Carbon 13 and D kinetic isotope effects in the reactions of CH₄ with O(D-1) and OH: New laboratory measurements and their implications for the isotopic composition of stratospheric methane, J. Geophys. Res.-Atmos., 106, 23127–23138, 2001. Schauble, E. A.: Applying Stable Isotope Fractionation \mboxTheory to New Systems, Rev. Mineral. Geochem., 55, 65–111, doi:10.2138/gsrmg.55.1.65, 2004. Thiemens, M. H.: History and applications of mass-independent isotope effects, Annu. Rev. Earth Pl. Sc., 34, 217–262, doi:10.1146/annurev.earth.34.031405.125026, 2006. Weston, R. E.: Oxygen isotope effects in the oxidation of methane to carbon monoxide, J. Phys. Chem. A, 105, 1656–1661, 2001. Yung, Y. L., Demore, W. B., and Pinto, J. P.: Isotopic exchange between carbon-dioxide and ozone via O(1D) in the stratosphere, Geophys. Res. Lett., 18, 13–16, 1991. Yung, Y. L., Lee, A. Y. T., Irion, F. W., DeMore, W. B., and Wen, J.: Carbon dioxide in the atmosphere: Isotopic exchange with ozone and its use as a tracer in the middle atmosphere, J. Geophys. Res.-Atmos., 102, 10857–10866, 1997. Zahn, A., Franz, P., Bechtel, C., Grooß, J.-U., and Röckmann, T.: Modelling the budget of middle atmospheric water vapour isotopes, Atmos. Chem. Phys., 6, 2073–2090, doi:10.5194/acp-6-2073-2006, 2006.