• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Tree sapling vitality and recovery following the unprecedented 2018 drought in central Europe

    2023-11-15 07:56:48MirelaBeloiuSchwenkeValeskaSchnlauCarlBeierkuhnlein
    Forest Ecosystems 2023年5期

    Mirela Beloiu Schwenke, Valeska Sch?nlau, Carl Beierkuhnlein

    a Institute of Terrestrial Ecosystems, ETH Zurich, Zurich, Switzerland

    b Department of Biogeography, University of Bayreuth, Universit?tsstra?e 30, 95447, Bayreuth, Germany

    c GIB Geographical Institute Bayreuth, University of Bayreuth, 95447, Bayreuth, Germany

    d BayCEER Bayreuth Center of Ecology and Environmental Research, University of Bayreuth, 95448, Bayreuth, Germany

    Keywords:Drought stress Temperate forest Deciduous trees Climate change Extreme events Soil depth

    ABSTRACT

    1.Introduction

    In the summer of 2018, Central Europe was affected by one of the most severe summer droughts since the beginning of the 21st century(Büntgen et al.,2021;Ionita et al.,2021).The frequency and intensity of extreme weather events,causing drought or heavy rain,are predicted to increase in the 21st century as a consequence of decreasing temperature gradients in the Northern Hemisphere and related stable Rossby Waves(Hoegh-Guldberg et al.,2018;IPCC,2019).In the face of global warming,vegetation is mostly affected by drought periods,particularly if these are combined with summer heat waves,as was the case in 2018.The increase in temperature and drought frequency, intensity, and duration are affecting the productivity, dynamics, regeneration, and mortality of forest ecosystems (Ciais et al., 2005; Hacket-Pain et al., 2016; Choat et al.,2018;Senf et al.,2018;Beloiu et al.,2022a).In 2018,and partially in 2019, temperate forests across Europe experienced prolonged droughts and heat waves that caused drought-related tree mortality(Brun et al., 2020; Schuldt et al., 2020; Hari et al., 2020; Beloiu et al.,2022b).Consequently, there is growing concern that large parts of forested areas will be affected due to the increasing frequency and intensity of extreme events(Albrich et al.,2020;Hammond et al.,2022).

    Temperate forests in Europe are characterized by a low diversity of tree species (Beierkuhnlein, 2007).Productivity is limited to a short vegetation period with moderate seasonal conditions in terms of temperature and moisture.It is therefore essential to understand the effects of repeated drought on the survival and distribution of the most important species.An increase in tree mortality is observed in Central Europe in recent years (Senf et al., 2021).The summer drought of 2018 severely affected temperate forests in Central Europe(Vitasse et al.,2019;Schuldt et al., 2020).In Germany, forest covers 32% of the land area with 54%being dominated by coniferous,31%by broadleaved and 13%by mixed forests (BMEL, 2012; Holzwarth et al., 2020).Of the coniferous forests,Picea abies and Pinus sylvestris are the most dominant species in Germany.They have been planted outside their natural range (in the lowlands)since the 18th century, where they replaced temperate broadleaved species such as Fagus sylvatica and Quercus spp.(Seim et al., 2022;Leuschner et al.,2022).Nowadays,however,these coniferous forests are severely affected by drought and subsequent bark beetle outbreaks(Seidl et al., 2016; Buras et al., 2018; Hlásny et al., 2021).Temperate broadleaved deciduous forests with natural regeneration dynamics of native tree species are expected to be more resilient to climate change.Many studies on broadleaved species focus on F.sylvatica, a dominant tree species in Central Europe.The 2018 drought challenged F.sylvatica's resilience under climate change(Frei et al.,2022;Langer and Buβkamp,2023), however the drivers underlying its sensitivity to drought across scales are not well known (Leuschner, 2020).Moreover, knowledge of the in-situ dynamics of native deciduous tree species populations is scarce(Beloiu et al.,2022b).

    The capacity of forests to withstand a disturbance,adapt,and recover their functional capacity is a key concept in ecology and is referred to as resilience (Holling, 1973; Turner, 2010).Forest recovery specifically refers to the process of returning the forest ecosystem to a functional state, where it can provide the ecological services and functions that it provided before the damage occurred(Ingrisch and Bahn,2018;Gessler et al.,2020).Forest resistance,on the other hand,refers to the capacity of a forest ecosystem to withstand and maintain its state and functions in the face of disturbance (Ingrisch and Bahn, 2018).Knowing which tree species are more resistant,versus those which are most affected,and how each recover from extreme droughts will help us better understand the processes that enhance forest stability and resilience.Furthermore, this plays a key role in understanding how to develop nature-based solutions to mitigate risk and prevent damages (Beierkuhnlein,2021).

    As young trees represent the potential forest of the future, the response of saplings to drought stress is important for assessing the potential for and limits to forest ecosystem adaptation(Beloiu,2022).It is not fully understood which tree species and species assemblages are most adaptable to increasing temperatures and summer drought (McDowell et al., 2008).Numerous research employs greenhouse experiments and model simulations to investigate the sensitivity of trees to climate extremes (Schall et al., 2012; Mette et al., 2013; Buhk et al., 2016).However, it is important to consider tree response to drought under natural conditions in the field including a variety of site conditions and demographic stages, as they have a major impact on growth and survival(Gimbel et al.,2015;Beloiu et al.,2020).

    Enhanced growth of understory vegetation cover in forests is often related to forest type,light availability,forest disturbances,topography,soil moisture, and nutrients.In the absence of forest encroachment, the structure and density of understory vegetation remain relatively stable over time(Cole et al.,2017).Litter cover and thickness are influenced by forest type, tree health, age, climate, and soil conditions (D?lle et al.,2017; Petritan et al., 2020).Deciduous trees produce more litter than conifers due to their annual leaf shedding,but their decomposition tends to be faster than that of conifers(Prescott et al.,2000;Barba et al.,2016).Under drought stress, litter cover can play both positive and negative roles in tree growth and survival.On the one hand,litter cover can act as a protective layer, reducing water loss from the soil and helping to maintain soil moisture levels(Knapp and Seastedt,1986;Deutsch et al.,2010).This can be especially beneficial during dry periods when water is scarce.On the other hand,a thick litter layer of especially deciduous trees has the potential to retain rainwater, hence limiting the water that reaches the soil and the roots (Walsh and Voigt, 1977).Understory vegetation increases soil water infiltration and reduces soil evapotranspiration(Dubbert et al.,2014).Thus,these factors,along with soil depth and water-holding capacity,influence soil drought(Miller and Poole, 1983).Even though understory vegetation and litter cover influence water availability both near the forest floor and in the topsoil, the influence of these factors on tree recovery after a drought is not known.

    Here we assess the response of tree sapling species to the 2018 drought and their subsequent recovery by 2021.We evaluate the influence of soil depth, understory vegetation, litter cover, and forest cover across different forest types on recovery patterns.Tree vitality was measured in 170 forest plots in central-eastern Germany during 2018-2022.We focus mainly on deciduous tree species, as these are naturally the dominant species in the temperate biome of Central Europe and are considered more suitable in the face of climate change.Specifically, we aim to address the following questions: (i) What are the differences between tree species in terms of their resistance to and recovery from the 2018 drought? (ii) What is the influence of soil depth, litter,understory vegetation, and forest cover across deciduous, mixed, and coniferous forests on tree resistance to drought and recovery from drought? Answering these questions will allow us to better understand the link between site conditions and tree species recovery which is critical for forest regeneration and climate change adaptation.

    2.Materials and methods

    2.1.Study area

    The study area is located in the Central Germany,spanning across the northern part of Bavaria (Fig.1).In this area, the greatest leaf loss and crown defoliation were observed in 2019 and 2020(StMELF,2020).The area is dominated by deciduous, mixed, and coniferous forests.The predominant tree species in the canopy are European beech (Fagus sylvatica L.), pedunculate oak (Quercus robur L.) and sessile oak (Quercus petraea (Matt.) Liebl.), Norway spruce (Picea abies (L.) Karst.) or Scots pine(Pinus sylvestris L.).The petrography of the bedrock is very diverse,including clay and siltstones,limestones,as well as sandstones(Doppler et al.,2004).The diverse geological conditions result in a broad spectrum of soil types in the study area.

    The climate of the study area is characterized by an annual average precipitation of 957 mm,and an annual average temperature of 7.1°C in the period 1970-2020 (DWD, 2021a, b; Beloiu et al., 2022b).Over the past 70 years, the study region experienced a temperature increase of about 1.2°C and a shift in the precipitation regime,with significantly less precipitation in the spring (Fig.2a and b).Between 2017 and 2021,substantial precipitation deficits occurred repeatedly (Fig.2c).The deficit of groundwater recharge(in the period 2011 to 2020 compared to the reference period 1971 to 2000)amounts to decreases of between 17%and 22%(LfU Bayern,2021).Historically,the most severe droughts were recorded in 1971-1974, 2003, and 2018.However, 2018 was an extremely dry year.In 2018,the study region experienced a combination of high temperatures and low precipitation in both spring and summer.Similarly, temperature and precipitation varied greatly in 2019.The precipitation and temperature data were extracted for the study sites from the raster dataset (1 km × 1 km) of the German Meteorological Service.

    2.2.Data acquisition

    Fig.1.Overview of the study area in Upper and Lower Franconia with the 170 sample plots.Forest cover based on digital landscape model (DLM) 1:250,000(BKG, 2019).

    Fig.2.Monthly precipitation and temperature patterns across the study region between the mean of the period 1961-1990 and 1991-2020(a,b)and between January 2017 and September 2021 based on the climatic data from the German Meteorological Service (c).

    Fieldwork was conducted in the years 2018,2019,2020,and 2021,in late summer from August to September.During this time, 170 plots between 261 and 722 m a.s.l.were monitored.The 170 plots were distributed over an area of approximately 15,680 ha(Fig.1).During the first data collection, at each plot, 10 young trees were marked (at least 1.3 m tall and DBH <10 cm)within a randomly selected area of 100 m2(Beloiu, 2022).In total, data were collected from 1,531 saplings.Tree vitality was recorded annually and assigned to a vitality class, and thereby assessed in terms of crown defoliation, as this is a visible indicator of drought stress.A total of 40 tree species were tagged,and individual statistics were performed using only the seven species with the highest sample size(n >50).Thus,after data cleaning,i.e.,removing the damaged or not re-found saplings, a total of 1,483 saplings were considered for analysis.Table 1 lists the species studied,their abundance,and their vitality during the study years.

    The species considered for the study represent widely distributed species in the study area and the species with the largest sample sizes.The species are sycamore maple (Acer pseudoplatanus L.), a widely distributed species, and F.sylvatica, the dominant species in Central Europe under natural conditions, which are both considered drought sensitive (Leuschner et al., 2001; San-Miguel-Ayanz et al., 2016).The others are mountain ash (Sorbus aucuparia L.) and silver birch (Betula pendula Roth.),as two pioneer tree species,common hornbeam(Carpinus betulus L.), whose distribution is limited by high summer temperatures,common hazel (Corylus avellana L.), which is adapted to warm climates(San-Miguel-Ayanz et al.,2016),and Quercus spp.,which are considered to be drought tolerant (H?rdtle et al., 2013).The two oak species Q.petraea and Q.robur were recorded separately but grouped under Quercus spp., as their characteristics show great plasticity and hybrid forms arise(Aas,1991;Jensen et al.,2008).

    The vitality classes are defined as follows:1=undamaged:no leaves are damaged by drought or show significant color changes;2=partially damaged: individual damaged leaves and branches of the young trees;but most of the leaves are undamaged; 3 = completely damaged: all leaves and branches of the plant show drought damage, this class includes specimens that died during or after the drought.

    The cover of understory vegetation,litter,and bare soil was recorded as a percentage for the 10 m × 10 m plot area in the summer of 2020(area fraction of litter,bare soil,and understory vegetation,respectively together summed to 100%).The percentage of these site-specific factors is an indicator of site microclimate and could influence species recovery in 2021.About 1-2 years following the drought, the die-off leaf mass decomposes and returns to normal levels (Prescott et al., 2000; Jacob et al., 2009).Hence, these measurements were done in the third year after the drought to avoid the potential impact of the 2018 drought on understory vegetation and litter cover.The understory vegetation is defined here by the herbaceous layer with forbs,ferns, non-grasses herbaceous species,and small shrubs from the family of Rosa spp.and Rubus spp., excluding tree saplings (D?lle et al., 2017).Besides mosses and ferns, a total of 41 plant families were represented within the plots.Forest canopy cover was recorded in percentage for each plot and the forest type(deciduous,coniferous and mixed)was determined.

    Table 1Number of saplings in each vitality class per year.No., the total number of tree saplings.Vitality classes: 1 = undamaged; 2 = partially damaged; 3 = completely damaged.

    Soil depth (horizon O, A and B in cm) defines the depth from the topsoil profile to the weathered bedrock and the depth of the loose material to which the roots can easily penetrate.The Pürckhauer gouge auger was used to determine the thickness of the soil.

    2.3.Data analysis

    The ability to recover after the severe 2018 drought of each tree species is expressed by the Recovery Index (RI) (Lloret et al., 2011).Recovery is estimated as the reverse ratio of vitality after and during the 2018 drought.The RI was calculated for the studied years, i.e.2018 to 2019 and 2018 to 2021.This resulted in three recovery classes (RI classes)with values 0,1,and 2,namely RI-0,RI-1,and RI-2.RI-0 implies an incomplete recovery in sapling performance and a decrease in vitality after the drought event; RI-1 implies a consistent vitality over the respective period; and RI-2 indicates that the performance of the respective tree species following drought is better in comparison to the drought period 2018 and thus that the species has improved its vitality.To statistically assess the drivers of recovery, ordinal logistic regression models from MASS package (Ripley, 2023) were used, where the response variable was sapling recovery from 2018 to 2019 and the explanatory variables were litter cover,understory vegetation cover,soil depth,forest type and forest canopy cover.

    Changes in sapling vitality between the years and between vitality classes were assessed using Fischer's exact test for count data and pairwise comparison(Fisher,1934;Herve,2023).To test the normality of the data distribution, the Shapiro-Wilk test was used.The influence of the environmental factors considered on recovery was examined using the Kruskal-Wallis test.Dunn's test for multiple comparisons was performed as a post hoc test(p-values were adjusted with the Bonferroni method)to determine the significantly different groups.The considered parameters were litter, understory vegetation, and soil depth.For all statistical calculations,the software RStudio(version 4.1.0),with additional packages FSA, RVAideMemoire (Herve, 2023), rcompanion (Mangiafico, 2019),ggpubr,psych,ggplot2,and car,was used.

    3.Results

    3.1.Sapling's resistance and recovery patterns

    The vitality of most tree species has significantly increased after the 2018 drought - showing a general recovery.Nevertheless, a significant but slow delayed decline in vitality was observed from 2019 to 2020,while vitality remained the same between 2020 and 2021(Fig.3).During the 2018 drought, out of a total of 1,483 saplings, only 35% showed strong resistance to drought,whereas 52%of all specimens were partially damaged and 13%were completely damaged(Fig.3).One year later,in 2019,34%of the partially damaged and 7%of the completely damaged saplings recovered.In the following years,2020 and 2021,>70%of the saplings maintained their high vitality.However,from 2019 to 2021,the percentage of partially and completely damaged saplings increased by 4% and 5%,respectively (Fig.3).The change within the vitality classes between the study years was significant(Fig.4a).

    Fig.3.Vitality of tree saplings from 2018 to 2021 and their percentage per class.

    3.2.Species-specific patterns

    Out of the seven species studied,F.sylvatica showed the highest resistance to the 2018 drought (47% of saplings undamaged), followed by B.pendula(40%),A.pseudoplatanus(38%),Quercus spp.(35%),C.avellana(35%),and C.betulus(15%)and S.aucuparia(8%)(Fig.4b-h).Fagus sylvatica showed a significant recovery from summer 2018 to summer 2019 andmaintenanceofthisconditioninthefollowingyears(Fig.4b).Similarly,the other species, B.pendula, A.pseudoplatanus, Quercus spp., and C.avellana, had a high recovery from 2018 to 2019 and rather stable maintenance of their vitality over the next years (Fig.4c-f).However,during summer of 2021, A.pseudoplatanus, and Quercus spp.recorded decreasing vitality compared to 2020(p <0.05).Carpinus betulus showed low resistance in 2018,but 80%of saplings recovered in 2019 and maintained their vitality over the years(Fig.4g).In contrast,the pioneer species S.aucuparia showed the lowest resistance in 2018, followed by a strong recovery in 2019 and a sharp decline in vitality in 2020 and 2021(Fig.4h).

    3.3.Drivers of sapling recovery

    Generally, saplings maintained their vitality when grew on deeper soils(p <0.05;Fig.5a).Coniferous forests had significantly deeper soils than deciduous or mixed forests(Appendix A, Fig.S1c).Saplings recovered well on shallow soils, but there were also saplings that did not recover on either soil depth, with the pattern tending towards shallow soils(Fig.5b).The significance of the relationship between soil depth and sapling recovery was mainly determined by plots with F.sylvatica (p <0.001), with some patterns present for the other species but not significant(p <0.05).

    Fig.4.Vitality status of the tree species' saplings during and after the 2018 drought in the years 2019, 2020, and 2021.Letters in the figures indicate significant differences between groups based on Fisher's exact test.

    Tree saplings growing at sites with high litter and low understory vegetation cover were most affected by drought in 2018, and they recovered better in sites with understory vegetation cover that was higher than litter cover(Fig.6a and b and Appendix B,Models 1 and 2).Specifically,the recovery of F.sylvatica,C.betulus,and A.pseudoplatanus after the drought was significantly influenced by understory vegetation (i.e.the herbaceous layer formed mainly by forbs and grasses, excluding tree saplings)and litter cover(Fig.6a-c),but no influence was found for the other species,B.pendula,Quercus spp.,C.avellana and S.aucuparia(p >0.05).We found no significant correlation between soil depth,precipitation,and understory vegetation,as well as between temperature and litter cover (p >0.05).Fagus sylvatica saplings were more resilient and recovered at sites with low litter but high vegetation cover(p <0.0001)(Fig.6d and e),but failed to recover at sites with high litter and low understory vegetation(Fig.6f).Acer pseudoplatanus saplings preserved their vitality from 2018 to 2021 in sites with low litter cover and high vegetation cover(Fig.6g-i).At sites with high litter and low vegetation cover,part of the A.pseudoplatanus saplings either did not recover at all or others recovered significantly.Hence,saplings at sites with high litter and low vegetation cover were most affected by drought in 2018.Carpinus betulus recovered better after the drought in sites with less litter cover and more vegetation cover (Fig.6j and k).However, sites with low litter and low understory vegetation cover presented both saplings that maintained vitality and saplings that were damaged but did not recover(Fig.6l).

    Fig.5.(a) Vitality of tree saplings in 2018 and (b) recovery of tree saplings from 2018 to 2021 as a function of soil depth.

    Fig.6.Differences in recovery from 2018 to 2021 given the litter and understory vegetation(i.e.ground layer vegetation)for all species and for Fagus sylvatica,Acer pseudoplatanus,Carpinus betulus;illustration of tree sapling recovery on sites with low(-)vegetation and(+)high litter and sites with high(+)vegetation and(-)low litter cover according to the statistical analysis.

    In addition, litter and understory vegetation cover showed clear patterns depending on the forest type.Litter cover,for example,increases from coniferous to mixed and deciduous forests, while understory vegetation cover and soil depth decrease(Appendix A,Fig.S1a,b,c).The 2018 summer drought had a greater impact on saplings in deciduous forests than on mixed and coniferous forests.Sapling recovery was higher in mixed and coniferous forests than in deciduous forests.Further, saplings growing in sites with high forest canopy cover had also reduced recovery compared to saplings growing in sites with low forest canopy cover(Appendix B,Models 3 and 4).

    4.Discussion

    4.1.Sapling resistance and recovery from the 2018 drought

    During the 2018 drought summer, considerable drought stress damage was observed across many forests in Central Europe (Schuldt et al.,2020).In this study,we assessed(1) the resistance and recovery of tree saplings growing in temperate forests during and after the 2018 summer drought and (2) the influence of soil depth, understory vegetation, and litter cover on sapling recovery.Our results show that 64% of all monitored tree saplings exhibited damage during 2018 summer.Nevertheless,by late summer 2019, the majority of saplings in each species had recovered(with only 21%still displaying damage).

    From 2019 to 2020, sapling vitality did not change substantially,except for S.aucuparia, which rapidly deteriorated in its vitality within one year.This species has an especially shallow root system, which therefore might benefit earliest from percolating rainwater.However,no such benefit was evident after 2019.The rainy growing season of 2021 thus did not have a significant positive effect on tree recovery visible in late summer 2021.This could be explained by the fact that, during the dry period from 2018 to 2020, tree species mainly invested in their belowground root biomass to increase water availability, which is why they showed more drought stress-induced damage aboveground,such as biomass shedding(Pretzsch et al.,2013).Furthermore,another study on mature trees showed increasing growth 1-5 years after the precipitation year (Jiang et al., 2019), while another one showed delayed consequences of drought,such as reduced growth and incomplete recovery for 3-4 years after severe drought(Anderegg et al.,2015).

    Further monitoring and the comparison of these results with results across different forest ecosystems could provide us a more general picture of the trends in saplings recovery subsequent drought stress under natural conditions.In this context,assessing forest responses to extreme climate events should be integrated into a forest monitoring network in order to allow continuous and comparable measurements between different forest ecosystems and across time (Zweifel et al., 2023).Moreover, more automatic approaches, should be considered for quantitative assessments of tree species health over larger areas(K?lin et al.,2019;Sandric et al.,2022;Beloiu et al.,2023).

    The two species C.betulus and C.avellana seem to have benefited directly from the higher precipitation in 2021.Similarities exist concerning their root system,which is medium-deep rooted(Kutschera and Lichtenegger,2002).Further,the recovery of tree species should not be considered in a generalizing manner,but rather as a function of the site and environmental conditions.Also, the performance over the years varies greatly between species, but also between individuals within a species.Although many studies attribute high drought sensitivity to F.sylvatica (Durrant et al., 2016; Martinez del Castillo et al., 2022),saplings in the study area have shown high resilience to drought.

    4.2.Drivers of sapling recovery

    Soil depth showed to be an important indicator of sapling resilience to drought, particularly for F.sylvatica, as it influences the availability of water and nutrients for plant roots (Miller and Poole, 1983).However,saplings growing in shallower soils recovered well after the 2018 drought stress.In general, deeper soils can store more water and nutrients,providing a buffer against drought stress.Further,coniferous forests with high understory vegetation are an indicator of more resilient forest sites compared to deciduous forests with high litter cover.In a previous study,Beloiu et al.(2020)showed that saplings with a larger diameter at breast height(dbh)and saplings growing in sites with a larger tree canopy cover were less resistant to drought.The high percentage of sapling recovery after the 2018 drought can be also explained by the increasing precipitation in spring 2019(Fig.2c)(Beloiu et al.,2022b).High recovery was possible even on shallow soils as a result of the sapling's capacity to take up water from the uppermost soil layer (Gessler et al., 2022).Other studies on mature trees and seedlings have shown a strong tree growth increase after drought events particularly for trees severely impacted(Pretzsch et al., 2013;Kunz et al.,2016).

    4.2.1.Fagus sylvatica and Carpinus betulus recovery

    Fagus sylvatica and C.betulus recovered better from the 2018 drought when there was both less litter cover and more understory vegetation at their site.First,this shows a correlation between damaged trees and high litter on the forest floor, confirming that leaf litter is a serious sign of drought stress.F.sylvatica, in turn, contributes to a widespread litter accumulation by an above-average litter production(about 900 g·.cm-2litter per year)(Kutschera and Lichtenegger,2002).The forest type and thus the origin of leaf biomass on the forest floor had a significant impact on the recovery of the saplings, as coniferous and deciduous trees have different decomposition rates and water retention.

    Even in the same environment, tree species have diverse root adaptations to cope with the stress imposed by drought (Maseda and Fernandez, 2006; Lübbe et al., 2017).Although it has been demonstrated that F.sylvatica and A.pseudoplatanus expand their water intake deeper as the soil dries up, the primary water uptake zone is just the top 0-0.7 m(Brinkmann et al.,2019).Another study revealed that F.sylvatica saplings used more than 50%of their water from the uppermost 5 cm soil layer at the beginning of the 2018 drought(Gessler et al.,2022).The deep-rooted Q.petraea preferentially takes up water from deeper soil layers most of the time,but shifts its water uptake to the percolated topsoil in response to rain.Whereas F.sylvatica takes up water from both the upper and deeper parts of the soil before, during, and after the precipitation event(Volkmann et al., 2016), using more than half of the precipitation that fell during the current growing season (Brinkmann et al., 2018).To quickly absorb rainwater,F.sylvatica concentrates its water intake on the top soil,especially during drought conditions.

    When the plot area has a high litter cover,it indicates that the humus forest soil is overlain by an additional layer through which rainwater must percolate before it reaches the topsoil,and thus the main root mass of F.sylvatica.On the one hand, this additional percolation duration represents time, especially under summer temperatures when much of the water is at risk of evaporating and is thus unable to reach the soil and roots.Walsh and Voight(1977),for example,showed that up to 25%of the rainwater from a low-intensity precipitation event can be completely retained by the soil through a thick layer of F.sylvatica leaves.As a result,a high litter cover,as is the case in deciduous forests,has negative effects due to water retention during precipitation events.In the situation of severe water shortage,as was the case from 2018 to 2020,the litter layer could determine the severity of tree damage and explain the significant differences among recovery classes.On the other hand, exposed soil(without litter cover) also means a higher risk of desiccation, allowing more direct evaporation from the upper substrate, but this was not the case in our study.

    More understory vegetation can contribute to a wetter microclimate near the soil as this layer of vegetation can serve as an evaporation barrier or shade provider just above the soil.Experiments concluded that the understory vegetation layer inhibits soil evapotranspiration and does not reduce the water-use efficiency of deciduous forests, even when evapotranspiration is high (Dubbert et al., 2015).For this reason,F.sylvatica may have recovered significantly better in areas with more understory vegetation, consistent with its stable performance at the highest vegetation cover (Fig.5).Thus, the species may be more dependent on the support of a moist-cool microclimate near the ground when water is scarce due to evaporation from the topsoil.

    The differences between RI classes are less significant for C.betulus than for F.sylvatica.Compared to F.sylvatica,C.betulus vitality remained the same in sites with high litter cover and low understory vegetation.Therefore, litter cover did not have a drastically negative effect on the performance of C.betulus.While litter cover can limit water percolation,it can also shade the soil, making the tree reliant on previous groundwater.Carpinus betulus does not differ from F.sylvatica in many characteristics on the surface, but there is a clear difference in terms of understory morphology:C.betulus roots deeper than F.sylvatica(K?stler et al.,1968;Warda,1998).Especially at a young age(up to 15 years),the considerable depth development of its taproot is striking,so that with a little dense substrate,depths of 1.2-1.3 m can be reached early(Gulder,1996).Thus,unlike F.sylvatica,it can generally draw water from deeper soil layers and rely less on percolating rainwater.This would be accompanied by the fact that a higher litter layer has a stabilizing effect on the vitality of C.betulus,whereas the vitality of F.sylvatica deteriorates with a higher one-C.betulus is not primarily dependent on percolating water in the topsoil and may instead benefit from the moisture-retaining effect of a litter layer.

    4.2.2.Acer pseudoplatanus recovery

    The isohydric species A.pseudoplatanus may benefit from additional nutrients from litter during drought(Pretzsch et al.,2012).However,the species showed more resistance with less litter cover.There are currently still some gaps in our knowledge about the exact nature of the root system of A.pseudoplatanus.A field study in the Franconian Forest,which is part of the study area,found the largest root proportion in the uppermost humic soil layers, but surprisingly high root depths were also observed(Nordmann, 2009).Kutschera and Lichtenegger (2002) also clearly attribute the largest proportion of its root mass to the upper soil layers.Furthermore, they state that it has the requirement of a balanced water supply,which is why very permeable soils should be saturated,at least in the deeper root zone.The problem is that groundwater levels have been falling constantly throughout Bavaria for several years, and lower soil layers thus tend to be less saturated with water.Because of this,it would make sense to focus root development on the upper soil.It is well known that A.pseudoplatanus responds to water deficiency by forming a very shallow and extraordinarily intensely branched root system (Dahmer,1997; Jensen et al., 2008).Like F.sylvatica, when the soil dries out, the area largely used for water uptake is limited to the upper 0-0.7 m(Brinkmann et al., 2019).A.pseudoplatanus may therefore rely more on percolating rainwater instead of subsoil water.

    Similar to F.sylvatica, the significant vitality deterioration of A.pseudoplatanus saplings on sites with less understory vegetation could be due to the degree of moisture and coolness removed from the topsoil by low litter cover being compensated by a higher cover of understory vegetation.In addition,understory vegetation can have a positive effect on the rootability of the soil (Chen et al., 2004).One of the few site criteria for A.pseudoplatanus to grow is that the soil be adequately rootable (Nordmann, 2009).While the extensive occurrence of certain grass species has often been associated with the poorer performance of tree species,the positive effects of a nitrate flora on its growth have been emphasized (Jensen et al., 2008).Our results showed that A.pseudoplatanus saplings that grew on sites with more understory vegetation were more damaged by the drought, possibly due to competition for resources,but also recovered better after the drought.

    5.Conclusions

    The 2018 drought was characterized by higher temperatures and significantly lower precipitation, which affected tree species across Central Europe.Our results showed that F.sylvatica had the highest resistance to drought,followed by B.pendula,A.pseudoplatanus,Quercus spp., C.avellana, C.betulus, and S.aucuparia.This information can be valuable for forest managers when selecting tree species for reforestation or afforestation projects in areas prone to drought.Generally, the deciduous tree species studied had high recovery capacity: despite being affected by the severe drought of 2018,they demonstrated the capacity to recover in 2019.The microclimate favoured by the local synergies between forest types(coniferous and mixed),soil depth and understorey vegetation, was most favourable for the resistance of the saplings and their recovery from drought.This finding highlights the resilience of deciduous tree species in response to extreme weather events, which is important to consider in forest management and conservation planning.Tree saplings growing in sites with deeper soils and more understory vegetation were more resistant to drought than those growing in sites with more litter cover.

    The slight decrease in vitality observed in several species by 2021,despite recovery from the 2018 drought, emphasizes the need for longterm monitoring of forest ecosystems.Continuously tracking the health of tree species can provide valuable insights into their long-term response to drought and other stressors.These findings can guide researchers and forest practitioners in implementing adaptive management strategies to enhance forest resilience to drought events.This may include selective tree species planting, soil management practices, and promoting understory vegetation to support tree health and survival during extreme climate events.Understanding how temperate deciduous tree species recover from drought stress is important for maintaining forest functions and adapting temperate forests to climate change.

    Funding

    This research received no external funding.

    Data availability

    The data used are already published on the TRY database (www..try-db.org) and by Beloiu et al.(2022)in the supplementary materials.

    Authors’contributions

    Mirela Beloiu: Conceptualization, Data acquisition, Investigation,Methodology,Visualization,Writing-Original Draft;Valeska Sch?nlau:Conceptualization, Data acquisition, Methodology, Formal analysis,Writing - Original Draft; Carl Beierkuhnlein: Conceptualization, Data acquisition,Supervision,Writing-Review&Editing.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    The authors would like to thank Reinhold Stahlmann and Sophy Kley for their support during the field measurements.

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.i.org/10.1016/j.fecs.2023.100140.

    少妇被粗大猛烈的视频| 日韩欧美在线乱码| 我要搜黄色片| 国产高清视频在线播放一区| 国产乱人偷精品视频| 黄色视频,在线免费观看| 天天躁夜夜躁狠狠久久av| 中文字幕人妻熟人妻熟丝袜美| 午夜视频国产福利| 午夜a级毛片| 在线观看66精品国产| 色尼玛亚洲综合影院| 精品久久久久久久末码| 欧美又色又爽又黄视频| 国产伦在线观看视频一区| a级毛色黄片| av在线播放精品| 午夜a级毛片| 成年版毛片免费区| 国产在视频线在精品| 我要搜黄色片| 国内揄拍国产精品人妻在线| 人妻制服诱惑在线中文字幕| 女的被弄到高潮叫床怎么办| 久久精品综合一区二区三区| 一个人看的www免费观看视频| 精品无人区乱码1区二区| 亚洲熟妇熟女久久| 亚洲中文字幕日韩| 午夜久久久久精精品| av卡一久久| 国产蜜桃级精品一区二区三区| 亚洲性久久影院| 最近的中文字幕免费完整| 少妇的逼好多水| 插逼视频在线观看| 亚洲av五月六月丁香网| 中文字幕久久专区| 最近手机中文字幕大全| 亚洲综合色惰| 亚洲精品一卡2卡三卡4卡5卡| 女的被弄到高潮叫床怎么办| 欧美激情在线99| 一a级毛片在线观看| av在线亚洲专区| 97超级碰碰碰精品色视频在线观看| 美女xxoo啪啪120秒动态图| 男人舔奶头视频| 日韩在线高清观看一区二区三区| 精品无人区乱码1区二区| 一个人免费在线观看电影| 日本色播在线视频| 欧美激情久久久久久爽电影| 18禁裸乳无遮挡免费网站照片| 免费在线观看影片大全网站| 欧美潮喷喷水| 一本久久中文字幕| 国产成人精品久久久久久| 如何舔出高潮| 美女xxoo啪啪120秒动态图| 久久精品久久久久久噜噜老黄 | 日韩一本色道免费dvd| 看片在线看免费视频| av女优亚洲男人天堂| 少妇的逼好多水| 此物有八面人人有两片| 97超碰精品成人国产| 美女 人体艺术 gogo| 亚洲最大成人av| 18禁黄网站禁片免费观看直播| 欧美激情久久久久久爽电影| 欧美绝顶高潮抽搐喷水| 天天躁日日操中文字幕| 少妇丰满av| 黄色配什么色好看| 欧美极品一区二区三区四区| 性插视频无遮挡在线免费观看| 日韩欧美精品v在线| 日日摸夜夜添夜夜添小说| 国产真实乱freesex| 成人午夜高清在线视频| 久久人人爽人人片av| 亚洲av美国av| 99久久久亚洲精品蜜臀av| 在线国产一区二区在线| 久久久久性生活片| 丝袜美腿在线中文| 久久久久久伊人网av| 一进一出抽搐gif免费好疼| 免费看光身美女| 午夜久久久久精精品| 简卡轻食公司| 中国国产av一级| 草草在线视频免费看| 国产高潮美女av| 97超碰精品成人国产| 久久99热6这里只有精品| 18禁在线播放成人免费| 国产真实伦视频高清在线观看| 久久久久国产精品人妻aⅴ院| 女的被弄到高潮叫床怎么办| 国产成人影院久久av| 国产男人的电影天堂91| 免费观看人在逋| 精品免费久久久久久久清纯| 丰满人妻一区二区三区视频av| 国产精品久久久久久亚洲av鲁大| 成人综合一区亚洲| videossex国产| 九九久久精品国产亚洲av麻豆| 免费在线观看影片大全网站| 国产伦精品一区二区三区四那| 久久精品夜色国产| 亚洲18禁久久av| 久久人人精品亚洲av| 久久久久国产网址| 可以在线观看毛片的网站| 特大巨黑吊av在线直播| 男人的好看免费观看在线视频| 亚洲自偷自拍三级| 欧美日本亚洲视频在线播放| 久久久久国产精品人妻aⅴ院| 日本三级黄在线观看| 1024手机看黄色片| 成人特级av手机在线观看| 少妇人妻一区二区三区视频| 国产视频一区二区在线看| 亚洲第一电影网av| 欧美丝袜亚洲另类| 99久久精品热视频| 欧美色欧美亚洲另类二区| 国内揄拍国产精品人妻在线| 中文字幕av成人在线电影| 天天躁日日操中文字幕| 五月伊人婷婷丁香| 蜜臀久久99精品久久宅男| 亚洲第一区二区三区不卡| 日韩,欧美,国产一区二区三区 | 日韩,欧美,国产一区二区三区 | 国产爱豆传媒在线观看| 波野结衣二区三区在线| 日韩大尺度精品在线看网址| 成年av动漫网址| 干丝袜人妻中文字幕| 可以在线观看的亚洲视频| 嫩草影院入口| 免费观看人在逋| 国产黄a三级三级三级人| 日韩制服骚丝袜av| www日本黄色视频网| 男女啪啪激烈高潮av片| 人人妻人人澡人人爽人人夜夜 | 九色成人免费人妻av| 亚洲成人中文字幕在线播放| 两个人的视频大全免费| 国产熟女欧美一区二区| 黄色配什么色好看| 国产成人一区二区在线| 免费人成在线观看视频色| 亚洲第一电影网av| 精品午夜福利视频在线观看一区| 97超级碰碰碰精品色视频在线观看| 午夜精品在线福利| 一级毛片我不卡| 亚洲国产精品成人综合色| 日韩欧美免费精品| 嫩草影院精品99| 国产精品日韩av在线免费观看| 亚州av有码| 精品一区二区免费观看| 国产午夜精品久久久久久一区二区三区 | 一区福利在线观看| 精品人妻视频免费看| 午夜福利在线在线| 免费在线观看影片大全网站| 女同久久另类99精品国产91| 国语自产精品视频在线第100页| 国产探花极品一区二区| 国产毛片a区久久久久| 一个人看视频在线观看www免费| 欧美高清成人免费视频www| 亚洲三级黄色毛片| 久久久久久久久久久丰满| 精品国内亚洲2022精品成人| 国产伦一二天堂av在线观看| 观看免费一级毛片| 国产精品永久免费网站| 免费黄网站久久成人精品| 成年免费大片在线观看| 亚洲四区av| 欧美成人一区二区免费高清观看| 嫩草影院新地址| 观看美女的网站| 99九九线精品视频在线观看视频| 亚洲国产欧美人成| 亚洲自偷自拍三级| 在线免费观看的www视频| 色综合亚洲欧美另类图片| 国产精品av视频在线免费观看| 美女免费视频网站| 国产真实伦视频高清在线观看| 亚洲乱码一区二区免费版| 蜜臀久久99精品久久宅男| 人人妻,人人澡人人爽秒播| 91午夜精品亚洲一区二区三区| av.在线天堂| 1000部很黄的大片| 99久久久亚洲精品蜜臀av| 悠悠久久av| 欧美一区二区精品小视频在线| 人妻夜夜爽99麻豆av| 国产乱人视频| 午夜影院日韩av| av在线观看视频网站免费| 18禁黄网站禁片免费观看直播| 午夜a级毛片| 国产精品国产三级国产av玫瑰| 深夜精品福利| 一级毛片电影观看 | 深爱激情五月婷婷| 欧美日韩乱码在线| 男女之事视频高清在线观看| 日韩欧美免费精品| 综合色丁香网| 欧美xxxx黑人xx丫x性爽| 18禁裸乳无遮挡免费网站照片| 久久久久九九精品影院| 国产亚洲av嫩草精品影院| 一进一出抽搐动态| 又爽又黄无遮挡网站| 女的被弄到高潮叫床怎么办| 极品教师在线视频| 级片在线观看| 看片在线看免费视频| 欧美绝顶高潮抽搐喷水| 三级男女做爰猛烈吃奶摸视频| 欧美三级亚洲精品| 男人舔奶头视频| 亚洲精品成人久久久久久| 久久久久久伊人网av| 国产精品三级大全| 国产精品女同一区二区软件| 日本a在线网址| 男人的好看免费观看在线视频| 秋霞在线观看毛片| 久久精品国产99精品国产亚洲性色| 乱码一卡2卡4卡精品| 日韩国内少妇激情av| 亚洲精华国产精华液的使用体验 | 麻豆乱淫一区二区| 在线看三级毛片| 成熟少妇高潮喷水视频| 99热这里只有是精品50| 欧美潮喷喷水| 国产成人一区二区在线| 成人亚洲欧美一区二区av| 91午夜精品亚洲一区二区三区| 18禁黄网站禁片免费观看直播| 菩萨蛮人人尽说江南好唐韦庄 | 我要看日韩黄色一级片| 我要搜黄色片| 俺也久久电影网| 女人被狂操c到高潮| 天堂影院成人在线观看| 国产av麻豆久久久久久久| 日本熟妇午夜| 女同久久另类99精品国产91| 22中文网久久字幕| 黄片wwwwww| АⅤ资源中文在线天堂| 男女之事视频高清在线观看| 老司机影院成人| 国产视频内射| 观看免费一级毛片| 亚洲无线观看免费| .国产精品久久| 亚洲电影在线观看av| 日本黄大片高清| 免费看a级黄色片| 黄片wwwwww| 国产精品免费一区二区三区在线| 一个人观看的视频www高清免费观看| 国产中年淑女户外野战色| 成人亚洲精品av一区二区| 国产精品一区二区免费欧美| 亚洲av中文字字幕乱码综合| 国产精品99久久久久久久久| 亚洲人成网站在线播放欧美日韩| 免费在线观看影片大全网站| 亚洲成人精品中文字幕电影| 丰满人妻一区二区三区视频av| 国产一区二区三区在线臀色熟女| 天天躁日日操中文字幕| 精品一区二区三区视频在线| 99热这里只有精品一区| 国产欧美日韩精品一区二区| 亚洲国产精品合色在线| 老司机福利观看| 精品国内亚洲2022精品成人| 亚洲av五月六月丁香网| 日本成人三级电影网站| 亚洲精品一区av在线观看| 色噜噜av男人的天堂激情| 国产一区二区激情短视频| 两个人视频免费观看高清| 在现免费观看毛片| 蜜桃亚洲精品一区二区三区| 成人美女网站在线观看视频| 久久久久久久久久成人| 亚洲欧美精品自产自拍| 国产真实伦视频高清在线观看| 亚洲av免费在线观看| 91麻豆精品激情在线观看国产| 色综合亚洲欧美另类图片| 久久精品久久久久久噜噜老黄 | 国产国拍精品亚洲av在线观看| 免费看美女性在线毛片视频| 久久九九热精品免费| 亚洲av不卡在线观看| 男女那种视频在线观看| 天天一区二区日本电影三级| 国产精品一区二区性色av| 男女啪啪激烈高潮av片| 久久精品夜夜夜夜夜久久蜜豆| a级毛片免费高清观看在线播放| 国产精品久久久久久av不卡| 午夜激情福利司机影院| 亚洲天堂国产精品一区在线| 国产精品99久久久久久久久| 91久久精品电影网| 99久久中文字幕三级久久日本| 日韩强制内射视频| 久久国内精品自在自线图片| 国产精品伦人一区二区| 啦啦啦观看免费观看视频高清| 日本熟妇午夜| 91久久精品电影网| 精品一区二区免费观看| 有码 亚洲区| 99久久无色码亚洲精品果冻| 午夜激情福利司机影院| 日本在线视频免费播放| 日韩成人伦理影院| 精品一区二区三区视频在线观看免费| 2021天堂中文幕一二区在线观| 亚洲精品色激情综合| 少妇猛男粗大的猛烈进出视频 | 99在线人妻在线中文字幕| 综合色丁香网| 三级毛片av免费| 老师上课跳d突然被开到最大视频| 91在线观看av| 亚洲一区二区三区色噜噜| 亚洲成人中文字幕在线播放| 有码 亚洲区| 午夜老司机福利剧场| 欧美3d第一页| 成人一区二区视频在线观看| 亚洲中文字幕日韩| 美女免费视频网站| 在线观看av片永久免费下载| 午夜视频国产福利| 亚洲av熟女| 国产亚洲精品av在线| 中文字幕精品亚洲无线码一区| 少妇人妻精品综合一区二区 | 丰满的人妻完整版| 两个人视频免费观看高清| 夜夜夜夜夜久久久久| 小蜜桃在线观看免费完整版高清| 国产午夜福利久久久久久| 色吧在线观看| 夜夜爽天天搞| 午夜精品国产一区二区电影 | 国产精品伦人一区二区| 国产aⅴ精品一区二区三区波| 老熟妇乱子伦视频在线观看| 无遮挡黄片免费观看| 狂野欧美激情性xxxx在线观看| 国产真实乱freesex| 久久久精品大字幕| 黄色日韩在线| 国产欧美日韩精品亚洲av| 天天躁日日操中文字幕| 老司机午夜福利在线观看视频| 国产午夜精品论理片| 国产精品久久久久久亚洲av鲁大| АⅤ资源中文在线天堂| 亚洲av免费高清在线观看| 久久久久精品国产欧美久久久| 伊人久久精品亚洲午夜| 亚洲成人精品中文字幕电影| 激情 狠狠 欧美| 99热全是精品| 看黄色毛片网站| 日韩精品青青久久久久久| 色哟哟·www| 国产亚洲精品久久久com| 亚洲美女搞黄在线观看 | 成人国产麻豆网| 噜噜噜噜噜久久久久久91| 久99久视频精品免费| 国产色婷婷99| 免费在线观看影片大全网站| 免费不卡的大黄色大毛片视频在线观看 | 欧美高清性xxxxhd video| 不卡视频在线观看欧美| 国产av在哪里看| 久久久久久久午夜电影| 久久久国产成人免费| 又黄又爽又刺激的免费视频.| 国产精品久久久久久久电影| 国产精品av视频在线免费观看| 久久6这里有精品| 韩国av在线不卡| 国产精品野战在线观看| 国产高清视频在线播放一区| 变态另类丝袜制服| 国产片特级美女逼逼视频| 国内精品久久久久精免费| 寂寞人妻少妇视频99o| 亚洲精品色激情综合| 久久午夜福利片| 日本在线视频免费播放| 国产亚洲精品综合一区在线观看| 亚洲色图av天堂| 热99在线观看视频| 久久九九热精品免费| 看黄色毛片网站| 亚洲成人中文字幕在线播放| 美女内射精品一级片tv| 国产精品一区二区免费欧美| 久久久久久九九精品二区国产| 欧美性感艳星| 成年av动漫网址| 精品熟女少妇av免费看| 淫秽高清视频在线观看| 国产v大片淫在线免费观看| 三级经典国产精品| 欧美精品国产亚洲| 夜夜夜夜夜久久久久| 色噜噜av男人的天堂激情| 亚洲精品粉嫩美女一区| 韩国av在线不卡| 少妇熟女aⅴ在线视频| 久久中文看片网| 简卡轻食公司| 免费看光身美女| 久久精品国产清高在天天线| 精品人妻熟女av久视频| 欧美绝顶高潮抽搐喷水| 日日撸夜夜添| 人妻久久中文字幕网| 精品午夜福利视频在线观看一区| 精品福利观看| 十八禁网站免费在线| 51国产日韩欧美| 欧美高清成人免费视频www| 久久久午夜欧美精品| 日日干狠狠操夜夜爽| 日韩欧美国产在线观看| 亚洲中文字幕一区二区三区有码在线看| 久久精品国产清高在天天线| 麻豆国产av国片精品| 免费一级毛片在线播放高清视频| 亚洲av中文av极速乱| 露出奶头的视频| 男人狂女人下面高潮的视频| 人人妻人人澡人人爽人人夜夜 | 日韩大尺度精品在线看网址| 色哟哟哟哟哟哟| 嫩草影院精品99| 老司机福利观看| 亚洲第一区二区三区不卡| 露出奶头的视频| 日本在线视频免费播放| 国产精品一区二区性色av| a级毛片a级免费在线| 精品久久久久久成人av| 精品一区二区三区视频在线观看免费| 亚洲av免费在线观看| 亚洲成人久久爱视频| 欧美人与善性xxx| 搡女人真爽免费视频火全软件 | 久久99热6这里只有精品| 国语自产精品视频在线第100页| 亚洲自偷自拍三级| 亚洲一区高清亚洲精品| h日本视频在线播放| 91久久精品国产一区二区三区| 国产av一区在线观看免费| 亚洲成人久久性| 观看美女的网站| 级片在线观看| 深夜a级毛片| 午夜精品在线福利| 中文资源天堂在线| 丰满的人妻完整版| 99热全是精品| 亚洲人成网站在线播放欧美日韩| 九九热线精品视视频播放| 免费黄网站久久成人精品| 午夜福利高清视频| а√天堂www在线а√下载| 真实男女啪啪啪动态图| 国产三级在线视频| 日本一二三区视频观看| 亚洲人成网站在线观看播放| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲一区高清亚洲精品| 欧美日韩一区二区视频在线观看视频在线 | 亚洲乱码一区二区免费版| 一区二区三区免费毛片| 国产精品人妻久久久久久| 九九久久精品国产亚洲av麻豆| 国产成人aa在线观看| 老熟妇乱子伦视频在线观看| 精品久久久久久久久av| 亚洲人成网站在线播| 日本-黄色视频高清免费观看| 久久久精品欧美日韩精品| 国产精品一区二区免费欧美| 国产 一区精品| 国产美女午夜福利| 在线观看66精品国产| 超碰av人人做人人爽久久| 免费观看精品视频网站| 男人舔女人下体高潮全视频| 久99久视频精品免费| .国产精品久久| av在线天堂中文字幕| 精品一区二区三区av网在线观看| 亚洲成人中文字幕在线播放| 可以在线观看的亚洲视频| 国产亚洲精品久久久久久毛片| 免费看美女性在线毛片视频| 日韩大尺度精品在线看网址| 色综合亚洲欧美另类图片| 成人鲁丝片一二三区免费| 国产爱豆传媒在线观看| 春色校园在线视频观看| 国内精品一区二区在线观看| 亚洲内射少妇av| 国产综合懂色| 欧美日韩精品成人综合77777| 能在线免费观看的黄片| 免费观看在线日韩| 夜夜爽天天搞| 搞女人的毛片| 欧美一区二区国产精品久久精品| 久久久a久久爽久久v久久| 男女下面进入的视频免费午夜| a级毛色黄片| 成人毛片a级毛片在线播放| 男人的好看免费观看在线视频| 免费无遮挡裸体视频| 观看美女的网站| 国产亚洲精品综合一区在线观看| АⅤ资源中文在线天堂| 日韩欧美国产在线观看| 午夜a级毛片| 国产单亲对白刺激| av天堂在线播放| 可以在线观看的亚洲视频| 欧美三级亚洲精品| 97热精品久久久久久| 99国产精品一区二区蜜桃av| 欧美潮喷喷水| 日本精品一区二区三区蜜桃| 久久久久久久久久久丰满| 3wmmmm亚洲av在线观看| 国产视频内射| 久久欧美精品欧美久久欧美| 久久久国产成人免费| 国产老妇女一区| 久久韩国三级中文字幕| 亚洲欧美日韩东京热| 久久久精品欧美日韩精品| 天堂av国产一区二区熟女人妻| 黑人高潮一二区| 国产 一区 欧美 日韩| 成人国产麻豆网| 国产精品久久电影中文字幕| 成人三级黄色视频| 又爽又黄无遮挡网站| 国产三级在线视频| 日本撒尿小便嘘嘘汇集6| 久久精品人妻少妇| 别揉我奶头~嗯~啊~动态视频| 18禁在线播放成人免费| 97超碰精品成人国产| 久久久精品94久久精品| 看片在线看免费视频| 久久精品综合一区二区三区| 欧美潮喷喷水| 精品久久国产蜜桃| 国产欧美日韩精品一区二区| 长腿黑丝高跟| 可以在线观看毛片的网站| АⅤ资源中文在线天堂| 精品久久久噜噜| 97超视频在线观看视频| 观看免费一级毛片| 精品乱码久久久久久99久播| 色尼玛亚洲综合影院| 男人舔女人下体高潮全视频| 欧美成人精品欧美一级黄| 97人妻精品一区二区三区麻豆| 夜夜夜夜夜久久久久| 无遮挡黄片免费观看| 免费看av在线观看网站| 久久人人精品亚洲av| 性欧美人与动物交配| 精华霜和精华液先用哪个| 国产精品伦人一区二区|