Measured and modelled source water δ18O based on tree-ring cellulose of larch and pine trees from the permafrost zone
iForest - Biogeosciences and Forestry, Volume 13, Issue 3, Pages 224-229 (2020)
doi: https://doi.org/10.3832/ifor3212-013
Published: Jun 19, 2020 - Copyright © 2020 SISEF
Research Articles
Abstract
To identify source water for trees growing on permafrost in Siberia, we applied mechanistic models that quantify physical and biochemical fractionation processes, leading to oxygen isotope variation (δ18O) in plant organic matter. These models allowed us to investigate the influence of a variety of climatic factors on tree-ring cellulose from two dominant species: Larix cajanderi Mayr. from northeastern Yakutia (69° 22′ N, 148° 25′ E, ~ 250 m a.s.l.) and Pinus sylvestris L. from Central Yakutia (62°14′ N, 129°37′ E, ~ 220 m a.s.l.). The climate of the region is highly continental with short growing seasons, low amount of precipitation and these forest ecosystems are growing on permafrost, which in turn impact the water cycle and climate variation in the δ18O of source water. We compared outputs of the Land surface Processes and eXchanges (LPX-Bern v. 1.3), and Roden-Lin-Ehleringer (RLE) models for the common period from 1945 to 2004. Based on our findings, trees from northeastern and central Yakutia may have access to additional thawed permafrost water during dry summer periods. Owing to differences in the soil structure, active thaw soil depth and root systems of trees at two Siberian sites, Larix cajanderi Mayr. trees can access water not more than from 50 cm depth, in contrast to Pinus sylvestris L. in Central Yakutia which can acquire water from up to 80 cm soil depth. The results enhance our understanding of the growth and survival of the trees in this extreme environment.
Keywords
Conifers, Climate, Drought, Permafrost Thaw Depth, Siberia, δ18O of Source Water
Introduction
Modern climatic changes are leading to significant air temperature increase in the Eurasian subarctic at even faster rates compared to the global average ([33], [14]). Siberian forests respond to both the timing and magnitude of changes in soil moisture as well as to soil temperature, permafrost thickness and distribution, which is in turn affected by snow and vegetation cover, soil texture and geothermal heat flux as well as atmospheric dryness ([5], [4], [8], [24]). Therefore, the fate of trees growing on permafrost undergoing climatic changes is of great interest, owing to the important role of permafrost in these ecosystems: (i) as additional water source ([43], [5], [6], [47]) for trees during droughts; and (ii) for the large amounts of carbon stored in these soils, which essentially contribute to the global carbon and water budget ([5]).
Increasing temperatures at higher latitudes, especially in areas with low precipitation inputs such as parts of Yakutia (200-250 mm year-1), enhance the risks of tree decline and mortality in boreal forests ([39], [4], [6], [44]). A number of studies report a pronounced increase of seasonal permafrost thaw depth in western Siberia ([43], [30], [34], [16]). Sugimoto et al. ([43]) reported, that the active soil layer depth at the central Yakutia can reach up to 1.4 m making thawed permafrost water available for Larix gmelinii Rupr. Moreover, the authors demonstrated that during wet summers in central Yakutia Larix gmelinii Rupr. can use rain water, while thawed permafrost water can be used by trees mainly during drought summers. Fyodorov-Davydov et al. ([16]) investigated the spatial and temporal trends in active soil layer for northern Yakutia and its relation to landscape and climate variables. The soil layer thaws during June-August and freezes again in autumn, defying the length of the short growing season. The few studies in Siberian regions show significant differences between sites in hydrological and geomorphological characteristics, which should be taken into account. To our knowledge, there are no direct seasonal or long-term measurements of the thawing depth of the active soil layer for northeastern Yakutia, which is the driest forested area in the world. Precipitation and thawed permafrost water are particularly crucial for trees growing in this region with severe temperature limitations and low amounts of precipitation ([2], [40], [39]).
Based on isotopic differences between precipitation and permafrost, variations of stable oxygen isotopes (δ18O) may help to reveal water relations in deep-rooted plants ([37]). Oxygen isotopes in plant organic matter are influenced by variations in the isotopic composition of the water source, which is closely related to the precipitation and soil water (after modification by evaporation at the soil surface and due to the shallow active layer depth).
In this study, we aim to determine the source water of Larix cajanderi Mayr. and Pinus sylvestris L. trees from the Siberian subarctic using measured and modeled δ18O variation in annual tree-ring cellulose at two contrasting study sites (northeastern and central Yakutia).
Materials and Methods
Study sites
Both study sites are located in the broad permafrost zone in northeastern Siberia: in northeastern Yakutia (NE - 69° 22′ N, 148° 25′ E, ~ 250 m a.s.l.) and Central Yakutia (CE - 62° 14′ N, 129° 37′ E, ~ 220 m a.s.l.) in Russia.
The mean winter temperatures (December-February) averaged for the period from 1948 to 2004 are -36 °C and -38 °C and in summer (June-July) +12 °C and +19 °C for NE and CE, respectively. However, minimum temperature range during winter can vary between -38 °C and -46 °C, while maximum day time summer temperatures range from +27 °C to +42 °C. Larix cajanderi Mayr. and Pinus sylvestris L. trees are the dominant tree species in NE and CE respectively. The growing period is rather short for trees at both sites and lasts around 60 days ([3]).
Low precipitation inputs are typical for both regions (200-250 mm year-1) and most of the precipitation falls in summer. The main difference between the two sites are soil conditions. Larix cajanderi Mayr. trees in NE are mainly growing on stones, rocks and mixtures of clay ([40]), while Pinus sylvestris L. in CE are growing on mixtures of sand and clay soils ([32]). The root system is also different between these species as they adapt to site-specific edaphic conditions. Larix cajanderi Mayr. roots are shallower than those of Pinus sylvestris L., which tend to be slightly deeper and more extensive in spread.
Local climate data
To assess the drivers of δ18O variation in tree-ring cellulose for measured and modeled datasets, we used daily maximum and minimum temperature, daily precipitation and daily relative humidity data from the Chokurdach weather station (70° 30′ N, 148° 08′ E, ~ 250-300 m a.s.l.) located close to the north-eastern site (NE) for the period from 1948 to 2004. For Central Yakutia (CE) we used data from the weather station in Yakutsk (62° 02′ N, 129° 44′ E, ~ 220 m a.s.l.) for the period from 1966 to 2004.
Annual δ18O in tree-ring cellulose chronologies
As a reference, we used δ18O values of Larix cajanderi Mayr. cellulose chronologies, which were constructed based on individually analysed annual tree rings from four trees from the NE site ([40]) for the period from 1880 to 2004 and from five trees for the period from 1900 to 2013 from the CE sites ([44]). These measured δ18O chronologies were used for comparative analyses with the LPX-Bern modeled data of the source water, δ18O of cellulose, and thaw depths for the common period (from 1901 to 2004) for both sites.
Cellulose samples from both study sites (NE and CE) were measured at the Stable Isotope Research Facility of the Paul Scherrer Institute (Villigen, Switzerland). For measurements of the oxygen isotope ratios from NE, cellulose samples were pyrolysed to CO at 1080 °C ([38]) in the continuous flow mode with an elemental analyzer (Carlo Erba EA-1110®, Milan, Italy) linked to an isotope-ratio mass-spectrometer (Delta S®, Finnigan MAT, Bremen, Germany). The CE samples were pyrolyzed using PYRO-cube® at 1420 °C (Elementar, Hanau, Germany) and analyzed with a Delta Plus XL® Mass spectrometer (Thermo Finnigan, Bremen, Germany). An earlier study showed that both systems yielded similar precisions (± 0.2‰.) and the values from the two instruments were highly correlated (better than the measurement precisions - [49]). The isotopic values were expressed in the δ notation relative to the international standard (eqn. 1):
where Rsample is the molar fraction of 18O/16O ratio of the sample and Rstandard of the Vienna Standard Mean Ocean Water for oxygen.
Estimation of the δ18O of source water using RLE model
Roden-Lin-Ehleringer (RLE) model was applied to determinate the δ18O of source water over time. For snowmelt and thawed permafrost water as well as rainwater we employed ranges between -28 and -25‰ and -6 to -11‰, respectively ([25], Prokushkin, personal communication). An earlier study by Saurer et al. ([37]) showed that soil water is a mixture of precipitation and thawed permafrost water, and that evaporation plays a minor role in this region. Barometric pressure was obtained as an estimate of 95 kPa from elevation tables in Pearcy et al. ([35]). Fan et al. ([12]) reported that values for stomatal conductance in leaves of Salix boganidensis from lowland Indigirka site in the northwestern Yakutia could range between 0.4 and 0.6 mol m-2 s-1. Streit et al. ([42]) reported that the stomatal conductance of Larix decidua needles was higher compared to Pinus mugo and varied from 0.2 to 0.6 mol m-2 s-1 during July at the high-elevated site in Stillberg, Davos (2180 m a.s.l.), Switzerland. To our knowledge, there are no direct measurements for Larix cajanderi Mayr. or Pinus sylvestris for Yakutia. For stomatal conductance we used value of 0.3 mol m-2 s-1, assuming a low stomatal conductance for conifer needles from a northern tree line in the cold Yakutia region.
The δ18O of xylem cellulose (δ18Ocx) was predicted as a function of both the isotopic ratio of the substrate sucrose and of the medium water at the site of cellulose synthesis, root, or xylem as part of the tree-ring ([36]). The δ18O in tree-ring cellulose is calculated as (eqn. 2):
where f0 is the fraction of exchange with leaf water (δ18Owx) and estimated to be 0.42 for oxygen from the best fit relationship derived from controlled experiments ([36]). The subscript cx is xylem cellulose, wx is xylem water and lw is leaf water at the site of sucrose synthesis. The parameter ε0 is a biochemical fractionation associated with the water-carbonyl oxygen interactions (27 ‰ - [41], [27]).
We modeled leaf water as (eqn. 3):
where δ18Olw is the isotope ratio at the site of evaporation, δ18Osource is a mixture of soil and precipitation water used by trees, εe = 1.0096, which is the temperature dependent equilibrium fractionation due to the phase change from liquid water to vapour ([29]), εk = 1.0285, which is the kinetic fractionation ([31]) due to diffusion of vapour into unsaturated air, δ18Ovapour is the isotopic composition of atmospheric water vapour, and ea/ei the ratio of the ambient to the leaf internal vapour pressure (kPa - [9], [10]).
The RLE model can predict source water δ18O values for both the NE and CE study sites by the following equation (eqn. 4):
where ƒ0 = 0.42, which is fraction of exchange with medium water; εoh and εoa = 27, which is heterotrophic and autotrophic fractionation factor during cellulose synthesis; δ18Olwm is the modeled leaf water ([46]).
Simulating δ18O in LPX Bern
The Land surface Processes and eXchanges (LPX-Bern v. 1.3) model ([22]) is a dynamic global vegetation model and features a process-based representation of the coupled terrestrial nitrogen, carbon and water cycle. Here only a natural vegetation, consisting of ten Plant Functional Types (PFTs) competing for resources and confined by bioclimatic limits, is considered. The δ18O isotope enabled version of the LPX-Bern model is introduced and discussed in detail by Keel et al. ([21]), which helps to reduce uncertainties in the interpretation of raw measured δ18O chronologies for investigating different temporal scales of the oxygen isotope patterns in tree rings. The δ18Ocell (stem cellulose isotopic composition) was calculated based on the eqn. 4.
The δ18Olw is calculated according to eqn. 3, but modified by including the Péclet effect ([13]) using the modeled rate of transpiration (eqn. 5):
where P = EL/CD; L is mixing path length (m); E is the rate of transpiration (mol m-2 s-1); C is the concentration of water 5.55 · 104 mol m-3; D is the diffusivity of H218O in water 2.66 · 10-9 m2 s-1; P is the Péclet number ([20]); δ18Olw Péclet is the Péclet corrected leaf water enrichment; δ18Olw from the eqn. 3 is the enrichment of the leaf water at the surface of evaporation.
In this study, the model runs for both study sites with a single grid cell configuration from 1901 to 2004. Monthly temperature, precipitation, cloud cover and number of wet days data are extracted from the 0.5° latitude/longitude grid cells closest to the study cells locations of the global CRU-TS3.22 climate dataset ([19]). Internally, the meteorological data are interpolated to a daily time-step, except for precipitation, where a stochastic weather generator is applied to compute daily precipitation according to Gerten et al. ([17]). The temperature and precipitation data correspond well (r = 0.96; p < 0.05) with the local weather stations (see above) from 1966 to 2004. Similarly, nitrogen deposition is supplied to the model ([26]). Furthermore, the oxygen module requires input values for the δ18O of soil water and water vapor, as well as relative humidity data, which are extracted from global transient isotope-enabled simulations of the coupled atmosphere-land surface model ECHAM5-JSBACH ([18]). Since the simulation only spans the period 1960-2012, earlier years are approximated by the values of the year 1960. Identical atmospheric CO2 concentrations are prescribed to both study sties from a global reconstruction ([11], [15], [28]).
The LPX-Bern δ18O simulations are useful for investigating different temporal scales of the oxygen isotope patterns in tree-ring, which help to reduce uncertainties in the interpretation of raw measured δ18O chronologies.
Climatological and statistical analyses
Pearson correlation coefficients were calculated between local weather station data and the CRU-TS3.22 climate dataset were calculated ([19]). Statistical calculations were performed in licensed version of the software Statistica® v. 13.3 (StatSoft, Tulsa, OK, USA).
Results
Modeled and measured δ18O tree-ring cellulose data
Using the LPX-Bern ([21]) and RLE ([36]) models, we obtained modeled δ18O tree-ring cellulose values for the NE (Fig. 1a) and CE (Fig. 1b) study sites in Yakutia.
Fig. 1 - The δ18O in tree-ring cellulose measured and modeled using LPX-Bern and RLE models. Modeled thaw depth as output parameter of the LPX model is presented for the (a) northeastern (NE) and (b) Central (CE) sites in Yakutia. Maximal seasonal thaw depth (LPX-TD_max) simulated by LPX-Bern model is presented in black lines for both sites (a, b).
Modeled LPX-Bern thaw soil depth varied much stronger for the NE site (in the range from 50-100 cm) compared to the modeled soil depth at the CE site (50-70 cm). Modeled values and data for δ18O in cellulose agree in so far that the tree-ring cellulose is more depleted in the heaver 18O isotope at the NE site than at the CE.
Measured δ18O data from the NE and CE are significantly correlated with RLE (r = 0.92 and r = 0.74, p < 0.05) and LPX-Bern (r = 0.59 and r = 0.46, p < 0.05), respectively for both studied regions. On average the offset between modeled RLE and measured δ18O data for both NE and CE is low (less than 1‰).
The LPX-Bern modeled thaw depth (maximum values) correlates significantly (r = -0.29, p < 0.05) with measured NE δ18O data only. The δ18O measured data from the NE and CE correlated significantly between each other (r = 0.20, p < 0.05) for the common period from 1901 to 2004, although the distance between the two sites is about 860 km. Modeled data for thaw permafrost soil depth at NE and CE also correlate significantly (r = 0.39, p < 0.05).
Modeling of source water
The δ18O leaf water in the RLE is calculated based on the isotopic composition of leaf water at the site of evaporation, which is not the same as bulk leaf water due to the contribution of unevaporative vein water and tissues ([45]) and the “Péclet effect” (see [13] which describe the convective / diffusive mixing between evaporatively enriched and un-enriched source water from the xylem).
A strong increase in estimated isotopic composition of the source water after the 1990s for NE was revealed, while an overall declining trend for CE was found (Fig. 2). Correlation coefficients between source water δ18O and July precipitation were significant at NE (r = -0.27, p < 0.05) and for May precipitation (r = -0.28, p < 0.05) at the CE site. However, the negative correlation with July precipitation was weak.
Fig. 2 - RLE modeled source water δ18O for (a) larch (red line) vs. July precipitation (black bars) and (b) pine trees (blue line) vs. May precipitation (blue bars).
A decrease in July precipitation at the NE site was detected after the year 2000, which is consistent with the modeled δ18O source water values.
Discussion and conclusion
Simulations of thawing depth and δ18O values in tree-ring cellulose, using the output of the RLE model as a proxy for source water over time, indicate that thawed permafrost water is available for trees from NE in addition to precipitation water. Based on this analysis we conclude that thawed permafrost water could be an important water source for trees, especially for periods when the amount of precipitation is lower than the evaporative water loss. Our model simulation is in line with earlier findings by Sugimoto et al. ([43]) for the central part of Yakutia, confirming that trees from the CE could use thawed permafrost water during dry events. Yet the question remains whether Larix cajanderi Mayr. from the NE site can utilize thawed permafrost water. There was also an indication of source water enrichment (18O) increased over time. This could be caused by rising temperatures resulting in higher evapotranspiration due to higher vapor pressure deficit (VPD) during drought conditions ([48]), but may also be the result of a diminishing influence of thawed permafrost water, which has a more negative isotope ratio than rain water ([37]).
Precipitation or snowmelt can be stored in soils or lost as run off, depending on permafrost depth. Source water for Pinus sylvestris L. trees from the CE site did not indicate a strong increase of thawed permafrost use as an additional water source as it was found at the NE. Opposite to the NE, trees from the CE have limited water availability caused by drought conditions in this area.
Our study showed that trees from NE and CE might have access to additional thawed permafrost water during dry summer periods, which is in line with the study by [43]. Due to differences in the soil structure, active thaw soil depth and root systems of trees at two Siberian sites, Larix cajanderi Mayr. trees cannot access water deeper than 50 cm, while Pinus sylvestris L. in Central Yakutia can access water as deep as 80 cm. This is opposite to the finding of Sugimoto et al. ([43]), who showed possible water access up to 1.4 m. These differences can be explained by different soil structures (sand in central Yakutia and rocks in northeastern Yakutia) as well as with species-specific rooting distributions, resulting in different rooting depths. Thus, trees have different access to water sources depending on their root distribution, i.e., surface rooting system (Larix cajanderi Mayr.) vs. roots proliferated to deeper soil layers (Pinus sylvestris L. - [40], [23], [7], [37], [24]).
We could show that the application of different oxygen isotope models ([21], [36]) helps to reveal the impact of important hydro-climatic factors such as the interplay between water sources; e.g., precipitation, soil moisture and frozen soil water and soil thawing depth, influencing tree-growth and physiology in an ecosystem under permafrost conditions. The discrepancies between modeled and measured data are in part due to different parameterizations, but more so due to different model designs for different scales, from whole plant to ecosystems. So far, direct seasonal measurement for the thawed permafrost depth or measurement for stomatal conductance at the remote study site, like NE are rare due to the harsh climatic conditions leading to costly and complex fieldwork at this region. Therefore, the modeling approach can be a powerful tool for the investigation of the influence of a variety of climatic factors on Siberian forest ecosystem water relations that impact on isotope fractionations in tree-ring cellulose from remote Siberian sites.
Furthermore, we must keep in mind that the heterogeneous soil structure will result in a great variation of δ18O values of the xylem water, especially with regard to time gaps between the isotopic signals of precipitation and what we find in xylem water. Depending on the depth of the active soil layer (where trees can access water) and its structure, these time gaps can be considerable ([1]). Such structural variations are difficult to take into account in any modeling approach. Further studies on adjustments and harmonization between LPX-Bern and RLE models and their parameterization based on a wider pool of measured data will provide the basis to enhance the relevance of these models to a powerful analytical tool. Although each model has its strengths and weaknesses, making use of their complementarity and strengths facilitates the enhancement of our understanding of the source water availability and its the complexity for trees in these boreal ecosystems under permafrost dynamics. For instance, the use of a process-based model, such as the LPX-Bern ([21]), allows investigating δ18O signals in tree rings as a proxy for hydrological and climate dynamics on a regional or global scale. Furthermore, factorial simulations (where each parameter is modified specifically) are a potential avenue to further improve our understanding of the linkage between the vegetation cover, climatic drivers and permafrost dynamics.
Acknowledgments
This work was supported by the Russian Science Foundation (RSF) (19-14188000 28), EU-FP7 Marie Curie Fellowship (EU-ISOTREC 235122, 909122). FJ and SL acknowledge support by the Swiss National Science Foundation (# 200020_172476). We thank Martin Werner for providing the soil water and humidity δ18O data from the ECHAM5-JSBACH model.
Author Contributions
OCh did the measurements and data interpretations for the NE Yakutia, and applied the RLE model to both studied sites. GT performed the experiments at the CE Yakutia site. SL modeled δ18O and thaw depth for both NE and CE sites. All authors significantly contributed to the writing of the manuscript.
Conflicts of Interest
The authors declare no conflict of interest.
References
Gscholar
CrossRef | Gscholar
CrossRef | Gscholar
CrossRef | Gscholar
Gscholar
Online | Gscholar
CrossRef | Gscholar
CrossRef | Gscholar
Online | Gscholar
Gscholar
CrossRef | Gscholar
Gscholar
CrossRef | Gscholar
Gscholar
Online | Gscholar
Gscholar
Gscholar
CrossRef | Gscholar
Authors’ Info
Authors’ Affiliation
Siberian Federal University, Krasnoyarsk, Svobodniy pr 79, RU-660041 (Russia)
Galina Timofeeva
Rolf Siegwolf 0000-0002-0249-0651
Matthias Saurer 0000-0002-3954-3534
Swiss Federal Institute for Forest, Snow and Landscape Research WSL, Zürcherstrasse 111 CH-8903 Birmensdorf (Switzerland)
Fortunat Joos 0000-0002-9483-6030
University of Bern, Climate and Environmental Physics, University of Bern, Sidlerstr. 5, CH-3012 Bern (Switzerland)
Fortunat Joos 0000-0002-9483-6030
Oeschger Centre for Climate Change Research, University of Bern, Falkenplatz 16, CH-3012 Bern (Switzerland)
Southern Oregon University, Biology Department, Ashland, OR 97520 (USA)
Corresponding author
Paper Info
Citation
Churakova-Sidorova OV, Lienert S, Timofeeva G, Siegwolf R, Roden J, Joos F, Saurer M (2020). Measured and modelled source water δ18O based on tree-ring cellulose of larch and pine trees from the permafrost zone. iForest 13: 224-229. - doi: 10.3832/ifor3212-013
Academic Editor
Rossella Guerrieri
Paper history
Received: Aug 07, 2019
Accepted: Apr 24, 2020
First online: Jun 19, 2020
Publication Date: Jun 30, 2020
Publication Time: 1.87 months
Copyright Information
© SISEF - The Italian Society of Silviculture and Forest Ecology 2020
Open Access
This article is distributed under the terms of the Creative Commons Attribution-Non Commercial 4.0 International (https://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Web Metrics
Breakdown by View Type
Article Usage
Total Article Views: 34483
(from publication date up to now)
Breakdown by View Type
HTML Page Views: 30357
Abstract Page Views: 1968
PDF Downloads: 1670
Citation/Reference Downloads: 9
XML Downloads: 479
Web Metrics
Days since publication: 1598
Overall contacts: 34483
Avg. contacts per week: 151.05
Article Citations
Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Feb 2023)
Total number of cites (since 2020): 3
Average cites per year: 0.75
Publication Metrics
by Dimensions ©
Articles citing this article
List of the papers citing this article based on CrossRef Cited-by.
Related Contents
iForest Similar Articles
Research Articles
Links between phenology and ecophysiology in a European beech forest
vol. 8, pp. 438-447 (online: 15 December 2014)
Research Articles
Effects of mild drought on the morphology of sun and shade needles in 20-year-old Norway spruce trees
vol. 12, pp. 27-34 (online: 10 January 2019)
Short Communications
Preliminary indications for diverging heat and drought sensitivities in Norway spruce and Scots pine in Central Europe
vol. 13, pp. 89-91 (online: 01 March 2020)
Research Articles
Spatial and temporal variation of drought impact on black locust (Robinia pseudoacacia L.) water status and growth
vol. 8, pp. 743-747 (online: 18 June 2015)
Research Articles
Phytopathogenic fungi in forest nurseries of Middle Siberia
vol. 13, pp. 507-512 (online: 05 November 2020)
Review Papers
Quantifying and modeling water availability in temperate forests: a review of drought and aridity indices
vol. 12, pp. 1-16 (online: 10 January 2019)
Research Articles
Effect of drought stress on some growth, morphological, physiological, and biochemical parameters of two different populations of Quercus brantii
vol. 11, pp. 212-220 (online: 01 March 2018)
Review Papers
Drought-induced mortality of Scots pines at the southern limits of its distribution in Europe: causes and consequences
vol. 3, pp. 95-97 (online: 15 July 2010)
Research Articles
Testing a dual isotope model to track carbon and water gas exchanges in a Mediterranean forest
vol. 2, pp. 59-66 (online: 18 March 2009)
Research Articles
Oak sprouts grow better than seedlings under drought stress
vol. 9, pp. 529-535 (online: 17 March 2016)
iForest Database Search
Search By Author
Search By Keyword
Google Scholar Search
Citing Articles
Search By Author
Search By Keywords
PubMed Search
Search By Author
Search By Keyword