The longitudinal variation of equatorial waves due to propagation on a varying zonal flow

The general 1-D theory of waves propagating on a zonally varying flow is developed from basic wave theory, and equations are derived for the variation of wavenumber and energy along ray paths. Different categories of behaviour are found, depending on the sign of the group velocity (cg) and a wave property, B. For B positive the wave energy and the wave number vary in the same sense, with maxima in relative easterlies or westerlies, depending on the sign of cg. Also the wave accumulation of Webster and Chang (1988) occurs where cg goes to zero. However for B negative they behave in opposite senses and wave accumulation does not occur. The zonal propagation of the gravest equatorial waves is analysed in detail using the theory. For non-dispersive Kelvin waves, B reduces to 2, and analytic solution is possible. B is positive for all the waves considered, except for the westward moving mixed Rossby-gravity (WMRG) wave which can have negative as well as positive B. Comparison is made between the observed climatologies of the individual equatorial waves and the result of pure propagation on the climatological upper tropospheric flow. The Kelvin wave distribution is in remarkable agreement, considering the approximations made. Some aspects of the WMRG and Rossby wave distributions are also in qualitative agreement. However the observed maxima in these waves in the winter westerlies in the eastern Pacific and Atlantic are not consistent with the theory. This is consistent with the importance of the sources of equatorial waves in these westerly duct regions due to higher latitude wave activity.

Dynamical trust and reputation computation model for B2C E-Commerce

Trust is one of the most important factors that influence the successful application of network service environments, such as e-commerce, wireless sensor networks, and online social networks. Computation models associated with trust and reputation have been paid special attention in both computer societies and service science in recent years. In this paper, a dynamical computation model of reputation for B2C e-commerce is proposed. Firstly, conceptions associated with trust and reputation are introduced, and the mathematical formula of trust for B2C e-commerce is given. Then a dynamical computation model of reputation is further proposed based on the conception of trust and the relationship between trust and reputation. In the proposed model, classical varying processes of reputation of B2C e-commerce are discussed. Furthermore, the iterative trust and reputation computation models are formulated via a set of difference equations based on the closed-loop feedback mechanism. Finally, a group of numerical simulation experiments are performed to illustrate the proposed model of trust and reputation. Experimental results show that the proposed model is effective in simulating the dynamical processes of trust and reputation for B2C e-commerce.

Crop yield reduction in the tropics under climate change: Processes and uncertainties

Many modelling studies examine the impacts of climate change on crop yield, but few explore either the underlying bio-physical processes, or the uncertainty inherent in the parameterisation of crop growth and development. We used a perturbed-parameter crop modelling method together with a regional climate model (PRECIS) driven by the 2071-2100 SRES A2 emissions scenario in order to examine processes and uncertainties in yield simulation. Crop simulations used the groundnut (i.e. peanut; Arachis hypogaea L.) version of the General Large-Area Model for annual crops (GLAM). Two sets of GLAM simulations were carried out: control simulations and fixed-duration simulations, where the impact of mean temperature on crop development rate was removed. Model results were compared to sensitivity tests using two other crop models of differing levels of complexity: CROPGRO, and the groundnut model of Hammer et al. [Hammer, G.L., Sinclair, T.R., Boote, K.J., Wright, G.C., Meinke, H., and Bell, M.J., 1995, A peanut simulation model: I. Model development and testing. Agron. J. 87, 1085-1093]. GLAM simulations were particularly sensitive to two processes. First, elevated vapour pressure deficit (VPD) consistently reduced yield. The same result was seen in some simulations using both other crop models. Second, GLAM crop duration was longer, and yield greater, when the optimal temperature for the rate of development was exceeded. Yield increases were also seen in one other crop model. Overall, the models differed in their response to super-optimal temperatures, and that difference increased with mean temperature; percentage changes in yield between current and future climates were as diverse as -50% and over +30% for the same input data. The first process has been observed in many crop experiments, whilst the second has not. Thus, we conclude that there is a need for: (i) more process-based modelling studies of the impact of VPD on assimilation, and (ii) more experimental studies at super-optimal temperatures. Using the GLAM results, central values and uncertainty ranges were projected for mean 2071-2100 crop yields in India. In the fixed-duration simulations, ensemble mean yields mostly rose by 10-30%. The full ensemble range was greater than this mean change (20-60% over most of India). In the control simulations, yield stimulation by elevated CO2 was more than offset by other processes-principally accelerated crop development rates at elevated, but sub-optimal, mean temperatures. Hence, the quantification of uncertainty can facilitate relatively robust indications of the likely sign of crop yield changes in future climates. (C) 2007 Elsevier B.V. All rights reserved.

On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy

In this paper, the available potential energy (APE) framework of Winters et al. (J. Fluid Mech., vol. 289, 1995, p. 115) is extended to the fully compressible Navier– Stokes equations, with the aims of clarifying (i) the nature of the energy conversions taking place in turbulent thermally stratified fluids; and (ii) the role of surface buoyancy fluxes in the Munk & Wunsch (Deep-Sea Res., vol. 45, 1998, p. 1977) constraint on the mechanical energy sources of stirring required to maintain diapycnal mixing in the oceans. The new framework reveals that the observed turbulent rate of increase in the background gravitational potential energy GPEr , commonly thought to occur at the expense of the diffusively dissipated APE, actually occurs at the expense of internal energy, as in the laminar case. The APE dissipated by molecular diffusion, on the other hand, is found to be converted into internal energy (IE), similar to the viscously dissipated kinetic energy KE. Turbulent stirring, therefore, does not introduce a new APE/GPEr mechanical-to-mechanical energy conversion, but simply enhances the existing IE/GPEr conversion rate, in addition to enhancing the viscous dissipation and the entropy production rates. This, in turn, implies that molecular diffusion contributes to the dissipation of the available mechanical energy ME =APE +KE, along with viscous dissipation. This result has important implications for the interpretation of the concepts of mixing efficiency γmixing and flux Richardson number Rf , for which new physically based definitions are proposed and contrasted with previous definitions. The new framework allows for a more rigorous and general re-derivation from the first principles of Munk & Wunsch (1998, hereafter MW98)’s constraint, also valid for a non-Boussinesq ocean: G(KE) ≈ 1 − ξ Rf ξ Rf Wr, forcing = 1 + (1 − ξ )γmixing ξ γmixing Wr, forcing , where G(KE) is the work rate done by the mechanical forcing, Wr, forcing is the rate of loss of GPEr due to high-latitude cooling and ξ is a nonlinearity parameter such that ξ =1 for a linear equation of state (as considered by MW98), but ξ <1 otherwise. The most important result is that G(APE), the work rate done by the surface buoyancy fluxes, must be numerically as large as Wr, forcing and, therefore, as important as the mechanical forcing in stirring and driving the oceans. As a consequence, the overall mixing efficiency of the oceans is likely to be larger than the value γmixing =0.2 presently used, thereby possibly eliminating the apparent shortfall in mechanical stirring energy that results from using γmixing =0.2 in the above formula.

Alternative methodologies - stretching the in vitro box

Current gas-based in vitro evaluation systems are extremely powerful research techniques. However they have the potential to generate a great deal more than simple fermentation dynamics. Details from four experiments are presented in which adaptation, and novel application, of an in vitro system allowed widely differing objectives to be examined. In the first two studies, complement methodologies were utilised. In such assays, an activity or outcome is inferred through the occurrence of a secondary event rather than by direct observation. Using an N-deficient incubation medium, the increase in starch fermentation, when supplemented with individual amino acids (i.e., known level of N) relative to that of urea (i.e., known quantity and N availability), provided an estimate of their microbial utilisation. Due to the low level of response observed with some arnino acids (notably methionine and lysine), it was concluded, that they may not need to be offered in a rumen-inert form to escape rumen microbial degradation. In another experiment, the extent to which degradation of plant cell wall components was inhibited by lipid supplementation was evaluated using fermentation gas release profiles of washed hay. The different responses due to lipid source and level of inclusion suggested that the degree of rumen protection required to ameliorate this depression was supplement dependent. That in vitro inocula differ in their microbial composition is of little interest per se, as long as the outcome is the same (i.e., that similar substrates are degraded at comparable rates and end-product release is equivalent). However where a microbial population is deficient in a particular activity, increasing the level of inoculation will have no benefit. Estimates of hydrolytic activity were obtained by examining fermentation kinetics of specific substrates. A number of studies identified a fundamental difference between rumen fluid and faecal inocula, with the latter having a lower fibrolytic activity, which could not be completely attributed to microbial numbers. The majority of forage maize is offered as an ensiled feed, however most of the information on which decisions such as choice of variety, crop management and harvesting date are made is based on fresh crop measurements. As such, an attempt was made to estimate ensiled maize quality from an in vitro analysis of the fresh crop. Fermentation profiles and chemical analysis confirmed changes in crop composition over the growing season, and loss of labile carbohydrates during ensiling. In addition, examination of degradation residues allowed metabolizable energy (ME) contents to be estimated. Due to difficulties associated with starch analysis, the observation that this parameter could be predicted by difference (together with an assumed degradability), allowed an estimate of ensiled maize ME to be developed from fresh material. In addition, the contribution of the main carbohydrates towards ME showed the importance of delaying harvest until maximum starch content has been achieved. (c) 2005 Elsevier B.V. All rights reserved.

A Riemann solver for compressible flow of a real gas in a Lagrangian frame

An approximate Riemann solver is presented for the compressible flow equations with a general (convex) equation of state in a Lagrangian frame of reference.

Cylindrically symmetric flows in a Lagrangian frame

An approximate Riemann solver, in a Lagrangian frame of reference, is presented for the compressible flow equations with cylindrical and spherical symmetry, including flow in a duct of variable cross section. The scheme is applied to a cylindrically symmetric problem involving the interaction of shocks.

Retrieval of global atmospheric electrical activity at a polluted urban site

The atmospheric electrical Potential Gradient (PG) arises from global thunderstorm activity, but surface measurements of the atmospheric Potential Gradient (PG) are influenced by global thunderstorms and local aerosol concentration changes. The local aerosol change can be monitored independently, and in some cases the concentration changes are closely related to PG changes. For these circumstances, a general theory to remove the local aerosol influence on PG measurements has been developed. Continuous measurements of PG and aerosol mass concentration were made during 24–31 Dec, 2005 within an urban environment at Reading, UK. The average diurnal variation of PG showed a double diurnal cycle, with maxima in the early morning and evening hours. The aerosol concentration has similar double maxima. Removing the aerosol using from the PG and aerosol correlation returns a single diurnal cycle, suggestive of the more global PG diurnal cycle.

Generalised Dirichelt-to-Neumann map in time dependent domains

We study the heat, linear Schrodinger and linear KdV equations in the domain l(t) < x < ∞, 0 < t < T, with prescribed initial and boundary conditions and with l(t) a given differentiable function. For the first two equations, we show that the unknown Neumann or Dirichlet boundary value can be computed as the solution of a linear Volterra integral equation with an explicit weakly singular kernel. This integral equation can be derived from the formal Fourier integral representation of the solution. For the linear KdV equation we show that the two unknown boundary values can be computed as the solution of a system of linear Volterra integral equations with explicit weakly singular kernels. The derivation in this case makes crucial use of analyticity and certain invariance properties in the complex spectral plane. The above Volterra equations are shown to admit a unique solution.

Drag produced by trapped lee waves and propagating mountain waves in a two-layer atmosphere

The surface drag force produced by trapped lee waves and upward propagating waves in non-hydrostatic stratified flow over a mountain ridge is explicitly calculated using linear theory for a two-layer atmosphere with piecewise-constant static stability and wind speed profiles. The behaviour of the drag normalized by its hydrostatic single-layer reference value is investigated as a function of the ratio of the Scorer parameters in the two layers l_2/l_1 and of the corresponding dimensionless interface height l_1 H, for selected values of the dimensionless ridge width l_1 a and ratio of wind speeds in the two layers. When l_2/l_1 → 1, the propagating wave drag approaches 1 in approximately hydrostatic conditions, and the trapped lee wave drag vanishes. As l_2/l_1 decreases, the propagating wave drag progressively displays an oscillatory behaviour with l_1 H, with maxima of increasing magnitude due to constructive interference of reflected waves in the lower layer. The trapped lee wave drag shows localized maxima associated with each resonant trapped lee wave mode, occurring for small l_2/l_1 and slightly higher values of l_1 H than the propagating wave drag maxima. As l1a decreases, i.e. the flow becomes more non-hydrostatic, the propagating wave drag decreases and the regions of non-zero trapped lee wave drag extend to higher l_2/l_1. These results are confirmed by numerical simulations for l_2/l_1 = 0.2. In parameter ranges of meteorological relevance, the trapped lee wave drag may have a magnitude comparable to that of propagating wave drag, and be larger than the reference single-layer drag. This may have implications for drag parametrization in global climate and weather-prediction models.

Limitations of small x resummation methods from F2data

We discuss several methods of calculating the DIS structure functions F2(x,Q2) based on BFKL-type small x resummations. Taking into account new HERA data ranging down to small xand low Q2, the pure leading order BFKL-based approach is excluded. Other methods based on high energy factorization are closer to conventional renormalization group equations. Despite several difficulties and ambiguities in combining the renormalization group equations with small x resummed terms, we find that a fit to the current data is hardly feasible, since the data in the low Q2 region are not as steep as the BFKL formalism predicts. Thus we conclude that deviations from the (successful) renormalization group approach towards summing up logarithms in 1/x are disfavoured by experiment.

Pulsating equilibria: stability through migration

This paper is to present a model of spatial equilibrium using a nonlinear generalization of Markov-chain type model, and to show the dynamic stability of a unique equilibrium. Even at an equilibrium, people continue to migrate among regions as well as among agent-types, and yet their overall distribution remain unchanged. The model is also adapted to suggest a theory of traffic distribution in a city.

A moving mesh approach for modelling avascular tumour growth

A key step in many numerical schemes for time-dependent partial differential equations with moving boundaries is to rescale the problem to a fixed numerical mesh. An alternative approach is to use a moving mesh that can be adapted to focus on specific features of the model. In this paper we present and discuss two different velocity-based moving mesh methods applied to a two-phase model of avascular tumour growth formulated by Breward et al. (2002) J. Math. Biol. 45(2), 125-152. Each method has one moving node which tracks the moving boundary. The first moving mesh method uses a mesh velocity proportional to the boundary velocity. The second moving mesh method uses local conservation of volume fraction of cells (masses). Our results demonstrate that these moving mesh methods produce accurate results, offering higher resolution where desired whilst preserving the balance of fluxes and sources in the governing equations.

Rapid changes in the redox conditions of the western Tethys Ocean during the early Aptian oceanic anoxic event

The early Aptian (125 to 121 Ma) records an episode of severe environmental change including a major perturbation of the carbon cycle, an oceanic anoxic event (OAE 1a, 122.5 Ma), a platform drowning episode and a biocalcification crisis. We propose to trace changes in the oxygenation state of the ocean during the early Aptian anoxic event using the redox-sensitive trace-element (RSTE) distribution, phosphorus accumulation rates (PARs) and organic-matter characterization in three different basins of the western Tethys. The following sections have been investigated: Gorgo a Cerbara (central Italy) in the Umbria Marche basin, Glaise (SE France) in the Vocontian basin and Cassis/La Bédoule (SE France) located in the Provencal basin. In the Gorgo a Cerbara section, RSTE distributions show a low background level along the main part of the section, contrasted by different maxima in concentrations within the Selli level. In the Glaise section, the Goguel level displays a weak increase in RSTE contents coeval with moderate TOC values. At Cassis/La Bédoule, no significant RSTE enrichments have been observed in sediments equivalent to the Selli level. These differences in the records of the geochemical proxies of the Selli level or its equivalent indicate the deposition under different redox conditions, probably related to the paleogeography. Our data indicate the development of anoxic–euxinic conditions in the deeper part of the Tethys during OAE 1a, whereas in the shallower environments, conditions were less reducing. Moreover, at Gorgo a Cerbara, the Selli level is characterized by rapid changes in the intensity of reducing conditions in the water column. Ocean eutrophication seems to be a major factor in the development and the persistence of anoxia as suggested by the PAR evolution. Higher PAR values at the onset of OAE 1a suggest an increase in nutrient input, whereas the return to lower values through the first part of the OAE 1a interval may be related to the weakened capacity to retain P in the sedimentary reservoir due to bottom-water oxygen depletion. This general pattern is contrasted by the data of Gorgo a Cerbara, where the sediments deposited during the OAE 1a interval show P-enrichments (mainly authigenic P). This is associated with maxima in TOC values and Corg:Ptot ratios, suggesting that a part of the remobilized P was trapped in the sediments and as such prevented from returning to the water column.

Re-examining the U.K.'s greatest tornado outbreak: forecasting the limited extent of tornadoes along a cold front

On 23 November 1981, a strong cold front swept across the U.K., producing tornadoes from the west to the east coasts. An extensive campaign to collect tornado reports by the Tornado and Storm Research Organisation (TORRO) resulted in 104 reports, the largest U.K. outbreak. The front was simulated with a convection-permitting numerical model down to 200-m horizontal grid spacing to better understand its evolution and meteorological environment. The event was typical of tornadoes in the U.K., with convective available potential energy (CAPE) less than 150 J kg-1, 0-1-km wind shear of 10-20 m s-1, and a narrow cold-frontal rainband forming precipitation cores and gaps. A line of cyclonic absolute vorticity existed along the front, with maxima as large as 0.04 s-1. Some hook-shaped misovortices bore kinematic similarity to supercells. The narrow swath along which the line was tornadic was bounded on the equatorward side by weak vorticity along the line and on the poleward side by zero CAPE, enclosing a region where the environment was otherwise favorable for tornadogenesis. To determine if the 104 tornado reports were plausible, first possible duplicate reports were eliminated, resulting in as few as 58 tornadoes to as many as 90. Second, the number of possible parent misovortices that may have spawned tornadoes is estimated from model output. The number of plausible tornado reports in the 200-m grid-spacing domain was 22 and as many as 44, whereas the model simulation was used to estimate 30 possible parent misovortices within this domain. These results suggest that 90 reports was plausible.