Photochemical Numerics for Global-Scale Modeling: Fidelity and GCM Testing

Elliott, Scott; Zhao, Xuepeng; Turco, R. P.; Kao, Chih-Yue Jim; Shen, Mei

doi:10.1175/1520-0450(1995)034<0694:pnfgsm>2.0.co;2

Cited by 10 publications

(9 citation statements)

References 0 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Implicit integration methods of varying sophistication were experimented with (206). Chemically intuitive groupings such as "families" were found to shorten computation times (207). The NO x species constitute an example.…”

Section: Modelsmentioning

confidence: 99%

ENERGY ANDMATERIALFLOWTHROUGH THEURBANECOSYSTEM

Decker

Elliott

Smith

et al. 2000

Annu. Rev. Energy. Environ.

Self Cite

338

226

View full text Add to dashboard Cite

Key Words megacity biogeochemistry, city succession, industrial metabolism, air pollution, system simulations s Abstract This paper reviews the available data and models on energy and material flows through the world's 25 largest cities. Throughput is categorized as stored, transformed, or passive for the major flow modes. The aggregate, fuel, food, water, and air cycles are all examined. Emphasis is placed on atmospheric pathways because the data are abundant. Relevant models of urban energy and material flows, demography, and atmospheric chemistry are discussed. Earth system-level loops from cities to neighboring ecosystems are identified. Megacities are somewhat independent of their immediate environment for food, fuel, and aggregate inputs, but all are constrained by their regional environment for supplying water and absorbing wastes. We elaborate on analogies with biological metabolism and ecosystem succession as useful conceptual frameworks for addressing urban ecological problems. We conclude that whereas data are numerous for some individual cities, cross-cutting compilations are lacking in biogeochemical analysis and modeling. Synthesis of the existing information will be a crucial first step. Cross-cutting field research and integrated, multidisciplinary simulations will be necessary.

show abstract

Section: Modelsmentioning

confidence: 99%

ENERGY ANDMATERIALFLOWTHROUGH THEURBANECOSYSTEM

Decker

Elliott

Smith

et al. 2000

Annu. Rev. Energy. Environ.

Self Cite

338

226

View full text Add to dashboard Cite

show abstract

“…That leads to the problem of determining how to partition each family among its member species. Our approach utilizes a unique multistep "projection" and partitioning scheme [Elliott et al, 1993[Elliott et al, , 1995Zhao, 1995] that ensures computational stability even with large time steps. Basically, the family continuity equation is solved explicitly using chemical rates averaged over the time step.…”

Section: Chemistry Componentmentioning

confidence: 99%

Aerosol‐induced chemical perturbations of stratospheric ozone: Three‐dimensional simulations and analysis of mechanisms

Zhao¹,

Turco²,

Kao³

et al. 1997

J. Geophys. Res.

Self Cite

View full text Add to dashboard Cite

Abstract. An atmospheric general circulation model is coupled with a stratospheric photochemical model to simulate the chemical/dynamical perturbations associated with background and volcanically perturbed aerosols in the lower stratosphere. The present work focuses on short-term anomalies at middle and high latitudes in the northern hemisphere, where large ozone depletions have been observed in late winter and early spring, particularly following the eruption of Mount Pinatubo. Five fully coupled simulations are analyzed, corresponding to a control case with only gas phase chemistry, and cases including heterogeneous chemistry on background aerosols, on E1 Chich6n-type, and on Pinatubo-type aerosols. It is found that heterogeneous reactions occurring on sulfate aerosols (background or postvolcanic) can strongly perturb the chemical partitioning in the lower stratosphere, leading to significant ozone depletion through enhanced chlorine, bromine, and odd-hydrogen catalytic cycles. In the Arctic lower stratosphere, the maximum zonal and March monthly mean local ozone reductions (with respect to the control case) can exceed 15% for the background aerosol case, 40% for the E1 Chich6n case, and 50% for the Pinatubo case. The corresponding zonal mean total column ozone decreases are roughly 5% and 15% for the background and volcanic aerosol cases, respectively. In the most extreme case tested (post-Pinatubo), a large ozone depletion below 30 mbar is offset to some extent by an ozone increase above that level. The results of a sensitivity study (in which the aerosols are distributed closer to the tropics, as might occur early after an eruption at low latitude) lead to relatively small total ozone depletions at northern high latitudes, and small ozone increases in the tropical lower stratosphere. The reduced impact on total ozone at high latitudes is associated both with local ozone increases above 30 mbar and with poleward transport of enhanced ozone from the tropical lower stratosphere. The ozone increase at low latitudes is the net result of compensating changes in the catalytic destruction cycles involving odd-nitrogen and chlorine species activated by heterogeneous processes at the low temperatures and abundant sunlight found near the tropical tropopause. Our simulations indicate that ozone variations triggered by volcanic injections of aerosols depend on the global distribution as well as the abundance of the particles and their evolution over time, and on the initial dynamical-radiative-chemical state of the atmosphere, which itself exhibits large seasonal and interannual variability.

show abstract

“…In one approach, researchers have sought to reduce the effective size of the JacobJan system matrix, thus reducing the number of computations needed to carry out the required matrix inversions. Elliott et al [1993Elliott et al [ , 1995Elliott et al [ , 1996 family. Owing to the longer lifetime of the family compared to its member species, the family rate equation can often be solved using standard algorithms, including low-order explicit numerical schemes, whereas similar solutions of the species rate equations would be unstable.…”

Section: Introductionmentioning

confidence: 99%

“…In one approach, researchers have sought to reduce the effective size of the JacobJan system matrix, thus reducing the number of computations needed to carry out the required matrix inversions. Elliott et al [1993Elliott et al [ , 1995Elliott et al [ , 1996 developed an approach employing "minimatrices" that represent submatrices of the Jacobian corresponding to specific families or groups of related species [e.g., Crutzen, 1971;Turco and Whitten, 1974;Cicerone et al, 1983;Brasseur and Solomon, 1984;Austin, 1991]. Families are selected so that their equivalent chemical rate equations, which are derived by summing the member species' rate equations [see Turco and Whitten, 1974], are much less stiff than those of the individual species.…”

Section: Introductionmentioning

confidence: 99%

“…The equivalent (summed) family rate equations can then be solved using relatively simple algorithms. Elliott et al [1993Elliott et al [ , 1995Elliott et al [ , 1996 demonstrated that this minimatrix approach can be highly efficient and accurate in many situations involving stratospheric and tropospheric chemical systcms. However, on a global grid the explicit family solutions can bccomc unstable in certain regimes, leading to substantial errors (see below).…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

A new family Jacobian solver for global three‐dimensional modeling of atmospheric chemistry

Zhao¹,

Turco²,

Shen³

1999

J. Geophys. Res.

Self Cite

View full text Add to dashboard Cite

Abstract. We present a new technique to solve complex sets of photochemical rate equations that is applicable to global modeling of the troposphere and stratosphere. The approach is based on the concept of "families" of species, whose chemical rate equations are tightly coupled. Variations of species concentrations within a family can be determined by inverting a linearized Jacobian matrix representing the family group. Since this group consists of a relatively small number of species the corresponding Jacobian has a low order (a minimatrix) compared to the Jacobian of the entire system. However, we go further and define a super-family that is the set of all families. The super-family is also solved by linearization and matrix inversion. The resulting Super-Family Matrix Inversion (SFMI) scheme is more stable and accurate than common family approaches. We discuss the numerical structure of the SFMI scheme and apply our algorithms to a comprehensive set of photochemical reactions. To evaluate performance, the SFMI scheme is compared with an optimized Gear solver. We find that the SFMI technique can be at least an order of magnitude more efficient than existing chemical solvers while maintaining relative errors in the calculations of 15% or less over a diurnal cycle. The largest SFMI errors arise at sunrise and sunset and during the evening when species concentrations may be very low. We show that sunrise/sunset errors can be minimized through a careful treatment of photodissociation during these periods; the nighttime deviations are negligible from the point of view of acceptable computational accuracy. The stability and flexibility of the SFMI algorithm should be sufficient for most modeling applications until major improvements in other modeling factors are achieved. In addition, because of its balanced computational design, SFMI can easily be adapted to parallel computing architectures. SFMI thus should allow practical long-term integrations of global chemistry coupled to general circulation and climate models, studies of interannual and interdecadal variability in atmospheric composition, simulations of past multidecadal trends owing to anthropogenic emissions, long-term forecasting associated with projected emissions, and sensitivity analyses for a wide range of physical and chemical parameters. To carry out practical calculations, the earliest and most straightforward explicit approaches for solving atmospheric chemical rate equations were quickly abandoned in favor of numerical implicitness (or semi-implicitness), with or without iterative procedures to stabilize the solutions [e.g., see Shimazaki and Laird, 1970; Turco and Whirten, 1974]. At a given time step, species concentrations could be repeatedly calculated using a simple semi implicit form of (1) by updating the chemical production and loss rates using species concentrations calculated at the previous iteration. While this gener-1767

show abstract

Photochemical Numerics for Global-Scale Modeling: Fidelity and GCM Testing

Cited by 10 publications

References 0 publications

ENERGY ANDMATERIALFLOWTHROUGH THEURBANECOSYSTEM

ENERGY ANDMATERIALFLOWTHROUGH THEURBANECOSYSTEM

Aerosol‐induced chemical perturbations of stratospheric ozone: Three‐dimensional simulations and analysis of mechanisms

A new family Jacobian solver for global three‐dimensional modeling of atmospheric chemistry

Contact Info

Product

Resources

About