Skip to main content

Spatiotemporal variations in productivity and water use efficiency across a temperate forest landscape of Northeast China

Abstract

Background

Forests are increasingly valued for non-timber ecosystem services in place of conventional wood fiber production. Biomass carbon sequestration is one of the key ecosystem services forests are relied upon for mitigating global climate change. However, planning for large-scale tree planting and managing established forest stands for carbon sequestration require careful consideration of the gain in biomass production and tradeoff for other regulatory services. How a tradeoff between forest production and conservation of water resources is shaped by the condition of forest stand and environmental factors remains a question of broad interest in sustainable forest ecosystem management.

Methods

We studied the spatiotemporal patterns of net primary productivity (NPP), evapotranspiration (ET), and water use efficiency (WUE), and their relationships with local climatic and forest stand factors over a temperate forest landscape in Changbai Mountain, Northeast China. The time series of spatial data on NPP and ET were extracted from the global remote sensing datasets for the MOD16A3 and MOD17A3 products for the period 2000–2014. The time series of spatial patterns of annual precipitation and annual mean temperature were obtained as grid maps for regional meteorological variables. Stand patches were categorized into the types of conifers, broadleaves, and mixed-wood, as well as age-classes of young, mid-age, near mature, mature, and old-growth stands, and by establishment into natural and planted. Information on stands and selective site variables were compiled from the Forest Inventory Datasets of China. Analyses were performed with Arc-GIS.

Results

Over the study period of 2000–2014, the landscape-level annual NPP varied between 311.7 and 573.6 gC∙m− 2∙a− 1, ET between 559.9 and 603.0 mm∙a− 1, and WUE between 0.54 and 1.01 gC∙m− 2∙mm− 1. Across the forest landscape, the mean annual NPP varied between 205.0 and 639.4 gC∙m− 2∙a− 1, ET between 441.5 and 784.0 mm∙a− 1, and WUE in the range of 0.46–1.10 gC∙m− 2∙mm− 1. The spatial variations in NPP, ET, and WUE were commonly attributable to forest type, stand age class and density, establishment mode, and temperature variables, with some effects of other selective factors on ET and WUE. The three forest types were significantly (p <  0.05) differentiated in the mean annual NPP, ET, and WUE: the coniferous forests were highest in NPP (505.3 ± 1.4 gC∙m− 2∙a− 1; n = 1041) and WUE (0.872 ± 0.004 gC∙m− 2∙mm− 1; n = 1041), and lowest in ET (584.1 ± 1.6 mm∙a− 1; n = 1041), followed by the mixed-wood forests (NPP: 500.5 ± 0.8 gC∙m− 2∙a− 1; WUE: 0.856 ± 0.02 gC∙m− 2∙mm− 1; ET: 589.3 ± 0.9 mm∙a− 1; n = 2156); whereas the broadleaved forests were lowest in NPP (491.6 ± 0.6 gC∙m− 2∙a− 1; n = 4428) and WUE (0.832 ± 0.02 gC∙m− 2∙mm− 1; n = 4428), and highest in ET (594.7 ± 0.6 mm∙a− 1; n = 4428). The mean annual NPP, ET and WUE increased with stand age typically in coniferous forests, and weakly in mixed-wood forests. The natural stands had significantly (p <  0.001) greater NPP and WUE than planted stands in conifers and mixed-woods.

Conclusions

Forest type and stand-age class are factors of key consideration in planning and managing forest landscape for achieving the compromise between regional carbon balance and water usage. Facilitating the establishment of coniferous trees and conservation of mature forest stands, and near-natural management of planted forests, can be favorable options for delivering ecosystem services of carbon sequestration and water conservation in temperate regions.

Background

Worldwide concerns on climate change and biodiversity decline have promoted global efforts in natural forest protection and expansion of planted forests, whilst the existing forests are increasingly used for non-timber ecosystem services as alternatives to conventional wood fiber production, leading to a reverse trend of changes in forest cover area and wood stocks in recent decades compared to earlier-stages of declines due to wood harvesting and clearance for agricultural production. While increases in forest cover area and improved forest productivity unquestionably benefit global forest carbon sequestration and biodiversity conservation, there are issues of consideration on their likely impacts on regional hydrological cycles (Jackson et al. 2005; Liu et al. 2016). Large-scale forest recovery and tree planting may result in increases in evapotranspiration (ET), hence reduction in stream flows (Farley et al. 2005; Sun et al. 2006; Sun et al. 2008), and in regions of water shortage, give rise to risk factors causing environment degradation (Wang and Cao 2011; Gao et al. 2014). Sustainable forest management and planning for establishment of planted forests need to take into consideration of forest productivity-ecosystem water usage relationships to compromise the balance between gain in forest productivity and conservation of water resources.

Many factors can affect the forest productivity-ecosystem water usage relationships (Sun et al. 2006; Liu et al. 2015a; Feng et al. 2016). Among the dominant factors and for given climatic conditions, forest types, tree community structure, and stand development are of direct relevance to forestry and forest management. The effects of forest types and tree community structure are largely a resultant of genetic determination of water relations of forest tree species, whereas stand development may affect tree water uptake and evapotranspiration due to temporal changes in rooting depth and canopy cover. An explicit knowledge on how forest types, including the factors of tree species and community structure, and stand development would shape the tradeoff between forest productivity and ecosystem water loss is a pre-requisite for assisting regional forestry planning and forest management in the context of carbon sequestration and water conservation.

The temperate forests of Northeast China were historically a major timber production base, with majority of the primary forests cleared or giving way to secondary forests due to intensive timber harvesting and conversation to larch plantations during the twentieth century (Zhu et al. 2008; Yu et al. 2011). It is now one of the target regions for forest protection and restoration in China in a national effort of conserving and improving forest resources. Proactive tree plantings are widely implemented for reforestation of the region for the benefits of non-timber forest ecosystem services. However, lack of clear understanding on the implications of tree species selection and planting regimes to long-term forest community stability and consideration of the tradeoff between productivity and water usage may undermine the delivery of desirable ecosystem services from the practices deployed in establishment and management of planted forests. In particular, analysis of recent climate change suggests increasing drought in Northeast China (Liu et al. 2011; Peng et al. 2012). This, coupled with the increasing regional demand for agricultural water usage, would raise an issue of consideration on the water economy of the planted forests.

Water use efficiency (WUE) is conceptually a measure of water economy of plants. It defines the effectiveness of water conservation strategy of plants at leaf or individual level, and quantifies the productivity-water loss tradeoffs at the ecosystem level. Previous studies have demonstrated that WUE varies greatly in space and time due to changes in climate, vegetation types, ontogeny of plants, and plant community development (e.g. Yu et al. 2008; Medrano et al. 2009; Poni et al. 2009; Tian et al. 2010; Tong et al. 2014; Liu et al. 2015a; Zhang et al. 2016). A recent study reveals an observed phenomenon that natural forests are better at carbon sequestration and water conservation than planted forests in China (Yu et al. 2019). It is yet unclear on the mechanisms for the poor performance of planted forests in productivity and water conservation. It also remains a question of research interest that how WUE would vary with forest types and stand development under a given climate. In this study, over a landscape of cool temperate forests in Northeast China, we examined the spatial variations in net primary productivity (NPP) and evapotranspiration (ET) based on remote sensing datasets MOD16 and MOD17, and determined their relationships with environmental and forest stand factors by incorporating the regional meteorological data and forest inventory data. Our primary objectives were to determine how the productivity and water usage would vary with forest type and stand development, and to identify the key site and climatic factors in affecting WUE. We hypothesize that with forest stand development, there is a stage when the dual ecosystem services of carbon sequestration and water conservation could be favorably compromised, i.e. achieving a synergy between carbon economy and water economy, and that the productivity-water usage relationship vary with forest types.

Materials and methods

Study area

The study area is located in Lushuihe, Fusong county of Jilin Province (latitude 42°24′–42°49′ N, longitude 127°29′–128°02′ E), China, covering a total land area of approximately 1200 km2 and elevation range from 400 to 1300 m a.s.l. The forest landscape is representative of the cool temperate forests in the Changbai Mountain of Northeast China, with natural forests accounting for more than 80% of the forested area. The general region belongs to northern temperate zone under the influence of continental monsoon climate (Qi et al. 2016). Figure 1 shows inter-annual changes in annual mean air temperature and annual precipitation in the study area for the period 1986–2015, based on data from a local meteorological station (latitude 42°29′35“ N, longitude 127°46’01” E, altitude 753.6 m a.s.l.).

Fig. 1
figure 1

Inter-annual variations in annual mean air temperature and annual precipitation for the period 1986–2015 based on data from a local meteorological station (latitude 42°29′35“, latitude 127°46’01”E, altitude 753.6 m a.s.l.) within the study area in Lushuihe of Jilin Province, Northeast China

The regional forests subjected to partial progressive logging starting in late 1960s, mostly by ways of clear-cutting, and some of the clear-cut areas were converted to larch or spruce plantations. For concern of fast depleting local forest resources and protection of forest cover, selective logging was implemented beginning in 1985 (Fang et al. 2016), then followed by enhanced conservation efforts under the “Natural Forest Protection Program” orchestrated by the State Administration of Forestry of P.R. China at the turn of twenty-first century. The main forest types include coniferous, broadleaved, and mixed-wood forests. The dominant trees species consists of Korean pine (Pinus koraiensis), mongolian oak (Quercus mongolica), mono maple (Acer mono), amur linden (Tilia amurensis), elm (Ulmus propinqua), manchurian ash (Fraxinus mandshurica), aspen (Populus davidiana), birch (Betula platyphylla) and manchurian fir (Abies holophylla).

Data collections and processing

The NPP datasets at spatial resolution of 1-km were obtained from the annual MOD17A3 products (version 55) of the NASA Earth Observation System (EOS) program, which are produced by the Numerical Terradynamic Simulation Group (NTSG) of University of Montana (UMT) (available from http://www.ntsg.umt.edu/project/modis/mod17.php) (Xia et al. 2015). The MOD17A3 products were generated using the MOD17 algorithm from the original NASA/USGS LPDAAC (Land Processes Distributed Active Archive Center) 8-day and annual v5 products, with data gaps in the 8-day temporal MODIS FPAR/LAI caused by cloudiness filled with information from accompanying quality assessment flags (Running et al. 2004; Zhao et al. 2005; Zhao and Running 2010).

The ET datasets were the annually MOD16A3 product with a 1-km resolution (available from http://www.ntsg.umt.edu/project/modis/mod16.php), which were calculated by using the improved ET algorithm based on the Penman-Monteith (Mu et al. 2011). The improved algorithm employed involves simplifying the calculation of vegetation cover fraction, calculating ET as the sum daytime and nighttime components, adding soil heat flux calculation, improving estimates of stomatal conductance, aerodynamic resistance and boundary layer resistance, separating dry canopy surface from the wet, and dividing soil surface into saturated wet and moist surfaces (Mu et al. 2011; Tang et al. 2016). The ecosystem WUE was calculated from MODIS NPP and ET.

Both NPP and ET data were not validated in this study due to inadequate ground-based measurements for our study area. However, many regional and global studies of NPP and ET have shown the reliability of the MODIS NPP and ET products (Turner et al. 2006; Zhao and Running 2010; Lathuillière et al. 2012; Mu et al. 2013). The two products have been widely validated around the world against measurements from eddy covariance flux tower data including AmeriFlux, Fluxnet-Canada, EuroFlux and ChinaFlux across multiple biomes (Zhao et al. 2005; Heinsch et al. 2006; Mu et al. 2011; Liu et al. 2015b; Huang et al. 2018), or from the forest inventory data (Neumann et al. 2015; Gulbeyaz et al. 2018). The use of MODIS products has also been reported to be a reliable approach in studies of regional NPP in Chinese ecosystems (Jia et al. 2012; Liu et al. 2015a; Zhang et al. 2015).

The meteorological data involving precipitation, mean, minimum and maximum temperature for the period 2000–2014 employed in this study are derived from the gridded daily dataset (CN05.1) with a spatial resolution of 0.25° × 0.25° (Wu and Gao 2013), which is constructed following the “anomaly approach” during the interpolation but with over 2400 observing stations in China (Wu et al. 2017). In the “anomaly approach”, a gridded climatology is first calculated, and then a gridded daily anomaly is added to the climatology to obtain the final dataset. Gridded daily meteorological values were converted to annual precipitation, and annual mean, annual mean maximum and annual mean minimum air temperatures, using spatial overlay analysis with map algebra method, respectively. Then a bilinear interpolation method was adopted to resample the gridded meteorological data on spatial resolution of 1-km.

Information on the attributes of forest stands were derived from the Category II database of Forest Inventory of Jilin Province for the period ending 2013, covering a total of 7763 compartments of varying sizes. Each compartment represented a forest stand of single type and relatively uniform structure, and contained information on spatial attributes, tree species identity, stand age class, stand density and establishment mode (i.e. established through natural regeneration or tree planting). The compartmentalized forest stands were firstly grouped into three broad types as coniferous (CF), broadleaved (BL), and mixed-wood (MW), and were then further classified by the mode of establishment (i.e. natural vs. planted) and stand age classes (young, mid-age, near-mature, mature, and over-mature) (Table 1). The stand age classes were made with reference to rotation in accordance with the standard survey and inventory protocols of the State Forestry Administration of China (Sun et al. 2016). The spatial distributions of the six classified forest types and stand age classes are illustrated in Fig. 2.

Table 1 Information on classified forest stands across the forest landscape within the study area in Lushuihe of Jilin Province, Northeast China
Fig. 2
figure 2

Spatial distributions of forests by classifications of type and regeneration mode in the study area of Lushuihe of Jilin Province, Northeast China. a1, natural broadleaved; a2, planted broadleaved; b1, natural coniferous; b2; planted coniferous; c1, natural mixed-wood; c2, planted mixed-wood

Terrain factors including elevation, slope, aspect and slope position, among which the elevation data were obtained from the Global Map data archives (https://globalmaps.github.io/), while others were derived from the Category II database of Forest Inventory of Jilin Province. Soil organic carbon data were extracted from the China soil map based on Harmonized World Soil Database (v1.1), which was made available by WESTDC (Environmental and Ecological Science Data Center for West China, National Natural Science Foundation of China, http://westdc.westgis.ac.cn).

Statistical analysis

To obtain the average values of annual NPP, ET and WUE of the whole study area, we firstly extracted the annual values of NPP, ET and WUE in each pixel, then the mean values were calculated as the mean of all pixels in a given year. The spatial maps of NPP, ET and WUE values were produced using the average values for the period 2000–2014 for each individual grid cell, with Arc-GIS10.1. A simple linear regression method was used to analyze the inter-annual variability and the variable trend of each pixel during 2000–2014. The slopes of linear regression represent the temporal and spatial trends: a positive value indicates an uptrend, and a negative value indicates a downtrend. A greater absolute value of slope corresponds to a more rapid variation. Spatial trend analysis was performed using ENVI5.0.

We used the multiple linear regression model to examine the influence of three categories of variables: vegetation factors (forest types, mode of establishment, age class, and stand density), site traits (elevation, slope, aspect, slope position and soil organic carbon (SOC)) and climate factors (annual precipitation, annual mean air temperature, annual mean maximum air temperature, and annual mean minimum air temperature) on the spatial variations of NPP, ET and WUE across a forest landscape. The analyses were performed using R3.4.3 software (R-Development Core Team 2018). Analysis of variance (ANOVA) was used to test the main and interactive effects of forest type, stand age class, and mode of establishment on NPP, ET and WUE, using GLM procedure of SPSS17.0 for Windows (SPSS Inc., Chicago, USA). The Duncan’s multi-range test was used to compare the differences among the main effects at a confidence level of p <  0.05.

Results

Over the period 2000–2014, the landscape-level annual NPP varied markedly inter-annually, ranging between 311.7 and 573.6 gC·m− 2·a− 1, with an average value of 453.8 gC·m− 2·a− 1 (Fig. 3). Against an uptrend from 2000 to 2005, there occurred a sharp decline in NPP from 2005 to 2007, followed by mostly a recovery trend thereafter until 2012. In contrast, the temporal changes in landscape-level ET were much less variable, keeping within a narrow range of 559.9–603.0 mm·a− 1 (Fig. 3). The pattern of inter-annual variations in landscape-level WUE followed that of NPP, varying over a wide range between 0.54 and 1.01 gC·m− 2·mm− 1 (Fig. 3). There was a significant uptrend in WUE for the period 2000–2006 (p <  0.05), and a tendency of increases in WUE for the later period starting in 2007. The temporal variations in NPP, ET and WUE by forest types were similar to those for the landscape-level.

Fig. 3
figure 3

Inter-annual changes in landscape-averaged NPP, ET and WUE during 2000–2014 in the study area of Lushuihe of Jilin Province, Northeast China

Spatially, the mean annual NPP over the whole study period varied between 205.0 and 639.4 gC·m− 2·a− 1 across the forest landscape, with 99% of the areas showing an increasing trend over the study period (Fig. 4a and b). The mean annual ET ranged between 441.5 and 784.0 mm·a− 1 with an average of 581.5 mm·a− 1 spatially, with most of the areas showing a decreasing trend over the period (Fig. 4c and d). The areas with higher mean annual NPP mostly corresponded to the lower mean annual ET, typically in old-growth coniferous or mixed-wood forest stands (Fig. 2). The spatial distribution of mean annual WUE was similar to that of mean annual NPP (Fig. 4e), with values in the range of 0.46–1.10 gC·m− 2·mm− 1.

Fig. 4
figure 4

Spatial distributions of (a) mean annual NPP, (b) trend of temporal changes in NPP, (c) mean annual ET, (d) trend of temporal changes in ET; (e) mean annual WUE, and (f) trend of temporal changes in WUE, over the period 2000–2014 in the study area of Lushuihe of Jilin Province, Northeast China

Multiple linear regression analysis identified forest type, stand age class, establishment mode, stand density, and three air temperature variables (annual mean, annual mean maximum, and annual mean minimum) as common factors significantly (p <  0.001) affecting the mean annual NPP, ET, and WUE (Table 2). The three air temperature variables differed in their effects on the mean annual NPP, ET, and WUE: while the annual mean air temperature had positive effects on the mean annual NPP and WUE, the effects of the annual mean maximum and minimum air temperatures were negative; the mean annual ET was negatively affected by the annual mean air temperature and positively by the annual mean maximum and minimum air temperatures. Apart from these common factors, mean annual NPP was significantly (p <  0.05) affected by slope and SOC, mean annual ET by slope position and annual precipitation, and mean annual WUE by slope position and SOC (Table 2).

Table 2 Summary of multiple linear regression for the effects of biotic and abiotic factors on the spatial variations of mean annual NPP, ET and WUE for the period 2000–2014 across the forest landscape in Lushuihe of Jilin Province, Northeast China. SOC, soil organic carbon; Tmean, mean annual air temperature; Tmax, mean annual maximum air temperature; Tmin, mean annual minimum air temperature

There were highly significant (p <  0.001) effects of forest type, stand age class, and an interaction between forest type and stand age class on the mean annual NPP, ET and WUE (Table 3). On average, the three forest types were significantly (p <  0.05) differentiated in the mean annual NPP, ET, and WUE: the coniferous forests were highest in NPP (505.3 ± 1.4 gC·m− 2·a− 1; n = 1041) and WUE (0.872 ± 0.004 gC·m− 2·mm− 1; n = 1041), and lowest in ET (584.1 ± 1.6 mm·a− 1; n = 1041), followed by the mixed-wood forests (NPP: 500.5 ± 0.8 gC·m− 2·a− 1; WUE: 0.856 ± 0.02 gC·m− 2·mm− 1; ET: 589.3 ± 0.9 mm·a− 1; n = 2156); whereas the broadleaved forests were lowest in NPP (491.6 ± 0.6 gC·m− 2·a− 1; n = 4428) and WUE (0.832 ± 0.02 gC·m− 2·mm− 1; n = 4428), and highest in ET (594.7 ± 0.6 mm·a− 1; n = 4428). In coniferous forests, there were clear trends of increasing NPP and WUE, and decreasing ET, with stand age class, whilst the mixed-wood forests displayed a weak but similar patterns of responses in NPP, WUE and ET with stand age class (Fig. 5). The broadleaved forests did not show a consistent pattern of change in NPP, ET, and WUE with stand age class (Fig. 5).

Table 3 Summary of F and p values in ANOVAs for testing the effects of forest type and stand age class on the mean annual NPP, ET and WUE for the period 2000–2014 across the forest landscape in Lushuihe of Jilin Province, Northeast China
Fig. 5
figure 5

Histogram of (A) mean annual NPP, (B) mean annual ET, and (C) mean annual WUE by forest type and stand age class over the period 2000–2014 in the study area of Lushuihe of Jilin Province, Northeast China. Vertical bars indicate standard errors. The lowercase letters above the bars designate significant differences among stand age classes within forest types at p <  0.05 (Duncan’s multi-range test)

There were also highly significant (p <  0.001) effects of establishment mode and an interaction between forest type and establishment mode on the mean annual NPP, ET and WUE (Table 4). Regardless of forest type and stand age class, the natural stands had highly significantly (p <  0.001) greater NPP (506.8 ± 0.6 gC·m− 2·a− 1; n = 6037) and WUE (0.877 ± 0.002 gC·m− 2·mm− 1; n = 6037), and smaller ET (582.8 ± 0.7 mm·a− 1; n = 6037), than the planted stands (NPP: 493.1 ± 0.8 gC·m− 2·a− 1; WUE: 0.835 ± 0.003 gC·m− 2·mm− 1; ET: 594.9 ± 1.0 mm·a− 1; n = 1587). However, by identifying forest types, the differences in NPP, ET and WUE between the natural and planted stands were most prominent in coniferous forests and significant (p <  0.05) in mixed-wood forests (Fig. 6). The broadleaved forests did not show any difference in NPP, ET and WUE between the natural and planted forests (Fig. 6).

Table 4 Summary of F and p values in ANOVAs for testing the effects of forest type and establishment mode on the mean annual NPP, ET and WUE for the period 2000–2014 across the forest landscape in Lushuihe of Jilin Province, Northeast China
Fig. 6
figure 6

Histogram of (A) mean annual NPP, (B) ET and (C) WUE by forest type and establishment mode over the period 2000–2014 in the study area of Lushuihe of Jilin Province, Northeast China. Vertical bars indicate standard errors. The lowercase letters above the bars designate significant differences between regeneration modes within forest types at p <  0.05 (Duncan’s multi-range test)

Discussion

In this study, we used the existing MOD17A3-NPP data, the MOD16A3-ET data, and the Category II Forest Inventory data to investigate the spatiotemporal variations and controls of productivity and water use efficiency over a temperate forest landscape in the Changbai Mountain of Northeast China. Our estimates of the means and spatial and temporal variations in NPP and ET based on MODIS products for the period 2000–2014 were in line with the ground-based measurements and model simulations for the general region (Luo 1996; Ni et al. 2001; Liu et al. 2013). Although some cautions may be needed in using the MODIS data which are not locally validated, many regional and global studies of NPP and ET have demonstrated the reliability of the MODIS NPP and ET products (Turner et al. 2006; Zhao and Running 2010; Lathuillière et al. 2012; Mu et al. 2013), including the studies of regional NPP in Chinese ecosystems (Jia et al. 2012; Liu et al. 2015a; Zhang et al. 2015). Nonetheless, there may occur some degree of uncertainty in our estimates of NPP, ET and WUE across the forest landscape.

Temporally, we found that the landscape-level NPP in our study area showed a sharp decline during 2005–2007, coinciding with an extreme drought event during 2005–2006, when an increase in temperature of up to 25.2% and a decline in precipitation of 12.6% were recorded (Zhang et al. 2016).

Previous studies have demonstrated that ecosystem WUE can be influenced by both biotic and abiotic factors (Hu et al. 2008; Liu et al. 2015a; Zhu et al. 2015). Here we found that, across the forest landscape in our study area, all the vegetation factors investigated, including forest type, stand age class and density, and establishment mode, and the three temperature variables, including mean annual mean, annual minimum were common factors significantly (p <  0.05) affecting the mean annual NPP, ET, and WUE, with other factors of site conditions and climate playing differential role in affecting NPP, ET and WUE. Clearly, forest type, stand density and developmental stage are key determinants of ecosystem NPP, ET and WUE in our study area, similar to the findings in other studies (Law et al. 2003; Irvine et al. 2004; Wang et al. 2011; Xiao et al. 2013; Neumann et al. 2015).

Also, climatic factors are important determinants of vegetation and play a crucial role in shaping spatial patterns and temporal dynamics of forest NPP, which have been investigated in many studies. Ni et al. (2001) reported that temperature and moisture are the dominant factors controlling the spatial distribution of NPP in China and the NPP of Chinese forests is positively correlated with annual mean temperature and precipitation. Mao et al. (2014) found that NPP was positively correlated with mean temperature, negatively with precipitation. Our result was consistent with the findings of Ni et al. (2001), although the effect of precipitation was not significant. This might be due to the relative narrow range of mean annual precipitation (905.8–995.0 mm·a− 1) across the forest landscape in our study area, means that precipitation was not the limited factor for NPP and partially lack of fine-grained data points for matching climatic variables and ecosystem NPP values. Interestingly, we found negative impacts of increasing mean annual minimum and maximum temperatures on NPP, possibly explainable by the effects of drought stress induced both by rising maximum temperature and insufficient winter-snows inferred by the rising minimum temperature.

The responses of WUE to climate factors are more complex as it depends upon the variations of NPP, ET and the extent to which the system is water limited (Huang et al. 2015). In this study, WUE was negatively affected by precipitation. Although both NPP and ET were positively affected by annual precipitation, an increase in precipitation would likely induce a greater increase in ET than NPP, consequently leading to decreased WUE (Huang et al. 2015). Moreover, increases in precipitation during the non-growing season could also aggravate nitrogen limitation due to enhanced nitrogen leaching and denitrification (Hovenden et al. 2014), resulting in a reduction in NPP and declining WUE. Contrary to the effect of precipitation, annual mean air temperature had a significant positive effect on WUE. This might be explained by that, as the length of growing season increases due to an increase in air temperature (Körner and Basler 2010; Gunderson et al. 2012), the photosynthetic rates are accelerated (Flanagan and Syed 2011), whilst the evapotranspiration rates tend to remain unchanged or decreased due to the stomatal regulation (Serrat-Capdevila et al. 2011). Consequently, the greater increase in NPP compared to ET leads to a positive trend of WUE.

In addition, the site conditions reflect the water and fertility of the forest stand, which indirectly affects the spatial pattern of NPP, ET and WUE. Our results showed that NPP was greatly influenced by site factors, including slope and SOC, such that NPP tends to be greater in the areas with steep slopes, where mostly primary forests remain intact due to access and harvesting difficulties in history.

Across the forest landscape of Lushuihe and over the study period 2000–2014, the mean annual NPP, ET and WUE were affected by forest types, stand age class, and how the forest stands were originated (i.e. establishment mode); the effects of stand age class and establishment mode appeared to be forest-type specific. Among the three forest types in this study, coniferous forests were generally highest in NPP and lowest in ET, hence highest in WUE, followed by mixed-wood forests. This contradict with some of the previous studies showing highest NPP values in mixed-wood forests compared with the coniferous and broadleaved forests in Northeast China (Zhao et al. 2011; Mao et al. 2014). In our study area, the coniferous forests were dominated by the Korean pine stands, which are climax forest type of the region and have been well protected and experienced less human disturbance. Typically in coniferous forests, with some tendency in mixed-wood forests, the mean annual NPP and WUE increased with stand age class and were significantly higher in naturally established stands than in planted stands. In well agreement with a recent study at the national scale by Yu et al. (2019), we found that naturally forests were indeed better at both productivity and water conservation than planted forests, but the effects appeared to be reserved for forest-type specific and demonstrated in stands with presence of coniferous trees at the landscape level. In broadleaved forests, there were no changes in NPP and WUE with stand age class and establishment mode, similar to the finding of Xiao et al. (2013) for a range of Chinese ecosystems.

Our findings of increasing NPP with stand age class up to old-growth forests were somewhat in contrary to commonly accepted view that NPP declines with age or peaks at middle-age then declines in the old age (Law et al. 2003; Pregitzer and Euskirchen 2004; Wang et al. 2011; Tang et al. 2014). This pattern of declining NPP in old-growth forests was based on the assumption that changes in productivity with age in complex, multi-aged, multispecies natural forests can be modelled simply as scaled-up versions of individual trees or even-aged stands (Carey et al. 2001), neglecting other factors such as climate, management history, successional stage, stand disturbance, stand density as well as soil fertility (Pregitzer and Euskirchen 2004; Vicca et al. 2012). There were studies showing that decrease in NPP with age was not general across ecoregions and with no decline in old stands (> 200 years old) in some ecoregions (Van Tuyl et al. 2005; Hudiburg et al. 2009). This signifies the importance of considering forest management history and historical disturbances in shaping the productivity-stand age relationship. Unlike in many other parts of the world, wood harvesting in Chinese forests had historically progressed with consideration of easy access and terrain conditions, translating to incomparable disturbances and stand age classes across the landscapes in forest regions.

Our study showed a significant decrease in ET with age class in coniferous forest, which is likely a consequence of decreases in transpiration with age as evidenced in many previous studies (Köstner 2001; Vertessy et al. 2001; Delzon and Loustau 2005), explainable by a decrease in stomatal conductance in old and taller trees (Irvine et al. 2004; Delzon and Loustau 2005).

Conclusions

In this study, across a temperate forest landscape in Northeast China, we depicted the patterns of spatial and temporal variations in ecosystem productivity and water usage and controlling factors. Forest type, stand development, and the regeneration scheme all played important roles in determining the ecosystem productivity and relationship with water usage. Overall, the ecosystem productivity and water usage relationship does not necessarily always occur as a simple tradeoff. Forest type and stand-age class are factors of key consideration in planning and managing forest landscape for achieving the compromise between regional carbon balance and water usage. Facilitating the establishment of coniferous trees and conservation of mature forest stands, and near-natural management of planted forests, can be favorable options for delivering ecosystem services of carbon sequestration and water conservation in temperate regions. Facilitating the establishment of coniferous trees and conservation of mature forest stands, and adoption of near-natural management of planted forests, can be viewed as favorable options for delivering ecosystem services of carbon sequestration and water conservation in temperate regions.

Abbreviations

ET:

Evapotranspiration

NPP:

Net primary productivity

Pre:

Mean annual precipitation

SOC:

Soil organic carbon

T max :

Annual mean maximum temperature

T mean :

Annual mean temperature

T min :

Annual mean minimum temperature

WUE:

Water use efficiency

References

  • Carey EV, Sala A, Keane R, Callaway RM (2001) Are old forests underestimated as global carbon sinks? Glob Change Biol 7:339–344

    Article  Google Scholar 

  • Delzon S, Loustau D (2005) Age-related decline in stand water use: sap flow and transpiration in a pine forest chronosequence. Agric For Meteorol 129:105–119

    Article  Google Scholar 

  • Fang XM, Yu DP, Zhou WM, Zhou L, Dai LM (2016) The effects of forest type on soil microbial activity in Changbai Mountain, Northeast China. Ann For Sci 73:473–482

    Article  Google Scholar 

  • Farley KA, Jobbágy EG, Jackson RB (2005) Effects of afforestation on water yield: a global synthesis with implications for policy. Glob Change Biol 11:1565–1576

    Article  Google Scholar 

  • Feng XM, Fu BJ, Piao SL, Wang S, Ciais P, Zeng ZZ, Lü YH, Zeng Y, Li Y, Jiang XH, Wu BF (2016) Revegetation in China’s loess plateau is approaching sustainable water resource limits. Nat Clim Chang 6:1019–1022

    Article  Google Scholar 

  • Flanagan LB, Syed KH (2011) Stimulation of both photosynthesis and respiration in response to warmer and drier conditions in a boreal peatland ecosystem. Glob Change Biol 17:2271–2287

    Article  Google Scholar 

  • Gao Y, Zhu XJ, Yu GR, He NP, Wang QF, Tian J (2014) Water use efficiency threshold for terrestrial ecosystem carbon sequestration in China under afforestation. Agric For Meteorol 195–196:32–37

    Article  Google Scholar 

  • Gulbeyaz O, Bond-Lamberty B, Akyurek Z, West TO (2018) A new approach to evaluate the MODIS annual NPP product (MOD17A3) using forest field data from Turkey. Int J Remote Sens 39:2560–2578

    Article  Google Scholar 

  • Gunderson CA, Edwards NT, Walker AV, O’Hara KH, Campion CM, Hanson PJ (2012) Forest phenology and a warmer climate - growing season extension in relation to climatic provenance. Glob Change Biol 18:2008–2025

    Article  Google Scholar 

  • Heinsch FA, Zhao MS, Running SW, Kimball JS, Nemani RR, Davis KJ, Bolstad PV, Cook BD, Desai AR, Ricciuto DM, Law BE, Oechel WC, Kwon H, Luo HY, Wofsy SC, Dunn AL, Munger JW, Baldocchi DD, Xu LK, Hollinger DY, Richardson AD, Stoy PC, Siqueira MBS, Monson RK, Burns SP, Flanagan LB (2006) Evaluation of remote sensing based terrestrial productivity from MODIS using regional tower eddy flux network observations. IEEE T Geosci Remote 44:1908–1925

    Article  Google Scholar 

  • Hovenden MJ, Newton PCD, Wills KE (2014) Seasonal not annual rainfall determines grassland biomass response to carbon dioxide. Nature 7511:583–586

    Article  Google Scholar 

  • Hu ZM, Yu GR, Fu YL, Sun XM, Li YN, Shi PL, Wang YF, Zheng ZM (2008) Effects of vegetation control on ecosystem water use efficiency within and among four grassland ecosystems in China. Glob Change Biol 14:1609–1619

    Article  Google Scholar 

  • Huang MT, Piao SL, Sun Y, Ciais P, Cheng L, Mao JF, Poulter B, Shi XY, Zeng ZZ, Wang YP (2015) Change in terrestrial ecosystem water-use efficiency over the last three decades. Glob Change Biol 21:2366–2378

    Article  Google Scholar 

  • Huang XJ, Ma MG, Wang XF, Tang XG, Yang H (2018) The uncertainty analysis of the MODIS GPP product in global maize croplands. Front Earth Sci 12:739–749

    Article  CAS  Google Scholar 

  • Hudiburg T, Law BE, Turner DP, Campbell J, Donato DC, Duane M (2009) Carbon dynamics of Oregon and northern California forests and potential land-based carbon storage. Ecol Appl 19:163–180

    Article  Google Scholar 

  • Irvine J, Law BE, Kurpius MR, Anthoni PM, Moore D, Schwarz PA (2004) Age-related changes in ecosystem structure and function and effects on water and carbon exchange in ponderosa pine. Tree Physiol 24:753–763

    Article  CAS  Google Scholar 

  • Jackson RB, Jobbágy EG, Avissar R, Roy SB, Barrett DJ, Cook CW, Farley KA, le Maitre DC, McCarl BA, Murray BC (2005) Trading water for carbon with biological carbon sequestration. Science 310:1944–1947

    Article  CAS  Google Scholar 

  • Jia ZZ, Liu SM, Xu ZW, Chen YJ, Zhu MJ (2012) Validation of remotely sensed evapotranspiration over the Hai River Basin, China. J Geophys Res-Atmos 117:D13113.

    Article  Google Scholar 

  • Körner C, Basler D (2010) Phenology under global warming. Science 327:1461–1462

    Article  Google Scholar 

  • Köstner B (2001) Evaporation and transpiration from forests in Central Europe-relevance of patch-level studies for spatial scaling. Meteorog Atmos Phys 76:69–82

    Article  Google Scholar 

  • Lathuillière MJ, Johnson MS, Donner SD (2012) Water use by terrestrial ecosystems: temporal variability in rainforest and agricultural contributions to evapotranspiration in Mato Grosso Brazil. Environ Res Lett 7:024024. https://doi.org/10.1088/1748-9326/7/2/024024

    Article  Google Scholar 

  • Law BE, Sun OJ, Campbell J, Van Tuyl S, Thornton PE (2003) Changes in carbon storage and fluxes in a chronosequence of ponderosa pine. Glob Change Biol 9:510–524

    Article  Google Scholar 

  • Liu XC, Xu ZX, Yu RH (2011) Trends of climate variability in China during the past decades. Climat Change 109:503–516

    Article  Google Scholar 

  • Liu YB, Xiao JF, Ju WM, Xu K, Zhou YL, Zhao YT (2016) Recent trends in vegetation greenness in China significantly altered annual evapotranspiration and water yield. Environ Res Lett 11:094010. https://doi.org/10.1088/1748-9326/11/9/094010

    Article  Google Scholar 

  • Liu YB, Xiao JF, Ju WM, Zhou YL, Wang SQ, Wu XC (2015a) Water use efficiency of China's terrestrial ecosystems and responses to drought. Sci Rep-Uk 5:13799. https://doi.org/10.1038/srep13799

    Article  Google Scholar 

  • Liu YB, Zhou YL, Ju WM, Chen JM, Wang SQ, He HL, Wang HM, Guan DX, Zhao FH, Li YL, Hao YB (2013) Evapotranspiration and water yield over China's landmass from 2000 to 2010. Hydrol Earth Syst Sc 17:4957–4980

    Article  Google Scholar 

  • Liu ZJ, Shao QQ, Liu JY (2015b) The performances of MODIS-GPP and -ET products in China and their sensitivity to input data (FPAR/LAI). Remote Sens 7:135–152

    Article  CAS  Google Scholar 

  • Luo TX (1996) Patterns of biological production and its mathematical models for main forest types of China. Ph.D. Dissertation, Committee of Synthesis Investigation of Natural Resources, Chinese Academy of Sciences, Beijing

    Google Scholar 

  • Mao DH, Wang ZM, Wu CS, Song KS, Ren CY (2014) Examining forest net primary productivity dynamics and driving forces in northeastern China during 1982–2010. Chinese Geogr Sci 24:631–646

    Article  Google Scholar 

  • Medrano H, Flexas J, Galmés J (2009) Variability in water use efficiency at the leaf level among Mediterranean plants with different growth forms. Plant Soil 317:17–29

    Article  CAS  Google Scholar 

  • Mu QZ, Zhao MS, Kimball JS, McDowell NG, Running SW (2013) A remotely sensed global terrestrial drought severity index. B Am Meteorol Soc 94:83–98

    Article  Google Scholar 

  • Mu QZ, Zhao MS, Running SW (2011) Improvements to a MODIS global terrestrial evapotranspiration algorithm. Remote Sens Environ 115:1781–1800

    Article  Google Scholar 

  • Neumann M, Zhao MS, Kindermann G, Hasenauer H (2015) Comparing MODIS net primary production estimates with terrestrial national forest inventory data in Austria. Remote Sens 7:3878–3906

    Article  Google Scholar 

  • Ni J, Zhang XS, Scurlock JMO (2001) Synthesis and analysis of biomass and net primary productivity in Chinese forests. Ann For Sci 58:351–384

    Article  Google Scholar 

  • Peng J, Dong WJ, Yuan WP, Zhang Y (2012) Responses of grassland and forest to temperature and precipitation changes in Northeast China. Adv Atmos Sci 29:1063–1077

    Article  Google Scholar 

  • Poni S, Bernizzoni F, Civardi S, Gatti M, Porro D, Camin F (2009) Performance and water-use efficiency (single-leaf vs. whole-canopy) of well-watered and half-stressed split-root Lambrusco grapevines grown in Po Valley (Italy). Agric Ecosyst Environ 129:97–106

    Article  Google Scholar 

  • Pregitzer KS, Euskirchen ES (2004) Carbon cycling and storage in world forests: biome patterns related to forest age. Glob Change Biol 10:2052–2077

    Article  Google Scholar 

  • Qi G, Chen H, Zhou L, Wang XC, Zhou WM, Qi L, Yang YH, Yang FL, Wang QL, Dai LM (2016) Carbon stock of larch plantations and its comparison with an old-growth forest in Northeast China. Chinese Geogr Sci 26:10–21

    Article  Google Scholar 

  • Running SW, Nemani RR, Heinsch FA, Zhao MS, Reeves M, Hashimoto H (2004) A continuous satellite-derived measure of global terrestrial primary production. BioScience 54:547–560

    Article  Google Scholar 

  • Serrat-Capdevila A, Scott RL, Shuttleworth WJ, Valdés JB (2011) Estimating evapotranspiration under warmer climates: insights from a semi-arid riparian system. J Hydrol 39:1–11

    Article  Google Scholar 

  • Sun G, Zhou GY, Zhang ZQ, Wei XH, McNulty SG, Vose JM (2006) Potential water yield reduction due to forestation across China. J Hydrol 328:548–558

    Article  Google Scholar 

  • Sun G, Zuo CQ, Liu SY, Liu ML, McNulty SG, Vose JM (2008) Watershed evapotranspiration increased due to changes in vegetation composition and structure under a subtropical climate. J Am Water Resour Assoc 44:1164–1175

    Article  Google Scholar 

  • Sun YH, Meng L, Tian L, Li GL, Sun OJ (2016) Assessing current stocks and future sequestration potential of forest biomass carbon in Daqing Mountain nature Reserve of Inner Mongolia, China. J Forestry Res 27:931–938

    Article  CAS  Google Scholar 

  • Tang JW, Luyssaert S, Richardson AD, Kutsch W, Janssens IA (2014) Steeper declines in forest photosynthesis than respiration explain age-driven decreases in forest growth. PNAS 111:8856–8860

    Article  CAS  Google Scholar 

  • Tang XG, Li HP, Xu XB, Luo JH, Li XY, Ding Z, Xie J (2016) Potential of MODIS data to track the variability in ecosystem water-use efficiency of temperate deciduous forests. Ecol Eng 91:381–391

    Article  Google Scholar 

  • Tian HQ, Chen GS, Liu ML, Zhang C, Sun G, Lu CQ, Xu XF, Ren W, Pan SF, Chappelka A (2010) Model estimates of net primary productivity, evapotranspiration, and water use efficiency in the terrestrial ecosystems of the southern United States during 1895–2007. Forest Ecol Manag 259:1311–1327

    Article  Google Scholar 

  • Tong XJ, Zhang JS, Meng P, Li J, Zheng N (2014) Ecosystem water use efficiency in a warm-temperate mixed plantation in the North China. J Hydrol 512:221–228

    Article  CAS  Google Scholar 

  • Turner DP, Ritts WD, Cohen WB, Gower ST, Running SW, Zhao MS, Costa MH, Kirschbaum AA, Ham JM, Saleska SR, Ahl DE (2006) Evaluation of MODIS NPP and GPP products across multiple biomes. Remote Sens Environ 102:282–292

    Article  Google Scholar 

  • Van Tuyl S, Law BE, Turner DP, Gitelman AI (2005) Variability in net primary production and carbon storage in biomass across Oregon forests - an assessment integrating data from forest inventories, intensive sites, and remote sensing. Forest Ecol Manag 209:273–291

    Article  Google Scholar 

  • Vertessy RA, Watson FGR, O’Sullivan SK (2001) Factors determining relations between stand age and catchment water balance in mountain ash forests. Forest Ecol Manag 143:13–26

    Article  Google Scholar 

  • Vicca S, Luyssaert S, Peñuelas J, Campioli M, Chapin FS, Ciais P, Heinemeyer A, Högberg P, Kutsch WL, Law BE, Malhi Y, Papale D, Piao SL, Reichstein M, Schulze ED, Janssens IA (2012) Fertile forests produce biomass more efficiently. Ecol Lett 15:520–526

    Article  CAS  Google Scholar 

  • Wang SQ, Zhou L, Chen JM, Ju WM, Feng XF, Wu WX (2011) Relationships between net primary productivity and stand age for several forest types and their influence on China’s carbon balance. J Environ Manag 92:1651–1662

    Article  Google Scholar 

  • Wang YF, Cao SX (2011) Carbon sequestration may have negative impacts on ecosystem health. Environ Sci Technol 45:1759–1760

    Article  CAS  Google Scholar 

  • Wu J, Gao XJ (2013) A gridded daily observation dataset over China region and comparison with the other datasets. Chin J Geophys 56:1102–1111

    Google Scholar 

  • Wu J, Gao XJ, Giorgi F, Chen DL (2017) Changes of effective temperature and cold/hot days in late decades over China based on a high resolution gridded observation dataset. Int J Climatol 37:788–800

    Article  CAS  Google Scholar 

  • Xia L, Wang F, Mu XM, Jin K, Sun WY, Gao P, Zhao GJ (2015) Water use efficiency of net primary production in global terrestrial ecosystems. J Earth Syst Sci 124:921–931

    Article  CAS  Google Scholar 

  • Xiao JF, Sun G, Chen JQ, Chen H, Chen SP, Dong G, Gao SH, Guo HQ, Guo JX, Han SJ, Kato T, Li YL, Lin GH, Lu WZ, Ma MG, McNulty S, Shao CL, Wang XF, Xie X, Zhang XD, Zhang ZQ, Zhao B, Zhou GS, Zhou J (2013) Carbon fluxes, evapotranspiration, and water use efficiency of terrestrial ecosystems in China. Agric For Meteorol 182–183:76–90

    Article  Google Scholar 

  • Yu DP, Zhou L, Zhou WM, Ding H, Wang QW, Wang Y, Wu XQ, Dai LM (2011) Forest management in Northeast China: history, problems, and challenges. Environ Manag 48:1122–1135

    Article  Google Scholar 

  • Yu GR, Song X, Wang QF, Liu YF, Guan DX, Yan JH, Sun XM, Zhang LM, Wen XF (2008) Water-use efficiency of forest ecosystems in eastern China and its relations to climatic variables. New Phytol 177:927–937

    Article  CAS  Google Scholar 

  • Yu Z, Liu SR, Wang JX, Wei XH, Schuler J, Sun PS, Harper R, Zegre N (2019) Natural forests exhibit higher carbon sequestration and lower water consumption than planted forests in China. Glob Change Biol 25:68–77

    Article  Google Scholar 

  • Zhang L, Tian J, He HL, Ren XL, Sun XM, Yu GR, Lu QQ, Lv LY (2015) Evaluation of water use efficiency derived from MODIS products against eddy variance measurements in China. Remote Sens 7:11183–11201

    Article  Google Scholar 

  • Zhang T, Peng J, Liang W, Yang YT, Liu YX (2016) Spatial-temporal patterns of water use efficiency and climate controls in China’s loess plateau during 2000-2010. Sci Total Environ 565:105–122

    Article  CAS  Google Scholar 

  • Zhao GS, Wang JB, Fan WY, Ying TY (2011) Vegetation net primary productivity in Northeast China in 2000-2008: simulation and seasonal change. Chin J Appl Ecol 22:621–630

    CAS  Google Scholar 

  • Zhao MS, Heinsch FA, Nemani RR, Running SW (2005) Improvements of the MODIS terrestrial gross and net primary production global dataset. Remote Sens Environ 95:164–176

    Article  Google Scholar 

  • Zhao MS, Running SW (2010) Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. Science 329:940–943

    Article  CAS  Google Scholar 

  • Zhu JJ, Liu ZG, Wang HX, Yan QL, Fang HY, Hu LL, Yu LZ (2008) Effects of site preparations on emergence and early establishment of Larix olgensis in montane regions of northeastern China. New Forest 36:247–260

    Article  Google Scholar 

  • Zhu XJ, Yu GR, Wang QF, Hu ZM, Zheng H, Li SG, Sun XM, Zhang YP, Yan JH, Wang HM, Zhao FH, Zhang JH, Shi PL, Li YN, Zhao L, Zhang FW, Hao YB (2015) Spatial variability of water use efficiency in China’s terrestrial ecosystems. Glob Planet Chang 129:37–44

    Article  Google Scholar 

Download references

Acknowledgements

We gratefully acknowledge the Lushuihe Forestry Bureau for providing the forest resource inventory data. We would like to thank Xiyue Zhang for her assistance in pre-processing the gridded meteorological data.

Funding

This research was funded by the Program of Forestry Research for the Public Benefits of Ministry of Finance of China (Grant No. 201404201).

Availability of data and materials

The original data are not owned by the authors and cannot be shared. The recomputed and processed data are available upon request.

Author information

Authors and Affiliations

Authors

Contributions

The conceptual framework for the study was developed by OJS and SL. XL collected and analyzed the data, and wrote the draft manuscript. TJAF and CJ contributed to data analysis and results interpretations. All authors contributed to discussions and writing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Osbert Jianxin Sun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, X., Farooqi, T.J.A., Jiang, C. et al. Spatiotemporal variations in productivity and water use efficiency across a temperate forest landscape of Northeast China. For. Ecosyst. 6, 22 (2019). https://doi.org/10.1186/s40663-019-0179-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40663-019-0179-x

Keywords