Estimation of Forest NPP and Carbon Sequestration in the Three Gorges Reservoir Area , Using the Biome-BGC Model

The Three Gorges Reservoir area is one of the most ecologically sensitive areas in China, and the forest landscape pattern in this region shows dramatic change due to the influence of the Three Gorges reservoir project. In this study, the locally parameterized Biome-BGC model, generated with long-term meteorological monitoring data, was used to simulate net primary productivity (NPP) and carbon density of the vegetation layer, the litter layer, and the soil layer for various forest types from 1992 to 2012 in this area. The total and unitary forest NPP presented obvious annual fluctuation under the combined influences of land use change and extreme weather events. Apart from the year 2006, from 1992 to 2012, the NPP values of each forest type showed an increasing trend, although the growth rates decreased. In 2006, due to abnormally high air temperatures and less precipitation, total and unit area forest NPP values decreased by 46.3% and 53.9%, respectively, compared to 2002. From 1992 to 2012, the carbon stocks of the forest vegetation layer, the litter layer, the soil layer, and the entire area gradually increased with decreasing growth rates. Additionally, forest carbon stocks were high in the east and the south and low in the west and the north. Generally, the forest productivity is greatly affected by the physiological and ecological characteristics of the plants themselves as well as the environmental factors, whereas total forest productivity is largely influenced by human activities. The increase in forest area and the optimization of the forest landscape pattern could improve the forest productivity and carbon sequestration.


Introduction
Climate warming is currently receiving considerable attention and has become a global research focus.Forests are the main vegetation cover in terrestrial ecosystems, and forest productivity and biomass play an important role in global climate change and circulation of substances [1].Via photosynthesis, forests fix atmospheric CO2 and accumulate carbon in their biomass, litter, and soil.In addition, due to biological respiration, litter decomposition, soil carbon degradation, and natural or artificial disturbances, forests emit carbon into the atmosphere in the form of CO2.Although the forest area only accounts for about one third of the total terrestrial ecosystem area, carbon storage in aboveground forest vegetation is equivalent to 80% of the aboveground total organic carbon of terrestrial ecosystems; belowground carbon storage accounts for 40% of that of the entire terrestrial ecosystem [2].
The Three Gorges Reservoir area (TGRA) of the Yangtze River is a new concept derived from the construction of the Three Gorges Dam Project since 1992.It refers to the region submerged due to the water level rise caused by the Three Gorges Dam (water level 175 m; height of the dam wall 185 m).The Three Gorges Reservoir Area (TGRA) of the Yangtze River is the most ecologically sensitive and fragile area in China with complex terrain and frequent natural disasters, which is located in the transition zone between a mountainous region and watershed region.As a special human activity disturbance, the Land Use/Cover Change (LUCC) was seriously affected by the construction of the reservoir, including the establishment of new towns and cities, migration of millions of people and the extended railway and highway nets [3].In recent decades, the built-up land increased the most, largely and rapidly, followed by water bodies.They were transformed from croplands, forests and grasslands.
Against the background of climate change and land use change, the present situation of the forest ecosystem in the Three Gorges Reservoir area, as well as its response to global climate change, has become the focus of a large number of studies both in China and around the world [4][5][6].This study focused on the temporal and spatial variations in forest productivity and carbon storage in the TGRA as well as the corresponding influencing factors.The research results are significant in revealing the extent and scope of human disturbance effects on forest productivity and forest carbon storage and to provide evidence for predicting the response of forest productivity to future climate change scenarios.Therefore, a locally parameterized Biome-BGC model was used to simulate net primary productivity (NPP) and carbon density of the vegetation layer, the litter layer, and the soil layer for various forest types from 1992 to 2012 in the TGRA in order to reveal the changes.

Study Area
The Three Gorges Reservoir Area is located at 106°16'-111°28'E, and 28°56'-31°44'N and covers a total of 26 district (counties), 22 of which are located in Chongqing Municipality directly under the Central Government while 4 are situated in Hubei Province.The reservoir catchment covers an area of about 5.8 × 10 4 km 2 , with a total population of about 22 million people.The climate is humid subtropical monsoon, with an annual average temperature of 16-18 °C and an annual precipitation of 1000-1300 mm, with early springs, hot summers, wet autumns, and warm winters [7].
The TGRA is situated on the eastern margin of the second step of the three geomorphic steps in China.It is at the junction of three tectonic units: the Daba Mountain fold belt, the parallel ridge valley in eastern Sichuan and the Sichuan-Hubei-Hunan-Guizhou uplift fold belt.Mountains account for about 74%, hills account for about 22%, plains and dams account for only 4%.The terrain of this area is high in the east and low in the west with an altitude range from 50 meters to 2900 meters.From the east to the west, the area can be divided into three parts: the eastern mountain part, the central parallel ridge valley part and the western hills (Figure 1).
The TGRA has rich plant resources.It contains 6088 species of vascular plants, belonging to 1428 genera and 208 families, and the species in this area account for about 20% of the total plant species in China [8].The main species are Pinaceae、Cupressaceae、Juglandaceae、Fagaceae、 Hamamelidaceae 、 Roseceae 、 Rutaceae 、 Anarcardiaceae 、 Theaceae 、 Cornaceae 、 Ericaceae 、 Betulaceae、Leguminosae、Gramineae and Pteridophyta [9,10].Based on the characteristics of the forest types and the main species composition, the forest vegetation in the TGRA can be divided into four major types, namely, coniferous forest, deciduous broadleaf forest, evergreen broadleaf forest, and shrub forest [11,12].Due to excessive exploitation and interference for decades, the natural forest area is extremely rare in this region, and most of the natural forests are secondary forests.The forest type structure is relatively simple, and Pinus massoniana and Cupressus funebris are the most widely distributed species with the highest abundance [13].The forests are mainly distributed in the eastern mountain part, the western Karst Mountains, and the low hilly areas in the south (Figure 2).

Spin-up and Run the Biome-BGC Model
The Biome-BGC model is one of the typical biogeochemical models.The Biome-BGC version 4.2 (Numerical Terradynamic Simulation Group, Missoula, MT, USA) was used in this study, which can predict the flow and storage of carbon, nitrogen, water, and radiation in terrestrial vegetation and soil, with a time-step of one day, based on daily meteorological data.In addition, the Biome-BGC model considers how global climate change and circulation of substances are affected by environmental conditions such as temperature, precipitation, solar radiation, soil texture, and atmospheric CO2 concentration [14][15][16][17].The Biome-BGC model has been widely used.For example, Ueyama M et al. applied this model to simulate the carbon flux of an Alaska black pine forest and verified the model by comparing the simulated values with the measured values by the eddy covariance method.Their results showed that gross primary productivity (GPP), net ecosystem carbon exchange capacity (NEE), and ecosystem respiration rates (RE) obtained by the two methods were highly correlated [18].Similarly, they also found in another study that the GPP, NEE, and RE values simulated by the Biome-BGC model were highly consistent with the observed values based on AsiaFlux [19].In China, this model has also been used to simulate the characteristics of NPP under current climatic conditions and to predict the NPP characteristics under future climate change scenarios, not only for a particular tree species, but also for forest ecosystems in different regions, such as Picea schrenkiana in Tianshan Mountains of Xinjiang [20], a Pinus elliottii plantation in a hilly Krasnozem area [21], a Quercus variabilis forest in the Beijing mountainous area [22], a Quercus variabilis forest [23], a Larix principis-rupprechtii forest in Beijing mountain area [24], a Robinia pseudoacacia forest [25], as well as the forest ecosystem in Fujian Province located in a subtropical zone and which has the highest forest coverage in China [26], and a forest ecosystem on the northern slope of Tianshan Mountains in the arid region of China [27].Compared with and verified by the measured values estimated according to the classical ecological survey method, the simulated values had a high correlation.Therefore, the Biome-BGC model can not only reproduce dynamic characteristics of short-term carbon fluxes but can also simulate long-term distribution patterns and dynamics of carbon storage as well as the responses of these processes to environmental changes, including climate change [28][29][30][31][32][33].
Firstly, the Biome-BGC model was spun-up based on a set of long-term meteorological data, site data, and ecophysiological data.As a hypothesis, the CO2 concentration values and nitrogen deposition values prior to the industrial revolution were used as a steady state ecosystem.After repeated and continuous simulation, the model reached an equilibrium state.In the state of equilibrium, the meteorological data, site data, and ecophysiological data during the study period were launched into the Biome-BGC in order to run the model and then obtain the final results [14,34].

Parameterization of the Biome-BGC Model
Three files are required to run the Biome-BGC model: initialization file, meteorological data file, ecophysiological constants file.
Initialization files are used to control the running process of the model.Firstly, the reading paths of the prepared meteorological and ecophysiological parameters files are specified.Then, the environmental information of each county needs to be described (Table 1), including latitude, altitude, effective depth of soil, soil texture, albedo, etc.Finally, the result variables of model simulation output are defined.The soil characteristic parameters in this study were based on the "National 1:1 million Digital Soil Map" developed by Nanjing Soil Research Institute of Chinese Academy of Sciences and Soil Environment Department of the Ministry of Agriculture of China.The meteorological data file is comprised of daily maximum temperature, daily minimum temperature, daily average temperature, daily precipitation, daily water vapor pressure difference, daily short wave radiation, and day length.Measured values from the national ground stations of China Meteorological Administration within the TGRA were used; each district (county) had a located meteorological station (Figures 1 and 2 shown the distribution of the stations), with a total of 34 years (from 1981 to 2014) of observations.
The ecophysiological parameter files of the Biome-BGC model are based on seven types: deciduous broadleaf forest (DBF), evergreen broadleaf forest (EBF), evergreen needle leaf forest (ENF), deciduous needle leaf forest (DNF), shrub forest (SHRUB), C3 grassland, and C4 grassland.The Numerical Terradynamic Simulation Group (NTSG) of the Montana University provided a set of default ecophysiological parameters, containing 42 parameters for each vegetation type.These default values are assigned based on the average value of parameters obtained from the literature research or the value determined by clustering analysis after literature retrieval and statistical analysis of each parameter for each community [35].According to the forest vegetation types and the main species in the TGAR, the ecophysiological constants of Pinus massoniana and oaks were used for the ENF and DBF, while the default parameters were used for the EBF and SHRUB.For the mixed coniferous broadleaved forest (MIX), the Biome-BGC model does not provide the corresponding ecophysiological parameters, making a direct assessment of these parameters impossible.However, previous studies have shown that we can effectively simulate the productivity of MIX forests by simulating the productivity of both vegetation types respectively and using a subsequent proportion conversion like 50% for each forest type [36,37].Some parameters, such as carbon allocation from roots to leaves, carbon allocation from stems to leaves, carbon allocation from thick roots to stems, carbon to nitrogen ratio of leaves, carbon to nitrogen ratio of leaf litter, carbon to nitrogen ratio of roots, and carbon to nitrogen ratio of woody parts, can be locally parameterized by actual measurement or literature findings [22,23,26,[38][39][40]; for other parameters that are difficult to be obtained, default values can be adopted, which are provided by the model (Table 2).

Validation of The Biome-BGC Model
In this study, the simulation results of the Biome-BGC model were validated through comparative analysis of observed values and simulated values, combined with relevant literature findings.The observed values of NPP were derived from the analytic-tree fitted biomass model, obtained from a sample-plot survey based on the sub-compartment data of forest resource inventory and planning projects in 2003 and 2013 [41,42].The carbon content in various plant organs, litters, and A and B layers of soil are determined by experiment, and the observed values of carbon density were estimated according to the carbon storage obtained from the biomass and the area data [43].Fu et al. used the CBM-CFS3 model as the basic platform and simulated annual NPP values (2.20-5.08Mg C hm −2 year −1 ) and carbon densities of various forest types in the Three Gorges Reservoir area by plotting age-harvest growth curves for the major forest types [44].Wang et al., based on forest resources inventory data from three periods (6, 7, and 8), used the forest stock volume (increment) extension method to estimate average NPP values of forests in the Three Gorges Reservoir area; NPP ranged between 2.85 and 6.19 Mg C hm −2 year −1 [45].The comparison results are shown in Tables 3  and 4. By comparison, the NPP values of the ENF and the Mixed coniferous and broadleaf forests, simulated by the Biome-BGC model, were slightly larger than the sample-plot measured values.The simulated carbon density values of the ENF and the Mixed coniferous and broad-leaf forest were close to the measured values, whereas the simulated carbon density values of the DBF, the EBF, and the Shrub forest were lower than the measured values.In spite of this, the values simulated via the Biome-BGC model were consistent with the values simulated by the CBM-CFS3 model and agreed with the estimation results based on the forest resource inventory.Considering that our research focused on the regional scale, the simulation results based on the Biome-BGC model can accurately reflect the real situation, verifying that this model is suitable for the Three Gorges Reservoir area.

Data Processing and Analysis
According to the different construction stages of the Three Gorges project, we selected five research time-points, namely, 1992, 1996, 2002, 2006, and 2012, to represent five stages of the project process, including the evaluation stage, main channel cutoff, formal impoundment, construction completion, and impoundment to 175 m for full operation, respectively.On this basis, the simulation results from the BIOME-BGC model together with the areas of different forest types derived from remote sensing images at research time-points with spatial geography information were used to analyze the dynamic changes in forest productivity and carbon storage in the reservoir area.
The total forest net primary productivity (NPP) for each respective year was obtained by calculating the NPP values from the Biome-BGC model and the corresponding areas of various forest types for each district (county).
where NPPtotal is the total forest NPP in a certain year (Tg C); NPPi,j and Si,j represent the NPP values (t hm −2 year −1 ) and area (hm −2 ) of the i-th (i = 1,2,3,...,5) forest type in the j-th (j = 1,2,3,...,20) district (county) in this year; NPPmean is the area-weighted average NPP value (t hm −2 year −1 ) of the forest in the reservoir area in the respective year.
The carbon stocks of forest vegetation, litter, and soil can be calculated by the carbon density values of the vegetation layer, the litter layer, and the soil layer obtained by the Biome-BGC and the corresponding area of various forest types for the district (county) in the respective year.Subsequently, the sum of the above three carbon stocks represents the total forest carbon stock.

Interannual Variation in NPP for Different Forest Types
The annual average values of NPP for various forest types in the reservoir area from 1981 to 2014 are shown in Figure 3.The simulation results indicated that the interannual variation of NPP was drastic, presenting five consecutive and alternatively fluctuating curves according to the different forest types.The higher and lower values exhibited consistent variation trends.Coniferous forests accounted for the largest area in the reservoir region and were mainly composed of Pinus massoniana forest, cedar forest, and Chinese fir forest, accounting for 35%-40% of the total forest area.The average NPP value was 550.29 g C m −2 year −1 , ranging between 262.93 and 807.71 g C m −2 year −1 , with a value difference of 544.78 g C m −2 year −1 between the maximum and the minimum values.Mixed coniferous and broadleaf forests were mainly dominated by Pinus massoniana, cedar, Chinese fir, Quercus variabilis, and oriental white oak, accounting for 15%-18% of the total forest area.The average NPP value was 504.49g C m −2 year −1 , varying between 307.02 and 654.63 g C m −2 year −1 .The deciduous broadleaf forest in the reservoir region is severely disturbed by human activities and presented a stripped distribution pattern.It accounted for 11% to 14% of the total forest area.The average NPP value was 458.69 g C m −2 year −1 , ranging between 351.12 and 531.62 g C m −2 year −1 .The evergreen broadleaf forest represented the zonal vegetation in the reservoir area, but excessive exploitation had led to serious disturbance and destruction of the original vegetation.The area with natural forest was extremely limited and mainly consists of secondary forests, accounting for only about 5% of the total forest area.The average NPP value was 456.40 g C m −2 year −1 , ranging between 248.69 and 607.67 g C m −2 year −1 .The shrub forest is dominated by Cotinus coggygria, Rhus javanica L., Coriaria nepalensis and Lespedeza bicolor, and its area accounted for about 30% of the total forest area [46].The average NPP value was 274.99 g C m −2 year −1 , ranging between 160.66 and 484.07 g C m −2 year −1 .From 1992 to 2012, the NPP contribution ratio of each forest type to the total forest NPP was relatively stable.To be specific, the NPP of the coniferous forest accounted for the largest proportion of the total NPP, reaching 43.9%-49.9%.The NPP of the mixed coniferous and the broadleaf forest accounted for the second largest proportion of the total NPP, namely, 16.4%-20.6%,while the proportions of shrub forest and deciduous broadleaf forest were similar, namely, 15.4%-18.4% and 13%-15.4%,respectively.The contribution of the evergreen broadleaf forest was the lowest and accounted for only 3.2%-4.7% of the total NPP.The NPP variation trends of different forest types were relatively consistent and similar to that of the total forest NPP.The period from 1992 to 1996 showed rapid forest growth, with significantly increasing NPP values for all forest types.From 1996 to 2002, the NPP values for all forest types basically showed an increasing trend.In 2006, the NPP values for all forest types significantly decreased.In 2012, apart from the increase in the NPP of deciduous broadleaf forest by 57.48% compared to that in 2006, NPP values of all forest types increased by more than 100% to 7.067 Tg C.

Forest Carbon Sequestration
The variations in the carbon sequestration of the forest vegetation layer, the litter layer, and the soil layer as well as the variations in carbon storage and carbon density of the entire area from 1992 to 2012 are shown in Table 6.
From 1992 to 2012, the carbon stocks of the forest vegetation layer, the litter layer, and the soil layer as well as the entire area have similar variation trends.Specifically speaking, the stocks of each layer and the entire area continuously increased from 1992 to 2012, while the growing rate began to slow down and stabilize after the rapid increase.From 1992 to 1996, the carbon stocks of the forest vegetation layer, the litter layer, and the soil layer, and the whole area increased by 2.04%-2.16%from 54.53 Tg, 6.678 Tg, 139.637Tg and 200.845Tg in 1992, respectively, to 55.707 Tg, 6.814 Tg, 142.607 Tg, and 204.128Tg in 1996.From 1996 to 2002, the growth rate was significantly higher, and the carbon stocks for each layer and the total area increased by 23.63%, 23.53%, 21.17%, and 21.92% to 68.872 Tg, 8.419 Tg, 172.795Tg, and 250.084Tg, respectively.From 2002 to 2006, the total forest carbon stock kept growing, although the growth rated decreased to 16.27%-17.46%;the carbon stocks of the three layers and the entire area were 80.899 Tg, 9.884 Tg, 200.916Tg, and 291.699Tg.The variation trends remained the same from 2006 to 2012, increased by 11.39%-13.13% to 90.117 T g, 11.016 Tg, 227.288Tg, and 328.421Tg, respectively, in 2012.

Spatial Evolution of Forest Productivity and Carbon Sequestration in the TGRA
As shown in Figure 4, regions with higher forest NPP values from 1992 to 2012 were mainly concentrated in the four counties in western Hubei in the eastern area of the TGRAT, and forest productivity gradually increased from west to east.The most significant change occurred in the counties Wuxi, Wushan, and Yunyang in the eastern hinterland.From 1992 to 1996 and further to 2002, forest NPP rose.Up to 2012, the forest productivities of Wulong County and Fengdu County, located in the Karst Mountains, also significantly improved from 80.046 GgT and 52.313 Gg T in 1992 to 179.673 Gg T and 118.145Gg T in 2012, respectively.The spatial variation of forest productivity was consistent with that of the forest and the fragmented forest patches.Therefore, increases in forest area and optimization of the forest landscape pattern can promote the improvement of forest productivity.As shown in Figure 5, the forest carbon stock in this area also gradually increased from west to east, and the regions with significant carbon stock variation were Wuxi County, Wushan County, and Yunyang County in the eastern hinterland, Wulong County, Fengdu County, and Shizhu County in the western hinterland, and Jiangjin District and Banan District in the western area of the reservoir area.From 1992 to 1996, the forest carbon stock in the eastern hinterland of the reservoir area increased significantly with continuously high values.From 2002 to 2006, the forest carbon stocks of the Jiangjin District and the Banan District, in the hilly areas at the paralleled ridge valley of the western part, increased significantly; while from 1992 to 2012, the forest carbon stocks of Wulong County, Fengdu County, and Shizhu County, in the southwestern hinterland of the reservoir area, increased gradually; this area became the green barrier of the southern reservoir area.Generally, the air temperature in the TGRA varied slightly, while precipitation presented an increasing trend from 1992 to 2002.With this influence, the unitary forest NPP in the reservoir area firstly increased significantly and then slightly decreased.Pearson´s correlation analysis of the NPP values of various forest types and the annual average temperature showed a significantly negative correlation between these factors (Table 7).Besides, Pearson´s correlation analysis on the NPP values and the annual precipitation showed a significantly positive correlation (Table 8).As a result of the high temperature and low precipitation in 2006, forest productivity for unit area was extremely low.We therefore assume that unitary forest productivity is greatly affected by the environmental factors such as solar radiation and precipitation in addition to the physiological and ecological characteristics of the plants themselves, whereas total forest productivity is largely influenced by human activities.

Influences of Land Use Change on Forest Carbon Stocks.
From 1992 to 2012, the areas of various landscape types changed differently [47].During these 20 years, crop land and forest land were the dominating land use types in this reservoir area, accounting for 93%−96% of the reservoir area.The area of crop land greatly decreased, and its decrease rate was accelerated with the progress of the Three Gorges project.Crop land was mainly transferred to construction land, forest land, water bodies, and the areas of forest land, construction land, or water increased over the years.Among these, construction land showed the most significant increase, increasing five-fold from 1992 to 2012, while water areas in 2012 were two times larger than that in 1992.Forest land increased by nearly 60% or 131.73 × 10 4 hm 2 compared to 1992.With the drastic change in land use since the implementation of reservoir project, the continuous increase in total NPP could be attributed to the significant increase in forest area.In addition, the average area of forest patches increased year by year from 1992 to 2012, while patch density decreased.Hereinto, the landscape dominance index of the coniferous forest increased significantly, and the forest landscape pattern changed from small patches to large patches and from a spatially scattered distribution to a continuous distribution, finally enhancing the internal and external connective of the landscape.
Accordingly, from 1992 to 2012, the carbon stocks of the forest vegetation layer, the litter layer, the soil layer, and the entire reservoir area, showed a continuously increasing trend, but the growth rates started to decrease and finally stabilized.The possible reasons could be attributed to the new forest plantation, forest landscape restoration caused by a number of important forestry ecological projects, including the Yangtze River Basin Shelter Forest System project, the Natural Forest Protection Project, and the Grain for Green Project, which were implemented after the severe floods in 1998.From 1992 to 2002, the total forest area in the reservoir region increased from 204.17 × 10 4 hm 2 to 254.93 × 10 4 hm 2 .

Conclusions
The Three Gorges Reservoir area is one of the most ecologically sensitive areas in China, and the forest landscape pattern in this region shows dramatic change due to the influence of the Three Gorges reservoir project.The forest NPP and carbon sequestration of various forest types in this area, which were estimated based on the Biome-BGC model, indicated that total and unitary forest NPP presented obvious annual fluctuation from 1992 to 2012 under the combined influences of land use change and extreme weather events.Apart from the year 2006, from 1992 to 2012, the NPP values of each forest type showed an increasing trend, although the growth rates decreased.In 2006, due to abnormally high air temperatures and less precipitation, total and unit area forest NPP values decreased by 46.3% and 53.9%, respectively, compared to 2002.From 1992 to 2012, the carbon stocks of the forest vegetation layer, the litter layer, the soil layer, and the entire area gradually increased with a decreasing growth rates.Forest carbon stocks were high in the east and the south and low in the west and the north.Regions with higher carbon stocks included the reservoir head, the Wushan-Wuxi section, the Shizhu-Wulong section, and the southern section of the Jiangjin District.From 1992 to 2012, the variation in the forest landscape pattern and the responses of forest productivity and carbon stocks indicated that the increase in forest area and the optimization of the forest landscape pattern could improve the forest productivity and carbon stocks.
In fact, the potential and historic land use changes influence not only the forest areas, but also the forest landscape patterns, such as patch size, fragmentation degree, and patches connectivity and so on.This paper only analyzed the changes in forest productivity and carbon storage caused by the changes of forest area, but did not deeply analyze the impacts of landscapes patterns.Therefore, it is an important issue for further study to clarify the contributions of landscape pattern to carbon storage based on the principle of edge effects.
Based on the Biome-BGC model, the simulated NPP values and carbon density values of the coniferous forest and the mixed coniferous and broadleaf forest were close to the measured values, while those of the deciduous broadleaf forest, the evergreen broadleaf forest, and the shrub forest were lower than the measured values.The possible reason could be the use of the default physiological and ecological parameters provided by the model, as well as the underestimate of the carbon density of forest soil.The effective approach to improve the precision of the Biome-BGC model is to provide more accurate spatial interpolation through the plot investigation data of meteorological stations and more physiological and ecological parameters of different vegetation types (e.g., bamboo forest, forest plantation, etc.).By combination of field plot measurement and model simulation, the conversion from point to region and a multi-scale simulation will be realized; also, the accuracy of scaling up or scaling down will be improved.

Figure 1 .
Figure 1.Location and administrative districts of the TGRA.

Figure 2 .
Figure 2. The distribution map of forest types in the TGRA in 2012.

Figure 3 .
Figure 3.Time series of NPP for various forests in TGRA from 1981 to 2014.

3. 2
.1.Forest Productivity The total forest NPP and the unit area NPP as well as the variation situations of NPP for various forest types during 1992-2012 are shown in Table 5.Both the total and unit area forest NPP in the TGRA presented an obvious interannual fluctuation.In 1992, the values were 7.827 Tg C and 3.833 Mg C hm −2 year −1 , respectively, while in 1996, total and unit area forest NPP values showed similar variations and increased by 44.8 and 41.8% compared to the values in 1992, reaching 11.338 Tg C and 5.437 Mg C hm −2 year −1 .In 2002, total forest NPP continued to increase, albeit with a declining growth rate, while it increased by 17.7% in 2002 compared with the value in 1996.In contrast, the unit area NPP value decreased slightly in 2002 by 3.72% compared with that in 1996.In 2006, total and unit area NPP values decreased by 46.3 and 53.9% compared with the values in 2002, reaching the minimum values of 7.171 Tg C and 2.415 Mg C hm −2 year −1 , respectively.In 2012, total forest NPP reached the maximum value of 14.999 Tg C, while the unit area NPP also increased to 4.465 Mg C hm −2 year −1 , which was 109.2 and 84.9% higher, respectively, compared to the values in 2006.

Figure 4 .
Figure 4. Spatial changes in forest NPP in the TGRA from 1992 to 2012.

Figure 5 .
Figure 5. Spatial changes in forest carbon storage in the TGRA from 1992 to 2012.

1 .
Influences of Air Temperature and Precipitation on Forest ProductivityForest productivity is not only controlled by the physiological and ecological characteristics of the plants themselves, but is also affected by solar radiation, precipitation, soil, and human activities.From 1981 to 2014, the annual average temperature in the TGRA showed an upward trend in the fluctuating process.During this period, the annual average temperature was 18.49 °C, with slight annual variations; a minimum temperature of 17.75 °C in 1989 and a maximum temperature of 19.59 °C in 2006.In 2006, the average temperature in China was abnormally high.For most regions within this area, the number of days with a maximum temperature of 35 °C or more was 20 times higher than that in the same period of previous years.In addition, the annual precipitation during 1981 to 2014 in the TGRA showed a larger interannual fluctuation and a slightly downward trend.Average annual precipitation was 1110.34 mm, and the interannual fluctuation was relatively larger, with a maximum annual precipitation of 1415.71mm (1998) and a minimum annual precipitation of 886.20 mm (2006).Maximum precipitation occurred in 1998, when the Yangtze River Basin suffered a catastrophic flood event; while in 2006, the high temperature and drought event in the Sichuan Province and the Chongqing City led to the decrease of annual precipitation by 259.6 mm compared with that in the same period of "normal" years, and annual precipitation in 2006 was the lowest since 1951 (Figure6).

Figure 6 .
Figure 6.Interannual variation of annual precipitation in Three Gorges Reservoir Area in 1981-2014.

Table 1 .
The environmental information of the counties in the TGRA.

Table 2 .
Value of ecophysiological parameters of forests used in Biome-BGC.

Table 3 .
Comparison between measured and simulated NPP values in the TGRA.

Table 4 .
Comparison of estimation results of carbon density for the main forest types in the TGRA.

Table 5 .
NPP of Forest Ecosystem in Three Gorges Reservoir Area in 1992~2012.

Table 6 .
Carbon storage of the Forest Ecosystem in TGR area from 1992 to 2012.

Table 7 .
Pearson correlation analysis between NPP of forest types and annual mean temperature in TGR.

Table 8 .
Pearson correlation analysis between NPP of five forest types and annual precipitation in TGR.