Global maps of soil temperature.

Research in global change ecology relies heavily on global climatic grids derived from estimates of air temperature in open areas at around 2 m above the ground. These climatic grids do not reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions occur and most terrestrial species reside. Here, we provide global maps of soil temperature and bioclimatic variables at a 1-km2 resolution for 0-5 and 5-15 cm soil depth. These maps were created by calculating the difference (i.e. offset) between in situ soil temperature measurements, based on time series from over 1200 1-km2 pixels (summarized from 8519 unique temperature sensors) across all the world's major terrestrial biomes, and coarse-grained air temperature estimates from ERA5-Land (an atmospheric reanalysis by the European Centre for Medium-Range Weather Forecasts). We show that mean annual soil temperature differs markedly from the corresponding gridded air temperature, by up to 10°C (mean = 3.0 ± 2.1°C), with substantial variation across biomes and seasons. Over the year, soils in cold and/or dry biomes are substantially warmer (+3.6 ± 2.3°C) than gridded air temperature, whereas soils in warm and humid environments are on average slightly cooler (-0.7 ± 2.3°C). The observed substantial and biome-specific offsets emphasize that the projected impacts of climate and climate change on near-surface biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments. The global soil-related bioclimatic variables provided here are an important step forward for any application in ecology and related disciplines. Nevertheless, we highlight the need to fill remaining geographic gaps by collecting more in situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications.

Drought-induced increase in water-use efficiency reduces secondary tree growth and tracheid wall thickness in a Mediterranean conifer.

In order to understand the impact of drought and intrinsic water-use efficiency (iWUE) on tree growth, we evaluated the relative importance of direct and indirect effects of water availability on secondary growth and xylem anatomy of Juniperus thurifera, a Mediterranean anisohydric conifer. Dendrochronological techniques, quantitative xylem anatomy, and (13)C/(12)C isotopic ratio were combined to develop standardized chronologies for iWUE, BAI (basal area increment), and anatomical variables on a 40-year-long annually resolved series for 20 trees. We tested the relationship between iWUE and secondary growth at short-term (annual) and long-term (decadal) temporal scales to evaluate whether gains in iWUE may lead to increases in secondary growth. We obtained a positive long-term correlation between iWUE and BAI, simultaneously with a negative short-term correlation between them. Furthermore, BAI and iWUE were correlated with anatomical traits related to carbon sink or storage (tracheid wall thickness and ray parenchyma amount), but no significant correlation with conductive traits (tracheid lumen) was found. Water availability during the growing season significantly modulated tree growth at the xylem level, where growth rates and wood anatomical traits were affected by June precipitation. Our results are consistent with a drought-induced limitation of tree growth response to rising CO2, despite the trend of rising iWUE being maintained. We also remark the usefulness of exploring this relationship at different temporal scales to fully understand the actual links between iWUE and secondary growth dynamics.

Thawing permafrost can mitigate warming-induced drought stress in boreal forest trees.

Perennially frozen soil, also known as permafrost, is important for the functioning and productivity of most of the boreal forest, the world's largest terrestrial biome. A better understanding of complex vegetation-permafrost interrelationships is needed to predict changes in local- to large-scale carbon, nutrient, and water cycle dynamics under future global warming. Here, we analyze tree-ring width and tree-ring stable isotope (C and O) measurements of Gmelin larch (Larix gmelinii (Rupr.) Rupr.) from six permafrost sites in the northern taiga of central Siberia. Our multi-parameter approach shows that changes in tree growth were predominantly controlled by the air and topsoil temperature and moisture content of the active soil and upper permafrost layers. The observed patterns range from strong growth limitations by early summer temperatures at higher elevations to significant growth controls by precipitation at warmer and well-drained lower-elevation sites. Enhanced radial tree growth is mainly found at sites with fast thawing upper mineral soil layers, and the comparison of tree-ring isotopes over five-year periods with different amounts of summer precipitation indicates that trees can prevent drought stress by accessing water from melted snow and seasonally frozen soil. Identifying the active soil and upper permafrost layers as central water resources for boreal tree growth during dry summers demonstrates the complexity of ecosystem responses to climatic changes.

New star on the stage: amount of ray parenchyma in tree rings shows a link to climate.

Tree-ring anatomy reflects the year-by-year impact of environmental factors on tree growth. Up to now, research in this field has mainly focused on the hydraulic architecture, with ray parenchyma neglected despite the growing recognition of its relevance for xylem function. Our aim was to address this gap by exploring the potential of the annual patterns of xylem parenchyma as a climate proxy. We constructed ring-width and ray-parenchyma chronologies from 1965 to 2004 for 20 Juniperus thurifera trees growing in a Mediterranean continental climate. Chronologies were related to climate records by means of correlation, multiple regression and partial correlation analyses. Ray parenchyma responded to climatic conditions at critical stages during the xylogenetic process; namely, at the end of the previous year's xylogenesis (October) and at the onset of earlywood (May) and latewood formation (August). Ray parenchyma-based chronologies have potential to complement ring-width chronologies as a tool for climate reconstructions. Furthermore, medium- and low-frequency signals in the variation of ray parenchyma may improve our understanding of how trees respond to environmental fluctuations and to global change.

VS-Cambium-Developer: A New Predictive Model of Cambium Functioning under the Influence of Environmental Factors.

Climate changes influence seasonal tree-ring formation. The result is a specific cell structure dependent on internal processes and external environmental factors. One way to investigate and analyze these relationships is to apply diverse simulation models of tree-ring growth. Here, we have proposed a new version of the VS-Cambium-Developer model (VS-CD model), which simulates the cambial activity process in conifers. The VS-CD model does not require the manual year-to-year calibration of parameters over a long-term cell production reconstruction or forecast. Instead, it estimates cell production and simulates the dynamics of radial cell development within the growing seasons. Thus, a new software based on R programming technology, able to efficiently adapt to the VS model online platform, has been developed. The model was tested on indirect observations of the cambium functioning in Larix sibirica trees from southern Siberia, namely on the measured annual cell production from 1963 to 2011. The VS-CD model proves to simulate cell production accurately. The results highlighted the efficiency of the presented model and contributed to filling the gap in the simulations of cambial activity, which is critical to predicting the potential impacts of changing environmental conditions on tree growth.

Sex determines xylem anatomy in a dioecious conifer: hydraulic consequences in a drier world.

Increased drought frequency and severity may reshape tree species distribution in arid environments. Dioecious tree species may be more sensitive to climate warming if sex-related vulnerability to drought occurs, since lower performance of one sex may drive differential stress tolerance, sex-related mortality rates and biased sex ratios. We explored the effect of sex and environment on branch hydraulic (hydraulic conductivity and vulnerability to embolism) and trunk anatomical traits in both sexes of the dioecious conifer Juniperus thurifera L. at two sites with contrasting water availability. Additionally, we tested for a trade-off between hydraulic safety (vulnerability to embolism) and efficiency (hydraulic conductivity). Vulnerability to embolism and hydraulic conductivity were unaffected by sex or site at branch level. In contrast, sex played a significant role in xylem anatomy. We found a trade-off between hydraulic safety and efficiency, with larger conductivities related to higher vulnerabilities to embolism. At the anatomical level, females' trunk showed xylem anatomical traits related to greater hydraulic efficiency (higher theoretical hydraulic conductivity) over safety (thinner tracheid walls, lower Mork's Index), whereas males' trunk anatomy followed a more conservative strategy, especially in the drier site. Reconciling the discrepancy between branch hydraulic function and trunk xylem anatomy would require a thorough and integrated understanding of the tree structure-function relationship at the whole-plant level. Nevertheless, lower construction costs and higher efficiency in females' xylem anatomy at trunk level might explain the previously observed higher growth rates in mesic habitats. However, prioritizing efficiency over safety in trunk construction might make females more sensitive to drought, endangering the species' persistence in a drier world.

Assessing Conifer Ray Parenchyma for Ecological Studies: Pitfalls and Guidelines.

Ray parenchyma is an essential tissue for tree functioning and survival. This living tissue plays a major role for storage and transport of water, nutrients, and non-structural carbohydrates (NSC), thus regulating xylem hydraulics and growth. However, despite the importance of rays for tree carbon and water relations, methodological challenges hamper knowledge about ray intra- and inter-tree variability and its ecological meaning. In this study we provide a methodological toolbox for soundly quantifying spatial and temporal variability of different ray features. Anatomical ray features were surveyed in different cutting planes (cross-sectional, tangential, and radial) using quantitative image analysis on stem-wood micro-sections sampled from 41 mature Scots pines (Pinus sylvestris). The percentage of ray surface (PERPAR), a proxy for ray volume, was compared among cutting planes and between early- and latewood to assess measurement-induced variability. Different tangential ray metrics were correlated to assess their similarities. The accuracy of cross-sectional and tangential measurements for PERPAR estimates as a function of number of samples and the measured wood surface was assessed using bootstrapping statistical technique. Tangential sections offered the best 3D insight of ray integration into the xylem and provided the most accurate estimates of PERPAR, with 10 samples of 4 mm(2) showing an estimate within ±6.0% of the true mean PERPAR (relative 95% confidence interval, CI95), and 20 samples of 4 mm(2) showing a CI95 of ±4.3%. Cross-sections were most efficient for establishment of time series, and facilitated comparisons with other widely used xylem anatomical features. Earlywood had significantly lower PERPAR (5.77 vs. 6.18%) and marginally fewer initiating rays than latewood. In comparison to tangential sections, PERPAR was systematically overestimated (6.50 vs. 4.92%) and required approximately twice the sample area for similar accuracy. Radial cuttings provided the least accurate PERPAR estimates. This evaluation of ray parenchyma in conifers and the presented guidelines regarding data accuracy as a function of measured wood surface and number of samples represent an important methodological reference for ray quantification, which will ultimately improve the understanding of the fundamental role of ray parenchyma tissue for the performance and survival of trees growing in stressed environments.