A Hamiltonian weak-wave model for shallow-water flow

A reduced dynamical model is derived which describes the interaction of weak inertia–gravity waves with nonlinear vortical motion in the context of rotating shallow–water flow. The formal scaling assumptions are (i) that there is a separation in timescales between the vortical motion and the inertia–gravity waves, and (ii) that the divergence is weak compared to the vorticity. The model is Hamiltonian, and possesses conservation laws analogous to those in the shallow–water equations. Unlike the shallow–water equations, the energy invariant is quadratic. Nonlinear stability theorems are derived for this system, and its linear eigenvalue properties are investigated in the context of some simple basic flows.

Non-ergodicity of inviscid two-dimensional flow on a beta-plane and on the surface of a rotating sphere

It is shown that, for a sufficiently large value of β, two-dimensional flow on a doubly-periodic beta-plane cannot be ergodic (phase-space filling) on the phase-space surface of constant energy and enstrophy. A corresponding result holds for flow on the surface of a rotating sphere, for a sufficiently rapid rotation rate Ω. This implies that the higher-order, non-quadratic invariants are exerting a significant influence on the statistical evolution of the flow. The proof relies on the existence of a finite-amplitude Liapunov stability theorem for zonally symmetric basic states with a non-vanishing absolute-vorticity gradient. When the domain size is much larger than the size of a typical eddy, then a sufficient condition for non-ergodicity is that the wave steepness ε < 1, where ε = 2[surd radical]2Z/βU in the planar case and $\epsilon = 2^{\frac{1}{4}} a^{\frac{5}{2}}Z^{\frac{7}{4}}/\Omega U^{\frac{5}{2}}$ in the spherical case, and where Z is the enstrophy, U the r.m.s. velocity, and a the radius of the sphere. This result may help to explain why numerical simulations of unforced beta-plane turbulence (in which ε decreases in time) seem to evolve into a non-ergodic regime at large scales.

A statistical model for sea surface diurnal warming driven by numercial weather predictions fluxes and winds

A statistical model is derived relating the diurnal variation of sea surface temperature (SST) to the net surface heat flux and surface wind speed from a numerical weather prediction (NWP) model. The model is derived using fluxes and winds from the European Centre for Medium-Range Weather Forecasting (ECMWF) NWP model and SSTs from the Spinning Enhanced Visible and Infrared Imager (SEVIRI). In the model, diurnal warming has a linear dependence on the net surface heat flux integrated since (approximately) dawn and an inverse quadratic dependence on the maximum of the surface wind speed in the same period. The model coefficients are found by matching, for a given integrated heat flux, the frequency distributions of the maximum wind speed and the observed warming. Diurnal cooling, where it occurs, is modelled as proportional to the integrated heat flux divided by the heat capacity of the seasonal mixed layer. The model reproduces the statistics (mean, standard deviation, and 95-percentile) of the diurnal variation of SST seen by SEVIRI and reproduces the geographical pattern of mean warming seen by the Advanced Microwave Scanning Radiometer (AMSR-E). We use the functional dependencies in the statistical model to test the behaviour of two physical model of diurnal warming that display contrasting systematic errors.

Retrievals of sea surface temperature from infrared imagery: origin and form of systematic errors

We show that retrievals of sea surface temperature from satellite infrared imagery are prone to two forms of systematic error: prior error (familiar from the theory of atmospheric sounding) and error arising from nonlinearity. These errors have different complex geographical variations, related to the differing geographical distributions of the main geophysical variables that determine clear-sky brightness-temperatures over the oceans. We show that such errors arise as an intrinsic consequence of the form of the retrieval (rather than as a consequence of sub-optimally specified retrieval coefficients, as is often assumed) and that the pattern of observed errors can be simulated in detail using radiative-transfer modelling. The prior error has the linear form familiar from atmospheric sounding. A quadratic equation for nonlinearity error is derived, and it is verified that the nonlinearity error exhibits predominantly quadratic behaviour in this case.

A finite element method for nonlinear elliptic problems

We present a Galerkin method with piecewise polynomial continuous elements for fully nonlinear elliptic equations. A key tool is the discretization proposed in Lakkis and Pryer, 2011, allowing us to work directly on the strong form of a linear PDE. An added benefit to making use of this discretization method is that a recovered (finite element) Hessian is a byproduct of the solution process. We build on the linear method and ultimately construct two different methodologies for the solution of second order fully nonlinear PDEs. Benchmark numerical results illustrate the convergence properties of the scheme for some test problems as well as the Monge–Amp`ere equation and the Pucci equation.

A modified PNN algorithm with optimal PD modeling using the orthogonal least squares method

In this paper a modified algorithm is suggested for developing polynomial neural network (PNN) models. Optimal partial description (PD) modeling is introduced at each layer of the PNN expansion, a task accomplished using the orthogonal least squares (OLS) method. Based on the initial PD models determined by the polynomial order and the number of PD inputs, OLS selects the most significant regressor terms reducing the output error variance. The method produces PNN models exhibiting a high level of accuracy and superior generalization capabilities. Additionally, parsimonious models are obtained comprising a considerably smaller number of parameters compared to the ones generated by means of the conventional PNN algorithm. Three benchmark examples are elaborated, including modeling of the gas furnace process as well as the iris and wine classification problems. Extensive simulation results and comparison with other methods in the literature, demonstrate the effectiveness of the suggested modeling approach.

Forecast encompassing tests and probability forecasts

We consider tests of forecast encompassing for probability forecasts, for both quadratic and logarithmic scoring rules. We propose test statistics for the null of forecast encompassing, present the limiting distributions of the test statistics, and investigate the impact of estimating the forecasting models' parameters on these distributions. The small-sample performance is investigated, in terms of small numbers of forecasts and model estimation sample sizes. We show the usefulness of the tests for the evaluation of recession probability forecasts from logit models with different leading indicators as explanatory variables, and for evaluating survey-based probability forecasts.

The Brownian traveller on manifolds

We study the inuence of the intrinsic curvature on the large time behaviour of the heat equation in a tubular neighbourhood of an unbounded geodesic in a two-dimensional Riemannian manifold. Since we consider killing boundary conditions, there is always an exponential-type decay for the heat semigroup. We show that this exponential-type decay is slower for positively curved manifolds comparing to the at case. As the main result, we establish a sharp extra polynomial-type decay for the heat semigroup on negatively curved manifolds comparing to the at case. The proof employs the existence of Hardy-type inequalities for the Dirichlet Laplacian in the tubular neighbourhoods on negatively curved manifolds and the method of self-similar variables and weighted Sobolev spaces for the heat equation.

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.

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr – Part 2: A new reconstruction of the interplanetary magnetic field

We present a new reconstruction of the interplanetary magnetic field (IMF, B) for 1846–2012 with a full analysis of errors, based on the homogeneously constructed IDV(1d)composite of geomagnetic activity presented in Part 1 (Lockwood et al., 2013a). Analysis of the dependence of the commonly used geomagnetic indices on solar wind parameters is presented which helps explain why annual means of interdiurnal range data, such as the new composite, depend only on the IMF with only a very weak influence of the solar wind flow speed. The best results are obtained using a polynomial (rather than a linear) fit of the form B = χ · (IDV(1d) − β)α with best-fit coefficients χ = 3.469, β = 1.393 nT, and α = 0.420. The results are contrasted with the reconstruction of the IMF since 1835 by Svalgaard and Cliver (2010).

Approximation by harmonic polynomials in star-shaped domains and exponential convergence of Trefftz hp-dGFEM

We study the approximation of harmonic functions by means of harmonic polynomials in two-dimensional, bounded, star-shaped domains. Assuming that the functions possess analytic extensions to a delta-neighbourhood of the domain, we prove exponential convergence of the approximation error with respect to the degree of the approximating harmonic polynomial. All the constants appearing in the bounds are explicit and depend only on the shape-regularity of the domain and on delta. We apply the obtained estimates to show exponential convergence with rate O(exp(−b square root N)), N being the number of degrees of freedom and b>0, of a hp-dGFEM discretisation of the Laplace equation based on piecewise harmonic polynomials. This result is an improvement over the classical rate O(exp(−b cubic root N )), and is due to the use of harmonic polynomial spaces, as opposed to complete polynomial spaces.

The effect of adding stinging nettle (Urtica dioica) haylage to a total mixed ration on performance and rumen function of lactating dairy cows

In vitro studies found that inclusion of dried stinging nettle (Urtica dioica) at 100 mg/g dry matter (DM) increased the pH of a rumen fluid inoculated fermentation buffer by 30% and the effect was persistent for 7 days. Our objective was to evaluate the effects of adding stinging nettle haylage to a total mixed ration on feed intake, eating and rumination activity, rumen pH, milk yield, and milk composition of lactating dairy cows. Six lactating Holstein-Friesian cows were used in a replicated 3 × 3 Latin Square design experiment with 3 treatments and 3 week periods. Treatments were a control (C) high-starch (311 g/kg DM) total mixed ration diet and two treatment diets containing 50 (N5) and 100 (N10) g nettle haylage (DM/kg) as a replacement for ryegrass silage (Lolium perenne). There was an increase (linear, P < 0.010) in the proportion of large particles and a reduction in medium (linear, P = 0.045) and fine particles (linear, P = 0.026) in the diet offered with increasing nettle inclusion. A numerical decrease (linear, P = 0.106) in DM intake (DMI) was observed as nettle inclusion in the diet increased. Milk yield averaged 20.3 kg/day and was not affected by diet. There was a decrease (quadratic, P = 0.01) in the time animals spent ruminating as nettle inclusion in the diet increased, in spite of an increase in the number of boli produced daily for the N5 diet (quadratic, P = 0.031). Animals fed the N10 diet spent less time with a rumen pH below 5.5 (P < 0.05) than cows fed the N5 diet. Averaged over an 8.5 h sampling period, there were no changes in the concentration or proportions of acetate or propionate in the rumen, but feeding nettle haylage reduced the concentrations of n-butyrate (quadratic, P < 0.001), i-butyrate (linear, P < 0.009) and n-caproate (linear, P < 0.003). Milk and fat and protein corrected milk yield were not affected when nettles replaced ryegrass silage in the diet of lactating dairy cows, despite a numerical reduction in feed intake. Rumination activity was reduced by the addition of nettle haylage to the diet, which may reflect differences in fibre structure between the nettle haylage and ryegrass silage fed. Changes observed in rumen pH suggest potential benefits of feeding nettle haylage for reducing rumen acidosis. However, the extent to which these effects were due to the fermentability and structure of the nettle haylage compared to the ryegrass silage fed, or a bioactive component of the nettles, is not certain

The effects of a soft drink tax in the UK

The majority of the UK population is either overweight or obese. Health economists, nutritionists and doctors are calling for the UK to follow the example of other European countries and introduce a tax on soft drinks as a result of the perception that high intakes contribute to diet-related disease. We use a demand model estimated with household-level data on beverage purchases in the UK to investigate the effects of a tax on soft drink consumption. The model is a Quadratic Almost Ideal Demand System, and censoring is handled by applying a double hurdle. Separate models are estimated for low, moderate and high consumers to allow for a differential impact on consumption between these groups. Applying different hypothetical tax rates, we conclude that understanding the nature of substitute/complement relationships is crucial in designing an effective policy as these relationships differ between consumers depending on their consumption level. The overall impact of a soft drink tax on calorie consumption is likely to be small.

A comparison of GOCE and drifter-based estimates of the North Atlantic steady-state surface circulation

Over the last decade, due to the Gravity Recovery And Climate Experiment (GRACE) mission and, more recently, the Gravity and steady state Ocean Circulation Explorer (GOCE) mission, our ability to measure the ocean’s mean dynamic topography (MDT) from space has improved dramatically. Here we use GOCE to measure surface current speeds in the North Atlantic and compare our results with a range of independent estimates that use drifter data to improve small scales. We find that, with filtering, GOCE can recover 70% of the Gulf Steam strength relative to the best drifter-based estimates. In the subpolar gyre the boundary currents obtained from GOCE are close to the drifter-based estimates. Crucial to this result is careful filtering which is required to remove small-scale errors, or noise, in the computed surface. We show that our heuristic noise metric, used to determine the degree of filtering, compares well with the quadratic sum of mean sea surface and formal geoid errors obtained from the error variance–covariance matrix associated with the GOCE gravity model. At a resolution of 100 km the North Atlantic mean GOCE MDT error before filtering is 5 cm with almost all of this coming from the GOCE gravity model.

How reward and emotional stimuli induce different reactions across the menstrual cycle

Despite widespread belief that moods are affected by the menstrual cycle, researchers on emotion and reward have not paid much attention to the menstrual cycle until recently. However, recent research has revealed different reactions to emotional stimuli and to rewarding stimuli across the different phases of the menstrual cycle. The current paper reviews the emerging literature on how ovarian hormone fluctuation during the menstrual cycle modulates reactions to emotional stimuli and to reward. Behavioral and neuroimaging studies in humans suggest that estrogen and progesterone have opposing influences. That is, it appears that estrogen enhances reactions to reward, but progesterone counters the facilitative effects of estrogen and decreases reactions to rewards. In contrast, reactions to emotionally arousing stimuli (particularly negative stimuli) appear to be decreased by estrogen but enhanced by progesterone. Potential factors that can modulate the effects of the ovarian hormones (e.g., an inverse quadratic function of hormones’ effects; the structural changes of the hippocampus across the menstrual cycle) are also discussed.