Evidence of non-stationary relationships between climate and forest responses: Increased sensitivity to climate change in Iberian forests

Long-term annual estimation of forest above ground biomass, canopy cover, and height from airborne and spaceborne sensors synergies in the Iberian Peninsula

The Mediterranean Basin has experienced substantial land use changes as traditional agriculture decreased and population migrated from rural to urban areas, which have resulted in a large forest cover increase. The combination of Landsat time series, providing spectral information, with lidar, offering three-dimensional insights, has emerged as a viable option for the large-scale cartography of forest structural attributes across large time spans. Here we develop and test a comprehensive framework to map forest above ground biomass, canopy cover and forest height in two regions spanning the most representative biomes in the peninsular Spain, Mediterranean (Madrid region) and temperate (Basque Country). As reference, we used lidar-based direct estimates of stand height and forest canopy cover. The reference biomass and volume were predicted from lidar metrics. Landsat time series predictors included annual temporal profiles of band reflectance and vegetation indices for the 1985-2023 period. Additional predictor variables including synthetic aperture radar, disturbance history, topography and forest type were also evaluated to optimize forest structural attributes retrieval. The estimates were independently validated at two temporal scales, i) the year of model calibration and ii) the year of the second lidar survey. The final models used as predictor variables only Landsat based metrics and topographic information, as the available SAR time-series were relatively short (1991-2011) and disturbance information did not decrease the estimation error. Model accuracies were higher in the Mediterranean forests when compared to the temperate forests (R2 = 0.6-0.8 vs. 0.4-0.5). Between the first (1985-1989) and the last (2020-2023) decades of the monitoring period the average forest cover increased from 21 ± 2% to 32 ± 1%, mean height increased from 6.6 ± 0.43 m to 7.9 ± 0.18 m and the mean biomass from 31.9 ± 3.6 t ha-1 to 50.4 ± 1 t ha-1 for the Mediterranean forests. In temperate forests, the average canopy cover increased from 55 ± 4% to 59 ± 3%, mean height increased from 15.8 ± 0.77 m to 17.3 ± 0.21m, while the growing stock volume increased from 137.8 ± 8.2 to 151.5 ± 3.8 m3 ha-1. Our results suggest that multispectral data can be successfully linked with lidar to provide continuous information on forest height, cover, and biomass trends.

Forest disturbance regimes and trends in continental Spain (1985-2023) using dense landsat time series

Forest disturbance regimes across biomes are being altered by interactive effects of global change. Establishing baselines for assessing change requires detailed quantitative data on past disturbance events, but such data are scarce and difficult to obtain over large spatial and temporal scales. The integration of remote sensing with dense time series analysis and cloud computing platforms is enhancing the ability to monitor historical disturbances, and especially non-stand replacing events along climatic gradients. Since the integration of such tools is still scarce in Mediterranean regions, here, we combine dense Landsat time series and the Continuous Change Detection and Classification - Spectral Mixture Analysis (CCDC-SMA) method to monitor forest disturbance in continental Spain from 1985 to 2023. We adapted the CCDC-SMA method for improved disturbance detection creating new spectral libraries representative of the study region, and quantified the year, month, severity, return interval, and type of disturbance (stand replacing, non-stand replacing) at a 30 m resolution. In addition, we characterised forest disturbance regimes and trends (patch size and severity, and frequency of events) of events larger than 0.5 ha at the national scale by biome (Mediterranean and temperate) and forest type (broadleaf, needleleaf and mixed). We quantified more than 2.9 million patches of disturbed forest, covering 4.6 Mha over the region and period studied. Forest disturbances were on average larger but less severe in the Mediterranean than in the temperate biome, and significantly larger and more severe in needleleaf than in mixed and broadleaf forests. Since the late 1980s, forest disturbances have decreased in size and severity while increasing in frequency across all biomes and forest types. These results have important implications as they confirm that disturbance regimes in continental Spain are changing and should therefore be considered in forest strategic planning for policy development and implementation.

Relative decline in density of Northern Hemisphere tree species in warm and arid regions of their climate niches

Although climate change is expected to drive tree species toward colder and wetter regions of their distribution, broadscale empirical evidence is lacking. One possibility is that past and present human activities in forests obscure or alter the effects of climate. Here, using data from more than two million monitored trees from 73 widely distributed species, we quantify changes in tree species density within their climatic niches across Northern Hemisphere forests. We observe a reduction in mean density across species, coupled with a tendency toward increasing tree size. However, the direction and magnitude of changes in density exhibit considerable variability between species, influenced by stand development that results from previous stand-level disturbances. Remarkably, when accounting for stand development, our findings show a significant change in density toward cold and wet climatic conditions for 43% of the species, compared to only 14% of species significantly changing their density toward warm and arid conditions in both early- and late-development stands. The observed changes in climate-driven density showed no clear association with species traits related to drought tolerance, recruitment and dispersal capacity, or resource use, nor with the temperature or aridity affiliation of the species, leaving the underlying mechanism uncertain. Forest conservation policies and associated management strategies might want to consider anticipated long-term species range shifts alongside the integration of contemporary within-distribution density changes.

Identifying drivers of non-stationary climate-growth relationships of European beech

The future performance of the widely abundant European beech (Fagus sylvatica L.) across its ecological amplitude is uncertain. Although beech is considered drought-sensitive and thus negatively affected by drought events, scientific evidence indicating increasing drought vulnerability under climate change on a cross-regional scale remains elusive. While evaluating changes in climate sensitivity of secondary growth offers a promising avenue, studies from productive, closed-canopy forests suffer from knowledge gaps, especially regarding the natural variability of climate sensitivity and how it relates to radial growth as an indicator of tree vitality. Since beech is sensitive to drought, we in this study use a drought index as a climate variable to account for the combined effects of temperature and water availability and explore how the drought sensitivity of secondary growth varies temporally in dependence on growth variability, growth trends, and climatic water availability across the species' ecological amplitude. Our results show that drought sensitivity is highly variable and non-stationary, though consistently higher at dry sites compared to moist sites. Increasing drought sensitivity can largely be explained by increasing climatic aridity, especially as it is exacerbated by climate change and trees' rank progression within forest communities, as (co-)dominant trees are more sensitive to extra-canopy climatic conditions than trees embedded in understories. However, during the driest periods of the 20th century, growth showed clear signs of being decoupled from climate. This may indicate fundamental changes in system behavior and be early-warning signals of decreasing drought tolerance. The multiple significant interaction terms in our model elucidate the complexity of European beech's drought sensitivity, which needs to be taken into consideration when assessing this species' response to climate change.

Structural diversity is better associated with forest productivity than species or functional diversity

Understanding the relationship between biodiversity and productivity can be advanced by improving metrics used to quantify biodiversity. Structural diversity, that is, variation of size and form of plant organs, is an emerging biodiversity metric. However, compared with the other biodiversity metrics, its relative importance in specific components of forest productivity, for example, recruitment of new individuals, biomass net change after accounting for mortality, is largely unknown, particularly across a large spatial scale with multiple influential gradients. To address the knowledge gap, we used USDA Forest Service Forest Inventory and Analysis (FIA) data across the southcentral USA from 2008 to 2017. We calculated forest biomass increments due to recruitment and growth and net change in biomass. Then, we quantified the effects of a range of abiotic and biotic variables on the biomass increments and net change. Our results showed that (1) Structural diversity was negatively associated with the two biomass increments and net change in biomass. The negative effects were supported by increased occurrences of insects and diseases with greater structural diversity. (2) Compared with species and functional diversity, structural diversity showed a better association with biomass increments and net change, suggested by its larger absolute values of standardized coefficients, and the effects of structural diversity were negative in contrast to species diversity. (3) The effects of structural diversity, stand age, and elevation differed between natural and planted forests that may stem from the differences in stand development and species composition between the two forest types. Together, structural diversity may represent an important dimension of biodiversity impacts on plant productivity, which could be related to the exacerbated disturbances with greater structural diversity.

A species' response to spatial climatic variation does not predict its response to climate change

The dominant paradigm for assessing ecological responses to climate change assumes that future states of individuals and populations can be predicted by current, species-wide performance variation across spatial climatic gradients. However, if the fates of ecological systems are better predicted by past responses to in situ climatic variation through time, this current analytical paradigm may be severely misleading. Empirically testing whether spatial or temporal climate responses better predict how species respond to climate change has been elusive, largely due to restrictive data requirements. Here, we leverage a newly collected network of ponderosa pine tree-ring time series to test whether statistically inferred responses to spatial versus temporal climatic variation better predict how trees have responded to recent climate change. When compared to observed tree growth responses to climate change since 1980, predictions derived from spatial climatic variation were wrong in both magnitude and direction. This was not the case for predictions derived from climatic variation through time, which were able to replicate observed responses well. Future climate scenarios through the end of the 21st century exacerbated these disparities. These results suggest that the currently dominant paradigm of forecasting the ecological impacts of climate change based on spatial climatic variation may be severely misleading over decadal to centennial timescales.

More than trees: Stand management can be used to improve ecosystem diversity, structure and functioning 20 years after forest restoration in drylands

In Mediterranean drylands, extensive areas have been restored by reforestation over the past decades to improve diversity, soil fertility, and tree natural regeneration, contributing to halting desertification and land degradation. However, evaluating reforestation success usually relies on tree survival, while holistic and long-term evaluations of reforestation success based on ecosystem diversity, structure and functioning are scarce. In this work, we provide the first assessment that combines the evaluation of planted trees and indicators of ecosystem diversity, structure, and functioning in established reforestations with three native Mediterranean species along a climatic gradient. We sampled 43 20-year-old stands with umbrella pine, holm oak and cork oak in Portugal, and tested the effects of tree species composition, stand management (i.e., differences in tree density and shrub cover), and edaphoclimatic conditions, on the size of planted trees, species diversity, structural complexity and indicators of ecosystem functioning related to productivity, soil nutrients and tree natural regeneration. Our results show that, after 20 years of reforestation, stand management was an essential driver of plant diversity and ecosystem functioning. Higher tree density, particularly of oaks, and higher shrub cover improved plant diversity, ecosystem productivity, and oak regeneration. The latter was also improved by structural complexity. Tree composition effects highlighted the importance of pine management to avoid competition. Since we evaluated these reforestations along a climatic gradient, we also conclude that climate influenced pine and holm oak size, ecosystem productivity, and soil C/N. Our research, by being based on assessing the long-term reforestation success in a more holistic way, highlighted the importance of stand management for improving ecosystem diversity and functioning in these restored systems. Practices such as increasing tree density up to ~800 trees/ha and allowing a shrub cover of ca. 30 %, may improve the ecological condition of future and currently reforested areas across the Mediterranean region.

Drought conditions, aridity and forest structure control the responses of Iberian holm oak woodlands to extreme droughts: A large-scale remote-sensing exploration in eastern Spain

Understanding how Mediterranean forests respond to the increasing frequency of extreme droughts and forest densification is crucial for effective land management in the present context of climate change and land abandonment. We study the responses of Iberian holm oak (Quercus ilex L.) woodlands to recent extreme droughts during 2000-2019 along broad gradients of climate aridity and forest structure. To this purpose, we apply large-scale remote-sensing using MODIS EVI as a primary production proxy in 5274 Q. ilex sites distributed within a 100,000 km2 region in eastern Spain. These woodlands were extensively affected by two extreme drought events in 2005 and 2012. Resistance, assessed as the capacity of the ecosystems to maintain primary production during drought, was significantly lower for semi-arid than for sub-humid and dry-transition conditions. Holm oak woodlands located in semi-arid areas of the region showed also poorer resilience to drought, characterized by low capacity to fully recover to their pre-drought production levels. Further, drought intensity and both pre- and post-drought hydric conditions controlled the variations of resistance, recovery and resilience between the two analyzed extreme drought events. Drought effects were particularly negative for dense Q. ilex stands under semi-arid climate conditions, where strong competition for scarce water resources reduced drought resistance. The observed drought vulnerability of semi-arid holm oak woodlands may affect the long-term stability of these dry forests. Adaptive management strategies, such as selective forest thinning, may be useful for improving drought responses in these more vulnerable semi-arid woodlands. Conversely, natural rewilding may more appropriately guide management actions for more humid areas, where densely developed Q. ilex woodlands show in general a high ability to maintain ecosystem primary production during drought.

Immediate and carry-over effects of late-spring frost and growing season drought on forest gross primary productivity capacity in the Northern Hemisphere

Forests are increasingly exposed to extreme global warming-induced climatic events. However, the immediate and carry-over effects of extreme events on forests are still poorly understood. Gross primary productivity (GPP) capacity is regarded as a good proxy of the ecosystem's functional stability, reflecting its physiological response to its surroundings. Using eddy covariance data from 34 forest sites in the Northern Hemisphere, we analyzed the immediate and carry-over effects of late-spring frost (LSF) and growing season drought on needle-leaf and broadleaf forests. Path analysis was applied to reveal the plausible reasons behind the varied responses of forests to extreme events. The results show that LSF had clear immediate effects on the GPP capacity of both needle-leaf and broadleaf forests. However, GPP capacity in needle-leaf forests was more sensitive to drought than in broadleaf forests. There was no interaction between LSF and drought in either needle-leaf or broadleaf forests. Drought effects were still visible when LSF and drought coexisted in needle-leaf forests. Path analysis further showed that the response of GPP capacity to drought differed between needle-leaf and broadleaf forests, mainly due to the difference in the sensitivity of canopy conductance. Moreover, LSF had a more severe and long-lasting carry-over effect on forests than drought. These results enrich our understanding of the mechanisms of forest response to extreme events across forest types.

Applying space-for-time substitution to infer the growth response to climate may lead to overestimation of tree maladaptation: Evidence from the North American White Spruce Network

Under climate change circumstances, increasing studies have reported the temporal instability of tree growth responses to climate, which poses a major challenge to linearly extrapolating past climate and future growth dynamics using tree-ring data. Space-for-time substitution (SFTS) is a potential solution to this problem that is widely used in the dendrochronology field to project past or future temporal growth response trajectories from contemporary spatial patterns. However, the projected accuracy of the SFTS in the climate effects on tree growth remains uncertain. Here, we empirically test the SFTS method by comparing the effect of spatial and temporal climate variations on climate responses of white spruce (Picea glauca), which has a transcontinental range in North America. We first applied a response surface regression model to capture the variations in growth responses along the spatial climate gradients. The results showed that the relationships between growth and June temperature varied along spatial climate gradients in a predictable way. And their relationships varied mainly along with local temperate condition. Then, the projected correlation coefficients between growth and climate using SFTS were compared against the observed. We found that the growth response changes caused by spatial versus temporal climate variations showed opposite trends. Moreover, the projected correlation coefficients using the SFTS were significantly lower than the observed. This finding suggests that applying the SFTS to project the growth response of white spruce might lead to an overestimation of the degree of tree maladaptation in future climate scenarios. And the overestimation is likely to get weaker from Alaska and Yukon Territory in the west to Quebec in the east. Although this is only a case study of the SFTS method for projecting tree growth response, our findings suggest that direct application of the SFTS method may not be applicable to all regions and all tree species.

Re-assessment of the climatic controls on the carbon and water fluxes of a boreal aspen forest over 1996-2016: Changing sensitivity to long-term climatic conditions

Recent evidence suggests that the relationships between climate and boreal tree growth are generally non-stationary; however, it remains uncertain whether the relationships between climate and carbon (C) fluxes of boreal forests are stationary or have changed over recent decades. In this study, we used continuous eddy-covariance and microclimate data over 21 years (1996-2016) from a 100-year-old trembling aspen stand in central Saskatchewan, Canada to assess the relationships between climate and ecosystem C and water fluxes. Over the study period, the most striking climatic event was a severe, 3-year drought (2001-2003). Gross ecosystem production (GEP) showed larger interannual variability than ecosystem respiration (R_e ) over 1996-2016, but R_e was the dominant component contributing to the interannual variation in net ecosystem production (NEP) during post-drought years. The interannual variations in evapotranspiration (ET) and C fluxes were primarily driven by temperature and secondarily by water availability. Two-factor linear models combining precipitation and temperature performed well in explaining the interannual variation in C and water fluxes (R²  > .5). The temperature sensitivities of all three C fluxes (NEP, GEP and R_e ) declined over the study period (p < .05), and, as a result, the phenological controls on annual NEP weakened. The decreasing temperature sensitivity of the C fluxes may reflect changes in forest structure, related to the over-maturity of the aspen stand at 100 years of age, and exacerbated by high tree mortality following the severe 2001-2003 drought. These results may provide an early warning signal of driver shift or even an abrupt status shift of aspen forest dynamics. They may also imply a universal weakening in the relationship between temperature and GEP as forests become over-mature, associated with the structural and compositional changes that accompany forest ageing.

Tree growth response to drought partially explains regional-scale growth and mortality patterns in Iberian forests

Tree-ring data has been widely used to inform about tree growth responses to drought at the individual scale, but less is known about how tree growth sensitivity to drought scales up driving changes in forest dynamics. Here, we related tree-ring growth chronologies and stand-level forest changes in basal area from two independent data sets to test if tree-ring responses to drought match stand forest dynamics (stand basal area growth, ingrowth, and mortality). We assessed if tree growth and changes in forest basal area covary as a function of spatial scale and tree taxa (gymnosperm or angiosperm). To this end, we compared a tree-ring network with stand data from the Spanish National Forest Inventory. We focused on the cumulative impact of drought on tree growth and demography in the period 1981-2005. Drought years were identified by the Standardized Precipitation Evapotranspiration Index, and their impacts on tree growth by quantifying tree-ring width reductions. We hypothesized that forests with greater drought impacts on tree growth will also show reduced stand basal area growth and ingrowth and enhanced mortality. This is expected to occur in forests dominated by gymnosperms on drought-prone regions. Cumulative growth reductions during dry years were higher in forests dominated by gymnosperms and presented a greater magnitude and spatial autocorrelation than for angiosperms. Cumulative drought-induced tree growth reductions and changes in forest basal area were related, but initial stand density and basal area were the main factors driving changes in basal area. In drought-prone gymnosperm forests, we observed that sites with greater growth reductions had lower stand basal area growth and greater mortality. Consequently, stand basal area, forest growth, and ingrowth in regions with large drought impacts was significantly lower than in regions less impacted by drought. Tree growth sensitivity to drought can be used as a predictor of gymnosperm demographic rates in terms of stand basal area growth and ingrowth at regional scales, but further studies may try to disentangle how initial stand density modulates such relationships. Drought-induced growth reductions and their cumulative impacts have strong potential to be used as early-warning indicators of regional forest vulnerability.

Divergent Tree Growth and the Response to Climate Warming and Humidification in the Tianshan Mountains, China

In recent decades, the global climate has changed significantly. The climate in Northwest China became warm-wet, especially in the Tianshan Mountains. In order to explore the response of tree growth to recent climate change, the two dominant trees species, Picea schrenkiana Fisch. et Mey. and Larix sibirica Ledeb., were studied with the dendrochronological method in the western Tianshan Mountains (WT) and the eastern Tianshan Mountains (ET). Our results showed that: (1) The tree growth of four sample sites in the WT significantly increased in recent decades, while the trees in the three sample sites in the ET significantly decreased. (2) In the WT, except for the Manas site, the tree-ring chronologies of the other three sites were significantly positively correlated with the mean annual minimum temperature. Tree-ring chronologies in the WT, except for Bangfanggou site, were significantly positively correlated with annual precipitation. In the ET, only the tree chronology of L. sibirica in the Balikun site was significantly negatively correlated with the annual temperatures, including the mean minimum, mean and mean maximum temperature. (3) The proportion of trees with a significant upward growth trend at each site decreased from west to east, and the proportion of trees with a significant downward growth trend at each site increased from west to east along the whole Tianshan Mountains. (4) The correlation of tree-ring chronologies with the annual temperature and annual precipitation was not stable during the study period. Warm-humidification promoted the growth of trees in the WT but inhibited tree growth in the ET, which may be exacerbated drought stress in the ET where the increase in precipitation was not enough to offset the increased evapotranspiration potential caused by warming.

Identifying Forest Structural Types along an Aridity Gradient in Peninsular Spain: Integrating Low-Density LiDAR, Forest Inventory, and Aridity Index

Forest structure is a key driver of forest functional processes. The characterization of forest structure across spatiotemporal scales is essential for forest monitoring and management. LiDAR data have proven particularly useful for cost-effectively estimating forest structural attributes. This paper evaluates the ability of combined forest inventory data and low-density discrete return airborne LiDAR data to discriminate main forest structural types in the Mediterranean-temperate transition ecotone. Firstly, we used six structural variables from the Spanish National Forest Inventory (SNFI) and an aridity index in a k-medoids algorithm to define the forest structural types. These variables were calculated for 2770 SNFI plots. We identified the main species for each structural type using the SNFI. Secondly, we developed a Random Forest model to predict the spatial distribution of structural types and create wall-to-wall maps from LiDAR data. The k-medoids clustering algorithm enabled the identification of four clusters of forest structures. A total of six out of forty-one potential LiDAR metrics were utilized in our Random Forest, after evaluating their importance in the Random Forest model. Selected metrics were, in decreasing order of importance, the percentage of all returns above 2 m, mean height of the canopy profile, the difference between the 90th and 50th height percentiles, the area under the canopy curve, and the 5th and the 95th percentile of the return heights. The model yielded an overall accuracy of 64.18%. The producer’s accuracy ranged between 36.11% and 88.93%. Our results confirm the potential of this approximation for the continuous monitoring of forest structures, which is key to guiding forest management in this region.

SilvAdapt.Net: A Site-Based Network of Adaptive Forest Management Related to Climate Change in Spain

Adaptive forest management (AFM) is an urgent need because of the uncertainty regarding how changes in the climate will affect the structure, composition and function of forests during the next decades. Current research initiatives for the long-term monitoring of impacts of silviculture are scattered and not integrated into research networks, with the consequent losses of opportunities and capacity for action. To increase the scientific and practical impacts of these experiences, it is necessary to establish logical frameworks that harmonize the information and help us to define the most appropriate treatments. In this context, a number of research groups in Spain have produced research achievements and know-how during the last decades that can allow for the improvement in AFM. These groups address the issue of AFM from different fields, such as ecophysiology, ecohydrology and forest ecology, thus resulting in valuable but dispersed expertise. The main objective of this work is to introduce a comprehensive strategy aimed to study the implementation of AFM in Spain. As a first step, a network of 34 experimental sites managed by 14 different research groups is proposed and justified. As a second step, the most important AFM impacts on Mediterranean pines, as one of the most extended natural and planted forest types in Spain, are presented. Finally, open questions dealing with key aspects when attempting to implement an AFM framework are discussed. This study is expected to contribute to better outlining the procedures and steps needed to implement regional frameworks for AFM.

High forest stand density exacerbates growth decline of conifers driven by warming but not broad-leaved trees in temperate mixed forest in northeast Asia

Increasing temperature over recent decades is expected to positively impact tree growth in humid regions. However, high stand density could increase the negative effects of warming-induced drought through inter-tree competition. How neighborhood competition impacts tree growth responding to climate change remains unclear. Here, we utilized the Changbai Mountain region in northeastern Asia as our study area. We quantified individual tree growth using tree-ring samples collected from three dominant tree species growing in three forest stand density levels. We estimated the effects of climate warming and forest stand density on growth processes and tested for a species-specific response to climate. Our results demonstrated that overall 25% of Korean pine, but only ~3% of Mongolian oak and ~ 4% of Manchurian ash experienced growth reduction. Increased forest density can also exacerbate growth reduction. We identified a climate turning point in 1984, where warming rapidly increased, and defined two groups, "enhance group" (EG) and "decline group" (DG), according to the individual tree growth trend after 1984. For the EG, climate warming increased temperature sensitivity, but the temperature sensitivity declined with increasing stand density for the whole study period. For the DG, tree growth sensitivity shifted from temperature to precipitation after 1984, driven by increased competition pressure under climate warming. Our study concludes that growth decline from warming-induced drought might be amplified by high forest stand density, was especially pronounced in conifer trees.

Linking functional traits with tree growth and forest productivity in Quercus ilex forests along a climatic gradient

Plant functional traits are highly plastic to changes in climatic factors and nutrient availability. However, the intraspecific plant response to abiotic factors and the overall effect on tree growth and productivity is still under debate. We studied forest productivity for 30 Quercus ilex subsp. ballota forests in Spain along a broad climatic gradient of aridity (mean annual precipitation from 321 to 858 mm). We used linear mixed models to quantify the effect of climatic and edaphic (soil nutrients, topography, and texture) factors on tree functional traits (leaf and branch traits), and subsequently, the effect of such functional traits and abiotic factors on the relative growth rate (RGR) of adult trees. We used piecewise structural equation models (SEMs) to determine the causal effect of intrinsic and extrinsic factors on forest productivity. Our results showed that tree functional traits were mainly explained by climatic and edaphic factors. Functional traits and tree biomass explained forest biomass and RGR, respectively, which ultimately explained forest productivity. In conclusion, intraspecific variability of functional traits has a significant effect on plant biomass and growth, which ultimately may explain forest productivity in Quercus ilex forests.

Long-Term Tree Regeneration of Fragmented Agroforestry Systems Under Varying Climatic Conditions

Iberian dehesas and montados are agroforestry systems protected by the European Habitats Directive due to high levels of biological diversity associated to their savannah-like structure. Tree scattering in dehesas, montados and other agroforestry systems is, however, known to compromise tree regeneration, although recent work suggests that it may protect tree populations from climate warming by alleviating plant-plant competition. We analyze how climatic conditions, tree isolation and their interactions influence the outcomes of regeneration stages, from flower production to early seedling establishment, using data gathered during the long-term monitoring (2001–2018) of ca. 300 Holm oak Quercus ilex trees located in central Spain. Holm oak reproductive effort, predispersal seed losses, and early seedling recruitment were sensitive to climate change, especially to year-round drought. Effort and early seedling recruitment decreased, while abortion and predispersal seed predation increased, with higher drought intensity. Spring warming increases pollination effectiveness, but had no further effect on acorn crops. Forest clearing seemed to have little scope to ameliorate these negative effects, as shown by weak or no interactive effects between the spatial configuration of trees (cover or isolation) and climate variables (spring temperature or drought intensity). Forest opening aimed at decreasing adult tree mortality under climate change scenarios would then have little or no effects on tree recruitment. Landscape-scale rotations alternating shrub encroachment and thinning along periods adapted to changing climate are proposed as the main management option to preserve both oak forests and dehesas in the long term.

Impacts of Climate Change on Hydroclimatic Conditions of U.S. National Forests and Grasslands

The conterminous United States includes national forests and grasslands that provide ecological, social, economic, recreational, and aesthetic services. Future climate change can alter long-term hydroclimatic conditions of national forests and grasslands and lead to negative consequences. This study characterizes shifts in hydroclimatology and basin characteristics of US National Forests (NFs) and National Grasslands (NGs) in response to climate change over the 21st century under the DRY, MIDDLE, and WET climate models with the representative concentration pathway (RCP) 8.5 emission scenario. Climatic projections for three climate models ranging from the driest to wettest conditions were obtained from the Multivariate Adaptive Constructed Analogs (MACA) dataset. Then, the variable infiltration capacity (VIC) model was used to model hydrological responses of the selected future climates. Changes in regional hydroclimatic conditions of NFs and NGs were assessed by the magnitude and direction of movements in the Budyko space. The Fu’s equation was applied to estimate changes in basin characteristics. The results indicate that NFs and NGs are likely to experience larger changes in basin characteristics compared to the average of the United States. In general, across the conterminous US, the NFs in mountainous regions are likely to have larger changes in hydroclimatic variables than NFs with lower elevation and NGs. Comparing Forest Service regions, Pacific Northwest, Intermountain, and Northern regions may have a less arid climate with lower freshwater availability. The Southwestern, Northern, Intermountain, and Rocky Mountain regions are likely to experience higher shifts in their basin characteristics. This study can help environmental scientists, and land and water managers improve future land management plans.