Stable isotope and modelling evidence for CO2 as a driver of glacial–interglacial vegetation shifts in southern Africa

Atmospheric CO2 concentration is hypothesized to influence vegetation distribution via tree–grass competition, with higher CO2 concentrations favouring trees. The stable carbon isotope (δ13C) signature of vegetation is influenced by the relative importance of C4 plants (including most tropical grasses) and C3 plants (including nearly all trees), and the degree of stomatal closure – a response to aridity – in C3 plants. Compound-specific δ13C analyses of leaf-wax biomarkers in sediment cores of an offshore South Atlantic transect are used here as a record of vegetation changes in subequatorial Africa. These data suggest a large increase in C3 relative to C4 plant dominance after the Last Glacial Maximum. Using a process-based biogeography model that explicitly simulates 13C discrimination, it is shown that precipitation and temperature changes cannot explain the observed shift in δ13C values. The physiological effect of increasing CO2 concentration is decisive, altering the C3/C4 balance and bringing the simulated and observed δ13C values into line. It is concluded that CO2 concentration itself was a key agent of vegetation change in tropical southern Africa during the last glacial–interglacial transition. Two additional inferences follow. First, long-term variations in terrestrial δ13Cvalues are not simply a proxy for regional rainfall, as has sometimes been assumed. Although precipitation and temperature changes have had major effects on vegetation in many regions of the world during the period between the Last Glacial Maximum and recent times, CO2 effects must also be taken into account, especially when reconstructing changes in climate between glacial and interglacial states. Second, rising CO2 concentration today is likely to be influencing tree–grass competition in a similar way, and thus contributing to the "woody thickening" observed in savannas worldwide. This second inference points to the importance of experiments to determine how vegetation composition in savannas is likely to be influenced by the continuing rise of CO2 concentration.

Data assimilation with correlated observation errors: experiments with a 1-D shallow water model

Remote sensing observations often have correlated errors, but the correlations are typically ignored in data assimilation for numerical weather prediction. The assumption of zero correlations is often used with data thinning methods, resulting in a loss of information. As operational centres move towards higher-resolution forecasting, there is a requirement to retain data providing detail on appropriate scales. Thus an alternative approach to dealing with observation error correlations is needed. In this article, we consider several approaches to approximating observation error correlation matrices: diagonal approximations, eigendecomposition approximations and Markov matrices. These approximations are applied in incremental variational assimilation experiments with a 1-D shallow water model using synthetic observations. Our experiments quantify analysis accuracy in comparison with a reference or ‘truth’ trajectory, as well as with analyses using the ‘true’ observation error covariance matrix. We show that it is often better to include an approximate correlation structure in the observation error covariance matrix than to incorrectly assume error independence. Furthermore, by choosing a suitable matrix approximation, it is feasible and computationally cheap to include error correlation structure in a variational data assimilation algorithm.

Statistical–dynamical downscaling of present day and future precipitation regimes in the Aksu river catchment in Central Asia

A statistical–dynamical downscaling (SDD) approach is applied to determine present day and future high-resolution rainfall distributions in the catchment of the river Aksu at the southern slopes of the Tienshan Mountains, Central Asia. First, a circulation weather type (CWT) classification is employed to define typical lower atmospheric flow regimes from ERA-40 reanalysis data. Selected representatives of each CWT are dynamically downscaled with the regional climate model COSMO-CLM 4.8 at a horizontal grid resolution of 0.0625°, using the ERA-40 reanalysis data as boundary conditions. Finally, the simulated representatives are recombined to obtain a high-resolution rainfall climatology for present day climate. The methodology is also applied to ensemble simulations of three different scenarios of the global climate model ECHAM5/MPI-OM1 to derive projections of rainfall changes until 2100. Comparisons of downscaled seasonal and annual rainfall with observational data suggest that the statistical–dynamical approach is appropriate to capture the observed present-day precipitation climatology over the low lands and the first elevations of the Tienshan Mountains. On the other hand, a strong bias is found at higher altitudes, where precipitation is clearly underestimated by SDD. The application of SDD to the ECHAM5/MPI-OM1 ensemble reveals that precipitation changes by the end of the 21st century depend on the season. While for autumn an increase of seasonal precipitation is found for all simulations, a decrease in precipitation is obtained during winter for most parts of the Aksu catchment. The spread between different ECHAM5/MPI-OM1 ensemble members is strongest in spring, where trends of opposite sign are found. The largest changes in rainfall are simulated for the summer season, which also shows the most pronounced spatial heterogeneity. Most ECHAM5/MPI-OM1 realizations indicate a decrease of annual precipitation over large parts of the Tienshan, and an increase restricted to the southeast of the study area. These results provide a good basis for downscaling present-day and future rainfall distributions for hydrological purposes.

Ensemble projections for wine production in the Douro Valley of Portugal

Wine production is largely governed by atmospheric conditions, such as air temperature and precipitation, together with soil management and viticultural/enological practices. Therefore, anthropogenic climate change is likely to have important impacts on the winemaking sector worldwide. An important winemaking region is the Portuguese Douro Valley, which is known by its world-famous Port Wine. The identification of robust relationships between atmospheric factors and wine parameters is of great relevance for the region. A multivariate linear regression analysis of a long wine production series (1932–2010) reveals that high rainfall and cool temperatures during budburst, shoot and inflorescence development (February-March) and warm temperatures during flowering and berry development (May) are generally favourable to high production. The probabilities of occurrence of three production categories (low, normal and high) are also modelled using multinomial logistic regression. Results show that both statistical models are valuable tools for predicting the production in a given year with a lead time of 3–4 months prior to harvest. These statistical models are applied to an ensemble of 16 regional climate model experiments following the SRES A1B scenario to estimate possible future changes. Wine production is projected to increase by about 10 % by the end of the 21st century, while the occurrence of high production years is expected to increase from 25 % to over 60 %. Nevertheless, further model development will be needed to include other aspects that may shape production in the future. In particular, the rising heat stress and/or changes in ripening conditions could limit the projected production increase in future decades.

Climate change scenarios for precipitation extremes in Portugal

Precipitation indices are commonly used as climate change indicators. Considering four Climate Variability and Predictability-recommended indices, this study assesses possible changes in their spatial patterns over Portugal under future climatic conditions. Precipitation data from the regional climate model Consortium for Small-Scale Modelling–Climate version of the Local Model (CCLM) ensemble simulations with ECHAM5/MPI-OM1 boundary conditions are used for this purpose. For recent–past, medians and probability density functions of the CCLM-based indices are validated against station-based and gridded observational dataset from ENSEMBLES-based (gridded daily precipitation data provided by the European Climate Assessment & Dataset project) indices. It is demonstrated that the model is able to realistically reproduce not only precipitation but also the corresponding extreme indices. Climate change projections for 2071–2100 (A1B and B1 SRES scenarios) reveal significant decreases in total precipitation, particularly in autumn over northwestern and southern Portugal, though changes exhibit distinct local and seasonal patterns and are typically stronger for A1B than for B1. The increase in winter precipitation over northeastern Portugal in A1B is the most important exception to the overall drying trend. Contributions of extreme precipitation events to total precipitation are also expected to increase, mainly in winter and spring over northeastern Portugal. Strong projected increases in the dry spell lengths in autumn and spring are also noteworthy, giving evidence for an extension of the dry season from summer to spring and autumn. Although no coupling analysis is undertaken, these changes are qualitatively related to modifications in the large-scale circulation over the Euro-Atlantic area, more specifically to shifts in the position of the Azores High and associated changes in the large-scale pressure gradient over the area.

Changes in storm track and cyclone activity in three SRES ensemble experiments with the ECHAM5/MPI-OM1 GCM

Synoptic activity over the Northern Hemisphere is evaluated in ensembles of ECHAM5/MPI-OM1 simulations for recent climate conditions (20C) and for three climate scenarios (following SRES A1B, A2, B1). A close agreement is found between the simulations for present day climate and the respective results from reanalysis. Significant changes in the winter mid-tropospheric storm tracks are detected in all three scenario simulations. Ensemble mean climate signals are rather similar, with particularly large activity increases downstream of the Atlantic storm track over Western Europe. The magnitude of this signal is largely dependent on the imposed change in forcing. However, differences between individual ensemble members may be large. With respect to the surface cyclones, the scenario runs produce a reduction in cyclonic track density over the mid-latitudes, even in the areas with increasing mid-tropospheric activity. The largest decrease in track densities occurs at subtropical latitudes, e.g., over the Mediterranean Basin. An increase of cyclone intensities is detected for limited areas (e.g., near Great Britain and Aleutian Isles) for the A1B and A2 experiments. The changes in synoptic activity are associated with alterations of the Northern Hemisphere circulation and background conditions (blocking frequencies, jet stream). The North Atlantic Oscillation index also shows increased values with enhanced forcing. With respect to the effects of changing synoptic activity, the regional change in cyclone intensities is accompanied by alterations of the extreme surface winds, with increasing values over Great Britain, North and Baltic Seas, as well as the areas with vanishing sea ice, and decreases over much of the subtropics.

Contrasting the responses of three different ground-based instruments to energetic electron precipitation

In order to make best use of the opportunities provided by space missions such as the Radiation Belt Storm Probes, we determine the response of complementary subionospheric radiowave propagation measurements (VLF), riometer absorption measurements (CNA), and GPS-produced total electron content (vTEC) to different energetic electron precipitation (EEP). We model the relative sensitivity and responses of these instruments to idealised monoenergetic beams of precipitating electrons, and more realistic EEP spectra chosen to represent radiation belts and substorm precipitation. In the monoenergetic beam case, we find riometers are more sensitive to the same EEP event occurring during the day than during the night, while subionospheric VLF shows the opposite relationship, and the change in vTEC is independent. In general, the subionospheric VLF measurements are much more sensitive than the other two techniques for EEP over 200 keV, responding to flux magnitudes two-three orders of magnitude smaller than detectable by a riometer. Detectable TEC changes only occur for extreme monoenergetic fluxes. For the radiation belt EEP case, clearly detectable subionospheric VLF responses are produced by daytime fluxes that are ~10 times lower than required for riometers, while nighttime fluxes can be 10,000 times lower. Riometers are likely to respond only to radiation belt fluxes during the largest EEP events and vTEC is unlikely to be significantly disturbed by radiation belt EEP. For the substorm EEP case both the riometer absorption and the subionospheric VLF technique respond significantly, as does the change in vTEC, which is likely to be detectable at ~3-4 TECu.

Dependence of winter precipitation over Portugal on NAO and baroclinic wave activity

The relationship between winter (DJF) rainfall over Portugal and the variable large scale circulation is addressed. It is shown that the poles of the sea level pressure (SLP) field variability associated with rainfall variability are shifted about 15° northward with respect to those used in standard definitions of the North Atlantic Oscillation (NAO). It is suggested that the influence of NAO on rainfall dominantly arises from the associated advection of humidity from the Atlantic Ocean. Rainfall is also related to different aspects of baroclinic wave activity, the variability of the latter quantity in turn being largely dependent on the NAO. A negative NAO index (leading to increased westerly surface geostrophic winds into Portugal) is associated with an increased number of deep (ps<980 hPa) surface lows over the central North Atlantic and of intermediate (980precipitation, but rather contribute to advection at their southern flanks. The other aspect of baroclinic wave activity varying with the NAO is the mid-tropospheric storm track (defined by the 500 hPa bandpass-filtered geopotential height variance). A possible local influence of the storm track due to vertical motions ahead of the upper air troughs cannot be unambiguously separated from the effect of advection. A separate influence of local surface cyclones over the Iberian peninsula which may, for instance, arise from the large scale ascent of air, is revealed by the statistics: for a given advection, rainfall amounts for months with local cyclone cores over the considered region tend to exceed those without. Copyright © 1999 Royal Meteorological Society

Sensitivity of Bolivian seasonally-dry tropical forest to precipitation and temperature changes over glacial-interglacial timescales

We used fossil pollen to investigate the response of the eastern Chiquitano seasonally-dry tropical forest (SDTF), lowland Bolivia, to high-amplitude climate change associated with glacial–interglacial cycles. Changes in the structure, composition and diversity of the past vegetation are compared with palaeoclimate data previously reconstructed from the same record, and these results shed light on the biogeographic history of today’s highly disjunct blocks of SDTF across South America. We demonstrate that lower glacial temperatures limited tropical forest in the Chiquitanía region, and suggest that SDTF was absent or restricted at latitudes below 17°S, the proposed location of the majority of the hypothesized ‘Pleistocene dry forest arc’ (PDFA). At 19500 yrs b.p., warming supported the establishment of a floristically-distinct SDTF, which showed little change throughout the glacial–Holocene transition, despite a shift to significantly wetter conditions beginning ca. 12500–12200 yrs b.p. Anadenanthera colubrina, a key SDTF taxon, arrived at 10000 yrs b.p., which coincides with the onset of drought associated with an extended dry season. Lasting until 3000 yrs b.p., Holocene drought caused a floristic shift to more drought-tolerant taxa and a reduction in α-diversity (shown by declining palynological richness), but closed-canopy forest was maintained throughout. In contrast to the PDFA, the modern distribution of SDTF most likely represents the greatest spatial coverage of these forests in southern South America since glacial times. We find that temperature is a key climatic control upon the distribution of lowland South American SDTF over glacial-interglacial timescales, and seasonality of rainfall exerts a strong control on their floristic composition.

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.

Localization of deep water formation: role of atmospheric moisture transport and geometrical constraints on ocean circulation

A series of coupled atmosphere–ocean–ice aquaplanet experiments is described in which topological constraints on ocean circulation are introduced to study the role of ocean circulation on the mean climate of the coupled system. It is imagined that the earth is completely covered by an ocean of uniform depth except for the presence or absence of narrow barriers that extend from the bottom of the ocean to the sea surface. The following four configurations are described: Aqua (no land), Ridge (one barrier extends from pole to pole), Drake (one barrier extends from the North Pole to 35°S), and DDrake (two such barriers are set 90° apart and join at the North Pole, separating the ocean into a large basin and a small basin, connected to the south). On moving from Aqua to Ridge to Drake to DDrake, the energy transports in the equilibrium solutions become increasingly “realistic,” culminating in DDrake, which has an uncanny resemblance to the present climate. Remarkably, the zonal-average climates of Drake and DDrake are strikingly similar, exhibiting almost identical heat and freshwater transports, and meridional overturning circulations. However, Drake and DDrake differ dramatically in their regional climates. The small and large basins of DDrake exhibit distinctive Atlantic-like and Pacific-like characteristics, respectively: the small basin is warmer, saltier, and denser at the surface than the large basin, and is the main site of deep water formation with a deep overturning circulation and strong northward ocean heat transport. A sensitivity experiment with DDrake demonstrates that the salinity contrast between the two basins, and hence the localization of deep convection, results from a deficit of precipitation, rather than an excess of evaporation, over the small basin. It is argued that the width of the small basin relative to the zonal fetch of atmospheric precipitation is the key to understanding this salinity contrast. Finally, it is argued that many gross features of the present climate are consequences of two topological asymmetries that have profound effects on ocean circulation: a meridional asymmetry (circumpolar flow in the Southern Hemisphere; blocked flow in the Northern Hemisphere) and a zonal asymmetry (a small basin and a large basin).

Controlling spurious diapycnal mixing in eddy-resolving height-coordinate ocean models: insights from virtual deliberate tracer release experiments

A perceived limitation of z-coordinate models associated with spurious diapycnal mixing in eddying, frontal flow, can be readily addressed through appropriate attention to the tracer advection schemes employed. It is demonstrated that tracer advection schemes developed by Prather and collaborators for application in the stratosphere, greatly improve the fidelity of eddying flows, reducing levels of spurious diapycnal mixing to below those directly measured in field experiments, ∼1 × 10−5 m2 s−1. This approach yields a model in which geostrophic eddies are quasi-adiabatic in the ocean interior, so that the residual-mean overturning circulation aligns almost perfectly with density contours. A reentrant channel configuration of the MIT General Circulation Model, that approximates the Antarctic Circumpolar Current, is used to examine these issues. Virtual analogs of ocean deliberate tracer release field experiments reinforce our conclusion, producing passive tracer solutions that parallel field experiments remarkably well.

Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.