Quantitative Assessment of the Human Appropriation of Net Primary Production (hanpp) in the Coastal Areas of Jiangsu, China

Global increases in population and consumption have raised concerns regarding the sustainability of the current and future use of natural resources. The human appropriation of net primary production (HANPP) provides a useful measure for determining human-derived alterations in the amount of biomass that is available in ecosystems each year. HANPP was calculated based on vegetation modelling, agricultural statistics, and remote sensing data on land use and land cover to assess the human impacts on ecosystems in the coastal areas of Jiangsu, China. The results showed that HANPP increased from 332 g¨C/m 2 /year in 2000 to 442 g¨C/m 2 /year in 2010, with an average annual increase of 2.9%. The proportion of appropriated net primary production increased from 50.3% to 71.0% of NPP pot , mainly driven by HANPP harv (harvested NPP) with an increase from 45.2% to 61.3% of NPP pot. Additionally, the spatial variation in average HANPP was striking among counties in the observed period with the lowest and highest values of 21.8% and 63.8% of NPP pot , respectively. Further analysis showed that observed levels of HANPP are high due to a high level of biomass harvest from cropland and the increases in fertilizer use, farmland irrigation rate and population and economic growth explain the trends in HANPP in the coastal areas of Jiangsu.


Introduction
Ecosystems are essential for human survival and development because they provide many types of services required by humans, e.g., food, fuel, and fiber.With the rapid development of economies and societies, huge impacts have been observed on the balance of ecosystems because of the over-consumption of natural resources [1].Human appropriation of net primary production (HANPP) is the amount of net primary production (NPP) appropriated by humans in a specific area and is extensively used as an ecological indicator to measure human interventions in the biosphere.Extensive attention has been given to HANPP, which is now used as an indicator that links natural processes with socio-economic processes, and it is used to measure the limitations of ecological environments in relation to the population and economy [2][3][4].However, most studies of HANPP have been conducted at global and national scales, and smaller scale studies are scarce.Smaller scale studies are necessary to better understand how HANPP can be used for land use policy and to demonstrate an important extension of HANPP.
The coastal areas of Jiangsu in China consist of three prefecture-level cities, namely Liangyungang, Yancheng and Nantong, which consists of three districts and 19 counties (Figure 1).The longitude and latitude of this region range from 118 ˝24 1 E to 121 ˝01 1 E and from 31 ˝38 1 N to 35 ˝08 1 N, with a total land area of 32,500 km 2 .This region is located at the junction of three productive areas: the China Coast, the Yangtze River and the Longhai-Lanxin Railway, which is also the head of the new Eurasian continental bridge.The climate in this area is classified as an intersection of the warm temperate zone and the northern subtropical zone, with an annual temperature of approximately 13-15 ˝C and annual rainfall of approximately 850-1080 mm.The vegetation in this area is dominated by cultivated crops, and the natural vegetation, including evergreen coniferous forest, evergreen broadleaved forest, and coastal salt vegetation, is distributed sparsely.According to the data of the Second Land Census of China, the percentages of cropland, forestland, construction land, water areas (including tidal flats), and unused land are 55.4%, 0.8%, 18.5%, 23.0%, and 2.3%, respectively.
Sustainability 2015, 7, 1-14 2 studies are necessary to better understand how HANPP can be used for land use policy and to demonstrate an important extension of HANPP.
The coastal areas of Jiangsu in China consist of three prefecture-level cities, namely Liangyungang, Yancheng and Nantong, which consists of three districts and 19 counties (Figure 1).The longitude and latitude of this region range from 118°24′ E to 121°01′ E and from 31°38′ N to 35°08′ N, with a total land area of 32,500 km 2 .This region is located at the junction of three productive areas: the China Coast, the Yangtze River and the Longhai-Lanxin Railway, which is also the head of the new Eurasian continental bridge.The climate in this area is classified as an intersection of the warm temperate zone and the northern subtropical zone, with an annual temperature of approximately 13-15 °C and annual rainfall of approximately 850-1080 mm.The vegetation in this area is dominated by cultivated crops, and the natural vegetation, including evergreen coniferous forest, evergreen broadleaved forest, and coastal salt vegetation, is distributed sparsely.According to the data of the Second Land Census of China, the percentages of cropland, forestland, construction land, water areas (including tidal flats), and unused land are 55.4%, 0.8%, 18.5%, 23.0%, and 2.3%, respectively.This area is a region of significant land use change; it is greatly affected by human activities and is also an ecologically fragile region where wetland resources and natural area reserves are densely distributed [5].Driven by the coastal development strategy in China, the ever-increasing use of ecosystems by humans has increased the pressure on the ecological environments in the coastal areas of Jiangsu [6,7].HANPP studies in China have rarely been reported [8] and are even lacking for the China Coast.The aims of this study were to estimate HANPP in the coastal areas of Jiangsu, to explore the temporal and spatial variation in HANPP, and to determine its driving forces in this area from 2000 to 2010.This study will be of value for revealing the intensity of human activities in terrestrial ecosystems and for providing important insights into sustainable development in the coastal areas of Jiangsu.Furthermore, this study is an attempt to extend HANPP research from global and national scales to a smaller scale.This area is a region of significant land use change; it is greatly affected by human activities and is also an ecologically fragile region where wetland resources and natural area reserves are densely distributed [5].Driven by the coastal development strategy in China, the ever-increasing use of ecosystems by humans has increased the pressure on the ecological environments in the coastal areas of Jiangsu [6,7].HANPP studies in China have rarely been reported [8] and are even lacking for the China Coast.The aims of this study were to estimate HANPP in the coastal areas of Jiangsu, to explore the temporal and spatial variation in HANPP, and to determine its driving forces in this area from 2000 to 2010.This study will be of value for revealing the intensity of human activities in terrestrial ecosystems and for providing important insights into sustainable development in the coastal areas of Jiangsu.Furthermore, this study is an attempt to extend HANPP research from global and national scales to a smaller scale.

Definition and Calculation of HANPP
Different studies use different definitions of HANPP, and these differences largely explain the range of reported results.Vitousek et al. calculated three estimated levels (low, intermediate, and high) of global HANPP according to the NPP used directly by humans and livestock, the entire NPP of human-dominated ecosystems and they further included NPP losses through land-conversions (by including a potential NPP estimate) [9].Imhoff et al. estimated the appropriation of NPP during harvesting and processing of forests and agriculture without accounting for the NPP losses through land use changes calculated by subtracting actual NPP from potential NPP [10].Haberl et al. [11] defined HANPP as the difference between the NPP of potential natural vegetation and the NPP remaining in ecosystems, and stated that HANPP is affected by two separate processes: (1) biomass harvesting and (2) land use and land cover changes.The Haberl method has become the standard for HANPP studies [12], and we chose to use the definition of HANPP outlined in Haberl, et al. [13].Accordingly, HANPP can be calculated as follows.
HANPP " HANPP luc `HANPP harv , HANPP luc " NPP pot ´NPP act (1) To calculate HANPP, some variables must first be determined, which include the NPP of the natural vegetation (NPP pot ), the NPP of the currently prevailing vegetation (NPP act ), and the harvested NPP (HANPP harv ).

Estimation of HANPP harv
More than one-half of the regions in the coastal areas of Jiangsu are covered by cultivated land, and the majority of harvests are related to agricultural activities.The wood harvest of the forestry industry was not considered in our study because the timber harvest data were unavailable and the proportion of the forestland in the study area was only approximately 0.8%.Thus, HANPP harv in our calculation focuses on biomass harvest from cropland.
We used crop yield data from the Jiangsu Agriculture Database (Statistics Bureau of Lianyungang, Yancheng and Nantong) to estimate the biomass harvested from croplands.Agricultural statistics were obtained for the period between 2000 and 2010.To analyze the spatial variability of HANPP, all harvest data were tabulated annually on a county scale (n = 21).Generally, only the edible portion of a crop was recorded in the agricultural harvest statistics.To incorporate information on the harvest residues, the reported production amounts were divided by the crop-specific ratios of the residue to economic yield, which were obtained from the literature [14][15][16][17].The harvest data were converted from fresh weight to dry weight using a literature-derived fraction of dry matter [18,19].
To estimate the agricultural harvest, production data were analyzed for 12 crops grown in the coastal areas of Jiangsu.The studied crops were wheat, barley, broad bean, pea, rice, maize, soybean, potato, peanut, rapeseed, sesame, and cotton and they represent an average of 82.7% of the total harvested acreage of the crops in the coastal areas of Jiangsu.Since the remaining cropped areas (17.3%) were primarily planted with vegetables for which related data sets were unavailable, HANPP estimates were not performed for these areas.The amount of HANPP harv in the coastal areas of Jiangsu was calculated by summing C y and C r as follows: Sustainability 2015, 7, 15857-15870 where C y is the total amount of NPP in the economic yield, C r is the total amount of NPP in the crop residues, Y i is the economic yield, F id is the dry matter fraction, F icy and F icr are the carbon fractions of the economic yield and residues, respectively, and R iry is the ratio of the residue to the economic yield.The details of these parameters are listed in Table 1.HANPP harv in our calculation includes the crop harvest and its residues (the straws and stover of wheat, rice, cotton, etc., and the leaves of potatoes and peanuts), whereas belowground NPP (with the exception of potatoes and peanuts), NPP losses during the growth period, and losses resulting from herbivory are not considered in the study.All of the results of this analysis are expressed in g C per year.

Estimation of NPP pot and NPP act
NPP pot was generally derived by the regression analysis of NPP data with climate or other factors.As an empirical model, the Miami model [20] determines NPP for a particular location as the minimum of the temperature and precipitation regression functions.It is widely used for its simplicity and relative accuracy and sometimes used as a baseline for global NPP evaluation [21][22][23][24][25].
The parameters of the Miami model were calibrated by fitting the NPP observations with temperature and precipitation data in the early years before the 1970s and must be re-parameterized with new sets of productivity and climate change observations [25].A new set of parameters was calibrated and optimized by Zaks, et al. in 2007 [25] based on the 3023 global NPP field observations and corresponding annual mean temperature and precipitation data.We chose the new calibrated Miami model to calculate NPP pot in this study.
NPP act in the coastal areas of Jiangsu from 2000 to 2010 was calculated using the Carnegie Ames Stanford Approach (CASA) model, which was proposed by Potter et al. in 1993 [26].NPP act calculated using the CASA model includes all biomass in the process of vegetation growth.The amount of NPP determined by absorbed photosynthetic active radiation (APAR) and light use efficiency (ε) was calculated as follows.
NPP px, tq =APAR px, tq ˆε px, tq APAR px, tq " SOL px, tq ˆFPAR px, tq ˆ0.5 ε px, tq " T 1 px, tq ˆT2 px, tq ˆW px, tq ˆε* (7) where t is the time; x is the spatial position; SOL(x,t) is the total amount of solar radiation (MJ/m 2 ); FPAR(x,t) is the fraction of photosynthetic active radiation, which was calculated using the method proposed by Los [27] that is based on the normalized difference vegetation index (NDVI) and the ratio vegetation index (SR); the constant 0.5 is the ratio of the amount of effective solar radiation absorbed by vegetation (wavelength range between 0.38 and 0.76 µm) to the total amount of solar radiation; T 1 (x,t) and T 2 (x,t) are the influence of temperature on the light use efficiency (i.e., the temperature scale) and calculated following Potter and Field,et al. [26,28]; W(x,t) is the influence of water on the light use efficiency (i.e., the water scale), calculated following Piao, et al. [29]; and ε* is the maximum light use efficiency, which was based on the results presented by Zhu et al. for typical vegetation in China [30].
The models are driven by NDVI, climate, and land use and land cover data.The monthly NDVI data were obtained from the MOD13Q1 16-day composite product that was downloaded from the Level 1 and Atmosphere Archive and Distribution System (LAADS) website (http://ladsweb.nascom.nasa.gov/data/).
The climate data, including temperature, precipitation and solar radiation, were obtained by interpolation from meteorological station data that were downloaded from the China Meteorological Data Sharing Service System website (http://cdc.nmic.cn/home.do).The solar radiation data were calculated following Zhang, et al. [31], using ground water vapour pressure and sunshine percentage, and the land use and land cover data were interpreted from TM images from 2000 and 2010.All data were projected to Albers projection and resampled to a resolution of 250 m ˆ250 m.

Trend Analysis
We used the following linear regression equation to analyse the HANPP trend.
where k is the trend slope and n is the number of years.F statistics were calculated as follows to test the significance of the trend, which reflects the confidence of the trend.
F " U ˆn ´2 Q U " where U is the sum of the squared errors, Q is the sum of the squares of the regression, y i is the observed value, ŷi is the regressed value, y is the average value, and n is the number of years (n = 11).
According to the F statistics for the significance test, we divided the trend into three classes: extremely significant (p < 0.01), significant (0.01 < p < 0.05) and not significant (p > 0.05).

Correlation Analysis
The relationship between HANPP and the driving factors was analysed using the Pearson correlation coefficient: where x i and y i are the observed values of the variables; x and y are the average values of the variables; and n is the number of samples.
The population, GDP, fertilizer use, and farmland irrigation area datasets used for the correlation analysis were obtained from the statistical yearbooks of Liangyungang, Yanchen, and Nantong (Statistics Bureau of Lianyungang, Yancheng and Nantong).

NPP pot and NPP act
NPP pot changed with fluctuations from 661 g¨C/m 2 /year to 623 g¨C/m 2 /year over the period from 2000-2010, and NPP act presented a basically similar trend to NPP pot , with values changing from 627 g¨C/m 2 /year to 574 g¨C/m 2 /year (Figure S1).The spatial distribution of the average NPP pot in the observed period indicated a decreasing trend from north to south in the study area, with values ranging between 603-832 g¨C/m 2 /year.The average NPP act over the period 2000-2010 was concentrated in the interval between 600-800 g¨C/m 2 /year, which mainly occurred in the cultivated land area of the counties.There was a gradient distribution of NPP act in the tidal flat reclamation zones, where NPP act increased gradually from the sea to the land.Land use and land cover had a significant impact on the spatial distribution of NPP act (See Figure S2).

HANPP harv
HANPP harv increased from 298 g¨C/m 2 /year to 393 g¨C/m 2 /year during the period from 2000 to 2010 (Figure 2a) but declined significantly in 2003 because of the reduced food crop yields, due in large part to the large decrease in the rice yields in Lianyungang and Yancheng caused by meteorological disasters, such as abnormal floods or long-term continuous rain.The proportions of HANPP harv increased from 45.2% to 61.3% of NPP pot , which exhibited a continually increasing trend after 2007.

NPPpot and NPPact
NPPpot changed with fluctuations from 661 g•C/m 2 /year to 623 g•C/m 2 /year over the period from 2000-2010, and NPPact presented a basically similar trend to NPPpot, with values changing from 627 g•C/m 2 /year to 574 g•C/m 2 /year (Figure S1).The spatial distribution of the average NPPpot in the observed period indicated a decreasing trend from north to south in the study area, with values ranging between 603-832 g•C/m 2 /year.The average NPPact over the period 2000-2010 was concentrated in the interval between 600-800 g•C/m 2 /year, which mainly occurred in the cultivated land area of the counties.There was a gradient distribution of NPPact in the tidal flat reclamation zones, where NPPact increased gradually from the sea to the land.Land use and land cover had a significant impact on the spatial distribution of NPPact (See Figure S2).

HANPPharv
HANPPharv increased from 298 g•C/m 2 /year to 393 g•C/m 2 /year during the period from 2000 to 2010 (Figure 2a) but declined significantly in 2003 because of the reduced food crop yields, due in large part to the large decrease in the rice yields in Lianyungang and Yancheng caused by meteorological disasters, such as abnormal floods or long-term continuous rain.The proportions of HANPPharv increased from 45.2% to 61.3% of NPPpot, which exhibited a continually increasing trend after 2007.Clear differences in HANPP harv trends among the cities of Lianyungang, Yancheng, and Nantong are evident in Figure 2b-d.In Lianyungang and Yancheng, HANPP harv increased significantly (p < 0.01) in contrast to Nantong (p > 0.05).HANPP harv in Lianyungang presented the greatest increase, i.e., from 34.9% to 61.0% of NPP pot between 2000 and 2010.In contrast, HANPP harv in Nantong exhibited the smallest increase, i.e., 45.6% to 47.4% of NPP pot .
To further reveal the spatial variations in HANPP harv trends, we performed linear regression analysis of the HANPP harv time series for the counties in the coastal areas of Jiangsu to obtain the slope (k) of HANPP harv for the period from 2000 and 2010, where k > 0 indicates an increasing and k < 0 indicates a decreasing HANPP harv trend.According to the F statistic for the significance test of the trend, we divided the trend into three classes: extremely significant (p < 0.01), significant (0.01 < p < 0.05) and not significant (p > 0.05).The results shown in Figure 3a illustrate that HANPP harv trends shifted from extremely significant increases to no significant trends in the counties from Lianyungang to Nantong, where the dominant crops changed from grain crops (including rice, wheat, barley, and maize) to economic crops (including peanut, rapeseed, sesame, and cotton) (Figure 3b, Table S1).The difference in the crop structure among the counties is one of the causes of the difference in HANPP harv trends.The stable increase in the yield of grain crops compared with economic crops may be another cause.To further reveal the spatial variations in HANPPharv trends, we performed linear regression analysis of the HANPPharv time series for the counties in the coastal areas of Jiangsu to obtain the slope (k) of HANPPharv for the period from 2000 and 2010, where k > 0 indicates an increasing and k < 0 indicates a decreasing HANPPharv trend.According to the F statistic for the significance test of the trend, we divided the trend into three classes: extremely significant (p < 0.01), significant (0.01 < p < 0.05) and not significant (p > 0.05).The results shown in Figure 3a illustrate that HANPPharv trends shifted from extremely significant increases to no significant trends in the counties from Lianyungang to Nantong, where the dominant crops changed from grain crops (including rice, wheat, barley, and maize) to economic crops (including peanut, rapeseed, sesame, and cotton) (Figure 3b, Table S1).The difference in the crop structure among the counties is one of the causes of the difference in HANPPharv trends.The stable increase in the yield of grain crops compared with economic crops may be another cause.

Comparison with Previous Studies
Land use classes in the study area were aggregated to five types, namely cropland, forestland, construction land, water areas, and unused land.The NPP pot of the natural waters (accounting for approximately 71% of the total water areas) including tidal flats, lakes, and Yangtze River was set to be equal to NPP act in our calculation.Although the proportion of construction land is relatively high, it did not result in a high HANPP luc proportion of NPP pot in the whole area because much of the construction land is villages (approximately 52% of the construction land), with little impervious surface and lots of trees or other vegetation.In our estimate, the average NPP pot and NPP act of construction land were 698 g¨C/m 2 /year and 463 g¨C/m 2 /year respectively for the years 2000-2010.Additionally, as stated in Section 2.2., HANPP harv in our calculation focus on biomass harvest from cropland.Thus, HANPP estimate in our study may be slightly lower than the actual value the coastal areas of Jiangsu for the period 2000-2010.
On a global scale, estimates of HANPP vary from 20% to 40% (Table 2) because of differences in the definition, method and conversion factors.On a country or regional scale, the spatial variation in HANPP is striking.Imhoff et al. generated a spatial balance sheet of the global NPP "supply" and "demand".This sheet shows that western Europe and south-central Asia appropriate more than 70% of their regional NPP, whereas South America appropriates approximately 6% of the NPP.As shown in Table 2, our results are similar to the estimates of HANPP in the United Kingdom, Philippines, Spain, Italy, Czech Republic, and China, where the proportion of appropriated NPP is greater than 55% because of the higher proportion of agricultural land.HANPP in Nova Scotia, Canada (25.5% of NPP pot ) is low because only 3%-8% of the area is covered by agricultural land.Land use is one of the important factors that influences HANPP.Generally, a high proportion of agricultural land leads to a high HANPP, which is consistent with the results obtained in the coastal areas of Jiangsu.Different HANPP trends can also be observed in countries and regions outlined in Table 2. Descending HANPP trends in the United Kingdom, Spain, and Italy, compared with increasing trends in China and the Philippines, which is consistent with our study area.The factors that influence HANPP trends will be discussed in the following sections.

Influence of Climate Conditions on HANPP
NPP pot and NPP act were estimated using the Miami model and CASA model, respectively, driven by the climate data, which exhibited annual fluctuations in the study area over the period from 2000 to 2010.Furthermore, NPP pot and NPP act varied significantly within individual years as a result of the high interannual variation in the climatic conditions.In particular, the 2007 annual mean temperature was the highest value recorded in Jiangsu province since 1950.
Calculations based on the data from 14 meteorological stations (19 stations in 2009 and 2010) and five radiation stations (Figure S3) indicated that NPP pot decreased (in 2004) or increased (in 2006 and 2007) significantly because of remarkable drought or high temperature and rainy climate, respectively, which resulted in the HANPP harv proportion of NPP pot being extremely high or low in the corresponding years (Figure 4a).Furthermore, a decrease in solar radiation caused by long-term continuous rain resulted in a decrease in NPP act which was calculated using the light use efficiency model (the CASA model).HANPP luc , which was determined as the difference between NPP pot and NPP act , increased significantly in 2006 (Figure 4a).
Additionally, HANPP luc , which accounted for a very small part of HANPP, exihibited a relatively strong sensitivity to the climate fluctuations.HANPP luc may drop below zero (in 2004 and 2008) under particular weather conditions, which result in that HANPP harv being higher than HANPP in these years.

Influence of Agricultural Production Conditions on HANPP
In an agricultural ecosystem, fertilizer use and irrigation area are important driving factors for the amount of harvested NPP [45].Here, fertilizer use per unit area (the amount of fertilizer use divided by the cropland area) and the farmland irrigation rate (the ratio of irrigation area to cropland area) were selected as indicators of the agricultural production conditions to analyze the response of HANPP harv to these factors in the coastal areas of Jiangsu.During the observed period, the fertilizer use increased from 643 kg/hm 2 to 710 kg/hm 2 , and the farmland irrigation rate increased from 74.2% to 80.6%.As shown in Figure 5, HANPP harv was significantly positively correlated with fertilizer use and farmland irrigation rate (p < 0.05).Overall, the increase in fertilizer use and irrigation rate on cropland resulted in an increase in crop yields in the study area.Generally, effective irrigation can weaken the effect of water stress and increase crop yields, and fertilizer use can increase the crop leaf nitrogen content, which accelerates the photosynthetic rate and also increases crop yields.4a).Furthermore, a decrease in solar radiation caused by long-term continuous rain resulted in a decrease in NPPact which was calculated using the light use efficiency model (the CASA model).HANPPluc, which was determined as the difference between NPPpot and NPPact, increased significantly in 2006 (Figure 4a).
Additionally, HANPPluc, which accounted for a very small part of HANPP, exihibited a relatively strong sensitivity to the climate fluctuations.HANPPluc may drop below zero (in 2004 and 2008) under particular weather conditions, which result in that HANPPharv being higher than HANPP in these years.

Influence of Agricultural Production Conditions on HANPP
In an agricultural ecosystem, fertilizer use and irrigation area are important driving factors for the amount of harvested NPP [45].Here, fertilizer use per unit area (the amount of fertilizer use divided by the cropland area) and the farmland irrigation rate (the ratio of irrigation area to cropland area) were selected as indicators of the agricultural production conditions to analyze the response of HANPPharv to these factors in the coastal areas of Jiangsu.During the observed period, the fertilizer use increased from 643 kg/hm 2 to 710 kg/hm 2 , and the farmland irrigation rate increased from 74.2% to 80.6%.As shown in Figure 5, HANPPharv was significantly positively correlated with fertilizer use and farmland irrigation rate (p < 0.05).Overall, the increase in fertilizer use and irrigation rate on cropland resulted in an increase in crop yields in the study area.Generally, effective irrigation can weaken the effect of water stress and increase crop yields, and fertilizer use can increase the crop leaf nitrogen content, which accelerates the photosynthetic rate and also increases crop yields.Similar studies on HANPP [36,38,39] indicated that agricultural intensification (including increasing fertilizer use) led to a decline in HANPP as areas of marginal productivity in agricultural land were abandoned or reforested, which is different from our results.In our study area HANPP increased with the increase in agricultural intensification because the cropland converted to construction land instead of being reforested.Similar studies on HANPP [36,38,39] indicated that agricultural intensification (including increasing fertilizer use) led to a decline in HANPP as areas of marginal productivity in agricultural land were abandoned or reforested, which is different from our results.In our study area HANPP increased with the increase in agricultural intensification because the cropland converted to construction land instead of being reforested.

Conclusions
In this study, HANPP was calculated in the coastal areas of Jiangsu from 2000 to 2010.It increased significantly over this period, and the proportion of appropriated NPP increased from

Conclusions
In this study, HANPP was calculated in the coastal areas of Jiangsu from 2000 to 2010.It increased significantly over this period, and the proportion of appropriated NPP increased from 50.3% to 71.0% of NPP pot .The HANPP harv proportion of HANPP ranged between 77.4% (in 2006) and 106.8% (in 2008) in the observed period, with a mean value of 92.7%.Similar work on HANPP indicated a relatively high HANPP, which was mainly determined by agricultural harvest.The results from our study are reasonable compared with those of previous studies and show that Haberl's method is suitable for use in the study area.Increased fertilizer use and irrigation are linked to improvements in crop yields, which result in increased biomass harvest in the study area dominated by cropland and also place further pressure on the local ecosystems.Additionally, this study showed that agricultural intensification did not lead to a decline in HANPP as cropland converted to construction land instead of being reforested in the study area, which provides essential information for agricultural policies.Further work is needed to improve the method by incorporating the local trade and by calibrating the conversion parameters according to field survey data.Finally, future studies on the relationships between HANPP trajectories at different economic development stages with the goal of achieving sustainable development will likely yield interesting results.

Figure 1 .
Figure 1.The location of the coastal areas of Jiangsu.

Figure 1 .
Figure 1.The location of the coastal areas of Jiangsu.

Figure 3 .
Figure 3.The significance of HANPPharv trends (a); and the average proportions of grain crops and economic crops (b) in the counties of the coastal areas of Jiangsu from 2000 to 2010.

Figure 3 .
Figure 3.The significance of HANPP harv trends (a); and the average proportions of grain crops and economic crops (b) in the counties of the coastal areas of Jiangsu from 2000 to 2010.

Sustainability 2015, 7 , 1 - 14 10 4 . 2 .
Influence of Climate Conditions on HANPP NPPpot and NPPact were estimated using the Miami model and CASA model, respectively, driven by the climate data, which exhibited annual fluctuations in the study area over the period from 2000 to 2010.Furthermore, NPPpot and NPPact varied significantly within individual years as a result of the high interannual variation in the climatic conditions.In particular, the 2007 annual mean temperature was the highest value recorded in Jiangsu province since 1950.Calculations based on the data from 14 meteorological stations (19 stations in 2009 and 2010) and five radiation stations (Figures S3) indicated that NPPpot decreased (in 2004) or increased (in 2006 and 2007) significantly because of remarkable drought or high temperature and rainy climate, respectively, which resulted in the HANPPharv proportion of NPPpot being extremely high or low in the corresponding years (Figure

Figure 5 .
Figure 5. Correlation between (a) HANPPharv and fertilizer use and (b) HANPPharv and farmland irrigation rate during the observed period.

Figure 5 .
Figure 5. Correlation between (a) HANPP harv and fertilizer use and (b) HANPP harv and farmland irrigation rate during the observed period.

Figure 6 .
Figure 6.Correlation between (a) HANPP and per capita GDP; (b) HANPP and population density during the observed period; (c) HANPPharv and per capita GDP, and (d) HANPPharv and population density in the 21 counties.

Figure 6 .
Figure 6.Correlation between (a) HANPP and per capita GDP; (b) HANPP and population density during the observed period; (c) HANPP harv and per capita GDP, and (d) HANPP harv and population density in the 21 counties.

Table 2 .
Comparisons between the estimates of HANPP in this paper and previous studies.