A stochastic parameterization for deep convection based on equilibrium statistics

A stochastic parameterization scheme for deep convection is described, suitable for use in both climate and NWP models. Theoretical arguments and the results of cloud-resolving models, are discussed in order to motivate the form of the scheme. In the deterministic limit, it tends to a spectrum of entraining/detraining plumes and is similar to other current parameterizations. The stochastic variability describes the local fluctuations about a large-scale equilibrium state. Plumes are drawn at random from a probability distribution function (pdf) that defines the chance of finding a plume of given cloud-base mass flux within each model grid box. The normalization of the pdf is given by the ensemble-mean mass flux, and this is computed with a CAPE closure method. The characteristics of each plume produced are determined using an adaptation of the plume model from the Kain-Fritsch parameterization. Initial tests in the single column version of the Unified Model verify that the scheme is effective in producing the desired distributions of convective variability without adversely affecting the mean state.

Adsorption, orientation and thermal decomposition of copper(II) hexafluoroacetylacetonate on rutile TiO2(110)

We have investigated the adsorption and thermal decomposition of copper hexafluoroacetylacetonate (Cu-11(hfaC)(2)) on single crystal rutile TiO2(110). Low energy electron diffraction shows that room temperature saturation coverage of the Cu-II(hfac)(2) adsorbate forms an ordered (2 x 1) over-layer. X-ray and ultra-violet photoemission spectroscopy of the saturated surface were recorded as the sample was annealed in a sequential manner to reveal decomposition pathways. The results show that the molecule dissociatively adsorbs by detachment of one of the two ligands to form hfac and Cu-1(hfac) which chemisorb to the substrate at 298 K. These ligands only begin to decompose once the surface temperature exceeds 473 K where Cu core level shifts indicate metallisation. This reduction from Cu(I) to Cu(0) takes place in the absence of an external reducing agent and without disproportionation and is accompanied by the onset of decomposition of the hfac ligands. Finally, C K-edge near edge X-ray absorption fine structure experiments indicate that both the ligands adsorb aligned in the < 001 > direction and we propose a model in which the hfac ligands adsorb on the 5-fold coordinated Ti atoms and the Cu-1(hfac) moiety attaches to the bridging O atoms in a square planar geometry. The calculated tilt angle for these combined geometries is approximately 10 degrees to the surface normal.

Layer-by-layer electrostatic entrapment of protein molecules on superparamagnetic nanoparticle: new strategy to enhance adsorption capacity and maintain biological activity

There is a recent interest to use inorganic-based magnetic nanoparticles as a vehicle to carry biomolecules for various biophysical applications, but direct attachment of the molecules is known to alter their conformation leading to attenuation in activity. In addition, surface immobilization has been limited to monolayer coverage. It is shown that alternate depositions of negatively charged protein molecules, typically bovine serum albumin (BSA) with a positively charged aminocarbohydrate template such as glycol chitosan (GC) on magnetic iron oxide nanoparticle surface as a colloid, are carried out under pH 7.4. Circular dichroism (CD) clearly reveals that the secondary structure of the entrapped BSA sequential depositions in this manner remains totally unaltered which is in sharp contrast to previous attempts. Probing the binding properties of the entrapped BSA using small molecules (Site I and Site II drug compounds) confirms for the first time the full retention of its biological activity as compared with native BSA, which also implies the ready accessibility of the entrapped protein molecules through the porous overlayers. This work clearly suggests a new method to immobilize and store protein molecules beyond monolayer adsorption on a magnetic nanoparticle surface without much structural alteration. This may find applications in magnetic recoverable enzymes or protein delivery.

Does reactive surface area depend on grain size? Results from pH 3, 25 circle C far-from-equilibrium flow-through dissolution experiments on anorthite and biotite

Laboratory determined mineral weathering rates need to be normalised to allow their extrapolation to natural systems. The principle normalisation terms used in the literature are mass, and geometric- and BET specific surface area (SSA). The purpose of this study was to determine how dissolution rates normalised to these terms vary with grain size. Different size fractions of anorthite and biotite ranging from 180-150 to 20-10 mu m were dissolved in pH 3, HCl at 25 degrees C in flow through reactors under far from equilibrium conditions. Steady state dissolution rates after 5376 h (anorthite) and 4992 h (biotite) were calculated from Si concentrations and were normalised to initial- and final- mass and geometric-, geometric edge- (biotite), and BET SSA. For anorthite, rates normalised to initial- and final-BET SSA ranged from 0.33 to 2.77 X 10(-10) mol(feldspar) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 5.74 to 8.88 X 10(-10) mol(feldspar) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.11 to 1.65 mol(feldspar) g(-1) s(-1). For biotite, rates normalised to initial- and final-BET SSA ranged from 1.02 to 2.03 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 3.26 to 16.21 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric edge SSA ranged from 59.46 to 111.32 x 10(-12) mol(biotite) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.81 to 6.93 X 10(-12) mol(biotite) g(-1) s(-1). For all normalising terms rates varied significantly (p <= 0.05) with grain size. The normalising terms which gave least variation in dissolution rate between grain sizes for anorthite were initial BET SSA and initial- and final-geometric SSA. This is consistent with: (1) dissolution being dominated by the slower dissolving but area dominant non-etched surfaces of the grains and, (2) the walls of etch pits and other dissolution features being relatively unreactive. These steady state normalised dissolution rates are likely to be constant with time. Normalisation to final BET SSA did not give constant ratios across grain size due to a non-uniform distribution of dissolution features. After dissolution coarser grains had a greater density of dissolution features with BET-measurable but unreactive wall surface area than the finer grains. The normalising term which gave the least variation in dissolution rates between grain sizes for biotite was initial BET SSA. Initial- and final-geometric edge SSA and final BET SSA gave the next least varied rates. The basal surfaces dissolved sufficiently rapidly to influence bulk dissolution rate and prevent geometric edge SSA normalised dissolution rates showing the least variation. Simple modelling indicated that biotite grain edges dissolved 71-132 times faster than basal surfaces. In this experiment, initial BET SSA best integrated the different areas and reactivities of the edge and basal surfaces of biotite. Steady state dissolution rates are likely to vary with time as dissolution alters the ratio of edge to basal surface area. Therefore they would be more properly termed pseudo-steady state rates, only appearing constant because the time period over which they were measured (1512 h) was less than the time period over wich they would change significantly. (c) 2006 Elsevier Inc. All rights reserved.

Salt marsh morphodynamics: An investigation of tidal flows and marsh channel equilibrium

Although numerous field studies have evaluated flow and transport processes in salt marsh channels, the overall role of channels in delivering and removing material from salt marsh platforms is still poorly characterised. In this paper, we consider this issue based on a numerical hydrodynamic model for a prototype marsh system and on a field survey of the cross-sectional geometry of a marsh channel network. Results of the numerical simulations indicate that the channel transfers approximately three times the volume of water that would be estimated from mass balance considerations alone. Marsh platform roughness exerts a significant influence on the partitioning of discharge between the channel and the marsh platform edge, alters flow patterns on the marsh platform due to its effects on channel-to-platform transfer and also controls the timing of peak discharge relative to marsh-edge overtopping. Although peak channel discharges and velocities are associated with the flood tide and marsh inundation, a larger volume of water is transferred by the channel during ebb flows, a portion of which transfer takes place after the tidal height is below the marsh platform. Detailed surveys of the marsh channels crossing a series of transects at Upper Stiffkey Marsh, north Norfolk, England, show that the total channel cross-sectional area increases linearly with catchment area in the inner part of the marsh, which is consistent with the increase in shoreward tidal prism removed by the channels. Toward the marsh edge, however, a deficit in the total cross-sectional area develops, suggesting that discharge partitioning between the marsh channels and the marsh platform edge may also be expressed in the morphology of marsh channel systems.

Adsorption of water vapor by bare soil in an olive grove in southern Spain

Data for water vapor adsorption and evaporation are presented for a bare soil (sandy loam, clay content 15%) in a southern Spanish olive grove. Water losses and gains were measured using eight high-precision minilysimeters, placed around an olive tree, which had been irrigated until the soil reached field capacity (similar to 0.22 m(3) m(-3)). They were subsequently left to dry for 10 days. A pair of lysimeters was situated at each of the main points of the compass (N, E, S, W), at a distance of 1 m (the inner set of lysimeters; ILS) and 2 m (the outer set of lysimeters; OLS), respectively, from the tree trunk. Distinct periods of moisture loss (evaporation) and moisture gain (vapor adsorption) could be distinguished for each day. Vapor adsorption often started just after noon and generally lasted until the (early) evening. Values of up to 0.7 mm of adsorbed water per day were measured. Adsorption was generally largest for the OLS (up to 100% more on a daily basis), and increased during the dry down. This was mainly the result of lower OLS surface soil moisture contents (period-average absolute difference similar to 0.005 m(3) m(-3)), as illustrated using various analyses employing a set of micrometeorological equations describing the exchange of water vapor between bare soil and the atmosphere. These analyses also showed that the amount of water vapor adsorbed by soils is very sensitive to changes in atmospheric forcing and surface variables. The use of empirical equations to estimate vapor adsorption is therefore not recommended.

Mesoscale simulations of organized convection: importance of convective equilibrium

The validity of convective parametrization breaks down at the resolution of mesoscale models, and the success of parametrized versus explicit treatments of convection is likely to depend on the large-scale environment. In this paper we examine the hypothesis that a key feature determining the sensitivity to the environment is whether the forcing of convection is sufficiently homogeneous and slowly varying that the convection can be considered to be in equilibrium. Two case studies of mesoscale convective systems over the UK, one where equilibrium conditions are expected and one where equilibrium is unlikely, are simulated using a mesoscale forecasting model. The time evolution of area-average convective available potential energy and the time evolution and magnitude of the timescale of convective adjustment are consistent with the hypothesis of equilibrium for case 1 and non-equilibrium for case 2. For each case, three experiments are performed with different partitionings between parametrized and explicit convection: fully parametrized convection, fully explicit convection and a simulation with significant amounts of both. In the equilibrium case, bulk properties of the convection such as area-integrated rain rates are insensitive to the treatment of convection. However, the detailed structure of the precipitation field changes; the simulation with parametrized convection behaves well and produces a smooth field that follows the forcing region, and the simulation with explicit convection has a small number of localized intense regions of precipitation that track with the mid-levelflow. For the non-equilibrium case, bulk properties of the convection such as area-integrated rain rates are sensitive to the treatment of convection. The simulation with explicit convection behaves similarly to the equilibrium case with a few localized precipitation regions. In contrast, the cumulus parametrization fails dramatically and develops intense propagating bows of precipitation that were not observed. The simulations with both parametrized and explicit convection follow the pattern seen in the other experiments, with a transition over the duration of the run from parametrized to explicit precipitation. The impact of convection on the large-scaleflow, as measured by upper-level wind and potential-vorticity perturbations, is very sensitive to the partitioning of convection for both cases. © Royal Meteorological Society, 2006. Contributions by P. A. Clark and M. E. B. Gray are Crown Copyright.

Interfacial structure of wild-type and mutant forms of puroindoline-b bound to DPPG monolayers

The interaction of wild-type puroindoline-b (Pin-b+) and two mutant forms having single residue substitutions (G46S or W44R) with L-alpha-dipalmitoylphosphatidyl-dl-glycerol (DPPG) as a Langmuir monolayer at the air/water interface was investigated by neutron reflectivity (NR) and Brewster angle microscopy (BAM). NR profiles were fitted using a three-layer model to enable differences in penetration of protein between the lipid headgroup and acyl regions to be determined. The data showed similar surface excesses for each of the three proteins at the interface; however, it was revealed that the depth of penetration of protein into the lipid region differed for each protein with Pin-b+ penetrating further into the acyl region of the lipid compared to the mutant forms of the protein that interacted with the headgroup region only. BAM images revealed that the domain structure of the DPPG monolayers was disrupted when Pin-b+ adsorption had reached equilibrium, suggesting protein penetration had led to compression of the lipid region. In contrast, the domain structure was unaffected by the W44R mutant, suggesting no change in compression of the lipid region and hence little or no penetration of protein into the lipid layer.

The equilibrium structure of HCN

The equilibrium structure of HCN has been determined from the previously published ground state rotational constants of eight isotopomers by using (B0‐Be) values obtained from a variational calculation of the vibration–rotation spectrum. The results are re(CH)=1.065 01(8) Å, and re(CN)=1.153 24(2) Å.

The infrared spectrum, equilibrium structure and harmonic and anharmonic force field of thioborine, HBS

Infrared spectra of the two stretching fundamentals of both HBS and DBS have been observed, using a continuous flow system through a multiple reflection long path cell at a pressure around 1 Torr and a Nicolet Fourier Transform spectrometer with a resolution of about 0•1 cm-1. The v3 BS stretching fundamental of DBS, near 1140 cm-1, is observed in strong Fermi resonance with the overtone of the bend 2v2. The bending fundamental v2 has not been observed and must be a very weak band. The analysis of the results in conjunction with earlier work gives the equilibrium structure (re(BH) = 1•1698(12) , re(BS) = 1•5978(3) ) and the harmonic and anharmonic force field.

The equilibrium structure of monofluoroacetylene

The equilibrium rotational constants Be of HCCF and DCCF have been determined from the ground state rotational constants B0, by determining the αr constants for all five fundamentals from the high-resolution vibrational—rotation spectrum making appropriate corrections for the effects of Fermi resonance. By combination with results from the 13C isotopomers and the recent ab initio calculations by Botschwina (Chem. Phys. Lett., 209 (1993) 117), the equilibrium structure is deduced to be: re(CH) = 1.0555(15) Å, re(CC) = 1.1955(8) Å and re(CF) = 1.2781(8) Å.

Influence of adsorption geometry in the heterogeneous enantioselective catalytic hydrogenation of a prototypical enone

Asymmetric catalysis is of paramount importance in organic synthesis and, in current practice, is achieved by means of homogeneous catalysts. The ability to catalyze such reactions heterogeneously would have a major impact both in the research laboratory and in the production of fine chemicals and pharmaceuticals, yet heterogeneous asymmetric hydrogenation of C═C bonds remains hardly explored. Very recently, we demonstrated how chiral ligands that anchor robustly to the surface of Pd nanoparticles promote asymmetric catalytic hydrogenation: ligand rigidity and stereochemistry emerged as key factors. Here, we address a complementary question: how does the enone reactant adsorb on the metal surface, and what implications does this have for the enantiodifferentiating interaction with the surface-tethered chiral modifiers? A reaction model is proposed, which correctly predicts the identity of the enantiomer experimentally observed in excess.