The influence of hilly terrain on canopy-atmosphere carbon dioxide exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO2 (N-ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO2 sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO2 exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO2. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N-ee = F-c, the vertical turbulent flux of CO2. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F-c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F-c-at canopy top differs by more than 20% compared to the flat-terrain case.

On defining base units in terms of fundamental constants

The definitions of the base units of the international system of units have been revised many times since the idea of such an international system was first conceived at the time of the French revolution. The objective today is to define all our units in terms of 'invariants of nature', i.e. by referencing our units to the fundamental constants of physics, or the properties of atoms, rather than the characteristics of our planet or of artefacts. This situation is reviewed, particularly in regard to finding a new definition of the kilogram to replace its present definition in terms of a prototype material artefact.

Observations of the inter-ocean exchange around South Africa

There is growing evidence that the interocean exchange south of Africa is an important link in the global overturning circulation of the ocean, the so‐called ocean conveyer belt. At this location, warm and salty Indian Ocean waters enter the South Atlantic and are pulled by currents that eventually reach the North Atlantic, where water cools and sinks. A major contributor to the exchange is the frequent shedding of ring eddies from the termination of the Agulhas Current south of the tip of Africa. This shedding is controlled by developments far upstream in the Indian Ocean, and variations in this ‘Agulhas Leakage’ can lead to changes in the rate and stability of the Atlantic overturning, with possible associated global climate variations [Weijer et al., 1999]. Regional climate variations in the tropical and subtropical Indian Ocean are known to affect the whole system of the Agulhas Current, including the interocean exchanges. This article reports on some of the seminal results of ongoing multinational, multidisciplinary projects that explore these issues.

Harmonic force fields from scaled SCF calculations: Program ASYM40

We report an extended version of our normal coordinate program ASYM40, which may be used to transform Cartesian force constants from ab initio calculations to a force field in nonredundant internal (symmetry) coordinates. When experimental data are available, scale factors for the theoretical force field may then be optimized by least-squares refinement. The alternative of refining an empirical force field to fit a wide variety of data, as with the previous version ASYM20, has been retained. We compare the results of least-squares refinement of the full harmonic force field with least-squares refinement of only the scale factors for an SCF calculated force field and conclude that the latter approach may be useful for large molecules where more sophisticated calculations are impractical. The refinement of scale factors for a theoretical force field is also useful when there are only limited spectroscopic data. The program will accept ab initio calculated force fields from any program that presents Cartesian force constants as output. The program is available through Quantum Chemistry Program Exchange.

Force constants and dipole‐moment derivatives of molecules from perturbed Hartree–Fock calculations

General expressions for the force constants and dipole‐moment derivatives of molecules are derived, and the problems arising in their practical application are reviewed. Great emphasis is placed on the use of the Hartree–Fock function as an approximate wavefunction, and a number of its properties are discussed and re‐emphasised. The main content of this paper is the development of a perturbed Hartree–Fock theory that makes possible the direct calculation of force constants and dipole‐moment derivatives from SCF–MO wavefunctions. Essentially the theory yields ∂ϕi / ∂RJα, the derivative of an MO with respect to a nuclear coordinate.

The calculation of force constants and normal coordinates IV: XH4 and XH3 molecules

Normal coordinate calculations of XH4 and XH3 molecules are reviewed and discussed. It is shown that for most of these molecules the true values of the force constants in the most General Harmonic Force Field can be uniquely determined only by making use of vibration-rotation interaction constants. It is emphasized that without these extra data the GFF is not determined. The results are compared with various model force fields for these molecules.

Analytical potentials for triatomic molecules IX. The prediction of anharmonic force constants from potential energy surfaces based on harmonic force fields and dissociation energies for SO2 and O3

Analytical potential energy functions which are valid at all dissociation limits have been derived for the ground states of SO2 and O3. The procedure involves minimizing the errors between the observed vibrational spectra and spectra calculated by a variational procedure. Good agreement is obtained between the observed and calculated spectra for both molecules. Comparisons are made between anharmonic force fields, previously determined from the spectral data, and the force fields obtained by differentiating the derived analytical functions at the equilibrium configurations.

Microwave spectra and centrifugal distortion constants of oxetane

The microwave spectra of oxetane (trimethylene oxide) and its three symmetrically deuterated isotopic species have been observed on a Hewlett-Packard microwave spectrometer from 26.5 to 40 GHz. For the parent species, the β-d2 and the αα′-d4 species, about 300 lines have been assigned for each molecule, and for the d6 species more than 600 lines have been assigned. The assignments range from v = 0 to v = 5 in the puckering vibration; although they are mostly Q transitions, either 3 or 4 R transitions have been observed for each vibrational state. The spectra have been interpreted using an effective rotational hamiltonian for each vibrational state, including five quartic distortion constants according to Watson's formulation, and a variable number of sextic distortion constants; in general, the lines are fitted to about ± 10 kHz. The distortion constants show an anomalous zig-zag dependence on the puckering vibrational quantum number, similar to that first observed for the rotational constants by Gwinn and coworkers. This is interpreted according to a simple modification of the standard theory of centrifugal distortion, involving the double minimum potential function in the puckering coordinate.

Force constants and dipole moment derivatives of molecules from perturbed Hartree-Fock calculations II. Applications to limited basis set SCF-MO wavefunctions

The perturbed Hartree–Fock theory developed in the preceding paper is applied to LiH, BH, and HF, using limited basis‐set SCF–MO wavefunctions derived by previous workers. The calculated values for the force constant ke and the dipole‐moment derivative μ(1) are (experimental values in parentheses): LiH, ke  =  1.618(1.026)mdyn/Å,μ(1)  =  −18.77(−2.0±0.3)D/ÅBH,ke  =  5.199(3.032)mdyn/Å,μ(1)  =  −1.03(−)D/Å;HF,ke  =  12.90(9.651)mdyn/Å,μ(1)  =  −2.15(+1.50)D/Å. The values of the force on the proton were calculated exactly and according to the Hellmann–Feynman theorem in each case, and the discrepancies show that none of the wavefunctions used are close to the Hartree–Fock limit, so that the large errors in ke and μ(1) are not surprising. However no difficulties arose in the perturbed Hartree–Fock calculation, so that the application of the theory to more accurate wavefunctions appears quite feasible.

The calculation of force constants and normal co-ordinates II: methyl fluoride

The mathematical difficulties which can arise in the force constant refinement procedure for calculating force constants and normal co-ordinates are described and discussed. The method has been applied to the methyl fluoride molecule, using an electronic computer. The best values of the twelve force constants in the most general harmonic potential field were obtained to fit twenty-two independently observed experimental data, these being the six vibration frequencies, three Coriolis zeta constants and two centrifugal stretching constants DJ and DJK, for both CH3F and CD3F. The calculations have been repeated both with and without anharmonicity corrections to the vibration frequencies. All the experimental data were weighted according to the reliability of the observations, and the corresponding standard errors and correlation coefficients of the force constants have been deduced. The final force constants are discussed briefly, and compared with previous treatments, particularly with a recent Urey-Bradley treatment for this molecule.

The two-dimensional anharmonic oscillator I: vibration-rotation constants of fulminic acid, HCNO

The lowest-wavenumber vibration of HCNO and DCNO, ν5, is known to involve a largeamplitude low-frequency anharmonic bending of the CH bond against the CNO frame. In this paper the anomalous vibrational dependence of the observed rotational constants B(v5, l5), and of the observed l-doubling interactions, is interpreted according to a simple effective vibration-rotation Hamiltonian in which the appropriate vibrational operators are averaged in an anharmonic potential surface over the normal coordinates (Q5x, Q5y). All of the data on both isotopes are interpreted according to a single potential surface having a minimum energy at a slightly bent configuration of the HCN angle ( 170°) with a maximum at the linear configuration about 2 cm−1 higher. The other coefficients in the Hamiltonian are also interpreted in terms of the structure and the harmonic and anharmonic force fields; the substitution structure at the “hypothetical linear configuration” determined in this way gives a CH bond length of 1.060 Å, in contrast to the value 1.027 Å determined from the ground-state rotational constants. We also discuss the difficulties in rationalizing our effective Hamiltonian in terms of more fundamental theory, as well as the success and limitations of its use in practice.

The calculation of force constants and normal coordinates III: constrained force fields for the methyl halides

A method is discussed for imposing any desired constraint on the force field obtained in a force constant refinement calculation. The application of this method to force constant refinement calculations for the methyl halide molecules is reported. All available data on the vibration frequencies, Coriolis interaction constants and centrifugal stretching constants of CH3X and CD3X molecules were used in the refinements, but despite this apparent abundance of data it was found that constraints were necessary in order to obtain a unique solution to the force field. The results of unconstrained calculations, and of three different constrained calculations, are reported in this paper. The constrained models reported are a Urey—Bradley force field, a modified valence force field, and a constraint based on orbital-following bond-hybridization arguments developed in the following paper. The results are discussed, and compared with previous results for these molecules. The third of the above models is found to reproduce the observed data better than either of the first two, and additional reasons are given for preferring this solution to the force field for the methyl halide molecules.

Signs of l-doubling constants in NF3

Computed infrared band contours are presented for the two degenerate fundamentals of the NF3 molecule, using the l-resonance band contour program described by Cartwright and Mills, with values of the vibration-rotation constants determined from the microwave spectrum by Morino and co-workers. Computed contours are presented for both possible signs of the l-doubling constant, and comparison with the observed band contours leads to the conclusion that q3 = −121.4 MHz and q4 = +51.4 MHz.