The counter-propagating Rossby-wave perspective on baroclinic instability. II: Application to the Charney model

The constant-density Charney model describes the simplest unstable basic state with a planetary-vorticity gradient, which is uniform and positive, and baroclinicity that is manifest as a negative contribution to the potential-vorticity (PV) gradient at the ground and positive vertical wind shear. Together, these ingredients satisfy the necessary conditions for baroclinic instability. In Part I it was shown how baroclinic growth on a general zonal basic state can be viewed as the interaction of pairs of ‘counter-propagating Rossby waves’ (CRWs) that can be constructed from a growing normal mode and its decaying complex conjugate. In this paper the normal-mode solutions for the Charney model are studied from the CRW perspective. Clear parallels can be drawn between the most unstable modes of the Charney model and the Eady model, in which the CRWs can be derived independently of the normal modes. However, the dispersion curves for the two models are very different; the Eady model has a short-wave cut-off, while the Charney model is unstable at short wavelengths. Beyond its maximum growth rate the Charney model has a neutral point at finite wavelength (r=1). Thereafter follows a succession of unstable branches, each with weaker growth than the last, separated by neutral points at integer r—the so-called ‘Green branches’. A separate branch of westward-propagating neutral modes also originates from each neutral point. By approximating the lower CRW as a Rossby edge wave and the upper CRW structure as a single PV peak with a spread proportional to the Rossby scale height, the main features of the ‘Charney branch’ (0neutral modes appearing at r=1. The behaviour of the other branches and neutral points is essentially the same when viewed from the CRW perspective, but with cancelling interior PV structures reducing the self and mutual interaction of the CRWs. The underlying dynamics determining the nature of all the solutions is the difference in the scale-dependence of PV inversion for boundary and interior PV anomalies, the Rossby-wave propagation mechanism and the CRW interaction. The behaviour of the Charney modes and the first neutral branch, which rely on tropospheric PV gradients, are arguably more applicable to the atmosphere than modes of the Eady model where the positive PV gradient exists only at the tropopause

DTM generation from LIDAR data using skewness balancing

LIght Detection And Ranging (LIDAR) data for terrain and land surveying has contributed to many environmental, engineering and civil applications. However, the analysis of Digital Surface Models (DSMs) from complex LIDAR data is still challenging. Commonly, the first task to investigate LIDAR data point clouds is to separate ground and object points as a preparatory step for further object classification. In this paper, the authors present a novel unsupervised segmentation algorithm-skewness balancing to separate object and ground points efficiently from high resolution LIDAR point clouds by exploiting statistical moments. The results presented in this paper have shown its robustness and its potential for commercial applications.

Linear amplification of marginally neutral baroclinic waves

Baroclinic wave development is investigated for unstable parallel shear flows in the limit of vanishing normal-mode growth rate. This development is described in terms of the propagation and interaction mechanisms of two coherent structures, called counter-propagating Rossby waves (CRWs). It is shown that, in this limit of vanishing normal-mode growth rate, arbitrary initial conditions produce sustained linear amplification of the marginally neutral normal mode (mNM). This linear excitation of the mNM is subsequently interpreted in terms of a resonance phenomenon. Moreover, while the mathematical character of the normal-mode problem changes abruptly as the bifurcation point in the dispersion diagram is encountered and crossed, it is shown that from an initial-value viewpoint, this transition is smooth. Consequently, the resonance interpretation remains relevant (albeit for a finite time) for wavenumbers slightly different from the ones defining cut-off points. The results are further applied to a two-layer version of the classic Eady model in which the upper rigid lid has been replaced by a simple stratosphere.

Comparison of EDTA and EDDS as potential soil amendments for enhanced phytoextraction of heavy metals

Phytoextraction has been proposed as an alternative remediation technology for soils polluted with heavy metals or radionuclides, but is generally conceived as too slow working. Enhancing the accumulation of trace pollutants in harvestable plant tissues is a prerequisite for the technology to be practically applicable. The chelating aminopolycarboxylic acid, ethylene diamine tetraacetate (EDTA), has been found to enhance shoot accumulation of heavy metals. However, the use of EDTA in phytoextraction may not be suitable due to its high environmental persistence, which may lead to groundwater contamination. This paper aims to assess whether ethylene diamine disuccinate (EDDS), a biodegradable chelator, can be used for enhanced phytoextraction purposes. A laboratory experiment was conducted to examine mobilisation of Cd, Cu, Cr, Ni, Pb and Zn into the soil solution upon application of EDTA or EDDS. The longevity of the induced mobilisation was monitored for a period of 40 days after application. Estimated effect half lives ranged between 3.8 and 7.5 days for EDDS, depending on the applied dose. The minimum observed effect half life of EDTA was 36 days, while for the highest applied dose no decrease was observed throughout the 40 day period of the mobilisation experiment. Performance of EDTA and EDDS for phytoextraction was evaluated by application to Helianthus annuus. Two other potential chelators, known for their biodegradability in comparison to EDTA, were tested in the plant experiment: nitrilo acetic acid (NTA) and citric acid. Uptake of heavy metals was higher in EDDS-treated pots than in EDTA-treated pots. The effects were still considered insufficiently high to consider efficient remediation. This may be partly due to the choice of timing for application of the soil amendment. Fixing the time of application at an earlier point before harvest may yield better results. NTA and citric acid induced no significant effects on heavy metal uptake. (C) 2004 Elsevier Ltd. All rights reserved.

Wetland restoration within agricultural watersheds: Balancing water quality protection with habitat conservation

Peat wetlands that have been restored from agricultural Land have the potential to act as Long term sources of phosphorus (P) and, therefore have to potenital to accelerate freshwater eutrophication. During a two-year study the water table in a eutrophic fen peat that was managed by pump drainage fluctuated annually between +20 cm and -60 cm relative to ground Level. This precise management was facilitated by the high hydraulic conductivity (K) of the humified peat (1.1 x 10(-5) m s(-1)) below around 60 cm depth. However, during one week of intermittent pumping, as much as 50 g ha(-1) dissolved P entered the pumped ditch. Summer. rainfall events and autumn reflooding also triggered P losses. The P Losses were attributed to the low P sorption capacity (217 mg kg(-1)) of the saturated peat below 60 cm, combined with its high K and the reductive dissolution of Fe bound P.

A potential energy surface for the ground state of acetylene, C2H2

A potential energy function has been derived for the ground state surface of C2H2 as a many-body expansion. The 2- and 3-body terms have been obtained by preliminary investigation of the ground state surfaces of CH2( 3B1) and C2H( 2Σ+). A 4-body term has been derived which reproduces the energy, geometry and harmonic force field of C2H2. The potential has a secondary minimum corresponding to the vinylidene structure and the geometry and energy of this are in close agreement with predictions from ab initio calculations. The saddle point for the HCCH-H2CC rearrangement is predicted to lie 2•530 eV above the acetylene minimum.

In situ degradability of dry matter and neutral-detergent fibre of thorn scrubland forage consumed by goats in the semi-arid region of north Mexico

Three goats provided with oesophageal and ruminal cannulae were used to determine variations in dry matter (DM) and neutral-detergent fibre (NDF) degradability of the forage consumed when grazing thorn scrubland in the semi-arid region of north Mexico, during two consecutive dry and wet periods. Ingesta samples were incubated intraruminally, the data were fitted to the exponential equation P = a + b (1-e(-ct)) and statistically analysed using a randomized-block design. Organic matter and crude protein (CP) contents were higher (P < 0.05) in the wet seasons. Values of NDF were similar in dry and wet season of both years whereas higher numerical values of acid-detergent fibre (ADF), lignin and cellulose were registered in the dry seasons. DM and NDF degradabilities after 24 and 48 h of ruminal incubation were higher (P < 0.05) in the wet seasons. Higher values (P < 0.05) in DM and NDF bag losses at zero time (A fraction) were registered in the two wet seasons. The insoluble but fermentable DM and NDF (B fractions) were higher (P < 0.05) in the 1999 wet season and variable in the rest of the studied period. Numerically higher values of DM and NDF c fraction were found in wet periods, whereas DM and NDF potential degradabilities were higher (P < 0.05) in the wet season in 1999 and similar across seasons in 2000. Lowest (P < 0.05) contents of CP in grazed forage, DM and NDF degradabilities after 48 h of ruminal incubation, and A, and B, and c fractions were observed in the dry seasons. Thus, these results may be related to both the lower feeding value of forage consumed by the animals and lower performance of livestock during this period. Then, the DM and NDF degradability after 48 h, together with the insoluble but fermentable matter and the c fraction permit the nutritive value of the forage consumed by grazing goats to be accurately described.

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcal/mol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively. (C) 2008 American Institute of Physics.

MULTIMODE quantum calculations of intramolecular vibrational energies of the water dimer and trimer using ab initio-based potential energy surfaces

We report vibrational configuration interaction calculations of the monomer fundamentals of (H2O)(2), (D2O)(2), (H2O)(3), and (D2O)(3) using the code MULTIMODE and full dimensional ab initio-based global potential energies surfaces (PESs). For the dimer the HBB PES [Huang , J. Chem. Phys 128, 034312 (2008)] is used and for the trimer a new PES, reported here, is used. The salient properties of the new trimer PES are presented and compared to previous single-point calculations and the vibrational energies are compared with experiments. (C) 2008 American Institute of Physics.

Case study: object and ground point separation on a river bank using airborne based LIDAR data

In the past decade, airborne based LIght Detection And Ranging (LIDAR) has been recognised by both the commercial and public sectors as a reliable and accurate source for land surveying in environmental, engineering and civil applications. Commonly, the first task to investigate LIDAR point clouds is to separate ground and object points. Skewness Balancing has been proven to be an efficient non-parametric unsupervised classification algorithm to address this challenge. Initially developed for moderate terrain, this algorithm needs to be adapted to handle sloped terrain. This paper addresses the difficulty of object and ground point separation in LIDAR data in hilly terrain. A case study on a diverse LIDAR data set in terms of data provider, resolution and LIDAR echo has been carried out. Several sites in urban and rural areas with man-made structure and vegetation in moderate and hilly terrain have been investigated and three categories have been identified. A deeper investigation on an urban scene with a river bank has been selected to extend the existing algorithm. The results show that an iterative use of Skewness Balancing is suitable for sloped terrain.

Bio-energy retains its mitigation potential under elevated CO2

Background If biofuels are to be a viable substitute for fossil fuels, it is essential that they retain their potential to mitigate climate change under future atmospheric conditions. Elevated atmospheric CO2 concentration [CO2] stimulates plant biomass production; however, the beneficial effects of increased production may be offset by higher energy costs in crop management. Methodology/Main findings We maintained full size poplar short rotation coppice (SRC) systems under both current ambient and future elevated [CO2] (550 ppm) and estimated their net energy and greenhouse gas balance. We show that a poplar SRC system is energy efficient and produces more energy than required for coppice management. Even more, elevated [CO2] will increase the net energy production and greenhouse gas balance of a SRC system with 18%. Managing the trees in shorter rotation cycles (i.e. 2 year cycles instead of 3 year cycles) will further enhance the benefits from elevated [CO2] on both the net energy and greenhouse gas balance. Conclusions/significance Adapting coppice management to the future atmospheric [CO2] is necessary to fully benefit from the climate mitigation potential of bio-energy systems. Further, a future increase in potential biomass production due to elevated [CO2] outweighs the increased production costs resulting in a northward extension of the area where SRC is greenhouse gas neutral. Currently, the main part of the European terrestrial carbon sink is found in forest biomass and attributed to harvesting less than the annual growth in wood. Because SRC is intensively managed, with a higher turnover in wood production than conventional forest, northward expansion of SRC is likely to erode the European terrestrial carbon sink.

Threshold-free object and ground point separation in LIDAR data

Light Detection And Ranging (LIDAR) is an important modality in terrain and land surveying for many environmental, engineering and civil applications. This paper presents the framework for a recently developed unsupervised classification algorithm called Skewness Balancing for object and ground point separation in airborne LIDAR data. The main advantages of the algorithm are threshold-freedom and independence from LIDAR data format and resolution, while preserving object and terrain details. The framework for Skewness Balancing has been built in this contribution with a prediction model in which unknown LIDAR tiles can be categorised as “hilly” or “moderate” terrains. Accuracy assessment of the model is carried out using cross-validation with an overall accuracy of 95%. An extension to the algorithm is developed to address the overclassification issue for hilly terrain. For moderate terrain, the results show that from the classified tiles detached objects (buildings and vegetation) and attached objects (bridges and motorway junctions) are separated from bare earth (ground, roads and yards) which makes Skewness Balancing ideal to be integrated into geographic information system (GIS) software packages.

Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Stratospheric heating by potential geoengineering aerosols

A fixed dynamical heating model is used to investigate the pattern of zonal-mean stratospheric temperature change resulting from geoengineering with aerosols composed of sulfate, titania, limestone and soot. Aerosol always heats the tropical lower stratosphere, but at the poles the response can be either heating, cooling, or neutral. The sign of the change in stratospheric Pole-Equator temperature difference depends on aerosol type, size and season. This has implications for modelling geoengineering impacts and the response of the stratospheric circulation.