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 surface winds over the Atlantic, Indian and Pacific Ocean sectors of the Southern Ocean in CMIP5 models: historical bias, forcing response, and state dependence

An assessment of the fifth Coupled Models Intercomparison Project (CMIP5) models’ simulation of the near-surface westerly wind jet position and strength over the Atlantic, Indian and Pacific sectors of the Southern Ocean is presented. Compared with reanalysis climatologies there is an equatorward bias of 3.7° (inter-model standard deviation of ± 2.2°) in the ensemble mean position of the zonal mean jet. The ensemble mean strength is biased slightly too weak, with the largest biases over the Pacific sector (-1.6±1.1 m/s, 27 -22%). An analysis of atmosphere-only (AMIP) experiments indicates that 41% of the zonal mean position bias comes from coupling of the ocean/ice models to the atmosphere. The response to future emissions scenarios (RCP4.5 and RCP8.5) is characterized by two phases: (i) the period of most rapid ozone recovery (2000-2049) during which there is insignificant change in summer; and (ii) the period 2050-2098 during which RCP4.5 simulations show no significant change but RCP8.5 simulations show poleward shifts (0.30, 0.19 and 0.28°/decade over the Atlantic, Indian and Pacific sectors respectively), and increases in strength (0.06, 0.08 and 0.15 m/s/decade respectively). The models with larger equatorward position biases generally show larger poleward shifts (i.e. state dependence). This inter-model relationship is strongest over the Pacific sector (r=-0.89) and insignificant over the Atlantic sector (r=-0.50). However, an assessment of jet structure shows that over the Atlantic sector jet shift is significantly correlated with jet width whereas over the Pacific sector the distance between the sub-polar and sub-tropical westerly jets appears to be more important.

A ‘Boscastle-type’ quasi-stationary convective system over the UK Southwest Peninsula

An investigation is presented of a quasi-stationary convective system (QSCS) which occurred over the UK Southwest Peninsula on 21 July 2010. This system was remarkably similar in its location and structure to one which caused devastating flash flooding in the coastal village of Boscastle, Cornwall on 16 August 2004. However, in the 2010 case rainfall accumulations were around four times smaller and no flooding was recorded. The more extreme nature of the Boscastle case is shown to be related to three factors: (1) higher rain rates, associated with a warmer and moister tropospheric column and deeper convective clouds; (2) a more stationary system, due to slower evolution of the large-scale flow; and (3) distribution of the heaviest precipitation over fewer river catchments. Overall, however, the synoptic setting of the two events was broadly similar, suggesting that such conditions favour the development of QSCSs over the Southwest Peninsula. A numerical simulation of the July 2010 event was performed using a 1.5-km grid length configuration of the Met Office Unified Model. This reveals that convection was repeatedly initiated through lifting of low-level air parcels along a quasi-stationary coastal convergence line. Sensitivity tests are used to show that this convergence line was a sea breeze front which temporarily stalled along the coastline due to the retarding influence of an offshore-directed background wind component. Several deficiencies are noted in the 1.5-km model’s representation of the storm system, including delayed convective initiation; however, significant improvements are observed when the grid length is reduced to 500 m. These result in part from an improved representation of the convergence line, which enhances the associated low-level ascent allowing air parcels to more readily reach their level of free convection. The implications of this finding for forecasting convective precipitation are discussed.

Characteristics of high-resolution versions of the Met Office unified model for forecasting convection over the United Kingdom

With many operational centers moving toward order 1-km-gridlength models for routine weather forecasting, this paper presents a systematic investigation of the properties of high-resolution versions of the Met Office Unified Model for short-range forecasting of convective rainfall events. The authors describe a suite of configurations of the Met Office Unified Model running with grid lengths of 12, 4, and 1 km and analyze results from these models for a number of convective cases from the summers of 2003, 2004, and 2005. The analysis includes subjective evaluation of the rainfall fields and comparisons of rainfall amounts, initiation, cell statistics, and a scale-selective verification technique. It is shown that the 4- and 1-km-gridlength models often give more realistic-looking precipitation fields because convection is represented explicitly rather than parameterized. However, the 4-km model representation suffers from large convective cells and delayed initiation because the grid length is too long to correctly reproduce the convection explicitly. These problems are not as evident in the 1-km model, although it does suffer from too numerous small cells in some situations. Both the 4- and 1-km models suffer from poor representation at the start of the forecast in the period when the high-resolution detail is spinning up from the lower-resolution (12 km) starting data used. A scale-selective precipitation verification technique implies that for later times in the forecasts (after the spinup period) the 1-km model performs better than the 12- and 4-km models for lower rainfall thresholds. For higher thresholds the 4-km model scores almost as well as the 1-km model, and both do better than the 12-km model.

Case study modeling of turbulent and mesoscale fluxes over the BOREAS region

Results from aircraft and surface observations provided evidence for the existence of mesoscale circulations over the Boreal Ecosystem-Atmosphere Study (BOREAS) domain. Using an integrated approach that included the use of analytical modeling, numerical modeling, and data analysis, we have found that there are substantial contributions to the total budgets of heat over the BOREAS domain generated by mesoscale circulations. This effect is largest when the synoptic flow is relatively weak, yet it is present under less favorable conditions, as shown by the case study presented here. While further analysis is warranted to document this effect, the existence of mesoscale flow is not surprising, since it is related to the presence of landscape patches, including lakes, which are of a size on the order of the local Rossby radius and which have spatial differences in maximum sensible heat flux of about 300 W m−2. We have also analyzed the vertical temperature profile simulated in our case study as well as high-resolution soundings and we have found vertical profiles of temperature change above the boundary layer height, which we attribute in part to mesoscale contributions. Our conclusion is that in regions with organized landscapes, such as BOREAS, even with relatively strong synoptic winds, dynamical scaling criteria should be used to assess whether mesoscale effects should be parameterized or explicitly resolved in numerical models of the atmosphere.

Will greenhouse gas-induced warming over the next 50 years lead to higher frequency and greater intensity of hurricanes?

The use of a high resolution atmospheric model at T106 resolution, for studying the influence of greenhouse warming on tropical storm climatology, is investigated. The same method for identifying the storms has been used as in a previous study by Bengtsson et al. The sea surface temperature anomalies have been taken from a previous transient climate change experiment, obtained with a low resolution ocean-atmosphere coupled model. The global distribution of the storms, at the time when the CO2 concentration in the atmosphere had doubled, agrees in geographical position and seasonal variability with that of the present climate, but the number of storms is significantly reduced, particularly at the Southern Hemisphere. The main reason to this, appear to be connected to changes in the large scale circulation, such as a weaker Hadley circulation and stronger upper air westerlies. The low level vorticity in the hurricane genesis regions is generally reduced compared to the present climate, while the vertical tropospheric wind shear is somewhat increased. Most tropical storm regions indicate reduced surface windspeeds and a slightly weaker hydrological cycle.

Multimodel assessment of the factors driving stratospheric ozone evolution over the 21st century

The evolution of stratospheric ozone from 1960 to 2100 is examined in simulations from 14 chemistry‐climate models, driven by prescribed levels of halogens and greenhouse gases. There is general agreement among the models that total column ozone reached a minimum around year 2000 at all latitudes, projected to be followed by an increase over the first half of the 21st century. In the second half of the 21st century, ozone is projected to continue increasing, level off, or even decrease depending on the latitude. Separation into partial columns above and below 20 hPa reveals that these latitudinal differences are almost completely caused by differences in the model projections of ozone in the lower stratosphere. At all latitudes, upper stratospheric ozone increases throughout the 21st century and is projected to return to 1960 levels well before the end of the century, although there is a spread among models in the dates that ozone returns to specific historical values. We find decreasing halogens and declining upper atmospheric temperatures, driven by increasing greenhouse gases, contribute almost equally to increases in upper stratospheric ozone. In the tropical lower stratosphere, an increase in upwelling causes a steady decrease in ozone through the 21st century, and total column ozone does not return to 1960 levels in most of the models. In contrast, lower stratospheric and total column ozone in middle and high latitudes increases during the 21st century, returning to 1960 levels well before the end of the century in most models.

Evidence for a climate signal in trends of global crop yield variability over the past 50 years

Low variability of crop production from year to year is desirable for many reasons, including reduced income risk and stability of supplies. Therefore, it is important to understand the nature of yield variability, whether it is changing through time, and how it varies between crops and regions. Previous studies have shown that national crop yield variability has changed in the past, with the direction and magnitude dependent on crop type and location. Whilst such studies acknowledge the importance of climate variability in determining yield variability, it has been assumed that its magnitude and its effect on crop production have not changed through time and, hence, that changes to yield variability have been due to non-climatic factors. We address this assumption by jointly examining yield and climate variability for three major crops (rice, wheat and maize) over the past 50 years. National yield time series and growing season temperature and precipitation were de-trended and related using multiple linear regression. Yield variability changed significantly in half of the crop–country combinations examined. For several crop–country combinations, changes in yield variability were related to changes in climate variability.

Geochemical changes recorded in Lynch's Crater, Northeastern Australia, over the past 50 ka

There are many geochemical reconstructions of environmental change in the mid and high latitudes but relatively few in the tropical latitudes, despite their considerable potential for reconstructing environmental processes that cannot be identified using more traditional proxies. Here we present one reconstruction of environmental change for the tropics. This reconstruction covers the past 50 ka using a suite of geochemical data from the high-resolution sequence of Lynch's Crater in northeast Queensland, Australia, a region highly sensitive to El Niño-Southern Oscillation (ENSO) activity. The 23 major oxides and trace elements measured could be summarised by extracting three axes using principal components analysis (accounting for 72% of the variability). The data indicate that the greatest variability in the geochemical data accounted for erosional activity within the catchment that was associated with past changes in the frequency of ENSO activity (though this was less sensitive during wetter periods, probably as a result of buffering by high vegetation cover). The remaining variability was largely explained by elements that form complexes with organic compounds (e.g., humic acids) and those that are important nutrients for specific vegetation types (and therefore a measure of vegetation distribution). For more detailed reconstructions, further work is required to disentangle the complex controls of elements within sedimentary sequences.

On the development of strong ridge episodes over the eastern North Atlantic

The occurrence of strong and persistent mid-latitude anticyclonic ridges over the Eastern North Atlantic is a major contributor to the occurrence of severe winter droughts over Western Iberia. We analyze the development of strong and persistent ridge episodes within 40–50°N; 40°W–5°E, which are defined as 300 hPa geopotential height anomalies above 50 gpm that persist for at least 10 consecutive days. Results suggest that the generation and maintenance of these episodes, with positive stratospheric geopotential anomalies over the North American continent and the adjacent North Pacific, are associated with an intensified polar jet. Such positive anomalies tend to detach from the main stratospheric anomaly and propagate eastwards and downwards as Rossby tropospheric waves. Furthermore, the Eastern North Atlantic ridge is generated and repeatedly reinforced until the stratospheric anomaly dissipates. Results also show evidence for waves breaking anticyclonically during the episodes, which is dynamically coherent with their persistency and quasi-stationarity.

Statistical uncertainty of changes in winter storms over the North Atlantic and Europe in an ensemble of transient climate simulations

Winter storms are among the most important natural hazards affecting Europe. We quantify changes in storm frequency and intensity over the North Atlantic and Europe under future climate scenarios in terms of return periods (RPs) considering uncertainties due to both sampling and methodology. RPs of North Atlantic storms' minimum central pressure (CP) and maximum vorticity (VOR) remain unchanged by 2100 for both the A1B and A2 scenarios compared to the present climate. Whereas shortened RPs for VOR of all intensities are detected for the area between British Isles/North-Sea/western Europe as early as 2040. However, the changes in storm VOR RP may be unrealistically large: a present day 50 (20) year event becomes approximately a 9 (5.5) year event in both A1B and A2 scenarios by 2100. The detected shortened RPs of storms implies a higher risk of occurrence of damaging wind events over Europe.

Evidence for changes in stratospheric transport and mixing over the past three decades based on multiple data sets and tropical leaky pipe analysis

Variability in the strength of the stratospheric Lagrangian mean meridional or Brewer-Dobson circulation and horizontal mixing into the tropics over the past three decades are examined using observations of stratospheric mean age of air and ozone. We use a simple representation of the stratosphere, the tropical leaky pipe (TLP) model, guided by mean meridional circulation and horizontal mixing changes in several reanalyses data sets and chemistry climate model (CCM) simulations, to help elucidate reasons for the observed changes in stratospheric mean age and ozone. We find that the TLP model is able to accurately simulate multiyear variability in ozone following recent major volcanic eruptions and the early 2000s sea surface temperature changes, as well as the lasting impact on mean age of relatively short-term circulation perturbations. We also find that the best quantitative agreement with the observed mean age and ozone trends over the past three decades is found assuming a small strengthening of the mean circulation in the lower stratosphere, a moderate weakening of the mean circulation in the middle and upper stratosphere, and a moderate increase in the horizontal mixing into the tropics. The mean age trends are strongly sensitive to trends in the horizontal mixing into the tropics, and the uncertainty in the mixing trends causes uncertainty in the mean circulation trends. Comparisons of the mean circulation and mixing changes suggested by the measurements with those from a recent suite of CCM runs reveal significant differences that may have important implications on the accurate simulation of future stratospheric climate.

High-latitude climate sensitivity to ice-sheet forcing over the last 120kyr

Interpretation of ice-core records is currently limited by paucity of modelling at adequate temporal and spatial resolutions. Several key questions relate to mechanisms of polar amplification and inter-hemispheric coupling on glacial/interglacial timescales. Here, we present the first results from a large set of global ocean–atmosphere climate model ‘snap-shot’ simulations covering the last 120 000 years using the Hadley Centre climate model (HadCM3) at up to 1 kyr temporal resolution. Two sets of simulations were performed in order to examine the roles of orbit and greenhouse gases versus ice-sheet forcing of orbital-scale climate change. A series of idealised Heinrich events were also simulated, but no changes to aerosols or vegetation were prescribed. This paper focuses on high latitudes and inter-hemispheric linkages. The simulations reproduce polar temperature trends well compared to ice-core reconstructions, although the magnitude is underestimated. Polar amplification varies with obliquity, but this variability is dampened by including variations in land ice coverage, while the overall amplification factor increases. The relatively constant amplification of Antarctic temperatures (with ice-sheet forcing included) suggests it is possible to use Antarctic temperature reconstructions to estimate global changes (which are roughly half the magnitude). Atlantic Ocean overturning circulation varies considerably only with the introduction of Northern Hemisphere ice sheets, but only weakens in the North Atlantic in the deep glacial, when ocean–sea-ice feedbacks result in the movement of the region of deep convection to lower latitudes and with the introduction of freshwater to the surface North Atlantic in order to simulate Heinrich events.