Recent trends in Inner Asian forest dynamics to temperature and precipitation indicate high sensitivity to climate change

Semi-arid ecosystems play an important role in regulating global climate with the fate of these ecosystems in the Anthropocene depending upon interactions among temperature, precipitation, and CO2. However, in cool-arid environments, precipitation is not the only limitation to forest productivity. Interactions between changes in precipitation and air temperature may enhance soil moisture stress while simultaneously extending growing season length, with unclear consequences for net carbon uptake. This study evaluates recent trends in productivity and phenology of Inner Asian forests (in Mongolia and Northern China) using satellite remote sensing, dendrochronology, and dynamic global vegetation model (DGVM) simulations to quantify the sensitivity of forest dynamics to decadal climate variability and trends. Trends in photosynthetically active radiation fraction (FPAR) between 1982 and 2010 show a greening of about 7% of the region in spring (March, April, May), and 3% of the area ‘browning’ during summertime (June, July, August). These satellite observations of FPAR are corroborated by trends in NPP simulated by the LPJ DGVM. Spring greening trends in FPAR are mainly explained by long-term trends in precipitation whereas summer browning trends are correlated with decreasing precipitation. Tree ring data from 25 sites confirm annual growth increments are mainly limited by summer precipitation (June, July, August) in Mongolia, and spring precipitation in northern China (March, April, May), with relatively weak prior-year lag effects. An ensemble of climate projections from the IPCC CMIP3 models indicates that warming temperatures (spring, summer) are expected to be associated with higher summer precipitation, which combined with CO2 causes large increases in NPP and possibly even greater forest cover in the Mongolian steppe. In the absence of a strong direct CO2 fertilization effect on plant growth (e.g., due to nutrient limitation), water stress or decreased carbon gain from higher autotrophic respiration results in decreased productivity and loss of forest cover. The fate of these semi-arid ecosystems thus appears to hinge upon the magnitude and subtleties of CO2 fertilization effects, for which experimental observations in arid systems are needed to test and refine vegetation models.

The impact of different glacial boundary conditions on atmospheric dynamics and precipitation in the North Atlantic region

Using a highly resolved atmospheric general circulation model, the impact of different glacial boundary conditions on precipitation and atmospheric dynamics in the North Atlantic region is investigated. Six 30-yr time slice experiments of the Last Glacial Maximum at 21 thousand years before the present (ka BP) and of a less pronounced glacial state – the Middle Weichselian (65 ka BP) – are compared to analyse the sensitivity to changes in the ice sheet distribution, in the radiative forcing and in the prescribed time-varying sea surface temperature and sea ice, which are taken from a lower-resolved, but fully coupled atmosphere-ocean general circulation model. The strongest differences are found for simulations with different heights of the Laurentide ice sheet. A high surface elevation of the Laurentide ice sheet leads to a southward displacement of the jet stream and the storm track in the North Atlantic region. These changes in the atmospheric dynamics generate a band of increased precipitation in the mid-latitudes across the Atlantic to southern Europe in winter, while the precipitation pattern in summer is only marginally affected. The impact of the radiative forcing differences between the two glacial periods and of the prescribed time-varying sea surface temperatures and sea ice are of second order importance compared to the one of the Laurentide ice sheet. They affect the atmospheric dynamics and precipitation in a similar but less pronounced manner compared with the topographic changes.

Late Pleistocene glaciation in the Central Andes: Temperature versus humidity control A case study from the eastern Bolivian Andes (17°S) and regional synthesis

A glacier–climate model was used to calculate climatic conditions in a test site on the east Andean slope around Cochabamba (17°S, Bolivia) for the time of the maximum Late Pleistocene glaciation. Results suggest a massive temperature reduction of about − 6.4 °C (+ 1.4/− 1.3 °C), combined with annual precipitation rates of about 1100 mm (+ 570 mm/− 280 mm). This implies no major change in annual precipitation compared with today. Summer precipitation was the source for the humidity in the past, as is the case today. This climate scenario argues for a maximum advance of the paleo-glaciers in the eastern cordillera during the global Last Glacial Maximum (LGM, 20 ka BP), which is confirmed by exposure age dates. In a synthesized view over the central Andes, the results point to an increased summer precipitation-driven Late Glacial (15–10 ka BP) maximum advance in the western part of the Altiplano (18°S–23°S), a temperature-driven maximum advance during full glacial times (LGM) in the eastern cordillera, and a pre- and post-LGM (32 ka BP/14 ka BP) maximum advance around 30°S related to increased precipitation and reduced temperature on the western slope of the Andes. The results indicate the importance of understanding the seasonality and details of the mass balance–climate interaction in order to disentangle drivers for the observed regionally asynchronous past glaciations in the central Andes.

Triple isotope (δD, δ17O, δ18O) study on precipitation, drip water and speleothem fluid inclusions for a Western Central European cave (NW Switzerland)

Deuterium (δD) and oxygen (δ18O) isotopes are powerful tracers of the hydrological cycle and have been extensively used for paleoclimate reconstructions as they can provide information on past precipitation, temperature and atmospheric circulation. More recently, the use of δ17O excess derived from precise measurement of δ17O and δ18O gives new and additional insights in tracing the hydrological cycle whereas uncertainties surround this proxy. However, 17O excess could provide additional information on the atmospheric conditions at the moisture source as well as about fractionations associated with transport and site processes. In this paper we trace water stable isotopes (δD,δ17O and δ18O) along their path from precipitation to cave drip water and finally to speleothem fluid inclusions for Milandre cave in northwestern Switzerland. A two year-long daily resolved precipitation isotope record close to the cave site is compared to collected cave drip water (3 months average resolution) and fluid inclusions of modern and Holocene stalagmites. Amount weighted mean δD,δ18O and δ17O are -71.0‰, -9.9‰, -5.2‰ for precipitation, -60.3‰, -8.7‰, -4.6‰ for cave drip water and -61.3‰, -8.3‰, -4.7‰ for recent fluid inclusions respectively. Second order parameters have also been derived in precipitation and drip water and present similar values with 18 per meg for 17O excess whereas d-excess is 1.5‰ more negative in drip water. Furthermore, the atmospheric signal is shifted towards enriched values in the drip water and fluid inclusions (Δ of ~ + 10‰ for δD). The isotopic composition of cave drip water exhibits a weak seasonal signal which is shifted by around 8 - 10 months (groundwater residence time) when compared to the precipitation. Moreover, we carried out the first δ17O measurement in speleothem fluid inclusions, as well as the first comparison of the δ17 O behaviour from the meteoric water to the fluid inclusions entrapment in speleothems. This study on precipitation, drip water and fluid inclusions will be used as a speleothem proxy calibration for Milandre cave in order to reconstruct paleotemperatures and moisture source variations for Western Central Europe.

The chironomid-temperature relationship: expression in nature and palaeoenvironmental implications

Fossils of chironomid larvae (non-biting midges) preserved in lake sediments are well-established palaeotemperature indicators which, with the aid of numerical chironomid-based inference models (transfer functions), can provide quantitative estimates of past temperature change. This approach to temperature reconstruction relies on the strong relationship between air and lake surface water temperature and the distribution of individual chironomid taxa (species, species groups, genera) that has been observed in different climate regions (arctic, subarctic, temperate and tropical) in both the Northern and Southern hemisphere. A major complicating factor for the use of chironomids for palaeoclimate reconstruction which increases the uncertainty associated with chironomid-based temperature estimates is that the exact nature of the mechanism responsible for the strong relationship between temperature and chironomid assemblages in lakes remains uncertain. While a number of authors have provided state of the art overviews of fossil chironomid palaeoecology and the use of chironomids for temperature reconstruction, few have focused on examining the ecological basis for this approach. Here, we review the nature of the relationship between chironomids and temperature based on the available ecological evidence. After discussing many of the surveys describing the distribution of chironomid taxa in lake surface sediments in relation to temperature, we also examine evidence from laboratory and field studies exploring the effects of temperature on chironomid physiology, life cycles and behaviour. We show that, even though a direct influence of water temperature on chironomid development, growth and survival is well described, chironomid palaeoclimatology is presently faced with the paradoxical situation that the relationship between chironomid distribution and temperature seems strongest in relatively deep, thermally stratified lakes in temperate and subarctic regions in which the benthic chironomid fauna lives largely decoupled from the direct influence of air and surface water temperature. This finding suggests that indirect effects of temperature on physical and chemical characteristics of lakes play an important role in determining the distribution of lake-living chironomid larvae. However, we also demonstrate that no single indirect mechanism has been identified that can explain the strong relationship between chironomid distribution and temperature in all regions and datasets presently available. This observation contrasts with the previously published hypothesis that climatic effects on lake nutrient status and productivity may be largely responsible for the apparent correlation between chironomid assemblage distribution and temperature. We conclude our review by summarizing the implications of our findings for chironomid-based palaeoclimatology and by pointing towards further avenues of research necessary to improve our mechanistic understanding of the chironomid-temperature relationship.

Forcing of stratospheric chemistry and dynamics during the Dalton Minimum

The response of atmospheric chemistry and dynamics to volcanic eruptions and to a decrease in solar activity during the Dalton Minimum is investigated with the fully coupled atmosphere–ocean chemistry general circulation model SOCOL-MPIOM (modeling tools for studies of SOlar Climate Ozone Links-Max Planck Institute Ocean Model) covering the time period 1780 to 1840 AD. We carried out several sensitivity ensemble experiments to separate the effects of (i) reduced solar ultra-violet (UV) irradiance, (ii) reduced solar visible and near infrared irradiance, (iii) enhanced galactic cosmic ray intensity as well as less intensive solar energetic proton events and auroral electron precipitation, and (iv) volcanic aerosols. The introduced changes of UV irradiance and volcanic aerosols significantly influence stratospheric dynamics in the early 19th century, whereas changes in the visible part of the spectrum and energetic particles have smaller effects. A reduction of UV irradiance by 15%, which represents the presently discussed highest estimate of UV irradiance change caused by solar activity changes, causes global ozone decrease below the stratopause reaching as much as 8% in the midlatitudes at 5 hPa and a significant stratospheric cooling of up to 2 °C in the mid-stratosphere and to 6 °C in the lower mesosphere. Changes in energetic particle precipitation lead only to minor changes in the yearly averaged temperature fields in the stratosphere. Volcanic aerosols heat the tropical lower stratosphere, allowing more water vapour to enter the tropical stratosphere, which, via HOx reactions, decreases upper stratospheric and mesospheric ozone by roughly 4%. Conversely, heterogeneous chemistry on aerosols reduces stratospheric NOx, leading to a 12% ozone increase in the tropics, whereas a decrease in ozone of up to 5% is found over Antarctica in boreal winter. The linear superposition of the different contributions is not equivalent to the response obtained in a simulation when all forcing factors are applied during the Dalton Minimum (DM) – this effect is especially well visible for NOx/NOy. Thus, this study also shows the non-linear behaviour of the coupled chemistry-climate system. Finally, we conclude that especially UV and volcanic eruptions dominate the changes in the ozone, temperature and dynamics while the NOx field is dominated by the energetic particle precipitation. Visible radiation changes have only very minor effects on both stratospheric dynamics and chemistry.

Present state of global wetland extent and wetland methane modelling: methodology of a model inter-comparison project (WETCHIMP)

The Wetland and Wetland CH4 Intercomparison of Models Project (WETCHIMP) was created to evaluate our present ability to simulate large-scale wetland characteristics and corresponding methane (CH4) emissions. A multi-model comparison is essential to evaluate the key uncertainties in the mechanisms and parameters leading to methane emissions. Ten modelling groups joined WETCHIMP to run eight global and two regional models with a common experimental protocol using the same climate and atmospheric carbon dioxide (CO2) forcing datasets. We reported the main conclusions from the intercomparison effort in a companion paper (Melton et al., 2013). Here we provide technical details for the six experiments, which included an equilibrium, a transient, and an optimized run plus three sensitivity experiments (temperature, precipitation, and atmospheric CO2 concentration). The diversity of approaches used by the models is summarized through a series of conceptual figures, and is used to evaluate the wide range of wetland extent and CH4 fluxes predicted by the models in the equilibrium run. We discuss relationships among the various approaches and patterns in consistencies of these model predictions. Within this group of models, there are three broad classes of methods used to estimate wetland extent: prescribed based on wetland distribution maps, prognostic relationships between hydrological states based on satellite observations, and explicit hydrological mass balances. A larger variety of approaches was used to estimate the net CH4 fluxes from wetland systems. Even though modelling of wetland extent and CH4 emissions has progressed significantly over recent decades, large uncertainties still exist when estimating CH4 emissions: there is little consensus on model structure or complexity due to knowledge gaps, different aims of the models, and the range of temporal and spatial resolutions of the models.

NGRIP CH4 concentration from 120 to 10 kyr before present and its relation to a δ15N temperature reconstruction from the same ice core

During the last glacial cycle, Greenland temperature showed many rapid temperature variations, the so-called Dansgaard–Oeschger (DO) events. The past atmospheric methane concentration closely followed these temperature variations, which implies that the warmings recorded in Greenland were probably hemispheric in extent. Here we substantially extend and complete the North Greenland Ice Core Project (NGRIP) methane record from the Preboreal Holocene (PB) back to the end of the last interglacial period with a mean time resolution of 54 yr. We relate the amplitudes of the methane increases associated with DO events to the amplitudes of the local Greenland NGRIP temperature increases derived from stable nitrogen isotope (δ15N) measurements, which have been performed along the same ice core (Kindler et al., 2014). We find the ratio to oscillate between 5 parts per billion (ppb) per °C and 18 ppb °C−1 with the approximate frequency of the precessional cycle. A remarkably high ratio of 25.5 ppb °C−1 is reached during the transition from the Younger Dryas (YD) to the PB. Analysis of the timing of the fast methane and temperature increases reveals significant lags of the methane increases relative to NGRIP temperature for DO events 5, 9, 10, 11, 13, 15, 19, and 20. These events generally have small methane increase rates and we hypothesize that the lag is caused by pronounced northward displacement of the source regions from stadial to interstadial. We further show that the relative interpolar concentration difference (rIPD) of methane is about 4.5% for the stadials between DO events 18 and 20, which is in the same order as in the stadials before and after DO event 2 around the Last Glacial Maximum. The rIPD of methane remains relatively stable throughout the full last glacial, with a tendency for elevated values during interstadial compared to stadial periods.

A climatological analysis of high-precipitation events in Dronning Maud Land, Antarctica, and associated large-scale atmospheric conditions

The link between high precipitation in Dronning Maud Land (DML), Antarctica, and the large-scale atmospheric circulation is investigated using ERA-Interim data for 1979–2009. High-precipitation events are analyzed at Halvfarryggen situated in the coastal region of DML and at Kohnen Station located in its interior. This study further includes a comprehensive comparison of high precipitation in ERA-Interim with precipitation data from the Antarctic Mesoscale Prediction System (AMPS) and snow accumulation measurements from automatic weather stations (AWSs), with the limitations of such a comparison being discussed. The ERA-Interim and AMPS precipitation data agree very well. However, the correspondence between high precipitation in ERA-Interim and high snow accumulation at the AWSs is relatively weak. High-precipitation events at both Halvfarryggen and Kohnen are typically associated with amplified upper level waves. This large-scale atmospheric flow pattern is preceded by the downstream development of a Rossby wave train from the eastern South Pacific several days before the precipitation event. At the surface, a cyclone located over the Weddell Sea is the main synoptic ingredient for high precipitation both at Halvfarryggen and at Kohnen. A blocking anticyclone downstream is not a requirement for high precipitation per se, but a larger share of blocking occurrences during the highest-precipitation days in DML suggests that these blocks strengthen the vertically integrated water vapor transport (IVT) into DML. A strong link between high precipitation and the IVT perpendicular to the local orography suggests that IVT could be used as a “proxy” for high precipitation, in particular over DML's interior.

Contrasting interannual and multidecadal NAO variability

Decadal and longer timescale variability in the winter North Atlantic Oscillation (NAO) has considerable impact on regional climate, yet it remains unclear what fraction of this variability is potentially predictable. This study takes a new approach to this question by demonstrating clear physical differences between NAO variability on interannual-decadal (<30 year) and multidecadal (>30 year) timescales. It is shown that on the shorter timescale the NAO is dominated by variations in the latitude of the North Atlantic jet and storm track, whereas on the longer timescale it represents changes in their strength instead. NAO variability on the two timescales is associated with different dynamical behaviour in terms of eddy-mean flow interaction, Rossby wave breaking and blocking. The two timescales also exhibit different regional impacts on temperature and precipitation and different relationships to sea surface temperatures. These results are derived from linear regression analysis of the Twentieth Century and NCEP-NCAR reanalyses and of a high-resolution HiGEM General Circulation Model control simulation, with additional analysis of a long sea level pressure reconstruction. Evidence is presented for an influence of the ocean circulation on the longer timescale variability of the NAO, which is particularly clear in the model data. As well as providing new evidence of potential predictability, these findings are shown to have implications for the reconstruction and interpretation of long climate records.

Parameterization of temperature sensitivity of spring phenology and its application in explaining diverse phenological responses to temperature change

Existing evidence of plant phenological change to temperature increase demonstrates that the phenological responsiveness is greater at warmer locations and in early-season plant species. Explanations of these findings are scarce and not settled. Some studies suggest considering phenology as one functional trait within a plant's life history strategy. In this study, we adapt an existing phenological model to derive a generalized sensitivity in space (SpaceSens) model for calculating temperature sensitivity of spring plant phenophases across species and locations. The SpaceSens model have three parameters, including the temperature at the onset date of phenophases (Tp), base temperature threshold (Tb) and the length of period (L) used to calculate the mean temperature when performing regression analysis between phenology and temperature. A case study on first leaf date of 20 plant species from eastern China shows that the change of Tp and Tb among different species accounts for interspecific difference in temperature sensitivity. Moreover, lower Tp at lower latitude is the main reason why spring phenological responsiveness is greater there. These results suggest that spring phenophases of more responsive, early-season plants (especially in low latitude) will probably continue to diverge from the other late-season plants with temperatures warming in the future.

Exploring Growth Variability and Crown Vitality of Sessile Oak (Quercus Petraea) in the Czech Republic

Unraveling climatic effects on growth of oak - Europe’s most ecologically and economically important forest species - has been the subject of many recent studies; however, more insight based on field data is necessary to better understand the relationship between climate and tree growth and to adapt forest management strategies to future climate change. In this report, we explore the influence of temperature, precipitation and drought variability on the productivity and vitality of oak stands in the Czech Highlands. We collected 180 cores from mature oaks (Quercus petraea) at four forest stands in the Czech Drahany Highlands. Standard dendromethods were used for sample preparation, ring width measurements, cross-dating, chronology development, and the assessment of growth-climate response patterns. Crown vitality was also evaluated, using the modified ICP Forests methodology. Late spring precipitation totals between May and June as well as the mean July temperature for the year of ring formation were found to be the most important factors for oak growth, whereas crown condition was significantly affected by spring and summer drought. This study is rep-resentative for similar bio-ecological habitats across Central Europe and can serve as a dendroclima-tological blueprint for earlier periods for which detailed meteorological information is missing .

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

A 600 years warm-season temperature record from varved sediments of Lago Plomo, Northern Patagonia, Chile (47°S)

High-resolution records of calibrated proxy data for the past millennium are fundamental to place current changes into the context of pre-industrial natural forced and unforced variability. Although the need for regional spatially-explicit comprehensive reconstructions is widely recognized, the proxy data sources are still scarce, particularly for the Southern Hemisphere and especially for South America. We present a 600-year long warm season temperature record from varved sediments of Lago Plomo, a proglacial lake of the Northern Patagonian Ice field in Southern Chile (46°59′S, 72°52′W, 203 m a.s.l.). The thickness of the bright summer sediment layer relative to the dark winter layer (measured as total brightness; % reflectance 400–730 nm) is calibrated against warm season SONDJF temperature (1900–2009; r = 0.58, p(aut) = 0.056, RE = 0.52; CE = 0.15, RMSEP = 0.28 °C; five-year triangular filtered data). In Lago Plomo, warm summer temperatures lead to enhanced glacier melt and suspended sediment transport, which results in thicker light summer layers and to brighter sediments. Although Patagonia shows pronounced regional differences in decadal temperature trends and variability, the 600 years temperature reconstruction from Lago Plomo compares favourably with other regional/continental temperature records, but also emphasizes significant regional differences for which no data and information existed so far. These regional differences seem to be real as they are also reflected in modern climate data sets (1900–2010). The reconstruction shows pronounced subdecadal – multidecadal variability with cold phases during parts of the Little Ice Age (16th and 18th centuries) and in the beginning of the 20th century. The most prominent warm phase is the 19th century which is as warm as the second half of the 20th century. The exceptional summer warmth AD 1780–1810 is also found in other archives of Northern Patagonia and Central Chile. Our record shows the delayed 20th century warming in the Southern Hemisphere. The comparison between winter precipitation and summer temperature (inter-seasonal coupling) from Lago Plomo reveals alternating phases with parallel and contrasting decadal trends of winter precipitation and summer temperature (positive and negative running correlations Rwinter PP; summer TT). This observation from the sediment proxy data is also confirmed by two sets of reanalysis data for the 20th century. Reanalysis data show that phases with negative correlations between winter precipitation and summer temperature (e.g., dry winters and warm summers) at Lago Plomo are characteristic for periods when circumpolar Westerly flow is displaced southward and enhanced around 60°S.

A Late-glacial and Holocene record of climatic change from a Swiss peat humification profile

Colorimetric measurements of alkaline extracts from two Swiss peat cores have provided a complete 14500-year-long record of peat humification, a proxy of effective precipitation. Peat from the cold Younger Dryas (11050–9550 cal. bc) was well preserved despite low levels of precipitation. A particularly dry period, peaking at c. 7100 cal. bc, is indicated by well-decomposed peat. Peat from c. 6750–4250 cal. bc shows a low degree of decomposition, indicating a wet bog surface despite relatively warm temperatures and therefore indicating high levels of precipitation. A sharp transition to higher levels of decomposition c. 4450–3750 cal. bc indicates a major transition to a drier bog surface. Subsequently, peat humification generally decreases towards the end of the deeper profile (c. cal. ad 1050), indicating a gradual return to wetter conditions. This gradual decrease is punctuated by periods of particularly low humification which appear to be due to shifts to higher levels of effective precipitation from c. 2500 to 1350 cal. bc, c. 1050 to 550 cal. bc, centered around 150 cal. bc, and from c. cal. ad 550 onwards. Anthropogenic influences appear to have affected peat humification at the site at least since the Middle Ages. This study indicates that humification in colder regions/time periods could be more affected by temperature than precipitation and vice versa.