How to parametrize urban-canopy drag to reproduce wind-direction effects within the canopy

The mean wind direction within an urban canopy changes with height when the incoming flow is not orthogonal to obstacle faces. This wind-turning effect is induced by complex processes and its modelling in urban-canopy (UC) parametrizations is difficult. Here we focus on the analysis of the spatially-averaged flow properties over an aligned array of cubes and their variation with incoming wind direction. For this purpose, Reynolds-averaged Navier–Stokes simulations previously compared, for a reduced number of incident wind directions, against direct numerical simulation results are used. The drag formulation of a UCparametrization ismodified and different drag coefficients are tested in order to reproduce the wind-turning effect within the canopy for oblique wind directions. The simulations carried out for a UC parametrization in one-dimensional mode indicate that a height-dependent drag coefficient is needed to capture this effect.

The effect of a new drag-law parameterization on ice shelf water plume dynamics

A drag law accounting for Ekman rotation adjacent to a flat, horizontal bou ndary is proposed for use in a plume model that is written in terms of the depth-mean velocity. The drag l aw contains a variable turning angle between the mean velocity and the drag imposed by the turbulent bound ary layer. The effect of the variable turning angle in the drag law is studied for a plume of ice shelf wat er (ISW) ascending and turning beneath an Antarctic ice shelf with draft decreasing away from the groundi ng line. As the ISW plume ascends the sloping ice shelf–ocean boundary, it can melt the ice shelf, wh ich alters the buoyancy forcing driving the plume motion. Under these conditions, the typical turning ang le is of order 10° over most of the plume area for a range of drag coefficients (the minus sign arises for th e Southern Hemisphere). The rotation of the drag with respect to the mean velocity is found to be signifi cant if the drag coefficient exceeds 0.003; in this case the plume body propagates farther along and across the b ase of the ice shelf than a plume with the standard quadratic drag law with no turning angle.

Impact of variable atmospheric and oceanic form drag on simulations of Arctic sea ice

Over Arctic sea ice, pressure ridges and floe andmelt pond edges all introduce discrete obstructions to the flow of air or water past the ice and are a source of form drag. In current climate models form drag is only accounted for by tuning the air–ice and ice–ocean drag coefficients, that is, by effectively altering the roughness length in a surface drag parameterization. The existing approach of the skin drag parameter tuning is poorly constrained by observations and fails to describe correctly the physics associated with the air–ice and ocean–ice drag. Here, the authors combine recent theoretical developments to deduce the total neutral form drag coefficients from properties of the ice cover such as ice concentration, vertical extent and area of the ridges, freeboard and floe draft, and the size of floes and melt ponds. The drag coefficients are incorporated into the Los Alamos Sea Ice Model (CICE) and show the influence of the new drag parameterization on the motion and state of the ice cover, with the most noticeable being a depletion of sea ice over the west boundary of the Arctic Ocean and over the Beaufort Sea. The new parameterization allows the drag coefficients to be coupled to the sea ice state and therefore to evolve spatially and temporally. It is found that the range of values predicted for the drag coefficients agree with the range of values measured in several regions of the Arctic. Finally, the implications of the new form drag formulation for the spinup or spindown of the Arctic Ocean are discussed.

The impact of variable sea ice roughness on changes in Arctic Ocean surface stress: a model study

The Arctic sea ice cover is thinning and retreating, causing changes in surface roughness that in turn modify the momentum flux from the atmosphere through the ice into the ocean. New model simulations comprising variable sea ice drag coefficients for both the air and water interface demonstrate that the heterogeneity in sea ice surface roughness significantly impacts the spatial distribution and trends of ocean surface stress during the last decades. Simulations with constant sea ice drag coefficients as used in most climate models show an increase in annual mean ocean surface stress (0.003 N/m2 per decade, 4.6%) due to the reduction of ice thickness leading to a weakening of the ice and accelerated ice drift. In contrast, with variable drag coefficients our simulations show annual mean ocean surface stress is declining at a rate of -0.002 N/m2 per decade (3.1%) over the period 1980-2013 because of a significant reduction in surface roughness associated with an increasingly thinner and younger sea ice cover. The effectiveness of sea ice in transferring momentum does not only depend on its resistive strength against the wind forcing but is also set by its top and bottom surface roughness varying with ice types and ice conditions. This reveals the need to account for sea ice surface roughness variations in climate simulations in order to correctly represent the implications of sea ice loss under global warming.

An analytical model of mountain wave drag for wind profiles with shear and curvature

An analytical model is developed to predict the surface drag exerted by internal gravity waves on an isolated axisymmetric mountain over which there is a stratified flow with a velocity profile that varies relatively slowly with height. The model is linear with respect to the perturbations induced by the mountain, and solves the Taylor–Goldstein equation with variable coefficients using a Wentzel–Kramers–Brillouin (WKB) approximation, formally valid for high Richardson numbers, Ri. The WKB solution is extended to a higher order than in previous studies, enabling a rigorous treatment of the effects of shear and curvature of the wind profile on the surface drag. In the hydrostatic approximation, closed formulas for the drag are derived for generic wind profiles, where the relative magnitude of the corrections to the leading-order drag (valid for a constant wind profile) does not depend on the detailed shape of the orography. The drag is found to vary proportionally to Ri21, decreasing as Ri decreases for a wind that varies linearly with height, and increasing as Ri decreases for a wind that rotates with height maintaining its magnitude. In these two cases the surface drag is predicted to be aligned with the surface wind. When one of the wind components varies linearly with height and the other is constant, the surface drag is misaligned with the surface wind, especially for relatively small Ri. All these results are shown to be in fairly good agreement with numerical simulations of mesoscale nonhydrostatic models, for high and even moderate values of Ri.

Reactive uptake coefficients for heterogeneous reaction of N2O5 with submicron aerosols of NaCl and natural sea salt

The kinetics of uptake of gaseous N2O5 on submicron aerosols containing NaCl and natural sea salt have been investigated in a flow reactor as a function of relative humidity (RH) in the range 30-80% at 295±2K and a total pressure of 1bar. The measured uptake coefficients, γ, were larger on the aerosols containing sea salt compared to those of pure NaCl, and in both cases increased with increasing RH. These observations are explained in terms of the variation in the size of the salt droplets, which leads to a limitation in the uptake rate into small particles. After correction for this effect the uptake coefficients are independent of relative humidity, and agree with those measured previously on larger droplets. A value of γ=0.025 is recommended for the reactive uptake coefficient for N2O5 on deliquesced sea salt droplets at 298K and RH>40%.

Parameter estimation using a particle method: inferring mixing coefficients from sea level observations

This paper presents a first attempt to estimate mixing parameters from sea level observations using a particle method based on importance sampling. The method is applied to an ensemble of 128 members of model simulations with a global ocean general circulation model of high complexity. Idealized twin experiments demonstrate that the method is able to accurately reconstruct mixing parameters from an observed mean sea level field when mixing is assumed to be spatially homogeneous. An experiment with inhomogeneous eddy coefficients fails because of the limited ensemble size. This is overcome by the introduction of local weighting, which is able to capture spatial variations in mixing qualitatively. As the sensitivity of sea level for variations in mixing is higher for low values of mixing coefficients, the method works relatively well in regions of low eddy activity.

CISK or WISHE as the mechanism for tropical cyclone intensification

Examination of conditional instability of the second kind (CISK) and wind-induced surface heat exchange (WISHE), two proposed mechanisms for tropical cyclone and polar low intensification, suggests that the sensitivity of the intensification rate of these disturbances to surface properties, such as surface friction and moisture supply, will be different for the two mechanisms. These sensitivities were examined by perturbing the surface characteristics in a numerical model with explicit convection. The intensification rate was found to have a strong positive dependence on the heat and moisture transfer coefficients, while remaining largely insensitive to the frictional drag coefficient. CISK does not predict the observed dependence of vortex intensification rate on the heat and moisture transfer coefficients, nor the insensitivity to the frictional drag coefficient since it anticipates that intensification rate is controlled by frictional convergence in the boundary layer. Since neither conditional instability nor boundary moisture content showed any significant sensitivity to the transfer coefficients, this is true of CISK using both the convective closures of Ooyama and of Charney and Eliassen. In comparison, the WISHE intensification mechanism does predict the observed increase in intensification rate with heat and moisture transfer coefficients, while not anticipating a direct influence from surface friction.

Time-varying coefficients in the almost ideal demand system: An empirical appraisal

This paper provides a generalisation of the structural time series version of the Almost Ideal Demand System (AIDS) that allows for time-varying coefficients (TVC/AIDS) in the presence of cross-equation constraints. An empirical appraisal of the TVC/AIDS is made using a dynamic AIDS with trending intercept as the baseline model with a data set from the Italian Household Budget Survey (1986-2001). The assessment is based on four criteria: adherence to theoretical constraints, statistical diagnostics on residuals, forecasting performance and economic meaningfulness. No clear evidence is found for superior performance of the TVC/AIDS, apart from improved short-term forecasts.

Octanol-water partition coefficients for predicting the effects of tannins in ruminant nutrition

Tannins can cause beneficial or harmful nutritional effects, but their great diversity has until now prevented a rational distinction between tannin structures and their nutritional responses. An attempt has been made to study this problem by examining the octanol-water solubilities of tannins. A relatively simple HPLC method has been developed for screening mixtures of plant tannins for their octanol-water partition coefficients (K-ow coefficients). Tannins were isolated from the fruits and leaves of different Acacia, Calliandra, Dichrostachys, and Piliostigma species, which are known to produce beneficial or harmful effects. The K-ow coefficients of these tannins ranged from 0.061 to 13.9, average coefficients of variation were 9.2% and recoveries were 107%. Acacia nilotica fruits and leaves had the highest K-ow coefficients, that is, 2.0 and 13.9, respectively. These A. nilotica products also have high concentrations of tannins. The combined effects of high octanol solubilities and high tannin concentrations may explain their negative effects on animal nutrition and health. It is known that compounds with high octanol solubilities are more easily absorbed into tissues, and it is, therefore, proposed that such compounds are more likely to cause toxicity problems especially if consumed in large quantities. According to the literature, tannins in human foods tend to have low K-ow coefficients, and this was confirmed for the tannins in Piliostigma thonningii fruits. Therefore, unconventional feeds or browse products should be screened not only for their tannin concentrations but also for low octanol-water partition coefficients in order to identify nutritionally safe feeds and to avoid potentially toxic feeds.

New magnetic nano-absorbent for the determination of n-octanol/water partition coefficients

A novel and generic miniaturization methodology for the determination of partition coefficient values of organic compounds in noctanol/water by using magnetic nanoparticles is, for the first time, described. We have successfully designed, synthesised and characterised new colloidal stable porous silica-encapsulated magnetic nanoparticles of controlled dimensions. These nanoparticles absorbing a tiny amount of n-octanol in their porous silica over-layer are homogeneously dispersed into a bulk aqueous phase (pH 7.40) containing an organic compound prior to magnetic separation. The small size of the particles and the efficient mixing allow a rapid establishment of the partition equilibrium of the organic compound between the solid supported n-octanol nano-droplets and the bulk aqueous phase. UV-vis spectrophotometry is then applied as a quantitative method to determine the concentration of the organic compound in the aqueous phase both before and after partitioning (after magnetic separation). log D values of organic compounds of pharmaceutical interest (0.65-3.50), determined by this novel methodology, were found to be in excellent agreement with the values measured by the shake-flask method in two independent laboratories, which are also consistent with the literature data. It was also found that this new technique gives a number of advantages such as providing an accurate measurement of log D value, a much shorter experimental time and a smaller sample size required. With this approach, the formation of a problematic emulsion, commonly encountered in shake-flask experiments, is eliminated. It is envisaged that this method could be applicable to the high throughput log D screening of drug candidates. (c) 2005 Elsevier B.V. All rights reserved.

Gas-phase rate coefficients for reactions of NO3, OH, O-3 and O(P-3) with unsaturated alcohols and ethers: correlations and structure-activity relations (SARs)

Experimental difficulties sometimes force modellers to use predicted rate coefficients for reactions of oxygenated volatile organic compounds (oVOCs). We examine here methods for making the predictions for reactions of atmospheric initiators of oxidation, NO3, OH, O-3 and O(P-3), with unsaturated alcohols and ethers. Logarithmic correlations are found between measured rate coefficients and calculated orbital energies, and these correlations may be used directly to estimate rate coefficients for compounds where measurements have not been performed. To provide a shortcut that obviates the need to calculate orbital energies, structure-activity relations (SARs) are developed. Our SARs are tested for predictive power against compounds for which experimental rate coefficients exist, and their accuracy is discussed. Estimated atmospheric lifetimes for oVOCs are presented. The SARs for alkenols successfully predict key rate coefficients, and thus can be used to enhance the scope of atmospheric models incorporating detailed chemistry. SARs for the ethers have more limited applicability, but can still be useful in improving tropospheric models. (C) 2008 Elsevier Ltd. All rights reserved.

Gas-phase rate coefficients for the reactions of NO3, OH and O-3 with alpha,beta-unsaturated esters and ketones: structure-activity relations (SARs)

Gas-phase rate coefficients for the atmospherically important reactions of NO3, OH and O-3 are predicted for 55 alpha,beta-unsaturated esters and ketones. The rate coefficients were calculated using a correlation described previously [Pfrang, C., King, M.D., C. E. Canosa-Mas, C.E., Wayne, R.P., 2006. Atmospheric Environment 40, 1170-1179]. These rate coefficients were used to extend structure-activity relations for predicting the rate coefficients for the reactions of NO3, OH or O-3 with alkenes to include alpha,beta-unsaturated esters and ketones. Conjugation of an alkene with an alpha,beta-keto or alpha,beta-ester group will reduce the value of a rate coefficient by a factor of similar to 110, similar to 2.5 and similar to 12 for reaction with NO3, OH or O-3, respectively. The actual identity of the alkyl group, R, in -C(O)R or -C(O)OR has only a small influence. An assessment of the reliability of the SAR is given that demonstrates that it is useful for reactions involving NO3 and OH, but less valuable for those of O-3 or peroxy nitrate esters. (c) 2006 Elsevier Ltd. All rights reserved.