Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Using high resolution tracer data to constrain water storage, flux and age estimates in a spatially distributed rainfall-runoff model

Acknowledgements The authors would like to thank Jonathan Dick, Josie Geris, Jason Lessels, and Claire Tunaley for data collection and Audrey Innes for lab sample preparation. We also thank Christian Birkel for discussions about the model structure and comments on an earlier draft of the paper. Climatic data were provided by Iain Malcolm and Marine Scotland Fisheries at the Freshwater Lab, Pitlochry. Additional precipitation data were provided by the UK Meteorological Office and the British Atmospheric Data Centre (BADC).We thank the European Research Council ERC (project GA 335910 VEWA) for funding the VeWa project.

340,000-year centennial-scale marine record of southern hemisphere climatic oscillation

In order to investigate rapid climatic changes at mid-southern latitudes, we have developed centennial-scale paleoceanographic records from the southwest Pacific that enable detailed comparison with Antarctic ice core records. These records suggest close coupling of mid-southern latitudes with Antarctic climate during deglacial and interglacial periods. Glacial sections display higher variability than is seen in Antarctic ice cores, which implies climatic decoupling between mid- and high southern latitudes due to enhanced circum-Antarctic circulation. Structural and temporal similarity with the Greenland ice core record is evident in glacial sections and suggests a degree of interhemispheric synchroneity not predicted from bipolar ice core correlations.

Age model and pollen analysis from site GIK16867 in the northern Angola Basin

Sporomorphs and dinoflagellate cysts from site GIK16867 in the northern Angola Basin record the vegetation history of the West African forest during the last 700 ka in relation to changes in salinity and productivity of the eastern Gulf of Guinea. During most cool and cold periods, the Afromontane forest, rather than the open grass-rich dry forest, expanded to lower altitudes partly replacing the lowland rain forest of the borderlands east of the Gulf of Guinea. Except in Stage 3, when oceanic productivity was high during a period of decreased atmospheric circulation, high oceanic productivity is correlated to strong winds. The response of marine productivity in the course of a climatic cycle, however, is earlier than that of wind vigour and makes wind-stress-induced oceanic upwelling in the area less likely. Monsoon variation is well illustrated by the pollen record of increased lowland rain forest that is paired to the dinoflagellate cyst record of decreased salinity forced by increased precipitation and run-off.

Age model and chemistry of sediment core STK7 from Steisslinger See

Stable isotope records for carbon and oxygen in bulk carbonates, carbon in bulk organic matter, and for total and chromium-reducible sulfur in a lacustrine sediment core from Lake Steisslingen (Southwest Germany) show several distinct and abrupt shifts during the last 15,000 years. Variations in the isotopic composition of authigenic carbonates indicate two major phases in the lake history. In the pre-Holocene, the hydrological budget of the lake was apparently stable. Variations of delta18O values of authigenic carbonates were, therefore, dominantly controlled by temperature changes. A decrease in the delta18Ocarb values of about 2 per mil at the Allerød/Younger Dryas transition is interpreted as a drop in mean annual air temperatures of approximately 5°C. An abrupt temperature increase of a similar magnitude is inferred at the Younger Dryas/Preboreal boundary. Throughout most of the Holocene, the isotopic composition of authigenic carbonates was influenced by marked changes in the hydrological budget of the lake. A major positive excursion in the delta13Ccarb and delta18Ocarb values at the beginning of the Atlantic and a smaller one in the Preboreal were related to evaporation effects, which indicate that dry climatic conditions must have prevailed at that time. A simultaneous increase in delta13C values of bulk organic matter at the beginning of the Atlantic suggests a high level of productivity in the lake. As a consequence, aqueous sulfate became limited as indicated by variations in the delta34S values of total and chromium-reducible sedimentary sulfur. Therefore, we conclude that the beginning of the Atlantic was characterized not only by dry but also by warm climatic conditions, which triggered a higher productivity in the lake. In the Subatlantic sediments, large variations in carbon, oxygen, and sulfur isotope ratios were observed as a result of human activities, causing considerable perturbations in the biogeochemical element cycling of Lake Steisslingen. Results obtained by the study of the continuous 15 ka record of Lake Steisslingen document clearly that isotopic proxy data from lacustrine sediments can provide useful information on environmental and climatic changes of local, regional, and in the case of the Younger Dryas event, of even hemispherical significance.

Mortality and climatic/environmental correlates of Widdringtonia cedargensis trees in the Cederberg, South Africa

Background: Conifer populations appear disproportionately threatened by global change. Most examples are, however, drawn from the northern hemisphere and long-term rates of population decline are not well documented as historical data are often lacking. We use a large and long-term (1931-2013) repeat photography dataset together with environmental data and fire records to account for the decline of the critically endangered Widdringtonia cedarbergensis. Eighty-seven historical and repeat photo-pairs were analysed to establish 20th century changes in W. cedarbergensis demography. A generalized linear mixed-effects model was fitted to determine the relative importance of environmental factors and fire-return interval on mortality for the species. Results: From an initial total of 1313 live trees in historical photographs, 74% had died and only 44 (3.4%) had recruited in the repeat photographs, leaving 387 live individuals. Juveniles (mature adults) had decreased (increased) from 27% (73%) to 8% (92%) over the intervening period. Our model demonstrates that mortality is related to greater fire frequency, higher temperatures, lower elevations, less rocky habitats and aspect (i.e. east-facing slopes had the least mortality). Conclusions: Our results show that W. cedarbergensis populations have declined significantly over the recorded period, with a pronounced decline in the last 30 years. Individuals that established in open habitats at lower, hotter elevations and experienced a greater fire frequency appear to be more vulnerable to mortality than individuals growing within protected, rocky environments at higher, cooler locations with less frequent fires. Climate models predict increasing temperatures for our study area (and likely increases in wildfires). If these predictions are realised, further declines in the species can be expected. Urgent management interventions, including seedling out-planting in fire-protected high elevation sites, reducing fire frequency in higher elevation populations, and assisted migration, should be considered.

Model simulation output data (CCSM3) related to uplift of Africa and Himalaya/Tibet Plateau and the influence on tropical African vegetation

Hominid evolution in the late Miocene has long been hypothesized to be linked to the retreat of the tropical rainforest in Africa. One cause for the climatic and vegetation change often considered was uplift of Africa, but also uplift of the Himalaya and the Tibetan Plateau was suggested to have impacted rainfall distribution over Africa. Recent proxy data suggest that in East Africa open grassland habitats were available to the common ancestors of hominins and apes long before their divergence and do not find evidence for a closed rainforest in the late Miocene. We used the coupled global general circulation model CCSM3 including an interactively coupled dynamic vegetation module to investigate the impact of topography on African hydro-climate and vegetation. We performed sensitivity experiments altering elevations of the Himalaya and the Tibetan Plateau as well as of East and Southern Africa. The simulations confirm the dominant impact of African topography for climate and vegetation development of the African tropics. Only a weak influence of prescribed Asian uplift on African climate could be detected. The model simulations show that rainforest coverage of Central Africa is strongly determined by the presence of elevated African topography. In East Africa, despite wetter conditions with lowered African topography, the conditions were not favorable enough to maintain a closed rainforest. A discussion of the results with respect to other model studies indicates a minor importance of vegetation-atmosphere or ocean-atmosphere feedbacks and a large dependence of the simulated vegetation response on the land surface/vegetation model.

A multi-model approach to evaluate the role of environmental variability and fishing pressure in sardine fisheries

Understanding the fluctuations in population abundance is a central question in fisheries. Sardine fisheries is of great importance to Portugal and is data-rich and of primary concern to fisheries managers. In Portugal, sub-stocks of Sardina pilchardus (sardine) are found in different regions: the Northwest (IXaCN), Southwest (IXaCS) and the South coast (IXaS-Algarve). Each of these sardine sub-stocks is affected differently by a unique set of climate and ocean conditions, mainly during larval development and recruitment, which will consequently affect sardine fisheries in the short term. Taking this hypothesis into consideration we examined the effects of hydrographic (river discharge), sea surface temperature, wind driven phenomena, upwelling, climatic (North Atlantic Oscillation) and fisheries variables (fishing effort) on S. pilchardus catch rates (landings per unit effort, LPUE, as a proxy for sardine biomass). A 20-year time series (1989-2009) was used, for the different subdivisions of the Portuguese coast (sardine sub-stocks). For the purpose of this analysis a multi-model approach was used, applying different time series models for data fitting (Dynamic Factor Analysis, Generalised Least Squares), forecasting (Autoregressive Integrated Moving Average), as well as Surplus Production stock assessment models. The different models were evaluated, compared and the most important variables explaining changes in LPUE were identified. The type of relationship between catch rates of sardine and environmental variables varied across regional scales due to region-specific recruitment responses. Seasonality plays an important role in sardine variability within the three study regions. In IXaCN autumn (season with minimum spawning activity, larvae and egg concentrations) SST, northerly wind and wind magnitude were negatively related with LPUE. In IXaCS none of the explanatory variables tested was clearly related with LPUE. In IXaS-Algarve (South Portugal) both spring (period when large abundances of larvae are found) northerly wind and wind magnitude were negatively related with LPUE, revealing that environmental effects match with the regional peak in spawning time. Overall, results suggest that management of small, short-lived pelagic species, such as sardine quotas/sustainable yields, should be adapted to a regional scale because of regional environmental variability.

An approximate method for the solution of an influence function foil rolling model

Fleck and Johnson (Int. J. Mech. Sci. 29 (1987) 507) and Fleck et al. (Proc. Inst. Mech. Eng. 206 (1992) 119) have developed foil rolling models which allow for large deformations in the roll profile, including the possibility that the rolls flatten completely. However, these models require computationally expensive iterative solution techniques. A new approach to the approximate solution of the Fleck et al. (1992) Influence Function Model has been developed using both analytic and approximation techniques. The numerical difficulties arising from solving an integral equation in the flattened region have been reduced by applying an Inverse Hilbert Transform to get an analytic expression for the pressure. The method described in this paper is applicable to cases where there is or there is not a flat region.