Groundwater Extraction in Floodplain Forests Reduces Radial Growth and Increases Summer Drought Sensitivity of Pedunculate Oak Trees (Quercus robur L.)

Skiadaresis, Georgios; Schwarz, Julia A.; Bauhus, Jürgen

doi:10.3389/ffgc.2019.00005

ORIGINAL RESEARCH article

Front. For. Glob. Change, 15 March 2019

Sec. Forest Disturbance

Volume 2 - 2019 | https://doi.org/10.3389/ffgc.2019.00005

Groundwater Extraction in Floodplain Forests Reduces Radial Growth and Increases Summer Drought Sensitivity of Pedunculate Oak Trees (Quercus robur L.)

$\r\nGeorgios Skiadaresis*$ Georgios Skiadaresis^*

Julia A. Schwarz

Jürgen Bauhus

Chair of Silviculture, Institute of Forest Sciences, University of Freiburg, Freiburg im Breisgau, Germany

In future, groundwater dependent ecosystems are likely to face increasing stress from direct and indirect effects of climate change. European floodplain forests already experience lower groundwater tables caused by river regulations and water extraction. Yet, the effects of more frequent and prolonged dry and hot periods in combination with groundwater extraction on floodplain forests are not well-understood. Here, we conducted a dendroecological study in floodplain forests dominated by pedunculate oak (Quercus robur), a widespread forest type in Europe. Since drought-related precipitation deficits are increasingly less compensated by available groundwater, we hypothesized that groundwater extraction would increase climate sensitivity of growth. Moreover, we expected that younger trees that established under the condition of lowered groundwater tables would be less sensitive to droughts in contrast to older trees that might not be able to adapt to abrupt changes in groundwater levels. In three regions in the Upper Rhine Valley of Germany, we analyzed retrospectively tree-ring growth in Q. robur stands where (a) water extraction had dramatically lowered groundwater levels and trees presumably lost direct access to it, or (b) trees have not been affected by water extraction. In addition, (c) nearby sessile oak (Q. petraea) stands at upland sites were analyzed to isolate the climate signal for each study region from a groundwater signal. In total, we analyzed ring-width chronologies from 18 forest stands and investigated their relationship with climate. Radial growth of oak trees was strongly correlated with summer drought (SPEI, Standardized Precipitation Evapotranspiration Index). Their sensitivity to summer drought increased significantly following the onset of groundwater extraction (88–49 years ago). However, sensitivity of oaks to summer drought at sites with groundwater contact remained relatively stable over time. Oaks at sites without water extraction took advantage of favorable growing conditions. Our results show that access to groundwater can buffer negative effects of summer drought on tree growth. Therefore, in the face of the ongoing climate change and more frequent occurrence of extreme drought events, we recommend that groundwater extraction from floodplain forests should be reduced rather than increased to maintain tree vitality in the long term.

Introduction

Groundwater is one of the most crucial and reliable freshwater resources of the planet. It provides water not only for human consumption, irrigation and industrial uses, but also for numerous ecosystems that depend either partly or exclusively on groundwater to cover their water needs (Bovolo et al., 2009; Orellana et al., 2012; Kløve et al., 2014). However, the quantity and quality of groundwater resources are globally at risk and are expected to be further affected by climate change in coming decades. The observed and projected changes in climate patterns will affect groundwater resources directly through temperature increases and alteration of precipitation regimes (i.e., amounts, forms and seasonal patterns) (Earman and Dettinger, 2011; Green et al., 2011). These direct effects of climate change on groundwater resources are highly uncertain (Taylor et al., 2013; Kløve et al., 2014). However, it is more certain that increasingly frequent and prolonged droughts will intensify groundwater appropriation for irrigation and industrial and domestic consumption (Bovolo et al., 2009; Earman and Dettinger, 2011; Green et al., 2011; Taylor et al., 2013; Kløve et al., 2014). This development will directly affect groundwater dependent ecosystems and amplify the negative effects of climate change on them (Kløve et al., 2014).

Floodplain forests are highly productive and biodiverse forest ecosystems that provide a variety of important ecological services including buffering of stream flow and improvement of water quality (Kozlowski, 2002; Schnitzler et al., 2005; Leuschner and Ellenberg, 2017). Trees in floodplain forests have access to both surface and groundwater to cover their water needs (Cunningham et al., 2011). This additional water resource can possibly explain the high productivity that has been commonly observed in floodplain forests (Cunningham et al., 2011). Alteration of groundwater regimes through water extraction and river regulations have been linked to reduced productivity and vitality and even to extensive dieback of floodplain forests (LeMaitre et al., 1999; Horner et al., 2009; Cunningham et al., 2011; Kath et al., 2014).

In Central Europe, the higher parts of hardwood floodplain forests are dominated by pedunculate oak (Quercus robur L.) (Leuschner and Ellenberg, 2017). Owing to their very fertile soils and the less frequent occurrence of floods, the majority of these forests have been converted to agricultural land in the last centuries (Hughes, 2003; Leuschner and Ellenberg, 2017). Remaining fragments of oak dominated floodplain forests in Europe are of crucial importance because of the multiple economic and ecological services they provide (Eaton et al., 2016; Leuschner and Ellenberg, 2017). However, extensive oak decline in recent decades across Europe has raised concerns about the future of these forests (Thomas et al., 2002; Thomas, 2008). A number of biotic and abiotic stressors are considered responsible for oak decline in Europe (Führer, 1998; Thomas, 2008). These include severe and repeated defoliation by insects, pathogen attacks, weather extremes such as droughts, severe winter and spring frosts, and decreased groundwater levels. Yet, extensive decline and mortality of oaks occurred in most cases only, when at least two stress factors interacted. However, as most of these factors are related to or controlled by climate, individual stress factors as well as their interactions should be examined (Führer, 1998; Thomas, 2008).

The study of tree-rings and hence radial growth of trees provides useful insights into changes of past growing conditions of forests and also facilitates investigations of short-term and long-term responses to abrupt and gradual changes in their environment. The influence of climatic variability is commonly examined as the relationship (i.e., correlation) between variation in annual stem-growth and climatic variables (hereafter referred to as “climate sensitivity”). The two most common oak species of Central Europe, pedunculate oak (Quercus robur L.) and sessile oak [Quercus petraea (Matt.) Liebl.] have been used extensively in dendrochronological studies (for a detailed review on the use of oaks in dendrochronology see the work by Haneca et al., 2009). Sensitivity of growth to drought has increased for both oak species in Central Europe over the last decades, but pedunculate oak appeared to be more sensitive to climatic variability than sessile oak (Friedrichs et al., 2008).

A complex interplay between groundwater availability and climatic variability controls tree growth in groundwater-fed forest ecosystems. However, the role of groundwater availability on tree responses to climatic fluctuations is poorly understood. To avoid or reduce negative impacts of increasing groundwater extraction in future, we need to understand how changes in groundwater availability affect the maintenance of forest structure and functioning (Zolfaghar et al., 2014). A number of studies have highlighted the importance of changes in groundwater regimes as a driver of oak decline and mortality (Thomas, 2008; Levanic et al., 2011; Stojanović et al., 2015). Other work has highlighted the role of groundwater availability for reducing the sensitivity of growth to climatic variability and especially to elevated summer temperatures for Abies balsamea and Pinus strobus growing in groundwater-fed wetlands in the North-eastern USA (Raney et al., 2016). However, most of these studies have focused on trees that were mature when groundwater levels dropped. To understand and predict responses of future oak generations to climate change and the associated decline in groundwater levels, a better understanding of tree-responses at different developmental stages is needed. It has been suggested that young tree-seedlings can adapt to a groundwater drawdown by extending their root systems to access water from deeper soil layers (Imada et al., 2008; Wang et al., 2015). However, the ability of trees to respond to changes in groundwater levels either by extending their roots to keep access to groundwater or by producing more roots in the upper soil layers is species-specific and might be related to the rate and extent of groundwater decline (Naumburg et al., 2005) and also to tree age. Previous studies on pedunculate oak indicated that old (more than 100 years old) in contrast to younger oaks (<100 years old) were not able to change their rooting pattern in response to lowered groundwater levels (Becker and Lévy, 1986; Thomas and Hartmann, 1998).

Although relationships between annual growth of pedunculate oak and climatic variables have been studied extensively, only few studies took into account additional changes in groundwater regimes. To our knowledge, no previous studies have investigated the effects of groundwater extraction on tree-growth sensitivity to climatic variability and especially to summer drought. Our overall aim for this study was to understand how lowered groundwater levels affect the sensitivity of growth to climate, in particular to summer droughts. We hypothesized that:

• Groundwater extraction increases climate sensitivity of tree growth because drought-related precipitation deficits are compensated by available groundwater to a lesser extent;

• The influence of lowered groundwater tables on climate-growth relationships is less pronounced in younger than in older trees, because young trees were established under the condition of water extraction and are therefore more likely to be adapted to the lowered water table.

Methods

Study Sites and Sampling Design

Tree-ring growth in oak dominated forest stands was analyzed in three regions in south-western Germany (Figure 1). The investigated regions are located in the Upper Rhine Valley and situated in former floodplain areas. River regulations and human activities along the river Rhine, such as expansion of settlements, industrialization and intensification of agriculture, have significantly altered its hydrology, morphology, and water quality during the past centuries and caused marked changes in the surrounding ecosystems (Dister et al., 1990; van Dijk et al., 1995). One study region is located close to the city of Freiburg, where groundwater levels have been lowered since 1930, initially due to river regulations and subsequently because of extraction for industrial uses. A second region is located 15 km north of Freiburg near the city of Emmendingen, where groundwater has been extracted for drinking water supply since 1961. The third region near Lampertheim is located approximately 200 km further north along the Upper Rhine Valley. There, massive extraction of groundwater from multiple wells within a forest area of approximately 10 000 hectares has taken place since 1969. This has resulted in dramatic decreases in groundwater levels and decline in forest vitality (Ahner and Schmidt, 2011).

FIGURE 1

Figure 1. Location of the study regions (A), the three study regions near Freiburg, Emmendingen, and Lampertheim (B), and locations of the investigated forest stands in each region (C, D and E for Lampertheim, Emmendingen and Freiburg, respectively). Stand IDs in (C–E) are acronyms denoting the study region and site type (first and second letter, respectively) of each investigated stand. F, E, L for Freiburg, Emmendingen, and Lampertheim, respectively. W indicates sites without groundwater extraction, C upland sites and E, S, and N are extraction sites (E for the extraction sites in Emmendingen and Freiburg, S and N for the Southern and Northern extraction sites, respectively, in Lampertheim).Numbers in Stand IDs indicate the average tree age of each stand. The maps were created using the OpenLayers Plugin (Sourcepole, 2017) in QGIS (QGIS Development Team, 2017).

In each of the three study regions, we selected oak stands at sites with contrasting groundwater regimes (Figure 2). First, we identified forest stands located adjacent to groundwater extraction wells (extraction sites). Second, reference sites were identified, where no extraction has taken place. In all regions, these sites were located in close proximity to the extraction sites but closer to rivers and streams. Trees at these wetter sites (no extraction), have still access to groundwater at least seasonally. In addition, we included a third type of stands on drier uplands sites with free draining soils (upland sites). In all regions these upland stands were dominated by sessile oak, while the other two site types where dominated by pedunculate oak.

FIGURE 2

Figure 2. Graphical representation of site-types in each study region. Oak stands from three different site-types with different groundwater conditions were sampled: Center, Quercus robur stands at sites adjacent to groundwater extraction wells, where trees have lost contact to groundwater (extraction sites), right, Q. robur stands at sites where trees have access to groundwater (no extraction sites) and left, upland Q. petraea stands, where trees never had contact to groundwater (upland sites).

To analyse trees of different age at the onset of groundwater extraction, we selected young and old forest stands for each site type and region. Average stand age (in 2016) ranged from 62 to 181 years for the extraction sites, from 40 to 160 years for sites without extraction and from 65 to 171 years for the upland sites. At Lampertheim, we selected two extraction sites because of the large extent of the affected forest area. In the study regions near Freiburg and Emmendingen, we selected one stand without extraction and one upland stand, respectively, due to the lack of oak stands from different age-classes. Our study design resulted in a total of 18 sampled stands.

Sampling and Dendrochronological Methods

We selected 20 dominant or co-dominant trees per stand except in the LC65 stand in Lampertheim, where only nine trees could be sampled (Table 1). A total of 349 trees were sampled for the purposes of this study. To avoid spatial correlation within a stand, we chose trees that were no direct competitors. The minimum distance between the selected trees within each stand varied among stands depending on their age and size. On average, minimum distance between older trees was 15 and 10 m between individuals in the younger stands. At these distances we cannot rule out below-ground connections between sampled trees through root grafts.

TABLE 1

Table 1. Stand characteristics of the 18 stands from 3 study regions sampled in this study.

First, we measured stem diameter at breast height (DBH at 1.3 m) and the height of each target tree. Additionally, we assessed tree vitality based on the crown conditions of target trees as proposed by Roloff (2001). Vital trees with full crowns and vital fine twigs were assigned to crown class 0; trees assigned to class (1). showed first signs of crown degradation and twig abscission; trees of class (2). were noticeably weakened already; and seriously weakened and dying trees with dead main branches and disintegration of the entire crown were assigned to crown class (3). Selected trees covered all vigor classes present in the studied stands.

We extracted two cores per tree from opposite directions at 1.3 m stem height, using an increment-borer with 0.5 mm diameter (Haglöf, Sweden). At upland sites, where trees were growing on a slope, the two cores were taken perpendicular to the slope direction to avoid sampling of reaction wood. Subsequently, cores were air-dried and smoothed using progressively finer sandpaper to increase visibility of ring boundaries. Increment cores were scanned at 2,400 dpi using an Epson Expression 1649XL flatbed color image scanner. We used the WinDendro software (Regent instruments, Quebec) to measure tree-ring widths from the acquired pictures. The resulting tree-ring width series were visually and statistically crossdated using the functions corr.rwl.seg() and ccf.series.rwl() (Bunn, 2010) from the Dendrochronology Program Library (dplR) (Bunn, 2008) in R (R Core Team, 2017). Prior to quantifying climate-growth relationships, age related trends were removed from the raw tree-ring series by using a 32-year spline with 50% frequency cut-off. This short spline removes the biological trends present in tree-ring series while simultaneously preserving annual to decadal variability in growth (Cook and Peters, 1981; Speer, 2001). Subsequently, the biweight robust mean was calculated to develop a residual chronology per stand. Stand-wise chronologies were truncated to have sample depth higher than five trees for each year to ensure adequate replicates for their whole length. Quality control, detrending and chronology development were performed using the dplR package in R.

Several descriptive statistics were calculated to assess the quality of raw (ar1) and detrended (Mglk, EPS, Rbar, ar1) chronologies for each stand. First order autocorrelation (ar1) was used to assess how current year's growth was affected by growth in the previous year (Fritts, 1976; Speer, 2010). Mean gleichläufigkeit (synchronicity) (Mglk) was used to assess similarity of detrended tree-ring series among trees from each stand (Schweingruber, 1988). The glk() function from the dplR package was used to compute glk (gleichläufigkeit) which performs pairwise comparison of all possible combinations between series within each stand. The expressed population signal (EPS) is an indicator of how well a chronology represents a theoretical infinite population (Wigley et al., 1984; Buras, 2017). It is used as an indicator of the common variability in detrended chronologies. Low values of EPS (commonly <0.85) indicate that the chronologies are dominated by individual tree signals rather than a consistent stand-level signal (Speer, 2010). Rbar is the mean correlation between tree-ring series within a stand. The Rbar was used as a measure of the common signal strength of the detrended chronologies.

Climate Data

Monthly temperature and precipitation data for each study region were obtained from the Climate Data Center of the German Weather Service (Deutscher Wetterdienst) using the meteorological stations closest to the study sites (Supplementary S.1). The distances between the selected meteorological stations and study sites were <20 km. For sites in the regions of Freiburg and Emmendingen, data from one station were used because the two sites are located within 10 km from the selected station. In Lampertheim, two weather stations were selected, one station for the extraction site in the northern part of the region and the no-extraction site and another one for the extraction site in the southern part and the upland site. For the sites in the regions of Freiburg and Emmendingen, climate data were obtained for the period between 1921 and 2016. For the study region near Lampertheim, climate data were obtained for the two sites in the northern part of the region for the period between 1896 and 2016 while for the sites in the southern part, data for the period between 1890 and 2016 were available. Due to the similarity (correlation coefficient for monthly precipitation and temperature series from the two stations were 0.9 and 1.0, respectively) and close proximity (25 km) of the two stations in the region of Lampertheim, climate data from the two stations were averaged over the common period (1896–2016) for all subsequent trend analysis. Temporal trends in seasonal precipitation and temperature were quantified and tested for significance by performing Mann-Kendall tests (Mann, 1945) using the Kendall package in R (McLeod, 2011). Rates of change per year were estimated by calculating non-parametric Sen's slopes (Sen, 1968) using the trend package in R (Pohlert, 2018).

Based on the obtained precipitation and temperature data, we calculated the Standardized Precipitation Evapotranspiration Index (SPEI) (Vicente-Serrano et al., 2010) using the SPEI package in R (Vicente-Serrano et al., 2010). Different periods (6, 12, 24, and 36 months) of the SPEI index were calculated for the months of the vegetation period (from March to September). Additionally, seasonal means for different scales of monthly SPEI values were calculated for the whole vegetation period and for individual seasons: spring (March, April, and May), summer (June, July, and August), autumn (September, October, and November), and winter (December of previous year, January, and February). We chose the multi-scalar SPEI index because previous studies have highlighted its superior performance in assessing drought impacts on ecological and hydrological systems and especially in predicting changes in forest growth (Vicente-Serrano et al., 2012, 2013; Bhuyan et al., 2017).

Statistical Analysis

Overall Climate-Growth Relationships

Overall climate-growth relationships were assessed by calculating bootstrapped Pearson's correlation coefficients from the detrended tree-ring chronologies and monthly and seasonal climatic data, using the package treeclim in R (Zang and Biondi, 2015). We first calculated correlation functions for the separate months of the vegetation period, from March to September, and then seasonal means for spring, summer, autumn, winter as well as for the whole vegetation period. Monthly and seasonal means were calculated for temperature, precipitation and the different periods of the standardized precipitation evapotranspiration index (SPEI). Coefficients of correlation between the stand-specific chronologies and climatic variables were averaged for each month and season, first including all stands, and subsequently separately for the different sites. To obtain average correlation coefficients, we transformed the correlation coefficients for the separate chronologies to Fisher's z scores prior to averaging and then back-transformed them to Pearson's correlation coefficients. This was done because the back-transformed average z score was found to be less biased than simply averaging correlation coefficients (Silver and Dunlap, 1987).

Temporal Development of Climate-Growth Relationships

To investigate the temporal variation of climate-growth relationships, we calculated bootstrapped moving correlation functions. These were calculated for consecutive 30-year windows, which were moved progressively by 1 year. Initially, moving correlation functions were calculated separately for all stand chronologies and monthly and seasonal climatic data as described above for the static correlation functions. To identify and better demonstrate changes related to groundwater extraction, correlation coefficients were presented relative to the date when the extraction of groundwater began (date when the extraction began = 0) for all sites within each of the three regions. Each correlation coefficient for a year x is the correlation coefficient between a tree-ring chronology and the specific climatic variable for the 30-year time window starting in the year x (time period starting in year x until year x+30).

Although we computed moving correlation functions for all climatic variables mentioned above (temperature, precipitation and different scales of SPEI), results presented and discussion are focused on the shorter time-scale of SPEI (6 months), since it was found to correlate best with growth of oaks at all the investigated sites. The calculation of both overall (static) and moving correlations for each of the 18 stands was performed for the maximum period during which climate and growth data overlapped (41–120 years).

Radial Stem-Growth Responses to Dry and Wet Summers

To properly interpret possible differences in sensitivity of annual stem growth to summer drought among trees from sites with differing groundwater regimes, we quantified growth responses for both dry and wet summer conditions. This was done by calculating absolute growth changes (radial growth in year t minus radial growth in year t-1) for individual trees excluding the first 20 years of juvenile growth and the time period before the onset of groundwater extraction. Years were characterized as dry or wet based on the mean SPEI 6 of the three summer months (June to August). Years with SPEI 6 values below −1 were characterized as dry, while years with values >1 were characterized as wet. We tested for differences in growth responses to dry and wet summer conditions among trees at sites with differing groundwater regimes through non-parametric Kruskal-Wallis tests. Trees from sites with similar groundwater conditions were pooled across the three study regions. Afterwards, we used Wilcoxon sign tests to perform pairwise comparisons among the three different sites (extraction, no extraction, and upland).

Results

Tree Ring Data, Chronology Characteristics, and Statistics

The descriptive statistics calculated for the different stands indicate an overall good quality of the developed chronologies (Table 2 and Supplementary S.2, S.3). Mean annual radial growth of the sampled trees ranged from 1.26 to 3.29 mm per year, where trees in the younger stands (<100 years old) had higher mean increments than older ones (more than 100 years old) at the majority of sites. This is in line with expected age-related decreasing growth trends. Annual growth of trees at the sites without groundwater extraction was higher than that of oaks from extraction sites over the entire range of cambial ages (Figure 3). Growth of sessile oaks at the upland sites was comparable to that of oaks at sites without extraction for cambial ages <100 years old. However, annual growth of trees at upland sites was lower than that of trees at the extraction and no extraction sites for cambial ages more than 100 years old. Series lengths of the developed chronologies were highly variable (ranging from 41 to 209 years) (Table 2 and Supplementary S.2). These differences were attributable to the different calendar years when groundwater extraction started in the three study regions.

TABLE 2

Table 2. Summary statistics of the 18 stand chronologies.

FIGURE 3

Figure 3. Tree ring width (TRW) vs. cambial age. Lines are fitted splines summarizing growth for the different sites. Shaded areas around lines denote 0.95 confidence intervals.

Climate Trends

Over the last 100 years, there was a steady and significant (p < 0.05) increase in mean temperature of all seasons and in all regions (on average 0.2°C per decade for the regions of Freiburg and Emmendingen, and 0.1°C per decade for the region of Lampertheim), except for winter temperature in the Lampertheim region, where the trend was not statistically significant (Figures 4A,B and Supplementary S.4). In contrast, the trends in precipitation were less consistent and varied among seasons and between the southern regions (Emmendingen and Freiburg) and the northern region (Lampertheim). For the southern regions, summer precipitation tended to decrease over the last 100 years (p = 0.06), while in the same region, mean precipitation in winter, spring and autumn showed no particular trends (Figure 4C and Supplementary S.4). Surprisingly, winter and spring precipitation in the region of Lampertheim showed a significantly positive trend during the examined period (p < 0.05, Figure 4D and Supplementary S.4).

FIGURE 4

Figure 4. Seasonal temperature and precipitation with trend lines (regression lines) for the investigated regions: (A) seasonal mean temperatures for Freiburg and Emmendingen, (B) seasonal mean temperatures for Lampertheim, (C) seasonal mean precipitation for Freiburg and Emmendingen, and (D) seasonal mean precipitation of Lampertheim. Winter includes previous year December to February of current year; spring = the months from March to May; summer season includes the months June to August and autumn the months from September to November. Gray areas denote 0.95 confidence intervals.

Climate-Growth Relationships

We assessed responses of tree growth to climatic variability by computing bootstrapped correlation functions between the stand-wise chronologies and the different monthly and seasonal climatic variables (precipitation, temperature and three different scales of the SPEI drought index). The SPEI series, both calculated for separate months and seasons, were better indicators of the growth responses of all investigated forest stands compared to precipitation or temperature series alone (lower growth was associated with more negative SPEI values) (Figure 5). Among all calculated SPEI-series, those referring to a shorter time-scale of 6 consecutive months, correlated best with inter-annual tree-growth variations across the different stands. Overall, growth of oaks in the investigated stands seemed to be strongly associated with summer drought as indicated by the highest positive correlation of growth data with SPEI 6 ending in July of the current year and mean SPEI 6 of the three summer months (June to August). These two SPEI indices showed significant correlations with growth series for 17 out of the 18 sampled stands (Figure 5). Further, a positive and significant correlation was found between radial growth and SPEI 6 of June of the current year for all 18 stands. The longer periods of SPEI showed in almost all cases fewer and weaker correlations with annual growth compared to those found for SPEI 6.

FIGURE 5

Figure 5. Averaged Pearson's correlation coefficients of all 18 detrended stand-wise growth chronologies with mean temperature, precipitation sums and three different scales of the Standardized Precipitation Evapotranspiration Index (SPEI) (6, 12, and 24 months). The vertical dashed line separates monthly from seasonal correlations. Seasonal means are calculated for spring (March to May), summer (June to August), autumn (September to November), winter (previous year December to February) as well as for the vegetation period (from March to September). The figures within bars indicate the numbers of stands with statistically significant correlations (maximum 18).

With a few exceptions, we found no significant relationship between tree growth and temperature across the studied stands (Figure 5). However, it appears that precipitation played a more important role for tree growth than temperature, as indicated by the stronger and higher number of significant correlations (Figure 5). Most importantly, mean monthly sums of precipitation during the vegetation period were strongly correlated with growth series for 12 out of the 18 sampled stands.

Climate-growth relationships for the individual stands and site-types (related to groundwater situation) were investigated for SPEI 6, as this climatic parameter was best associated with growth in the sampled stands regardless of age and site conditions (Figure 5). Correlations of tree-ring chronologies with SPEI 6 series ending in the months May to September and with mean seasonal SPEI 6 series for summer (June to August) and the whole vegetation period (May to September) were positive and significant for all site types (Figure 6). Drought during the summer months seemed to play the most important role in controlling tree growth, in particular for the lowland sites (sites with and without extraction), where correlations between annual growth and SPEI 6 series were stronger than at the upland sites (Figure 6).

FIGURE 6

Figure 6. Averaged Pearson's correlation coefficients of the detrended tree-ring chronologies with monthly and seasonal climatic parameters by site-type (A for the extraction sites, B for no extraction sites, and C for upland sites). Horizontal dashed lines indicate levels of significance at 0.05. Vertical dashed lines separate monthly from seasonal correlations.

Remarkably, the strongest correlations between tree-growth and SPEI 6 were observed for trees growing at sites without groundwater extraction followed by trees from the extraction sites while the lowest correlations were found for trees from the upland sites (Figure 6C). For trees from sites without extraction, the strength of the relationship between growth and SPEI 6 during the summer months as well as during the vegetation period tended to decrease with average stand age (Figure 7B). On the other hand, no such trends could be observed at the extraction (Figure 7A) and upland sites (Figure 7C).

FIGURE 7

Figure 7. Pearson's correlation coefficients of the separate detrended tree-ring chronologies with seasonal averages of SPEI 6 (A for the extraction sites, B for no extraction sites, and C for upland sites). Seasonal means of SPEI 6 are calculated for spring (March to May), summer (June to August) as well as for the vegetation period (from March to September). Triangles and dots indicate significant and non-significant correlations respectively, at p < 0.05.

In addition, trees at sites without groundwater extraction showed significantly higher growth increases in years with favorable summer conditions (SPEI 6 of summer months > 1) (p < 0.01) compared to trees at the extraction and upland sites (Figure 8). However, in years with adverse summer conditions (SPEI 6 of summer months <-1), trees from the different site types showed more similar growth responses. Here, trees at the no extraction sites had slightly (although statistically not significant) stronger growth depressions than trees at the extraction and upland sites (Figure 8).

FIGURE 8

Figure 8. Absolute radial growth changes (radial growth in year t minus radial growth in year t-1) in dry (mean SPEI 6 of June, July and August <-1) and wet (mean SPEI 6 of June, July and August >1). summers. Stars indicate statistically significant differences among pairs connected with brackets (**p < 0.01, ***p < 0.001, n.s, not significant).

Temporal Variability in Climate-Growth Relationships

Temporal changes in sensitivity of radial growth to climate were assessed by computing moving-window correlation functions between the detrended stand-wise chronologies and climatic variables (precipitation, temperature and three different scales of SPEI index). Moving correlation analyses revealed a general instability and an overall strengthening of the relationship between growth and SPEI 6 over time for the investigated period (Figures 9, 10). Interestingly, this increasing sensitivity of tree-growth to drought was interrupted for the 30-year long windows starting between 1920 and 1950, at all sites (Supplementary S.5).

FIGURE 9

Figure 9. Moving correlations between monthly and mean seasonal values of SPEI 6 and growth for the different oak stands. (A for the extraction sites, B for no extraction sites, and C for upland sites). The window length is 30 years and the window offset is 1 year. Black stars indicate significance at p < 0.05. The continuous vertical lines indicate the starting date of the extraction. The dashed lines are placed 30 years before the extraction started. Note that correlation coefficients are presented here relative to the date when the extraction of groundwater began (date when the extraction began = 0). Time since extraction indicates the beginning of each 30-year window. For instance, correlation coefficient in year x is the correlation coefficient between growth and the specific climate parameter for the 30-year period starting in year x. Stand IDs are acronyms denoting the study region and site type (first and second letter, respectively) of each investigated stand. F, E, L for Freiburg, Emmendingen and Lampertheim, respectively. W indicates no extraction sites, C upland sites and E, S, and N extraction sites (E for the extraction sites in Emmendingen and Freiburg, S and N for the Southern and Northern extraction sites, respectively, in Lampertheim). Numbers in Stand IDs indicate the average tree age of each stand.

FIGURE 10

Figure 10. Moving bootstrapped correlations between annual growth and mean SPEI 6 of summer months (June to August) for the three different site-types [extraction (A), no extraction (B), and upland (C)]. Red triangles and gray dots indicate significant and non-significant correlations respectively, at p < 0.05. Each point represents the correlation coefficient of a single stand calculated for a 30-year-long window. Time on the x-axis represents the starting year of each 30-year-long window, which is relative to the beginning of the extraction (extraction date set to 0) due to different onsets of extraction among regions. Solid vertical lines indicate the dates when the extraction of groundwater began, while the dashed lines mark 30 years before the onset of groundwater extraction.

Temporal Variability in Growth Sensitivity to Drought at the Extraction Sites

The strength of the relationship between tree-growth and monthly and seasonal values of SPEI 6 during the vegetation period at the extraction sites has been increasing over the last decades, indicating a higher growth sensitivity to seasonal and in particular to summer droughts (Figures 9, 10). This is highlighted not only by the more positive values of correlation coefficients, but also by the larger number of significant correlations over time (Figures 9, 10). For each of the stands at extraction sites, the increase in growth sensitivity to monthly and seasonal values of SPEI 6 becomes especially apparent soon after the onset of groundwater extraction (i.e., as soon as the 30-year long windows approach the extraction date “0”). This is indicated by the much higher proportion of significant correlation coefficients in the post- compared to pre-extraction period (Figures 9, 10).

Temporal Variability in Growth Sensitivity to Drought at the Upland Sites

The correlation between growth and SPEI 6 was weaker for sessile oaks growing in the drier upland stands compared to pendunculate oaks at lowland sites (extraction and no extraction) for the regions of Emmendingen and Lampertheim. In contrast, correlation coefficients of trees growing in upland stands in Freiburg were more similar to those of trees from the lowland sites in this region (Figure 9). Remarkably, in the first half of the investigated period, there were only very few significant drought-growth relationships for upland stands in the regions of Lampertheim and Emmendingen (Figure 9). However, similar to the results of the extraction sites, sensitivity of growth to drought during the vegetation period increased over the last decades when combining data from all upland sites (Figure 10C).

Temporal Variability in Growth Sensitivity to Drought at the No Extraction Sites

In contrast to extraction and upland sites following the period of lower sensitivity of tree growth to drought (1920–1950), correlation coefficients at sites without extraction reached the previous levels and did not exceed them as it did at the extraction and upland sites (Figures 9, 10 and Supplementary S.5). Interestingly, there was no drop in the correlation coefficients at the site without extraction in Emmendingen (stand ID = EW77, Figure 9B) during this period. Here the overall correlation coefficients remained stable and higher than in stands from the other two regions (Freiburg and Lampertheim). These differences might be attributable to the younger age of stands at the no extraction site in the region of Emmendingen (see also Table 1).

Discussion

Our study design (Figure 2) facilitated the investigation and direct comparison of climate-growth relationships of trees growing under differing groundwater regimes. We were able to quantify the development of climate-growth relationships over approximately the past 100 years including for trees that have been exposed to lowered groundwater levels for several decades (49–88 years). The findings of this study support the idea that changes in groundwater levels can modify climate-growth relationships in groundwater-dependent forest ecosystems. Contrary to our expectation, there appeared to be no consistent relationship between tree age and radial growth sensitivity to climatic variability across sites with differing in soil water regime.

Climatic Sensitivity of Radial Growth in Oaks

The positive correlations between growth of oaks and the seasonal SPEI 6 for summer (i.e., the average SPEI 6 of June, July, and August) (Figures 5, 6) suggest that drought severity during summer but also during the first months of the current vegetation period control growth. This was found for all three investigated regions and regardless of distance to groundwater table or stand age. The observed higher sensitivity of both oak species to short-term drought stress and especially during summer months is consistent with other findings from studies across Europe (Drobyshev et al., 2008; Mérian et al., 2011; Bhuyan et al., 2017; Árvai et al., 2018) including observations from Central-West Germany (Friedrichs et al., 2008). However, sessile oaks in the drier uplands showed an overall lower sensitivity to climatic parameters compared to pedunculate oaks at the lowland sites (with and without extraction) (Figure 6C). In addition, trees at upland sites were indeed less responsive to both favorable and unfavorable summer conditions (Figure 8). This difference between upland and lowland sites might be attributable to species-specific climatic sensitivity, as earlier studies report higher climatic sensitivity for pedunculate than sessile oaks (Friedrichs et al., 2008).

Several studies suggest that inter-annual variation observed in tree-ring widths of pedunculate and sessile oaks is largely related to the variability of latewood width and to a lesser extent to that of earlywood width (Zhang, 1997; González and Eckstein, 2003; Kern et al., 2013). The strong relationship between climatic conditions in summer (June to August) and radial growth that we observed here for both oak species may be explained by the strong dependence of latewood formation on water availability during these months (González and Eckstein, 2003). Similar relationships can be expected in other ring-porous species (see for example Zhang, 1997; González and Eckstein, 2003 and references therein). In addition, our results indicate that precipitation during the vegetation period (March to September) plays an important role in controlling oak growth while temperature seems to have a comparably lower impact on growth of the studied trees. Similarly, higher growth sensitivity to drought (Palmer Drought Severity Index, PDSI) and precipitation than to temperature has been reported for oaks sampled within a large network of sites in the same region (Friedrichs et al., 2008) and for north-eastern Hungary (Kern et al., 2013).

The Influence of Tree Age on Growth Sensitivity in Oaks

Sensitivity of oaks to climatic variability at sites without groundwater extraction appeared to decrease with stand age, whereas no clear trend occurred at the other site types (Figure 7). Previous studies have yielded contradictory results concerning the effect of age on climatic sensitivity of trees (e.g., Rozas, 2005; Rozas et al., 2009; Copenheaver et al., 2011; Voelker, 2011; Konter et al., 2016). Any age-related effects have been attributed to stand dynamics and/or physiological and morphological changes with tree development (Day et al., 2002). With increasing age and height, water is transferred over longer distances within trees, which requires adjustments of water conducting elements (Ryan et al., 2006). For example, vessel size increased with cambial age in ring-porous oak species (e.g., Rita et al., 2016). However, increases in conduit diameter can render trees more susceptible to drought-induced cavitation, potentially increasing the susceptibility of older trees to drought (Ryan et al., 2006; Fonti et al., 2010; Voelker, 2011). From our data, we cannot explain why climate sensitivity of tree-growth at sites without water extraction was lower in older than in younger trees. A possible explanation could be that large trees which can store more non-structural C may be able to support radial growth to a larger extent from these reserves and are therefore less reliant on the actual climatic conditions (Galiano et al., 2011; Sala et al., 2012).

Interestingly, growth sensitivity to climatic variables in trees at sites without groundwater extraction was on average higher than in trees at the extraction and upland sites (Figures 6, 7). It may be that the comparably high groundwater tables during most months of the year cause oak trees to develop a relatively shallow root system, which in turn could have led to an increased sensitivity to summer drought as water availability in the upper soil becomes limited (Oosterbaan and Nabuurs, 1991; Leuschner and Ellenberg, 2017). Indeed, growth of trees at the sites without water extraction was most responsive to both wet and dry summers, demonstrating a higher ability of these trees to take advantage of favorable summer conditions but also slightly stronger growth depressions during periods of reduced soil water availability (Figure 8). However, annual growth of trees at sites without groundwater extraction was superior to that of trees at the extraction sites (Figure 3 and Table 2). This indicates that negative effects of summer droughts are outweighed by growth benefits from access to groundwater under favorable growing conditions (Leuschner and Ellenberg, 2017). In addition, this indicates that sensitivity in terms of the strength of climate-growth correlations on its own is not a sufficient indicator of growth stress. Only in combination with the analysis of growth deviations under favorable and unfavorable conditions can this indicator be interpreted properly. If sensitivity of radial growth increases as a result of positive growth responses to favorable growing conditions, this should not be read as a negative sign, as has been done in many previous studies of climate-growth relationships (e.g., von Oheimb et al., 2014; Caminero et al., 2018).

In contrast to sites without water extraction, sensitivity of trees to drought at the groundwater extraction sites showed no particular trend with increasing stand age. The lower sensitivity of trees between 60 and 100 years old at the extraction sites, compared to trees of similar age growing at the wetter sites, might be a consequence of their lack of positive growth responses under favorable conditions (Figure 8). This lack of response in radial growth under favorable conditions could be a consequence of higher investment of C into the root system (e.g., Wang et al., 2015) or lower amounts of stored non-structural C. However, we have no quantitative information on rooting pattern and belowground biomass in relation to groundwater levels.

Effect of Groundwater Regime on Temporal Development of Climate-Drought Relationships

Moving correlation analyses revealed non-stationary relationships between annual growth of oak trees and short-term drought indices (SPEI 6) (Figures 9, 10 and Supplementary S5). Similar changes in climate-growth relationships over time for pedunculate and sessile oaks have been observed in several studies conducted in Central-Europe (e.g., Friedrichs et al., 2008, 2009; Mérian et al., 2011; Árvai et al., 2018). We observed a substantial decrease in the strength of drought-growth relationships for the 30-year long windows between 1920 and 1950 in comparison to the period before. Likewise, Friedrichs et al. (2008) reported a similar reduction of the relationship between growth and drought for exactly the same period for both pedunculate and sessile oaks in Central-West Germany. The observed reduced sensitivity of annual-growth to drought during this period demonstrates clearly the complexity when studying climate-growth relationships. Surprisingly, despite low precipitation and higher than average temperatures during this period, growth levels increased (Friedrichs et al., 2008). Following this period, trees growing at the extraction and upland sites demonstrated a greater growth sensitivity to drought, while no such effect was found for trees at sites without extraction.

The influence of summer drought on growth of oaks at the extraction sites increased following the beginning of groundwater extraction in all regions as indicated by the stronger climate-growth relationships over time. As we did not observe such an increase in strength of the relationship between summer-drought and growth at adjacent sites without groundwater extraction over the same period, the stronger climate-growth relationships at extraction sites can be regarded as a result of lowered groundwater tables. This increased sensitivity to summer drought was accompanied by overall lower growth levels compared to trees at sites without extraction over the period following the onset of groundwater extraction (Figure 3). Additionally, oak trees at the extraction sites had a lower ability to respond positively to years with favorable summer conditions (Figure 8). The differences observed in the response of trees to summer-drought between extraction and no extraction sites corroborate findings of Raney et al. (2016) and further support the idea that negative effects of elevated summer temperatures and reduced precipitation on tree-growth can be mitigated by access to groundwater. In contrast to our results for trees at lowland sites without extraction, we observed an increasing sensitivity of radial growth to summer drought for sessile oak trees at the drier uplands over the investigated period (Figures 9, 10 and Supplementary S.5). Our results are consistent with results of other studies, which report a stronger relationship between growth of sessile oaks and drought during recent decades compared to the beginning of the twentieth century (Friedrichs et al., 2008, 2009). It is very likely that the observed trends of increasing drought sensitivity in sessile oaks are indeed driven by decreases in summer precipitation that were observed in the region of Freiburg and Emmendingen but most importantly by increases in summer temperature, which were observed for all studied regions (Figure 4).

Several additional non-climatic factors can influence radial growth of trees, potentially affecting the climatic signal laid-down in tree-ring records. These include massive defoliation through periodical outbreaks of herbivorous insects (Büntgen et al., 2009; Black et al., 2016) and masting events that divert the available C from radial growth to seed production (Di Filippo et al., 2007; Genet et al., 2010; Hacket-Pain et al., 2015). The influence of such factors on climate-growth relationships could not be tested for the time periods examined here, although it might have explained some of the noise in our data.

Conclusions

In this study, we investigated the effects of groundwater extraction on climate-growth relationships for pedunculate and sessile oak trees covering a wide range of tree ages growing under differing groundwater regimes. The evidence presented in this study suggests that groundwater extraction impairs radial growth and increases growth sensitivity of trees to summer drought while access to groundwater appears to mitigate the negative effects of increasing summer water-deficits on tree growth. In future, floodplain forests will likely experience more stress from the combined effects of increasing groundwater extraction rates and more frequent and prolonged warm and dry conditions. Therefore, the location of groundwater extraction sites as well as the rates and timing of water extraction should be carefully considered to ensure that tree vitality and ecosystem health are not impaired.

Data Availability

The datasets generated for this study are available on request to the corresponding author.

Author Contributions

JB conceived the study. JS and JB acquired funding for the project. GS, JB, and JS designed the study. GS carried out data collection, performed the analysis and drafted the manuscript. All authors discussed and interpreted the results and contributed to the writing of the final manuscript.

Funding

This work was supported by the DRIeR project funded under the framework of the Water Research Network by the Ministry of Science, Research, and Arts of the German Federal State of Baden Wuerttemberg (AZ. 7532.21/2.1.6).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank the forestry districts of the City of Freiburg, Emmendingen and Lampertheim for their support and excellent collaboration. We also like to thank Jan Geyer for his assistance in fieldwork and Dr. Willy Tegel for his valuable comments on the manuscript. The article processing charge was funded by the German Research Foundation (DFG) and the University of Freiburg in the funding programme Open Access Publishing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/ffgc.2019.00005/full#supplementary-material

References

Ahner, J., and Schmidt, M. (2011). “Modellierung der Einzelbaummortalität der Eiche unter besonderer Berücksichtigung von Grundwasserabsenkungen – Eine Untersuchung im Hessischen Ried,” in Presented at the Sektion Ertragskunde Jahrestagung 2011. Cottbus, 157.

Google Scholar

Árvai, M., Morgós, A., and Kern, Z. (2018). Growth-climate relations and the enhancement of drought signals in pedunculate oak (Quercus robur L.) tree-ring chronology in Eastern Hungary. IForest Biogeosci. For. 11, 267–274. doi: 10.3832/ifor2348-011