Runge-Kutta IMEX schemes for the Horizontally Explicit/Vertically Implicit (HEVI) solution of wave equations

Many operational weather forecasting centres use semi-implicit time-stepping schemes because of their good efficiency. However, as computers become ever more parallel, horizontally explicit solutions of the equations of atmospheric motion might become an attractive alternative due to the additional inter-processor communication of implicit methods. Implicit and explicit (IMEX) time-stepping schemes have long been combined in models of the atmosphere using semi-implicit, split-explicit or HEVI splitting. However, most studies of the accuracy and stability of IMEX schemes have been limited to the parabolic case of advection–diffusion equations. We demonstrate how a number of Runge–Kutta IMEX schemes can be used to solve hyperbolic wave equations either semi-implicitly or HEVI. A new form of HEVI splitting is proposed, UfPreb, which dramatically improves accuracy and stability of simulations of gravity waves in stratified flow. As a consequence it is found that there are HEVI schemes that do not lose accuracy in comparison to semi-implicit ones. The stability limits of a number of variations of trapezoidal implicit and some Runge–Kutta IMEX schemes are found and the schemes are tested on two vertical slice cases using the compressible Boussinesq equations split into various combinations of implicit and explicit terms. Some of the Runge–Kutta schemes are found to be beneficial over trapezoidal, especially since they damp high frequencies without dropping to first-order accuracy. We test schemes that are not formally accurate for stiff systems but in stiff limits (nearly incompressible) and find that they can perform well. The scheme ARK2(2,3,2) performs the best in the tests.

Optimising the FAMOUS climate model: inclusion of global carbon cycling

FAMOUS fills an important role in the hierarchy of climate models, both explicitly resolving atmospheric and oceanic dynamics yet being sufficiently computationally efficient that either very long simulations or large ensembles are possible. An improved set of carbon cycle parameters for this model has been found using a perturbed physics ensemble technique. This is an important step towards building the "Earth System" modelling capability of FAMOUS, which is a reduced resolution, and hence faster running, version of the Hadley Centre Climate model, HadCM3. Two separate 100 member perturbed parameter ensembles were performed; one for the land surface and one for the ocean. The land surface scheme was tested against present-day and past representations of vegetation and the ocean ensemble was tested against observations of nitrate. An advantage of using a relatively fast climate model is that a large number of simulations can be run and hence the model parameter space (a large source of climate model uncertainty) can be more thoroughly sampled. This has the associated benefit of being able to assess the sensitivity of model results to changes in each parameter. The climatologies of surface and tropospheric air temperature and precipitation are improved relative to previous versions of FAMOUS. The improved representation of upper atmosphere temperatures is driven by improved ozone concentrations near the tropopause and better upper level winds.

Implicit learning and imperceptible influence: Syncretic literacy of multilingual Chinese children.

This article reports on an ethnographic study involving the literacy practices of two multilingual Chinese children from two similar yet different cultural and linguistic contexts: Montreal and Singapore. Using syncretism as a theoretical tool, this inquiry examines how family environment and support facilitate children’s process of becoming literate in multiple languages. Informed by sociocultural theory, the inquiry looks in particular at the role of grandparents in the syncretic literacy practices of children. Through comparative analysis, the study reveals similarities and differences that, when considered together, contribute to our understanding of multilingual children’s creative forms of learning with regard to their rich literacy resources in multiple languages, the imperceptible influences of mediators, various learning styles and syncretic literacy practices.

Punishing the foreigner: implicit discrimination in the Premier League based on oppositional identity

We present the first empirical study to reveal the presence of implicit discrimination in a non-experimental setting. By using a large dataset of in-match data in the English Premier League, we show that white referees award significantly more yellow cards against non-white players of oppositional identity. We argue that this is the result of implicit discrimination by showing that this discriminatory behaviour: (i) increases in how rushed the referee is before making a decision, and (ii) it increases in the level of ambiguity of the decision. The variation in (i) and (ii) cannot be explained by any form of conscious discrimination such as taste-based or statistical discrimination. Moreover, we show that oppositional identity players do not differ in their behaviour from other players along several dimensions related to aggressiveness and style of play providing further evidence that this is not statistical discrimination.

Explicit and implicit subject bias in the ABS journal quality guide

This paper addresses issues raised in two recent papers published in this journal about the UK Association of Business Schools' Journal Quality Guide (ABS Guide). While much of the debate about journal rankings in general, and the ABS Guide in particular, has focused on the construction, power and (mis)use of these rankings, this paper differs in that it explains and provides evidence about explicit and implicit biases in the ABS Guide. In so doing, it poses potentially difficult questions that the editors of the ABS Guide need to address and urgently rectify if the ABS Guide seeks to build and retain legitimacy. In particular, the evidence in this paper shows explicit bias in the ABS Guide against several subject areas, including accounting and finance. It also shows implicit bias against accounting and finance when comparing journal rankings in sub-areas shared between accounting and finance and the broader business management subject areas.

Effect of inclusion of lactulose and Lactobacillus plantarum on the intestinal environment and performance of piglets at weaning

The prebiotic lactulose, a probiotic strain of Lactobacillus plantarum (L. plantarum) and a synbiotic combination of these two agents were evaluated as growth promoters in 25–39-day old commercial weaning pigs. Ninety-six weaning pigs were allocated into 32 pens, taking initial weight into account, and distributed into four groups as follows: a control diet (CTR), the same diet supplemented daily with L. plantarum (109 CFU/mL sprayed on top; 20 mL/pig) (LPN); 10 g/kg lactulose (LAC) or a combination of both treatments (SYN). At day 14, eight piglets from each group were euthanized and proximal colon digesta was sampled for luminal pH, short-chain fatty acids (SCFA) and lactic acid concentrations. Deoxyribonucleic acid was extracted from colonic digesta and the microbial community was profiled by terminal restriction fragment length polymorphism analysis (T-RFLP) and qPCR. Blood urea nitrogen (BUN) and acute-phase proteins (Pig-MAP) were measured. Lactulose treatment (LAC) improved feed intake (P<0.05), average daily gain (P<0.01), feed:gain ratio (P<0.05) and reduced BUN (P<0.01). Both, LAC and LPN treatment, decreased the Enterobacteriaceae:Lactobacillus spp. ratio in the colonic luminal contents (P<0.05). Moreover LPN treatment promoted a decrease in the percentage of branched fatty acids (P<0.01) suggesting a reduction in proteolytic microbial activity. Microbial profiling of colonic luminal contents by T-RFLP revealed changes in some microbial species. Terminal restriction fragments (TRFs) compatible with Bifidobacterium thermoacidophilum were more frequently detected in experimental diets compared to CTR (P<0.05). Pigs receiving SYN diet demonstrated the combined positive effects of individual LAC and LPN treatment although we were not able to show a specific increase in the probiotic strain with the inclusion of lactulose. Collectively, these data suggest the combination of lactulose and L. plantarum acts as a complementary synbiotic, but not as a synergistic combination.

Differential effect of beetroot bread on postprandial DBP according to Glu298Asp polymorphism in the eNOS gene: a pilot study.

Our objective was to investigate whether the presence of Glu298Asp polymorphism in the endothelial NO synthase (eNOS) gene differentially affects the postprandial blood pressure response to dietary nitrate-rich beetroot bread. A randomised, single-blind, controlled, crossover acute pilot study was performed in 14 healthy men (mean age: 34±9 years) who were retrospectively genotyped for Glu298Asp polymorphism (7GG; T carriers 7). Volunteers were randomised to receive 200 g beetroot-enriched bread (1.1 mmol nitrate) or control bread (no beetroot; 0.01 mmol nitrate) on two separate occasions 10 days apart. Baseline and incremental area under the curve of blood pressure and NOx (nitrate/nitrite) were measured for a 6-h postprandial period. A treatment × genotype interaction was observed for diastolic blood pressure (P<0.02), which was significantly lower in T carriers (P<0.01) after consumption of beetroot bread compared with control bread. No significant differences were observed in the GG group. The beneficial diastolic blood pressure reduction was observed only in the T carriers of the Glu298Asp polymorphism in the eNOS gene after consumption of nitrate-rich beetroot bread. These data require confirmation in a larger population group.

The composite-tendency RAW filter in semi-implicit integrations

The time discretization in weather and climate models introduces truncation errors that limit the accuracy of the simulations. Recent work has yielded a method for reducing the amplitude errors in leapfrog integrations from first-order to fifth-order. This improvement is achieved by replacing the Robert--Asselin filter with the RAW filter and using a linear combination of the unfiltered and filtered states to compute the tendency term. The purpose of the present paper is to apply the composite-tendency RAW-filtered leapfrog scheme to semi-implicit integrations. A theoretical analysis shows that the stability and accuracy are unaffected by the introduction of the implicitly treated mode. The scheme is tested in semi-implicit numerical integrations in both a simple nonlinear stiff system and a medium-complexity atmospheric general circulation model, and yields substantial improvements in both cases. We conclude that the composite-tendency RAW-filtered leapfrog scheme is suitable for use in semi-implicit integrations.

Differential interference effects of negative emotional states on subsequent semantic and perceptual processing

Past studies have revealed that encountering negative events interferes with cognitive processing of subsequent stimuli. The present study investigates whether negative events affect semantic and perceptual processing differently. Presentation of negative pictures produced slower reaction times than neutral or positive pictures in tasks that require semantic processing, such as natural or man-made judgments about drawings of objects, commonness judgments about objects, and categorical judgments about pairs of words. In contrast, negative picture presentation did not slow down judgments in subsequent perceptual processing (e.g., color judgments about words, size judgments about objects). The subjective arousal level of negative pictures did not modulate the interference effects on semantic or perceptual processing. These findings indicate that encountering negative emotional events interferes with semantic processing of subsequent stimuli more strongly than perceptual processing, and that not all types of subsequent cognitive processing are impaired by negative events.

Differential brain activity during emotional versus nonemotional reversal learning

The ability to change an established stimulus–behavior association based on feedback is critical for adaptive social behaviors. This ability has been examined in reversal learning tasks, where participants first learn a stimulus–response association (e.g., select a particular object to get a reward) and then need to alter their response when reinforcement contingencies change. Although substantial evidence demonstrates that the OFC is a critical region for reversal learning, previous studies have not distinguished reversal learning for emotional associations from neutral associations. The current study examined whether OFC plays similar roles in emotional versus neutral reversal learning. The OFC showed greater activity during reversals of stimulus–outcome associations for negative outcomes than for neutral outcomes. Similar OFC activity was also observed during reversals involving positive outcomes. Furthermore, OFC activity is more inversely correlated with amygdala activity during negative reversals than during neutral reversals. Overall, our results indicate that the OFC is more activated by emotional than neutral reversal learning and that OFC's interactions with the amygdala are greater for negative than neutral reversal learning.

Coalescence of particles by differential sedimentation

We consider a three dimensional system consisting of a large number of small spherical particles, distributed in a range of sizes and heights (with uniform distribution in the horizontal direction). Particles move vertically at a size-dependent terminal velocity. They are either allowed to merge whenever they cross or there is a size ratio criterion enforced to account for collision efficiency. Such a system may be described, in mean field approximation, by the Smoluchowski kinetic equation with a differential sedimentation kernel. We obtain self-similar steady-state and time-dependent solutions to the kinetic equation, using methods borrowed from weak turbulence theory. Analytical results are compared with direct numerical simulations (DNS) of moving and merging particles, and a good agreement is found.

Numerical analyses of Runge–Kutta implicit–explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models

With the prospect of exascale computing, computational methods requiring only local data become especially attractive. Consequently, the typical domain decomposition of atmospheric models means horizontally-explicit vertically-implicit (HEVI) time-stepping schemes warrant further attention. In this analysis, Runge-Kutta implicit-explicit schemes from the literature are analysed for their stability and accuracy using a von Neumann stability analysis of two linear systems. Attention is paid to the numerical phase to indicate the behaviour of phase and group velocities. Where the analysis is tractable, analytically derived expressions are considered. For more complicated cases, amplification factors have been numerically generated and the associated amplitudes and phase diagnosed. Analysis of a system describing acoustic waves has necessitated attributing the three resultant eigenvalues to the three physical modes of the system. To do so, a series of algorithms has been devised to track the eigenvalues across the frequency space. The result enables analysis of whether the schemes exactly preserve the non-divergent mode; and whether there is evidence of spurious reversal in the direction of group velocities or asymmetry in the damping for the pair of acoustic modes. Frequency ranges that span next-generation high-resolution weather models to coarse-resolution climate models are considered; and a comparison is made of errors accumulated from multiple stability-constrained shorter time-steps from the HEVI scheme with a single integration from a fully implicit scheme over the same time interval. Two schemes, “Trap2(2,3,2)” and “UJ3(1,3,2)”, both already used in atmospheric models, are identified as offering consistently good stability and representation of phase across all the analyses. Furthermore, according to a simple measure of computational cost, “Trap2(2,3,2)” is the least expensive.

Curl free pressure gradients over orography in a solution of the fully compressible Euler equations with implicit treatment of acoustic and gravity waves

Steep orography can cause noisy solutions and instability in models of the atmosphere. A new technique for modelling flow over orography is introduced which guarantees curl free gradients on arbitrary grids, implying that the pressure gradient term is not a spurious source of vorticity. This mimetic property leads to better hydrostatic balance and better energy conservation on test cases using terrain following grids. Curl-free gradients are achieved by using the co-variant components of velocity over orography rather than the usual horizontal and vertical components. In addition, gravity and acoustic waves are treated implicitly without the need for mean and perturbation variables or a hydrostatic reference profile. This enables a straightforward description of the implicit treatment of gravity waves. Results are presented of a resting atmosphere over orography and the curl-free pressure gradient formulation is advantageous. Results of gravity waves over orography are insensitive to the placement of terrain-following layers. The model with implicit gravity waves is stable in strongly stratified conditions, with N∆t up to at least 10 (where N is the Brunt-V ̈ais ̈al ̈a frequency). A warm bubble rising over orography is simulated and the curl free pressure gradient formulation gives much more accurate results for this test case than a model without this mimetic property.

The inclusion of forage mixtures in the diet of growing dairy heifers: impacts on digestion, energy utilisation, and methane emissions

Intensive farming focusing on monoculture grass species to maximise forage production has led to a reduction in the extent and diversity of species-rich grasslands. However, plant communities with higher species number (richness) are a potential strategy for more sustainable production and mitigation of greenhouse gas (GHG) emissions. Research has indicated the need to understand opportunities that forage mixtures can offer sustainable ruminant production systems. The objective of the two experiments reported here were to evaluate multiple species forage mixtures in comparison to ryegrass-dominant pasture, when conserved or grazed, on digestion, energy utilisation, N excretion, and methane emissions by growing 10–15 month old heifers. Experiment 1 was a 4 × 4 Latin square design with five week periods. Four forage treatments of: (1) ryegrass (control); permanent pasture with perennial ryegrass (Lolium perenne); (2) clover; a ryegrass:red clover (Trifolium pratense) mixture; (3) trefoil; a ryegrass:birdsfoot trefoil (Lotus corniculatus) mixture; and (4) flowers; a ryegrass:wild flower mixture of predominately sorrel (Rumex acetosa), ox-eye daisy (Leucanthemum vulgare), yarrow (Achillea millefolium), knapweed (Centaurea nigra) and ribwort plantain (Plantago lanceolata), were fed as haylages to four dairy heifers. Measurements included digestibility, N excretion, and energy utilisation (including methane emissions measured in respiration chambers). Experiment 2 used 12 different dairy heifers grazing three of the same forage treatments used to make haylage in experiment 1 (ryegrass, clover and flowers) and methane emissions were estimated using the sulphur hexafluoride (SF6) tracer technique. Distribution of ryegrass to other species (dry matter (DM) basis) was approximately 70:30 (clover), 80:20 (trefoil), and 40:60 (flowers) for experiment 1. During the first and second grazing rotations (respectively) in experiment 2, perennial ryegrass accounted for 95 and 98% of DM in ryegrass, and 84 and 52% of DM in clover, with red clover accounting for almost all of the remainder. In the flowers mixture, perennial ryegrass was 52% of the DM in the first grazing rotation and only 30% in the second, with a variety of other flower species occupying the remainder. Across both experiments, compared to the forage mixtures (clover, trefoil and flowers), ryegrass had a higher crude protein (CP) content (P < 0.001, 187 vs. 115 g kg −1 DM) and DM intake (P < 0.05, 9.0 vs. 8.1 kg day −1). Heifers in experiment 1 fed ryegrass, compared to the forage mixtures, had greater total tract digestibility (g kg −1) of DM (DMD; P < 0.008, 713 vs. 641) and CP (CPD, P < 0.001, 699 vs. 475), and used more intake energy (%) for body tissue deposition (P < 0.05, 2.6 vs. −4.9). For both experiments, heifers fed flowers differed the most compared to the ryegrass control for a number of measurements. Compared to ryegrass, flowers had 40% lower CP content (P < 0.001, 113 vs. 187 g kg −1), 18% lower DMD (P < 0.01, 585 vs. 713 g kg −1), 42% lower CPD (P < 0.001, 407 vs. 699 g kg −1), and 10% lower methane yield (P < 0.05, 22.6 vs. 25.1 g kg −1 DM intake). This study has shown inclusion of flowers in forage mixtures resulted in a lower CP concentration, digestibility and intake. These differences were due in part to sward management and maturity at harvest. Further research is needed to determine how best to exploit the potential environmental benefits of forage mixtures in sustainable ruminant production systems.