Characteristics of Atmospheric Transport Using Three Numerical Formulations for Atmospheric Dynamics in a Single GCM Framework

Philip J. Rasch NCAR,* Boulder, Colorado

Search for other papers by Philip J. Rasch in
Current site
Google Scholar
PubMed
Close
,
Danielle B. Coleman NCAR,* Boulder, Colorado

Search for other papers by Danielle B. Coleman in
Current site
Google Scholar
PubMed
Close
,
Natalie Mahowald NCAR,* Boulder, Colorado

Search for other papers by Natalie Mahowald in
Current site
Google Scholar
PubMed
Close
,
David L. Williamson NCAR,* Boulder, Colorado

Search for other papers by David L. Williamson in
Current site
Google Scholar
PubMed
Close
,
Shian-Jiann Lin Geophysical Fluid Dynamics Laboratory, Princeton University, Princeton, New Jersey

Search for other papers by Shian-Jiann Lin in
Current site
Google Scholar
PubMed
Close
,
Byron A. Boville NCAR,* Boulder, Colorado

Search for other papers by Byron A. Boville in
Current site
Google Scholar
PubMed
Close
, and
Peter Hess NCAR,* Boulder, Colorado

Search for other papers by Peter Hess in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

This study examines the sensitivity of a number of important archetypical tracer problems to the numerical method used to solve the equations of tracer transport and atmospheric dynamics. The tracers' scenarios were constructed to exercise the model for a variety of problems relevant to understanding and modeling the physical, dynamical, and chemical aspects of the climate system. The use of spectral, semi-Lagrangian, and finite volume (FV) numerical methods for the equations is explored. All subgrid-scale physical parameterizations were the same in all model simulations.

The model behavior with a few short simulations with passive tracers is explored, and with much longer simulations of radon, SF6, ozone, a tracer designed to mimic some aspects of a biospheric source/sink of CO2, and a suite of tracers designed around the conservation laws for thermodynamics and mass in the model.

Large differences were seen near the tropopause in the model, where the FV core shows a much reduced level of vertical and meridional mixing. There was also evidence that the subtropical subsidence regions are more isolated from Tropics and midlatitudes in the FV core than seen in the other model simulations. There are also big differences in the stratosphere, particularly for age of air in the stratosphere and ozone. A comparison with estimated age of air from CO2 and SF6 measurements in the stratosphere suggest that the FV core is behaving most realistically.

A neutral biosphere (NB) test case is used to explore issues of diurnal and seasonal rectification of a tracer with sources and sinks at the surface. The sources and sinks have a zero annual average, and the rectification is associated with temporal correlations between the sources and sinks, and transport. The test suggests that the rectification is strongly influenced by the resolved-scale dynamics (i.e., the dynamical core) and that the numerical formulation for dynamics and transport still plays a critical role in the distribution of NB-like species. Since the distribution of species driven by these processes have a strong influence on the interpretation of the “missing sink” for CO2 and the interpretation of climate change associated with anthropogenic forcing herein, these issues should not be neglected.

The spectral core showed the largest departures from the predicted nonlinear relationship required by the equations for thermodynamics and mass conservations. The FV and semi-Lagrangian dynamics (SLD) models both produced errors a factor of 2 lower. The SLD model shows a small but systematic bias in its ability to maintain this relationship that was not present in the FV simulation.

The results of the study indicate that for virtually all of these problems, the model numerics still have a large role in influencing the model solutions. It was frequently the case that the differences in solutions resulting from varying the numerics still exceed the differences in the simulations resulting from significant physical perturbations (like changes in greenhouse gas forcing). This does not mean that the response of the system to physical changes is not correct. When results are consistent using different numerical formulations for dynamics and transport it lends confidence to one's conclusions, but it does indicate that some caution is required in interpreting the results.

The results from this study favor use of the FV core for tracer transport and model dynamics. The FV core is, unlike the others, conservative, less diffusive (e.g., maintains strong gradients better), and maintains the nonlinear relationships among variables required by thermodynamic and mass conservation constraints more accurately.

* The National Center for Atmospheric Research is sponsored by the National Science Foundation

Corresponding author address: Dr. Philip J. Rasch, NCAR, P.O. Box 3000, Boulder, CO 80307-3000. Email: pjr@ncar.ucar.edu

Abstract

This study examines the sensitivity of a number of important archetypical tracer problems to the numerical method used to solve the equations of tracer transport and atmospheric dynamics. The tracers' scenarios were constructed to exercise the model for a variety of problems relevant to understanding and modeling the physical, dynamical, and chemical aspects of the climate system. The use of spectral, semi-Lagrangian, and finite volume (FV) numerical methods for the equations is explored. All subgrid-scale physical parameterizations were the same in all model simulations.

The model behavior with a few short simulations with passive tracers is explored, and with much longer simulations of radon, SF6, ozone, a tracer designed to mimic some aspects of a biospheric source/sink of CO2, and a suite of tracers designed around the conservation laws for thermodynamics and mass in the model.

Large differences were seen near the tropopause in the model, where the FV core shows a much reduced level of vertical and meridional mixing. There was also evidence that the subtropical subsidence regions are more isolated from Tropics and midlatitudes in the FV core than seen in the other model simulations. There are also big differences in the stratosphere, particularly for age of air in the stratosphere and ozone. A comparison with estimated age of air from CO2 and SF6 measurements in the stratosphere suggest that the FV core is behaving most realistically.

A neutral biosphere (NB) test case is used to explore issues of diurnal and seasonal rectification of a tracer with sources and sinks at the surface. The sources and sinks have a zero annual average, and the rectification is associated with temporal correlations between the sources and sinks, and transport. The test suggests that the rectification is strongly influenced by the resolved-scale dynamics (i.e., the dynamical core) and that the numerical formulation for dynamics and transport still plays a critical role in the distribution of NB-like species. Since the distribution of species driven by these processes have a strong influence on the interpretation of the “missing sink” for CO2 and the interpretation of climate change associated with anthropogenic forcing herein, these issues should not be neglected.

The spectral core showed the largest departures from the predicted nonlinear relationship required by the equations for thermodynamics and mass conservations. The FV and semi-Lagrangian dynamics (SLD) models both produced errors a factor of 2 lower. The SLD model shows a small but systematic bias in its ability to maintain this relationship that was not present in the FV simulation.

The results of the study indicate that for virtually all of these problems, the model numerics still have a large role in influencing the model solutions. It was frequently the case that the differences in solutions resulting from varying the numerics still exceed the differences in the simulations resulting from significant physical perturbations (like changes in greenhouse gas forcing). This does not mean that the response of the system to physical changes is not correct. When results are consistent using different numerical formulations for dynamics and transport it lends confidence to one's conclusions, but it does indicate that some caution is required in interpreting the results.

The results from this study favor use of the FV core for tracer transport and model dynamics. The FV core is, unlike the others, conservative, less diffusive (e.g., maintains strong gradients better), and maintains the nonlinear relationships among variables required by thermodynamic and mass conservation constraints more accurately.

* The National Center for Atmospheric Research is sponsored by the National Science Foundation

Corresponding author address: Dr. Philip J. Rasch, NCAR, P.O. Box 3000, Boulder, CO 80307-3000. Email: pjr@ncar.ucar.edu

Save
  • Andrews, A. E., and Coauthors, 2001: Mean ages of stratospheric air derived from in situ observations of CO2, CH4, and N2O. J. Geophys. Res, 106 , 3229532314.

    • Search Google Scholar
    • Export Citation
  • Babiker, M., J. Reilly, M. Mayer, R. Eckaus, I. S. Wing, and R. Hyman, 2001: The MIT Emissions Prediction and Policy Analysis (EPPA) model: Revisions, sensitivities, and comparison of results. MIT Joint Program on the Policy and Science of Global Change Tech. Rep. 71, Massachusetts Institute of Technology, Cambridge, MA, 90 pp.

  • Bey, I., and Coauthors, 2001: Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation. J. Geophys. Res, 106 , 2307323096.

    • Search Google Scholar
    • Export Citation
  • Boering, K. A., S. C. Wofsy, B. C. Daube, H. R. Schneider, M. Loewenstein, J. R. Podolske, and T. J. Conway, 1996: Stratospheric transport rates and mean age distribution derived from observations of atmospheric CO2 and N2O. Science, 274 , 13401343.

    • Search Google Scholar
    • Export Citation
  • Boville, B. A., and X. Cheng, 1988: Upper boundary effects in a general circulation model. J. Atmos. Sci, 45 , 25922606.

  • Boville, B. A., and C. S. Bretherton, 2003: Heating and kinetic energy dissipation in the NCAR Community Atmosphere Model. J. Climate, 16 , 38773887.

    • Search Google Scholar
    • Export Citation
  • Boville, B. A., D. Jacob, and M. Prather, 1997: Global tracer transport models, proceedings of a WCRP workshop on the parameterization of subgrid-scale tracer transport. Tech. Doc., WMO/TD 823 26, World Meteorological Organization, 100 pp.

  • Boville, B. A., P. J. Rasch, J. J. Hack, and J. R. McCaa, 2006: Representation of clouds and precipitation processes in the Community Atmosphere Model version 3 (CAM3). J. Climate, 19 , 21842198.

    • Search Google Scholar
    • Export Citation
  • Butchart, N., and A. A. Sciafe, 2001: Removal of chlorofluorocarbons by increased mass exchange between the stratosphere and troposphere in a changing climate. Nature, 410 , 799802.

    • Search Google Scholar
    • Export Citation
  • Collins, W. D., and Coauthors, 2004: Description of the NCAR Community Atmosphere Model: CAM3.0. Tech. Rep. NCAR/TN-464+STR, National Center for Atmospheric Research, Boulder, CO, 226 pp. [Available online at http://www.ccsm.ucar.edu/models/atm-cam.].

  • Collins, W. D., and Coauthors, 2006a: The Community Climate System Model, version 3 (CCSM3). J. Climate, 19 , 21222143.

  • Collins, W. D., and Coauthors, 2006b: The formulation and atmospheric simulation of the Community Atmosphere Model version 3: CAM3. J. Climate, 19 , 21402157.

    • Search Google Scholar
    • Export Citation
  • Collins, W. J., R. G. Derwent, B. Garnier, C. E. Johnson, M. G. Sanderson, and S. S. Stevenson, 2003: Effect of stratosphere-troposphere exchange on the future tropospheric ozone trend. J. Geophys. Res, 108 .8528, doi:10.1029/2002JD002617.

    • Search Google Scholar
    • Export Citation
  • Crutzen, P. J., M. G. Lawrence, and U. Poschl, 1999: On the background photochemistry of tropospheric ozone. Tellus, 51A , 123146.

  • Denning, A., I. Y. Fung, and D. Randall, 1995: Latitudinal gradient of atmospheric CO2 due to seasonal exchange with land biota. Nature, 376 , 240243.

    • Search Google Scholar
    • Export Citation
  • Denning, A., G. Collatz, C. Zhang, D. Randall, J. A. Berry, P. J. Sellers, G. A. Colello, and D. A. Daazlich, 1996a: Simulations of terrestrial carbon metabolism and atmospheric CO2 in a general circulation model. Part 1: Surface carbon fluxes. Tellus, 48B , 521542.

    • Search Google Scholar
    • Export Citation
  • Denning, A., D. Randall, G. Collatz, and P. Sellers, 1996b: Simulations of terrestrial carbon metabolism and atmospheric CO2 in a general circulation model. Part 2: Spatial and temporal variations of atmospheric CO2. Tellus, 48B , 543567.

    • Search Google Scholar
    • Export Citation
  • Denning, A., and Coauthors, 1999: Three-dimensional transport and concentration of SF6. A model intercomparison study (TransCom 2). Tellus, 51B , 266297.

    • Search Google Scholar
    • Export Citation
  • Douglass, A. R., R. S. Stolarski, S. E. Strahan, and P. S. Connell, 2004: Radicals and reservoirs in the GMI chemistry and transport model: Comparison to measurements. J. Geophys. Res, 109 .D16302, doi:10.1029/2004JD004632.

    • Search Google Scholar
    • Export Citation
  • Dutton, J. A., 1976: The Ceaseless Wind. McGraw Hill, 579 pp.

  • Elkins, J. W., and Coauthors, 1996: Airborne gas chromatograph for in situ measurement of long-lived species in the upper troposphere and lower stratosphere. Geophys. Res. Lett, 23 , 347350.

    • Search Google Scholar
    • Export Citation
  • Eluszkiewicz, J., R. S. Hemler, J. D. Mahlman, L. Bruhwiler, and L. L. Takacs, 2000: Sensitivity of age-of-air calculations to the choice of advection scheme. J. Atmos. Sci, 57 , 31853201.

    • Search Google Scholar
    • Export Citation
  • Feichter, J., and P. J. Crutzen, 1990: Parameterization of vertical tracer transport due to deep cumulus convection in a global transport model and its evaluation with 222Radon measurements. Tellus, 42B , 100117.

    • Search Google Scholar
    • Export Citation
  • Genthon, C., and A. Armengaud, 1995: Radon 222 as a comparative tracer of transport and mixing in two general circulation models the atmosphere. J. Geophys. Res, 100 , 28492866.

    • Search Google Scholar
    • Export Citation
  • Gurney, K. R., and Coauthors, 2002: Towards robust regional estimates of CO2 sources and sinks using atmospheric transport models. Nature, 415 , 626630.

    • Search Google Scholar
    • Export Citation
  • Gurney, K. R., and Coauthors, 2003: TransCom 3 CO2 inversion intercomparison: 1. Annual mean control results and sensitivity to transport and prior flux information. Tellus, 55B , 555579.

    • Search Google Scholar
    • Export Citation
  • Hall, T. M., and R. A. Plumb, 1994: Age as a diagnostic of stratospheric transport. J. Geophys. Res, 99 , 10591070.

  • Hall, T. M., K. Boering, and R. A. Plumb, 1999: Evaluation of transport in stratospheric models. J. Geophys. Res, 104 , 1881518839.

  • Harnisch, J. R., R. Borchers, P. Fabian, and M. Maiss, 1996: Tropospheric trends for CF4 and C2F6 since 1982 derived from SF6 dated stratospheric air. Geophys. Res. Lett, 23 , 10991102.

    • Search Google Scholar
    • Export Citation
  • Jacob, D. J., and M. J. Prather, 1990: Radon-222 as a test of convective transport in a general circulation model. Tellus, 42B , 118134.

    • Search Google Scholar
    • Export Citation
  • Jacob, D. J., and Coauthors, 1997: Intercomparison of convective and synoptic transport in global models using 222Rn and other tracers. J. Geophys. Res, 102 , 59535970.

    • Search Google Scholar
    • Export Citation
  • Johnson, D. R., A. J. Lenzen, T. H. Zapotocny, and T. K. Schaack, 2000: Numerical uncertainties in simulation of reversible isentropic processes and entropy conservation. J. Climate, 13 , 38603884.

    • Search Google Scholar
    • Export Citation
  • Johnson, D. R., A. J. Lenzen, T. H. Zapotocny, and T. K. Schaack, 2002: Numerical uncertainties in simulation of reversible isentropic processes and entropy conservation: Part II. J. Climate, 15 , 17771804.

    • Search Google Scholar
    • Export Citation
  • Kiehl, J. T., C. Shields, J. J. Hack, and W. D. Collins, 2006: The climate sensitivity of the Community Climate System Model version 3 (CCSM3). J. Climate, 19 , 25842596.

    • Search Google Scholar
    • Export Citation
  • Lin, S-J., 2004: A “vertically Lagrangian” finite-volume dynamical core for global atmospheric models. Mon. Wea. Rev, 132 , 22932307.

    • Search Google Scholar
    • Export Citation
  • Liu, S., J. McAfee, and R. Cicerone, 1984: Radon222 and tropospheric vertical transport. J. Geophys. Res, 89 , 72917297.

  • McLinden, C. A., S. C. Olsen, B. Hannegan, O. Wild, and M. J. Prather, 2000: Stratospheric ozone in 3-d models: A simple chemistry and the cross-tropopause flux. J. Geophys. Res, 105 , 1465314665.

    • Search Google Scholar
    • Export Citation
  • Olsen, S., and J. T. Randerson, 2004: Differences between surface and column atmospheric CO2 and implications for carbon cycle research. J. Geophys. Res, 109 .D02301, doi:10.1029/2003JD003968.

    • Search Google Scholar
    • Export Citation
  • Otto-Bliesner, B. L., R. A. Tomas, E. C. Brady, C. M. Ammann, and Z. Kothavala, 2006: Climate sensitivity of moderate- and low-resolution versions of CCSM3 to preindustrial forcings. J. Climate, 19 , 25672583.

    • Search Google Scholar
    • Export Citation
  • Pearman, G. I., and P. Hyson, 1980: Activities of the global biosphere as reflected in atmospheric CO2 records. J. Geophys. Res, 85 , 44684474.

    • Search Google Scholar
    • Export Citation
  • Potter, C. S., J. T. Randerson, C. B. Field, P. A. Matson, P. M. Vitousek, H. A. Mooney, and S. A. Klooster, 1993: Terrestrial ecosystem production: A process model based on global satellite and surface data. Global Biogeochem. Cycles, 8 , 811841.

    • Search Google Scholar
    • Export Citation
  • Prather, M. J., and E. E. Remsberg, 1992: The atmospheric effects of stratospheric aircraft: Report of the 1992 models and measurements workshop. NASA Tech. Rep. 1292, 3 vols., 764 pp.

  • Pyle, J., and M. Prather, 1996: Global tracer transport models, report of a scientific symposium, Bermuda, 1990. WMO Tech. Doc. 24, World Meteorological Organization, 186 pp.

  • Randerson, J. T., M. V. Thompson, T. J. Conway, I. Y. Fung, and C. B. Field, 1997: The contribution of terrestrial sources and sinks to trends in the seasonal cycle of atmospheric carbon dioxide. Global Biogeochem. Cycles, 11 , 535560.

    • Search Google Scholar
    • Export Citation
  • Rasch, P. J., 1986: Toward atmospheres without tops: Absorbing upper boundary conditions for numerical models. Quart. J. Roy. Meteor. Soc, 112 , 11951218.

    • Search Google Scholar
    • Export Citation
  • Rasch, P. J., and D. L. Williamson, 1990: Computational aspects of moisture transport in global models of the atmosphere. Quart. J. Roy. Meteor. Soc, 116 , 10711090.

    • Search Google Scholar
    • Export Citation
  • Rasch, P. J., and D. L. Williamson, 1991: Sensitivity of a general circulation model climate to the moisture transport formulation. J. Geophys. Res, 96 , 1312313137.

    • Search Google Scholar
    • Export Citation
  • Rasch, P. J., B. A. Boville, and G. P. Brasseur, 1995: A three-dimensional general circulation model with coupled chemistry for the middle atmosphere. J. Geophys. Res, 100 , 90419071.

    • Search Google Scholar
    • Export Citation
  • Rasch, P. J., and Coauthors, 2000: A comparison of scavenging and deposition processes in global models; Results from the WCRP Cambridge workshop of 1995. Tellus, 52B , 10251056.

    • Search Google Scholar
    • Export Citation
  • Rasch, P., and Coauthors, 2006: Characterization of tropical transient activity in the CAM3 atmospheric hydrologic cycle. J. Climate, 19 , 22222242.

    • Search Google Scholar
    • Export Citation
  • Ray, E. A., F. L. Moore, J. W. Elkins, G. S. Dutton, D. W. Fahey, H. Vmel, S. J. Oltmans, and K. H. Rosenlof, 1999: Transport into the northern hemisphere lowermost stratosphere revealed by in situ tracer measurements. J. Geophys. Res, 104 , 2656526580.

    • Search Google Scholar
    • Export Citation
  • Rind, D., and J. Lerner, 1996: The use of online tracers as a diagnostic tool in GCM model development. J. Geophys. Res, 101 , 1266712683.

    • Search Google Scholar
    • Export Citation
  • Rind, D., J. Lerner, K. Shah, and R. Suozzo, 1999: Use of online tracers as a diagnostic tool in general circulation model development, 2. transport between the troposphere and stratosphere. J. Geophys. Res, 104 , 91519167.

    • Search Google Scholar
    • Export Citation
  • Rind, D., J. Lerner, and C. McLinden, 2001: Changes of tracer distributions in the doubled CO2 climate. J. Geophys. Res, 106 , 2806128079.

    • Search Google Scholar
    • Export Citation
  • Waugh, D. W., and T. M. Hall, 2002: Age of stratospheric air: Theory, observations and models. Rev. Geophys, 40 .1010, doi:10.1029/2000RG000101.

    • Search Google Scholar
    • Export Citation
  • Williamson, D. L., 1983: Description of NCAR Community Climate Model (CCM0B). Tech. Rep. NCAR/TN-210+STR, National Center for Atmospheric Research, Boulder, CO, 88 pp.

  • Williamson, D. L., and J. G. Olson, 1994: Climate simulations with a semi-Lagrangian version of the NCAR Community Climate Model. Mon. Wea. Rev, 122 , 15941610.

    • Search Google Scholar
    • Export Citation
  • Williamson, D. L., and P. J. Rasch, 1994: Water vapor transport in the NCAR CCM2. Tellus, 46A , 3451.

  • Williamson, D. L., and J. G. Olson, 1998: A comparison of semi-Lagrangian and Eulerian polar climate simulations. Mon. Wea. Rev, 126 , 9911000.

    • Search Google Scholar
    • Export Citation
  • Zeng, G., and J. A. Pyle, 2003: Changes in tropospheric ozone between 2000 and 2100 modeled in a chemistry-climate model. Geophys. Res. Lett, 30 .1392, doi:10.1029/2002GL016708.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 459 92 13
PDF Downloads 231 58 3