GSFC Code 916: Atmospheric Chemistry and
Dynamics Branch

Parameterized Chemistry and Transport Model


New development | PCTM Home | Gallery | Experiment summaries | Archived datasets | Publications list | Contact


The parameterized chemistry and transport model (PCTM) is one of several numerical tools which the Atmospheric Chemistry and Dynamics Branch uses to study the evolution of trace compounds in the Earth's atmosphere. Please browse through these pages to become familiar with the PCTM and the scope of issues to which it is applied.


Welcome to the PCTM

The PCTM numerically solves the constituent continuity equation
Constituent continuity equation
over the entire globe. In Equation 1, x is the constituent mixing ratio, t is time, V is the velocity vector, and P and L(x) are, respectively, the production and loss rates for x.

Each step of the integration is accomplished with two processes. Advection is done first and produces an "intermediate" consituent configuration, x*. In prediction form, Equation 1 is written in terms of x* as

Constituent continuity equation in finite difference form
The PCTM uses the flux-form semi-Lagrangian advection scheme of Lin and Rood [1996] which is often applied with a Zelasek-type multidimensional flux correction.

Second, the parameterization of the photochemistry is applied to x* by

Prediction equation for consituent mixing ratio
L has dimensions of s-1 and is called the loss frequency. It arises from a Taylor's series expansion of L(x). The left hand side of Equation 2 is the predicted constituent mixing ratio which serves as x(t) for the next iteration.

The coefficients P and L are specified in a variety of ways. The simplest and most computationally efficient method involves reading tables of P and L where each is a function of, say, latitude, pressure, and time. For surface source gases, the surface exchange in kg m-2 s-1 is converted to ground-level production rates based on the volume of each respective surface layer grid box. For many short-lived tropospheric gases, the loss frequencies are determined by multiplying published reaction rates by the local number density of the particular reactant.

The subgrid-scale transport due to convection process and turbulence in the planetary boundary layer (PBL) is included in the PCTM. The parameterizations of these subgrid-scale processes are done in the GCMs. The output of the cloud mass flux and turbulence diffusion coefficients from the data assimilation models or GCMs are used by the CTM to compute the subgrid-scale transport of the tracers.

The PCTM is run in the so-called offline mode with the dynamics supplied by archived meteorological datasets. The winds, temperatures, surface pressures, and cloud mass flux and detrainment are most often provided by the Goddard Earth Observing System (GEOS) Data Assimilation System (DAS). This reflects our emphasis on model validation through the use of observed constituent mixing ratios.

Recently, the model has been improved in the transport as well as in the subgrid parameterization in order to enhance the accuracy and to ensure mass conservation (see New Development).

References

Lin, S.-J., and R. B. Rood, Multidimensional flux-form semi-Lagrangian transport schemes, Mon. Wea. Rev., 124, 2046-2070, 1996.


Updated: Apr 30, 2002
Author: Zhengxin Zhu (SSAI) zhu@code916.gsfc.nasa.gov

Responsible NASA organization/official: Dr. A. R. Douglass, Atmospheric Chemist ry and Dynamics Branch