El Nino-mean state-seasonal cycle interactions in a multi-model ensemble

The modelled El Nino-mean state-seasonal cycle interactions in 23 coupled ocean-atmosphere GCMs, including the recent IPCC AR4 models, are assessed and compared to observations and theory. The models show a clear improvement over previous generations in simulating the tropical Pacific climatology. Systematic biases still include too strong mean and seasonal cycle of trade winds. El Nino amplitude is shown to be an inverse function of the mean trade winds in agreement with the observed shift of 1976 and with theoretical studies. El Nino amplitude is further shown to be an inverse function of the relative strength of the seasonal cycle. When most of the energy is within the seasonal cycle, little is left for inter-annual signals and vice versa. An interannual coupling strength (ICS) is defined and its relation with the modelled El Nino frequency is compared to that predicted by theoretical models. An assessment of the modelled El Nino in term of SST mode (S-mode) or thermocline mode (T-mode) shows that most models are locked into a S-mode and that only a few models exhibit a hybrid mode, like in observations. It is concluded that several basic El Nino-mean state-seasonal cycle relationships proposed by either theory or analysis of observations seem to be reproduced by CGCMs. This is especially true for the amplitude of El Nino and is less clear for its frequency. Most of these relationships, first established for the pre-industrial control simulations, hold for the double and quadruple CO2 stabilized scenarios. The models that exhibit the largest El Nino amplitude change in these greenhouse gas (GHG) increase scenarios are those that exhibit a mode change towards a T-mode (either from S-mode to hybrid or hybrid to T-mode). This follows the observed 1976 climate shift in the tropical Pacific, and supports the-still debated-finding of studies that associated this shift to increased GHGs. In many respects, these models are also among those that best simulate the tropical Pacific climatology (ECHAM5/MPI-OM, GFDL-CM2.0, GFDL-CM2.1, MRI-CGM2.3.2, UKMO-HadCM3). Results from this large subset of models suggest the likelihood of increased El Nino amplitude in a warmer climate, though there is considerable spread of El Nino behaviour among the models and the changes in the subsurface thermocline properties that may be important for El Nino change could not be assessed. There are no clear indications of an El Nino frequency change with increased GHG.

Characteristics and variability of storm tracks in the north Pacific, Bering Sea, and Alaska

The North Pacific and Bering Sea regions represent loci of cyclogenesis and storm track activity. In this paper climatological properties of extratropical storms in the North Pacific/Bering Sea are presented based upon aggregate statistics of individual storm tracks calculated by means of a feature-tracking algorithm run using NCEP–NCAR reanalysis data from 1948/49 to 2008, provided by the NOAA/Earth System Research Laboratory and the Cooperative Institute for Research in Environmental Sciences, Climate Diagnostics Center. Storm identification is based on the 850-hPa relative vorticity field (ζ) instead of the often-used mean sea level pressure; ζ is a prognostic field, a good indicator of synoptic-scale dynamics, and is directly related to the wind speed. Emphasis extends beyond winter to provide detailed consideration of all seasons. Results show that the interseasonal variability is not as large during the spring and autumn seasons. Most of the storm variables—genesis, intensity, track density—exhibited a maxima pattern that was oriented along a zonal axis. From season to season this axis underwent a north–south shift and, in some cases, a rotation to the northeast. This was determined to be a result of zonal heating variations and midtropospheric moisture patterns. Barotropic processes have an influence in shaping the downstream end of storm tracks and, together with the blocking influence of the coastal orography of northwest North America, result in high lysis concentrations, effectively making the Gulf of Alaska the “graveyard” of Pacific storms. Summer storms tended to be longest in duration. Temporal trends tended to be weak over the study area. SST did not emerge as a major cyclogenesis control in the Gulf of Alaska.

Effects of fluctuating daily surface fluxes on the time-mean oceanic circulation

The effect of fluctuating daily surface fluxes on the time-mean oceanic circulation is studied using an empirical flux model. The model produces fluctuating fluxes resulting from atmospheric variability and includes oceanic feedbacks on the fluxes. Numerical experiments were carried out by driving an ocean general circulation model with three different versions of the empirical model. It is found that fluctuating daily fluxes lead to an increase in the meridional overturning circulation (MOC) of the Atlantic of about 1 Sv and a decrease in the Antarctic circumpolar current (ACC) of about 32 Sv. The changes are approximately 7% of the MOC and 16% of the ACC obtained without fluctuating daily fluxes. The fluctuating fluxes change the intensity and the depth of vertical mixing. This, in turn, changes the density field and thus the circulation. Fluctuating buoyancy fluxes change the vertical mixing in a non-linear way: they tend to increase the convective mixing in mostly stable regions and to decrease the convective mixing in mostly unstable regions. The ACC changes are related to the enhanced mixing in the subtropical and the mid-latitude Southern Ocean and reduced mixing in the high-latitude Southern Ocean. The enhanced mixing is related to an increase in the frequency and the depth of convective events. As these events bring more dense water downward, the mixing changes lead to a reduction in meridional gradient of the depth-integrated density in the Southern Ocean and hence the strength of the ACC. The MOC changes are related to more subtle density changes. It is found that the vertical mixing in a latitudinal strip in the northern North Atlantic is more strongly enhanced due to fluctuating fluxes than the mixing in a latitudinal strip in the South Atlantic. This leads to an increase in the density difference between the two strips, which can be responsible for the increase in the Atlantic MOC.

Evaluation of forest snow processes models (SnowMIP2)

Thirty‐three snowpack models of varying complexity and purpose were evaluated across a wide range of hydrometeorological and forest canopy conditions at five Northern Hemisphere locations, for up to two winter snow seasons. Modeled estimates of snow water equivalent (SWE) or depth were compared to observations at forest and open sites at each location. Precipitation phase and duration of above‐freezing air temperatures are shown to be major influences on divergence and convergence of modeled estimates of the subcanopy snowpack. When models are considered collectively at all locations, comparisons with observations show that it is harder to model SWE at forested sites than open sites. There is no universal “best” model for all sites or locations, but comparison of the consistency of individual model performances relative to one another at different sites shows that there is less consistency at forest sites than open sites, and even less consistency between forest and open sites in the same year. A good performance by a model at a forest site is therefore unlikely to mean a good model performance by the same model at an open site (and vice versa). Calibration of models at forest sites provides lower errors than uncalibrated models at three out of four locations. However, benefits of calibration do not translate to subsequent years, and benefits gained by models calibrated for forest snow processes are not translated to open conditions.

Recent oppositely directed trends in solar climate forcings and the global mean surface air temperature. II. Different reconstructions of the total solar irradiance variation and dependence on response time scale

We have previously placed the solar contribution to recent global warming in context using observations and without recourse to climate models. It was shown that all solar forcings of climate have declined since 1987. The present paper extends that analysis to include the effects of the various time constants with which the Earth’s climate system might react to solar forcing. The solar input waveform over the past 100 years is defined using observed and inferred galactic cosmic ray fluxes, valid for either a direct effect of cosmic rays on climate or an effect via their known correlation with total solar irradiance (TSI), or for a combination of the two. The implications, and the relative merits, of the various TSI composite data series are discussed and independent tests reveal that the PMOD composite used in our previous paper is the most realistic. Use of the ACRIM composite, which shows a rise in TSI over recent decades, is shown to be inconsistent with most published evidence for solar influences on pre-industrial climate. The conclusions of our previous paper, that solar forcing has declined over the past 20 years while surface air temperatures have continued to rise, are shown to apply for the full range of potential time constants for the climate response to the variations in the solar forcings.

Recent changes in solar outputs and the global mean surface temperature. III. Analysis of contributions to global mean air surface temperature rise

A multivariate fit to the variation in global mean surface air temperature anomaly over the past half century is presented. The fit procedure allows for the effect of response time on the waveform, amplitude and lag of each radiative forcing input, and each is allowed to have its own time constant. It is shown that the contribution of solar variability to the temperature trend since 1987 is small and downward; the best estimate is -1.3% and the 2sigma confidence level sets the uncertainty range of -0.7 to -1.9%. The result is the same if one quantifies the solar variation using galactic cosmic ray fluxes (for which the analysis can be extended back to 1953) or the most accurate total solar irradiance data composite. The rise in the global mean air surface temperatures is predominantly associated with a linear increase that represents the combined effects of changes in anthropogenic well-mixed greenhouse gases and aerosols, although, in recent decades, there is also a considerable contribution by a relative lack of major volcanic eruptions. The best estimate is that the anthropogenic factors contribute 75% of the rise since 1987, with an uncertainty range (set by the 2sigma confidence level using an AR(1) noise model) of 49–160%; thus, the uncertainty is large, but we can state that at least half of the temperature trend comes from the linear term and that this term could explain the entire rise. The results are consistent with the intergovernmental panel on climate change (IPCC) estimates of the changes in radiative forcing (given for 1961–1995) and are here combined with those estimates to find the response times, equilibrium climate sensitivities and pertinent heat capacities (i.e. the depth into the oceans to which a given radiative forcing variation penetrates) of the quasi-periodic (decadal-scale) input forcing variations. As shown by previous studies, the decadal-scale variations do not penetrate as deeply into the oceans as the longer term drifts and have shorter response times. Hence, conclusions about the response to century-scale forcing changes (and hence the associated equilibrium climate sensitivity and the temperature rise commitment) cannot be made from studies of the response to shorter period forcing changes.

Association between ovarian follicle development and pregnancy rates in dairy cows undergoing spontaneous oestrous cycles

Ovarian follicle development continues in a wave-like manner during the bovine oestrous cycle giving rise to variation in the duration of ovulatory follicle development. The objectives of the present study were to determine whether a relationship exists between the duration of ovulatory follicle development and pregnancy rates following artificial insemination (AI) in dairy cows undergoing spontaneous oestrous cycles, and to identify factors influencing follicle turnover and pregnancy rate and the relationship between these two variables. Follicle development was monitored by daily transrectal ultrasonography from 10 days after oestrus until the subsequent oestrus in 158 lactating dairy cows. The cows were artificially inseminated following the second observed oestrus and pregnancy was diagnosed 35 days later. The predominant pattern of follicle development was two follicle waves (74.7%) with three follicle waves in 22.1% of oestrous cycles and four or more follicle waves in 3.2% of oestrous cycles. The interval from ovulatory follicle emergence to oestrus (EOI) was 3 days longer (P < 0.0001) in cows with two follicle waves than in those with three waves. Ovulatory follicles from two-wave oestrous cycles grew more slowly but were approximately 2 mm larger (P < 0.0001) on the day of oestrus. Twin ovulations were observed in 14.2% of oestrous cycles and occurred more frequently (P < 0.001) in three-wave oestrous cycles; consequently EOI was shorter in cows with twin ovulations. Overall, 57.0% of the cows were diagnosed pregnant 35 days after AI. Linear logistic regression analysis revealed an inverse relationship between EOI and the proportion of cows diagnosed pregnant, among all cows (n = 158; P < 0.01) and amongst those with single ovulations (n = 145; P < 0.05). Mean EOI was approximately I day shorter (P < 0.01) in cows that became pregnant than in non-pregnant cows; however, pregnancy rates did not differ significantly among cows with different patterns of follicle development. These findings confirm and extend previous observations in pharmacologically manipulated cattle and show, for the first time, that in dairy cows undergoing spontaneous oestrous cycles, natural variation in the duration of post-emergence ovulatory follicle development has a significant effect on pregnancy rate, presumably reflecting variation in oocyte developmental competence.

Simulation of the impact of high temperature stress on annual crop yields

Brief periods of high temperature which occur near flowering can severely reduce the yield of annual crops such as wheat and groundnut. A parameterisation of this well-documented effect is presented for groundnut (i.e. peanut; Arachis hypogaeaL.). This parameterisation was combined with an existing crop model, allowing the impact of season-mean temperature, and of brief high-temperature episodes at various times near flowering, to be both independently and jointly examined. The extended crop model was tested with independent data from controlled environment experiments and field experiments. The impact of total crop duration was captured, with simulated duration being within 5% of observations for the range of season-mean temperatures used (20-28 degrees C). In simulations across nine differently timed high temperature events, eight of the absolute differences between observed and simulated yield were less than 10% of the control (no-stress) yield. The parameterisation of high temperature stress also allows the simulation of heat tolerance across different genotypes. Three parameter sets, representing tolerant, moderately sensitive and sensitive genotypes were developed and assessed. The new parameterisation can be used in climate change studies to estimate the impact of heat stress on yield. It can also be used to assess the potential for adaptation of cropping systems to increased temperature threshold exceedance via the choice of genotype characteristics. (c) 2005 Elsevier B.V. All rights reserved.

Simulation of crop yields using ERA-40: Limits to skill and nonstationarity in weather-yield relationships

Reanalysis data provide an excellent test bed for impacts prediction systems. because they represent an upper limit on the skill of climate models. Indian groundnut (Arachis hypogaea L.) yields have been simulated using the General Large-Area Model (GLAM) for annual crops and the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-yr reanalysis (ERA-40). The ability of ERA-40 to represent the Indian summer monsoon has been examined. The ability of GLAM. when driven with daily ERA-40 data, to model both observed yields and observed relationships between subseasonal weather and yield has been assessed. Mean yields "were simulated well across much of India. Correlations between observed and modeled yields, where these are significant. are comparable to correlations between observed yields and ERA-40 rainfall. Uncertainties due to the input planting window, crop duration, and weather data have been examined. A reduction in the root-mean-square error of simulated yields was achieved by applying bias correction techniques to the precipitation. The stability of the relationship between weather and yield over time has been examined. Weather-yield correlations vary on decadal time scales. and this has direct implications for the accuracy of yield simulations. Analysis of the skewness of both detrended yields and precipitation suggest that nonclimatic factors are partly responsible for this nonstationarity. Evidence from other studies, including data on cereal and pulse yields, indicates that this result is not particular to groundnut yield. The detection and modeling of nonstationary weather-yield relationships emerges from this study as an important part of the process of understanding and predicting the impacts of climate variability and change on crop yields.

Influence of foliar diseases and their control by fungicides on grain yield and quality in wheat

Twenty-eight field experiments on sandy-loam soils in the UK (1982-2003) are reviewed by relating the extension of the green area duration of the flag leaf (GLADF) by fungicides to effects on yield and quality of winter wheat. Over all experiments mean grain yield = 8.85t ha(-1) at 85% DM. With regards quality, mean values were: thousand grain weight (TGW) = 44.5 g; specific weight (SWT) = 76.9 kg hl(-1); crude protein concentration (CP (N x 5.7)) = 12.5 % DM; Hagberg falling number (HFN) = 285 s; and sodium dodecyl sulphate (SDS)-sedimentation volume = 69ml. For each day (d) that fungicides increased GLADF there were associated average increases in yield (0.144 1 ha(-1) d(-1), se 0.0049, df = 333), TGW (0.56 gd(-1), se = 0.017) and SWT (0.22 kg hl(-1) d(-1), se 0.011). Some curvature was evident in all these relationships. When GLADF was delayed beyond 700 degrees Cd after anthesis, as was possible in cool wet seasons, responses were curtailed, or less reliable. Despite this apparent terminal sink limitation, fungicide effects on sink size, eg endosperm cell numbers or maximum water mass per grain, were not prerequisites for large effects on grain yield, TGW or SWT. Fungicide effects on CP were variable. Although the average response of CP was negative (-0.029%DM/d; se = 0.00338), this depended on cultivar and disease controlled. Controlling biotrophs such as rusts, (Puccinia spp.) tended to increase CP, whereas controlling a more necrotrophic pathogen (Septoria tritici) usually reducedCP. Irrespective of pathogen controlled, delaying senescence of the flag leaf was associated with increased nitrogen yields in the grain (averaging 2.24 kg N ha-1 d(-1), se = 0.0848) due to both increased N uptake into the above ground crop, and also more efficient remobilisation of N from leaf laminas. When sulphur availability appeared to be adequate, fungicide x cultivar interactions were similar on S as for CP, although N:S ratios tended to decline (i.e. improve for bread making) when S. tritici was controlled. On average, SDS-sedimentation volume declined (-0. 18 ml/d, se = 0.027) with increased GLADF, broadly commensurate with the average effect on CP. Hagberg falling number decreased as fungicide increased GLADF (-2.73 s/d, se = 0.178), indicating an increase in alpha-amylase activity.

Effect of dietary vitamin A restriction on marbling and conjugated-linoleic acid (CLA) content in Holstein steers

To determine the effect of duration of dietary vitamin A restriction on site of fat deposition in growing cattle, 60 Holstein steers (BW = 218.4 ± 6.55 kg) were fed a diet based on high-moisture corn with 2,200 IU supplemental vitamin A/kg DM (C) or no supplemental vitamin A for a long (243 d; LR) or short (131 d; SR) restriction prior to harvest at 243 d. The SR steers were fed the C diet for the first 112 d. Steers were penned individually and fed for ad libitum intake. Jugular vein blood samples for serum retinol analysis were collected on d 1, 112, and 243. Carcass samples were collected for composition analysis. Subcutaneous fat samples were collected for fatty acid composition. Fat samples from the i.m. and s.c. depot were collected to measure adipocyte size and density. Feedlot performance (ADG, DMI, and G:F) was not affected (P > 0.05) by vitamin A restriction. On d 243, the i.m. fat content of the LM was 33% greater (P < 0.05) for LR than for SR and C steers (5.6 vs. 3.9 and 4.2% ether extract, respectively). Depth of back fat and KPH percentage were not affected (P = 0.44 and 0.80, respectively) by vitamin A restriction. Carcass weight, composition of edible carcass, and yield grade were similar among treatments (P > 0.10). Liver retinol (LR = 6.1, SR = 6.5, and C = 44.7 µg/g; P < 0.01) was reduced in LR and SR vs. C steers. On d 243, LR and SR steers had similar serum retinol concentrations, and these were lower (P < 0.01) than those of C steers (LR = 21.2, SR = 25.2, and C = 36.9 µg/dL). Intramuscular adipose cellularity (adipocyte/mm2 and mean adipocyte diameter) on d 112 and d 243 was not affected (P > 0.10) by vitamin A restriction. Restricting vitamin A intake for 243 d increased i.m fat percentage without affecting s.c. or visceral fat deposition, feedlot performance, or carcass weight. Restricting vitamin A intake for 131 d at the end of the finishing period appears to be insufficient to affect the site of fat deposition in Holstein steers.

Digestion, rumen fermentation and circulating concentrations of insulin, growth hormone and IGF-1 in steers given maize silages harvested at three stages of maturity

Advancing maturity of forage maize is associated with increases in the proportion of dry matter (DM) and starch and decreases in the proportions of structural carbohydrates in the ensiled crop. Three maize silages (286 (low, L), 329 (medium, M) and 379 (high, H) g DM per kg fresh weight) plus a concentrate formulated to give isonitrogenous intakes were offered to Holstein-Friesian steers fitted with a cannula in the dorsal sac of the rumen and a 'T' piece cannula in the proximal duodenum in an experiment with a cross-over design that allowed four collection periods. Nutrient flow to the duodenum was estimated using chromium-EDTA. Steers consumed approximately 0(.)6 kg DM per day less of diet L compared with the other two diets (P=0(.)026), resulting in less DM being digested (P=0(.)005) but digestibility did not differ between diets. Similar results were obtained for organic matter. There were no differences between diets in the intake or digestibility of neutral-detergent fibre. Intake, duodenal flow and faecal output of starch were greater for steers offered diets M and H compared with those given diet L (P < 0(.)05). In all diets rumen digestion contributed to over 90% of total digestion of starch, although rumen digestibility declined significantly with advancing maize maturity (P=0(.)002). Molar proportions of acetic acid were higher in diet H (P < 0(.)05) whilst proportions of propionic acid and n-butyric acid were higher in diets M and L. There were no significant differences between diets in mean rumen pH or ammonia concentrations. Mean circulating concentrations of insulin were higher (P=0(.)009) in cattle given diets L and M compared with diet H. There were no differences between diets in the mean circulating concentration of growth hormone, or the frequency, amplitude and duration of growth hormone pulses, or the mean circulating concentrations of IGF-1. Changes in forage composition that accompany advancing maize maturity affect overall silage digestion and circulating concentrations of insulin.