Decadal climate prediction: an update from the trenches

This paper provides an update on research in the relatively new and fast-moving field of decadal climate prediction, and addresses the use of decadal climate predictions not only for potential users of such information but also for improving our understanding of processes in the climate system. External forcing influences the predictions throughout, but their contributions to predictive skill become dominant after most of the improved skill from initialization with observations vanishes after about six to nine years. Recent multi-model results suggest that there is relatively more decadal predictive skill in the North Atlantic, western Pacific, and Indian Oceans than in other regions of the world oceans. Aspects of decadal variability of SSTs, like the mid-1970s shift in the Pacific, the mid-1990s shift in the northern North Atlantic and western Pacific, and the early-2000s hiatus, are better represented in initialized hindcasts compared to uninitialized simulations. There is evidence of higher skill in initialized multi-model ensemble decadal hindcasts than in single model results, with multi-model initialized predictions for near term climate showing somewhat less global warming than uninitialized simulations. Some decadal hindcasts have shown statistically reliable predictions of surface temperature over various land and ocean regions for lead times of up to 6-9 years, but this needs to be investigated in a wider set of models. As in the early days of El Niño-Southern Oscillation (ENSO) prediction, improvements to models will reduce the need for bias adjustment, and increase the reliability, and thus usefulness, of decadal climate predictions in the future.

Inhomogeneous two-dimensional turbulence in the atmosphere

Faced with the strongly nonlinear and apparently random behaviour of the energy-containing scales in the atmosphere, geophysical fluid dynamicists have attempted to understand the synoptic-scale atmospheric flow within the context of two-dimensional homogeneous turbulence theory (e.g. FJØRTOFT [1]; LEITH [2]). However atmospheric observations (BOER and SHEPHERD [3] and Fig.1) show that the synoptic-scale transient flow evolves in the presence of a planetary-scale, quasi-stationary background flow which is approximately zonal (east-west). Classical homogeneous 2-D turbulence theory is therefore not strictly applicable to the transient flow. One is led instead to study 2-D turbulence in the presence of a large-scale (barotropically stable) zonal jet inhomogeneity.

Large-scale two-dimensional turbulence in the atmosphere

Global FGGE data are used to investigate several aspects of large-scale turbulence in the atmosphere. The approach follows that for two-dimensional, nondivergent turbulent flows which are homogeneous and isotropic on the sphere. Spectra of kinetic energy, enstrophy and available potential energy are obtained for both the stationary and transient parts of the flow. Nonlinear interaction terms and fluxes of energy and enstrophy through wavenumber space are calculated and compared with the theory. A possible method of parameterizing the interactions with unresolved scales is considered. Two rather different flow regimes are found in wavenumber space. The high-wavenumber regime is dominated by the transient components of the flow and exhibits, at least approximately, several of the conditions characterizing homogeneous and isotropic turbulence. This region of wavenumber space also displays some of the features of an enstrophy-cascading inertial subrange. The low-wavenumber region, on the other hand, is dominated by the stationary component of the flow, exhibits marked anisotropy and, in contrast to the high-wavenumber regime, displays a marked change between January and July.

Using rain-use efficiency to explore livestock production trends in rangelands in the Transkei, South Africa

We qualitatively describe the condition of communally managed rangelands in the Transkei, South Africa, using GIS and high resolution near-infrared imagery. Using livestock census data from 28 magisterial districts in the Transkei, we explored the trends in livestock biomass from 1923–1998. The area had been subjected to intensive herbivory by domestic livestock during that period, and the high livestock biomass had been blamed for the perceived degradation or ‘overgrazing’ of the region. Our assessment used the concept rain-use efficiency (RUE) (kg dry matter ha–1 mm–1) to determine whether there is evidence of change in the efficiency of the system to produce domestic livestock. We calculated RUE from annual livestock numbers and the mean annual rainfall for each district. We found no evidence of a decline in rain-use efficiency between the two assessment periods (1923–1944, 1945–1998). There was evidence of a shift in the ratio of sheep to goats between 1923 and 1998, with goat numbers increasing (greater than twofold) relative to sheep in eight districts. This trend may be associated with changes in the structure of vegetation. We conclude that this region is not showing evidence of system run down that affects domestic livestock production.

Northern Hemisphere atmospheric blocking as simulated by 15 atmospheric general circulation models in the period 1979-1988

As a part of the Atmospheric Model Intercomparison Project (AMIP), the behaviour of 15 general circulation models has been analysed in order to diagnose and compare the ability of the different models in simulating Northern Hemisphere midlatitude atmospheric blocking. In accordance with the established AMIP procedure, the 10-year model integrations were performed using prescribed, time-evolving monthly mean observed SSTs spanning the period January 1979–December 1988. Atmospheric observational data (ECMWF analyses) over the same period have been also used to verify the models results. The models involved in this comparison represent a wide spectrum of model complexity, with different horizontal and vertical resolution, numerical techniques and physical parametrizations, and exhibit large differences in blocking behaviour. Nevertheless, a few common features can be found, such as the general tendency to underestimate both blocking frequency and the average duration of blocks. The problem of the possible relationship between model blocking and model systematic errors has also been assessed, although without resorting to ad-hoc numerical experimentation it is impossible to relate with certainty particular model deficiencies in representing blocking to precise parts of the model formulation.

Some results from an intercomparison of the climates simulated by 14 atmospheric general circulation models

Some climatological information from 14 atmospheric general circulation models is presented and compared in order to assess the ability of a broad group of models to simulate current climate. The quantities considered are cross sections of temperature, zonal wind, and meridional stream function together with latitudinal distributions of mean sea level pressure and precipitation rate. The nature of the deficiencies in the simulated climates that are common to all models and those which differ among models is investigated; the general improvement in the ability of models to simulate certain aspects of the climate is shown; consideration is given to the effect of increasing resolution on simulated climate; and approaches to understanding and reducing model deficiencies are discussed. The information presented here is a subset of a more voluminous compilation which is available in report form (Boer et al., 1991). This report contains essentially the same text, but results from all 14 models are presented together with additional results in the form of geographical distributions of surface variables and certain difference statistics.

An intercomparison of the climates simulated by 14 atmospheric general circulation models

Climatological information from fourteen atmospheric general circulation models is presented and compared in order to assess the ability of a broad group of models to simulate current climate. The quantities considered are cross sections of temperature, zonal wind and meridional stream function together with latitudinal distributions of mean sea-level pressure and precipitation rate. The nature of the deficiencies in the simulated climates that are common to all models and those which differ among models is investigated, general improvement in the ability of models to simulate certain aspects of the climate is shown, consideration is given to the effect of increasing resolution on simulated climate and approaches to the understanding and reduction of model deficiencies are discussed.

Northern Hemisphere atmospheric blocking as simulated by 15 atmospheric general circulation models in the period 1979-1988 (Results from an AMIP diagnostic subproject)

As a part of the Atmospheric Model Intercomparison Project (AMIP), the behaviour of 15 general circulation models has been analysed in order to diagnose and compare the ability of the different models in simulating Northern Hemisphere midlatitude atmospheric blocking. In accordance with the established AMIP procedure, the 10-year model integrations were performed using prescribed, time-evolving monthly mean observed SSTs spanning the period January 1979–December 1988. Atmospheric observational data (ECMWF analyses) over the same period have been also used to verify the models results. The models involved in this comparison represent a wide spectrum of model complexity, with different horizontal and vertical resolution, numerical techniques and physical parametrizations, and exhibit large differences in blocking behaviour. Nevertheless, a few common features can be found, such as the general tendency to underestimate both blocking frequency and the average duration of blocks. The problem of the possible relationship between model blocking and model systematic errors has also been assessed, although without resorting to ad-hoc numerical experimentation it is impossible to relate with certainty particular model deficiencies in representing blocking to precise parts of the model formulation.

Variation in genes related to hepatic lipid metabolism and changes in waist circumference and body weight.

We analysed single nucleotide polymorphisms (SNPs) tagging the genetic variability of six candidate genes (ATF6, FABP1, LPIN2, LPIN3, MLXIPL and MTTP) involved in the regulation of hepatic lipid metabolism, an important regulatory site of energy balance for associations with body mass index (BMI) and changes in weight and waist circumference. We also investigated effect modification by sex and dietary intake. Data of 6,287 individuals participating in the European prospective investigation into cancer and nutrition were included in the analyses. Data on weight and waist circumference were followed up for 6.9 ± 2.5 years. Association of 69 tagSNPs with baseline BMI and annual changes in weight as well as waist circumference were investigated using linear regression analysis. Interactions with sex, GI and intake of carbohydrates, fat as well as saturated, monounsaturated and polyunsaturated fatty acids were examined by including multiplicative SNP-covariate terms into the regression model. Neither baseline BMI nor annual weight or waist circumference changes were significantly associated with variation in the selected genes in the entire study population after correction for multiple testing. One SNP (rs1164) in LPIN2 appeared to be significantly interacting with sex (p = 0.0003) and was associated with greater annual weight gain in men (56.8 ± 23.7 g/year per allele, p = 0.02) than in women (-25.5 ± 19.8 g/year per allele, p = 0.2). With respect to gene-nutrient interaction, we could not detect any significant interactions when accounting for multiple testing. Therefore, out of our six candidate genes, LPIN2 may be considered as a candidate for further studies.

Association of vitamin D status with arterial blood pressure and hypertension risk: a mendelian randomisation study

BACKGROUND: Low plasma 25-hydroxyvitamin D (25[OH]D) concentration is associated with high arterial blood pressure and hypertension risk, but whether this association is causal is unknown. We used a mendelian randomisation approach to test whether 25(OH)D concentration is causally associated with blood pressure and hypertension risk. METHODS: In this mendelian randomisation study, we generated an allele score (25[OH]D synthesis score) based on variants of genes that affect 25(OH)D synthesis or substrate availability (CYP2R1 and DHCR7), which we used as a proxy for 25(OH)D concentration. We meta-analysed data for up to 108 173 individuals from 35 studies in the D-CarDia collaboration to investigate associations between the allele score and blood pressure measurements. We complemented these analyses with previously published summary statistics from the International Consortium on Blood Pressure (ICBP), the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and the Global Blood Pressure Genetics (Global BPGen) consortium. FINDINGS: In phenotypic analyses (up to n=49 363), increased 25(OH)D concentration was associated with decreased systolic blood pressure (β per 10% increase, -0·12 mm Hg, 95% CI -0·20 to -0·04; p=0·003) and reduced odds of hypertension (odds ratio [OR] 0·98, 95% CI 0·97-0·99; p=0·0003), but not with decreased diastolic blood pressure (β per 10% increase, -0·02 mm Hg, -0·08 to 0·03; p=0·37). In meta-analyses in which we combined data from D-CarDia and the ICBP (n=146 581, after exclusion of overlapping studies), each 25(OH)D-increasing allele of the synthesis score was associated with a change of -0·10 mm Hg in systolic blood pressure (-0·21 to -0·0001; p=0·0498) and a change of -0·08 mm Hg in diastolic blood pressure (-0·15 to -0·02; p=0·01). When D-CarDia and consortia data for hypertension were meta-analysed together (n=142 255), the synthesis score was associated with a reduced odds of hypertension (OR per allele, 0·98, 0·96-0·99; p=0·001). In instrumental variable analysis, each 10% increase in genetically instrumented 25(OH)D concentration was associated with a change of -0·29 mm Hg in diastolic blood pressure (-0·52 to -0·07; p=0·01), a change of -0·37 mm Hg in systolic blood pressure (-0·73 to 0·003; p=0·052), and an 8·1% decreased odds of hypertension (OR 0·92, 0·87-0·97; p=0·002). INTERPRETATION: Increased plasma concentrations of 25(OH)D might reduce the risk of hypertension. This finding warrants further investigation in an independent, similarly powered study.

Effect of fall-grazed sericea lespedeza (Lespedeza cuneata) on gastrointestinal nematode infections of growing goats

High prevalence of anthelmintic-resistant gastrointestinal nematodes (GIN) in goats has increased pressure to find effective, alternative non-synthetic control methods, one of which is adding forage of the high condensed tannin (CT) legume sericea lespedeza (SL; Lespedeza cuneata) to the animal's diet. Previous work has demonstrated good efficacy of dried SL (hay, pellets) against small ruminant GIN, but information is lacking on consumption of fresh SL, particularly during the late summer–autumn period in the southern USA when perennial warm-season grass pastures are often low in quality. A study was designed to determine the effects of autumn (September–November) consumption of fresh SL forage, grass pasture (predominantly bermudagrass, BG; Cynodon dactylon), or a combination of SL + BG forage by young goats [intact male Spanish kids, 9 months old (20.7 ± 1.1 kg), n = 10/treatment group] on their GIN infection status. Three forage paddocks (0.40 ha) were set up at the Fort Valley State University Agricultural Research Station (Fort Valley, GA) for an 8-week trial. The goats in each paddock were supplemented with a commercial feed pellet at 0.45 kg/head/d for the first 4 weeks of the trial, and 0.27 kg/head/d for the final 4 weeks. Forage samples taken at the start of the trial were analyzed for crude protein (CP), neutral detergent fiber (NDF), and acid detergent fiber (ADF) content, and a separate set of SL samples was analyzed for CT in leaves, stems, and whole plant using the benzyl mercaptan thiolysis method. Animal weights were taken at the start and end of the trial, and fecal and blood samples were collected weekly for determination of fecal egg counts (FEC) and packed cell volume (PCV), respectively. Adult GIN was recovered from the abomasum and small intestines of all goats at the end of the experiment for counting and speciation. The CP levels were highest for SL forage, intermediate for SL + BG, and lowest for BG forage samples, while NDF and ADF values were the opposite, with highest levels in BG and lowest in SL forage samples. Sericea lespedeza leaves had more CT than stems (16.0 g vs. 3.3 g/100 g dry weight), a slightly higher percentage of PDs (98% vs. 94%, respectively) and polymers of larger mean degrees of polymerization (42 vs. 18, respectively). There were no differences in average daily gain or blood PCV between the treatment groups, but SL goats had lower FEC (P < 0.05) than the BG or SL + BG forage goats throughout most of the trial. The SL + BG goats had lower FEC than the BG forage animals by the end of the trial (week 8, P < 0.05). The SL goats had lower numbers (P < 0.05) of male Haemonchus contortus and tended to have fewer female (P < 0.10) and total (P < 0.07) H. contortus compared with the BG goats. The predominant GIN in all the goats was Trichostrongylus colubriformis (73% of total GIN). As a low-input forage with activity against pathogenic GIN (H. contortus), SL has a potential to reduce producers’ dependence upon synthetic anthelmintics and also to fill the autumn ‘window’ in good-quality fresh forages for goat grazing in the southern USA.

Association between Mediterranean and Nordic diet scores and changes in weight and waist circumference: influence of FTO and TCF7L2 loci

BACKGROUND: Several studies have shown that adherence to the Mediterranean Diet measured by using the Mediterranean diet score (MDS) is associated with lower obesity risk. The newly proposed Nordic Diet could hold similar beneficial effects. Because of the increasing focus on the interaction between diet and genetic predisposition to adiposity, studies should consider both diet and genetics. OBJECTIVE: We investigated whether FTO rs9939609 and TCF7L2 rs7903146 modified the association between the MDS and Nordic diet score (NDS) and changes in weight (Δweight), waist circumference (ΔWC), and waist circumference adjusted for body mass index (BMI) (ΔWCBMI). DESIGN: We conducted a case-cohort study with a median follow-up of 6.8 y that included 11,048 participants from 5 European countries; 5552 of these subjects were cases defined as individuals with the greatest degree of unexplained weight gain during follow-up. A randomly selected subcohort included 6548 participants, including 5496 noncases. Cases and noncases were compared in analyses by using logistic regression. Continuous traits (ie, Δweight, ΔWC, and ΔWCBMI) were analyzed by using linear regression models in the random subcohort. Interactions were tested by including interaction terms in models. RESULTS: A higher MDS was significantly inversely associated with case status (OR: 0.98; 95% CI: 0.96, 1.00), ΔWC (β = -0.010 cm/y; 95% CI: -0.020, -0.001 cm/y), and ΔWCBMI (β = -0.008; 95% CI:-0.015, -0.001) per 1-point increment but not Δweight (P = 0.53). The NDS was not significantly associated with any outcome. There was a borderline significant interaction between the MDS and TCF7L2 rs7903146 on weight gain (P = 0.05), which suggested a beneficial effect of the MDS only in subjects who carried 1 or 2 risk alleles. FTO did not modify observed associations. CONCLUSIONS: A high MDS is associated with a lower ΔWC and ΔWCBMI, regardless of FTO and TCF7L2 risk alleles. For Δweight, findings were less clear, but the effect may depend on the TCF7L2 rs7903146 variant. The NDS was not associated with anthropometric changes during follow-up.