Calibration and bias correction of climate projections for crop modelling: an idealised case study over Europe

Producing projections of future crop yields requires careful thought about the appropriate use of atmosphere-ocean global climate model (AOGCM) simulations. Here we describe and demonstrate multiple methods for ‘calibrating’ climate projections using an ensemble of AOGCM simulations in a ‘perfect sibling’ framework. Crucially, this type of analysis assesses the ability of each calibration methodology to produce reliable estimates of future climate, which is not possible just using historical observations. This type of approach could be more widely adopted for assessing calibration methodologies for crop modelling. The calibration methods assessed include the commonly used ‘delta’ (change factor) and ‘nudging’ (bias correction) approaches. We focus on daily maximum temperature in summer over Europe for this idealised case study, but the methods can be generalised to other variables and other regions. The calibration methods, which are relatively easy to implement given appropriate observations, produce more robust projections of future daily maximum temperatures and heat stress than using raw model output. The choice over which calibration method to use will likely depend on the situation, but change factor approaches tend to perform best in our examples. Finally, we demonstrate that the uncertainty due to the choice of calibration methodology is a significant contributor to the total uncertainty in future climate projections for impact studies. We conclude that utilising a variety of calibration methods on output from a wide range of AOGCMs is essential to produce climate data that will ensure robust and reliable crop yield projections.

Effluent quality and ammonia emissions from out-wintering pads in England, Wales and Ireland

Out-wintering pads offer a reduced cost system for wintering cattle, minimising damage to pasture, providing animal welfare and production benefits, and generate, potentially, a more manageable effluent and lower ammonia emissions. The objectives of the present study were (i) to contribute to improved understanding of the factors impacting on effluent quality, ammonia emissions and animal welfare via observations on four farm-based out-wintering pads (ComOWPs) in England, Wales and Ireland and more detailed studies undertaken on four experimental OWPs (ExpOWPs) constructed at Rothamsted Research North Wyke, Devon, England and (ii) to corroborate the effluent quality data from both the ComOWPs and the ExpOWPs, with findings in the literature. Woodchip size, feeding management and area allowance were the treatment factors applied on the ExpOWPs. These three factors were randomised across the four ExpOWPs, over four 6–7 week periods. Effluent quality from the ExpOWPs was sampled frequently in a flow proportional way and analysed for total N (TN); total P (TP); total solids (TS); ammonium-N (NH4+-N); nitrate-N (NO3−-N). Beef cattle were periodically weighed for determination of live weight gain (LWG). An approximate nitrogen balance was calculated as a means of understanding its partitioning and fate during and after the ExpOWPs use. Effluent quality from the ComOWPs was sampled frequently, also in a flow-proportional way, and analysed for TN, TP, TS, NH4+-N, NO3−-N, total K and COD. Effluent quality data from the ExpOWPs showed no significant differences (P > 0.05) between treatments, with average concentrations of 1095 mg l−1, and 806 mg l−1, for TN and NH4+-N, respectively. Average effluent concentrations from the ComOWPs were 356 mg l−1 TN and 124 mg l−1 NH4+-N. Ammonia emissions from the ExpOWPs showed no significant differences (P > 0.05) between the treatments, with average mean emission rates of 2.5 g m−2 d−1 NH3-N, respectively. A positive correlation was established between NH3-N emission rate and wind speed. Emission rates from the ComOWPs ranged from 0.7 to 1.6 g m−2 d−1 NH3-N. Average daily LWG on the ExpOWPs was 1.33 kg steer−1 d−1. The effluent from both the ComOWPs and ExpOWPs were more similar with dirty water and of consistently lower strength than beef cattle slurry, as supported by findings in the literature, and therefore, it is suggested to be subject to the regulatory requirements of dirty water rather than slurry.

Reconstruction of gusty wind speed time series from autonomous data logger records

The collection of wind speed time series by means of digital data loggers occurs in many domains, including civil engineering, environmental sciences and wind turbine technology. Since averaging intervals are often significantly larger than typical system time scales, the information lost has to be recovered in order to reconstruct the true dynamics of the system. In the present work we present a simple algorithm capable of generating a real-time wind speed time series from data logger records containing the average, maximum, and minimum values of the wind speed in a fixed interval, as well as the standard deviation. The signal is generated from a generalized random Fourier series. The spectrum can be matched to any desired theoretical or measured frequency distribution. Extreme values are specified through a postprocessing step based on the concept of constrained simulation. Applications of the algorithm to 10-min wind speed records logged at a test site at 60 m height above the ground show that the recorded 10-min values can be reproduced by the simulated time series to a high degree of accuracy.

Enhanced levels of cold shock proteins in Listeria monocytogenes LO28 upon exposure to low temperature and high hydrostatic pressure

Listeria monocytogenes is a psychrotrophic food-borne pathogen that is problematic for the food industry because of its ubiquitous distribution in nature and its ability to grow at low temperatures and in the presence of high salt concentrations. Here we demonstrate that the process of adaptation to low temperature after cold shock includes elevated levels of cold shock proteins (CSPs) and that the levels of CSPs are also elevated after treatment with high hydrostatic pressure (HHP). Two-dimensional gel electrophoresis combined with Western blotting performed with anti-CspB of Bacillus subtilis was used to identify four 7-kDa proteins, designated Csp1, Csp2, Csp3, and Csp4. In addition, Southern blotting revealed four chromosomal DNA fragments that reacted with a csp probe, which also indicated that a CSP family is present in L. monocytogenes LO28. After a cold shock in which the temperature was decreased from 37°C to 10°C the levels of Csp1 and Csp3 increased 10- and 3.5-fold, respectively, but the levels of Csp2 and Csp4 were not elevated. Pressurization of L. monocytogenes LO28 cells resulted in 3.5- and 2-fold increases in the levels of Csp1 and Csp2, respectively. Strikingly, the level of survival after pressurization of cold-shocked cells was 100-fold higher than that of cells growing exponentially at 37°C. These findings imply that cold-shocked cells are protected from HHP treatment, which may affect the efficiency of combined preservation techniques.

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.

Detection and attribution of climate change: a regional perspective

The Intergovernmental Panel on Climate Change fourth assessment report, published in 2007 came to a more confident assessment of the causes of global temperature change than previous reports and concluded that ‘it is likely that there has been significant anthropogenic warming over the past 50 years averaged over each continent except Antarctica.’ Since then, warming over Antarctica has also been attributed to human influence, and further evidence has accumulated attributing a much wider range of climate changes to human activities. Such changes are broadly consistent with theoretical understanding, and climate model simulations, of how the planet is expected to respond. This paper reviews this evidence from a regional perspective to reflect a growing interest in understanding the regional effects of climate change, which can differ markedly across the globe. We set out the methodological basis for detection and attribution and discuss the spatial scales on which it is possible to make robust attribution statements. We review the evidence showing significant human-induced changes in regional temperatures, and for the effects of external forcings on changes in the hydrological cycle, the cryosphere, circulation changes, oceanic changes, and changes in extremes. We then discuss future challenges for the science of attribution. To better assess the pace of change, and to understand more about the regional changes to which societies need to adapt, we will need to refine our understanding of the effects of external forcing and internal variability

Climate and human influences on global biomass burning over the past two millennia

Large, well-documented wildfires have recently generated worldwide attention, and raised concerns about the impacts of humans and climate change on wildfire regimes. However, comparatively little is known about the patterns and driving forces of global fire activity before the twentieth century. Here we compile sedimentary charcoal records spanning six continents to document trends in both natural and anthropogenic biomass burning for the past two millennia. We find that global biomass burning declined from AD 1 to 1750, before rising sharply between 1750 and 1870. Global burning then declined abruptly after 1870. The early decline in biomass burning occurred in concert with a global cooling trend and despite a rise in the human population. We suggest the subsequent rise was linked to increasing human influences, such as population growth and land-use changes. Our compilation suggests that the final decline occurred despite increasing air temperatures and population. We attribute this reduction in the amount of biomass burned over the past 150 years to the global expansion of intensive grazing, agriculture and fire management.

Assessment of temperature, trace species, and ozone in chemistry-climate model simulations of the recent past

Simulations of the stratosphere from thirteen coupled chemistry-climate models (CCMs) are evaluated to provide guidance for the interpretation of ozone predictions made by the same CCMs. The focus of the evaluation is on how well the fields and processes that are important for determining the ozone distribution are represented in the simulations of the recent past. The core period of the evaluation is from 1980 to 1999 but long-term trends are compared for an extended period (1960–2004). Comparisons of polar high-latitude temperatures show that most CCMs have only small biases in the Northern Hemisphere in winter and spring, but still have cold biases in the Southern Hemisphere spring below 10 hPa. Most CCMs display the correct stratospheric response of polar temperatures to wave forcing in the Northern, but not in the Southern Hemisphere. Global long-term stratospheric temperature trends are in reasonable agreement with satellite and radiosonde observations. Comparisons of simulations of methane, mean age of air, and propagation of the annual cycle in water vapor show a wide spread in the results, indicating differences in transport. However, for around half the models there is reasonable agreement with observations. In these models the mean age of air and the water vapor tape recorder signal are generally better than reported in previous model intercomparisons. Comparisons of the water vapor and inorganic chlorine (Cly) fields also show a large intermodel spread. Differences in tropical water vapor mixing ratios in the lower stratosphere are primarily related to biases in the simulated tropical tropopause temperatures and not transport. The spread in Cly, which is largest in the polar lower stratosphere, appears to be primarily related to transport differences. In general the amplitude and phase of the annual cycle in total ozone is well simulated apart from the southern high latitudes. Most CCMs show reasonable agreement with observed total ozone trends and variability on a global scale, but a greater spread in the ozone trends in polar regions in spring, especially in the Arctic. In conclusion, despite the wide range of skills in representing different processes assessed here, there is sufficient agreement between the majority of the CCMs and the observations that some confidence can be placed in their predictions.

Global biomass burning: a synthesis and review of Holocene paleofire records and their controls.

We synthesize existing sedimentary charcoal records to reconstruct Holocene fire history at regional, continental and global scales. The reconstructions are compared with the two potential controls of burning at these broad scales – changes in climate and human activities – to assess their relative importance on trends in biomass burning. Here we consider several hypotheses that have been advanced to explain the Holocene record of fire, including climate, human activities and synergies between the two. Our results suggest that 1) episodes of high fire activity were relatively common in the early Holocene and were consistent with climate changes despite low global temperatures and low levels of biomass burning globally; 2) there is little evidence from the paleofire record to support the Early Anthropocene Hypothesis of human modification of the global carbon cycle; 3) there was a nearly-global increase in fire activity from 3 to 2 ka that is difficult to explain with either climate or humans, but the widespread and synchronous nature of the increase suggests at least a partial climate forcing; and 4) burning during the past century generally decreased but was spatially variable; it declined sharply in many areas, but there were also large increases (e.g., Australia and parts of Europe). Our analysis does not exclude an important role for human activities on global biomass burning during the Holocene, but instead provides evidence for a pervasive influence of climate across multiple spatial and temporal scales.

Arctic ozone depletion in 2002–2003 measured by ASUR and comparison with POAM observations

We present ozone loss estimated from airborne measurements taken during January–February and March in the Arctic winter 2002/2003. The first half of the winter was characterized by unusually cold temperatures and the second half by a major stratospheric sudden warming around 15–18 January 2003. The potential vorticity maps show a vortex split in the lower stratosphere during the major warming (MW) in late January and during the minor warming in mid-February due to wave 1 amplification. However, the warming can be termed as a vortex displacement event as there was no vortex split during the MW period at 10 hPa. Very low temperatures, large areas of polar stratospheric clouds (PSCs), and high chlorine activation triggered significant ozone loss in the early winter, as the vortex moved to the midlatitude regions. The ozone depletion derived from the ASUR measurements sampled inside the vortex, in conjunction with the Mimosa-Chim model tracer, shows a maximum of 1.3 ± 0.2 ppmv at 450–500 K by late March. The partial column loss derived from the ASUR ozone profiles reaches up to 61 ± 4 DU in 400–550 K in the same period. The evolution of ozone and ozone loss assessed from the ASUR measurements is in very good agreement with POAM observations. The reduction in ozone estimated from the POAM measurements shows a similar maximum of 1.3 ± 0.2 ppmv at 400–500 K or 63 ± 4 DU in 400–550 K in late March. Our study reveals that the Arctic winter 2002/2003 was unique as it had three minor warmings and a MW, yet showed large loss in ozone. No such feature was observed in any other Arctic winter in the 1989–2010 period. In addition, an unusually large ozone loss in December, around 0.5 ± 0.2 ppmv at 450–500 K or 12 ± 1 DU in 400–550 K, was estimated for the first time in the Arctic. A careful and detailed diagnosis with all available published results for this winter exhibits an average ozone loss of 1.5 ± 0.3 ppmv at 450–500 K or 65 ± 5 DU in 400–550 K by the end of March, which exactly matches the ozone depletion derived from the ASUR, POAM and model data. The early ozone loss together with considerable loss afterwards put the warm Arctic winter 2002/2003 amongst the moderately cold winters in terms of the significance of the ozone loss.

A combined statistical and dynamical approach for downscaling large-scale footprints of European windstorms

The occurrence of mid-latitude windstorms is related to strong socio-economic effects. For detailed and reliable regional impact studies, large datasets of high-resolution wind fields are required. In this study, a statistical downscaling approach in combination with dynamical downscaling is introduced to derive storm related gust speeds on a high-resolution grid over Europe. Multiple linear regression models are trained using reanalysis data and wind gusts from regional climate model simulations for a sample of 100 top ranking windstorm events. The method is computationally inexpensive and reproduces individual windstorm footprints adequately. Compared to observations, the results for Germany are at least as good as pure dynamical downscaling. This new tool can be easily applied to large ensembles of general circulation model simulations and thus contribute to a better understanding of the regional impact of windstorms based on decadal and climate change projections.

Statistical modelling of grapevine yield in the Port Wine region under present and future climate conditions

The impact of projected climate change on wine production was analysed for the Demarcated Region of Douro, Portugal. A statistical grapevine yield model (GYM) was developed using climate parameters as predictors. Statistically significant correlations were identified between annual yield and monthly mean temperatures and monthly precipitation totals during the growing cycle. These atmospheric factors control grapevine yield in the region, with the GYM explaining 50.4% of the total variance in the yield time series in recent decades. Anomalously high March rainfall (during budburst, shoot and inflorescence development) favours yield, as well as anomalously high temperatures and low precipitation amounts in May and June (May: flowering and June: berry development). The GYM was applied to a regional climate model output, which was shown to realistically reproduce the GYM predictors. Finally, using ensemble simulations under the A1B emission scenario, projections for GYM-derived yield in the Douro Region, and for the whole of the twenty-first century, were analysed. A slight upward trend in yield is projected to occur until about 2050, followed by a steep and continuous increase until the end of the twenty-first century, when yield is projected to be about 800 kg/ha above current values. While this estimate is based on meteorological parameters alone, changes due to elevated CO2 may further enhance this effect. In spite of the associated uncertainties, it can be stated that projected climate change may significantly benefit wine yield in the Douro Valley.

European storminess and associated circulation weather types: future changes deduced from a multi-model ensemble of GCM simulations

A range of possible changes in the frequency and characteristics of European wind storms under future climate conditions was investigated on the basis of a multi-model ensemble of 9 coupled global climate model (GCM) simulations for the 20th and 21st centuries following the IPCC SRES A1B scenario. A multi-model approach allowed an estimation of the (un)certainties of the climate change signals. General changes in large-scale atmospheric flow were analysed, the occurrence of wind storms was quantified, and atmospheric features associated with wind storm events were considered. Identified storm days were investigated according to atmospheric circulation, associated pressure patterns, cyclone tracks and wind speed patterns. Validation against reanalysis data revealed that the GCMs are in general capable of realistically reproducing characteristics of European circulation weather types (CWTs) and wind storms. Results are given with respect to frequency of occurrence, storm-associated flow conditions, cyclone tracks and specific wind speed patterns. Under anthropogenic climate change conditions (SRES A1B scenario), increased frequency of westerly flow during winter is detected over the central European investigation area. In the ensemble mean, the number of detected wind storm days increases between 19 and 33% for 2 different measures of storminess, only 1 GCM revealed less storm days. The increased number of storm days detected in most models is disproportionately high compared to the related CWT changes. The mean intensity of cyclones associated with storm days in the ensemble mean increases by about 10 (±10)% in the Eastern Atlantic, near the British Isles and in the North Sea. Accordingly, wind speeds associated with storm events increase significantly by about 5 (±5)% over large parts of central Europe, mainly on days with westerly flow. The basic conclusions of this work remain valid if different ensemble contructions are considered, leaving out an outlier model or including multiple runs of one particular model.

Quantitative estimates of glacial and Holocene temperature and precipitation change in lowland Amazonian Bolivia

Quantitative estimates of temperature and precipitation change during the late Pleistocene and Holocene have been difficult to obtain for much of the lowland Neotropics. Using two published lacustrine pollen records and a climate-vegetation model based on the modern abundance distributions of 154 Neotropical plant families, we demonstrate how family-level counts of fossil pollen can be used to quantitatively reconstruct tropical paleoclimate and provide needed information on historic patterns of climatic change. With this family-level analysis, we show that one area of the lowland tropics, northeastern Bolivia, experienced cooling (1–3 °C) and drying (400 mm/yr), relative to present, during the late Pleistocene (50,000–12,000 calendar years before present [cal. yr B.P.]). Immediately prior to the Last Glacial Maximum (LGM, ca. 21,000 cal. yr B.P.), we observe a distinct transition from cooler temperatures and variable precipitation to a period of warmer temperatures and relative dryness that extends to the middle Holocene (5000–3000 cal. yr B.P.). This prolonged reduction in precipitation occurs against the backdrop of increasing atmospheric CO2 concentrations, indicating that the presence of mixed savanna and dry-forest communities in northeastern Bolivia durng the LGM was not solely the result of low CO2 levels, as suggested previously, but also lower precipitation. The results of our analysis demonstrate the potential for using the distribution and abundance structure of modern Neotropical plant families to infer paleoclimate from the fossil pollen record.

DK-1D: a drift-kinetic simulation tool for modelling the shear Alfvén wave and its interaction with collisionless plasma

We present a highly accurate tool for the simulation of shear Alfven waves (SAW) in collisionless plasma. SAW are important in space plasma environments because for small perpendicular scale lengths they can support an electric field parallel to the ambient magnetic field. Electrons can be accelerated by the parallel electric field and these waves have been implicated as the source of vibrant auroral displays. However, the parallel electric field carried by SAW is small in comparison to the perpendicular electric field of the wave, making it difficult to measure directly in the laboratory, or by satellites in the near-Earth plasma environment. In this paper, we present a simulation code that provides a means to study in detail the SAW-particle interaction in both space and laboratory plasma. Using idealised, small-amplitude propagating waves with a single perpendicular wavenumber, the simulation code accurately reproduces the damping rates and parallel electric field amplitudes predicted by linear theory for varying temperatures and perpendicular scale lengths. We present a rigorous kinetic derivation of the parallel electric field strength for small-amplitude SAW and show that commonly-used inertial and kinetic approximations are valid except for where the ratio of thermal to Alfv\'{e}n speed is between 0.7 and 1.0. We also present nonlinear simulations of large-amplitude waves and show that in cases of strong damping, the damping rates and parallel electric field strength deviate from linear predictions when wave energies are greater than only a few percent of the plasma kinetic energy, a situation which is often observed in the magnetosphere. The drift-kinetic code provides reliable, testable predictions of the parallel electric field strength which can be investigated directly in the laboratory, and will help to bridge the gap between studies of SAW in man-made and naturally occuring plasma.