Next Article in Journal
Traceability System for Improved Utilization of Solid Biofuel from Agricultural Prunings
Next Article in Special Issue
Analysis of Blue and Green Water Consumption at the Irrigation District Scale
Previous Article in Journal
Czech Machinery Cluster and Its Role in Sustainable Development of Moravian-Silesian Enterprises during the Post-Transformation Era
Previous Article in Special Issue
Sustainability Index Evaluation of the Rainwater Harvesting System in Six US Urban Cities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Potential Climate Change Effects on the Rice Yield and Water Footprint in the Nanliujiang Catchment, China

1
State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing 100038, China
2
College of Civil Engineering and Architecture, Guangxi University, Nanning 530004, China
*
Authors to whom correspondence should be addressed.
Sustainability 2018, 10(2), 242; https://0-doi-org.brum.beds.ac.uk/10.3390/su10020242
Submission received: 24 November 2017 / Revised: 6 January 2018 / Accepted: 8 January 2018 / Published: 23 January 2018
(This article belongs to the Special Issue Impacts of Climate Change on Hydrology, Water Quality and Ecology)

Abstract

:
The Nanliujiang catchment is one of major rice production bases of South China. Irrigation districts play an important role in rice production which requires a large quantity of water. There are potential risks on future climate change in response to rice production, agricultural irrigation water use and pollution control locally. The SWAT model was used to quantify the yield and water footprint (WF) of rice in this catchment. A combined method of automatic and manual sub-basin delineation was used for the model setup in this work to reflect the differences between irrigation districts in yield and water use of rice. We validated our simulations against observed leaf area index, biomass and yield of rice, evapotranspiration and runoff. The outputs of three GCMs (GFDL-ESM2M, IPSL-CM5A-LR and HadGEM2-ES) under three RCPs (RCP2.6, 4.5, 8.5) were fed to the SWAT model. The results showed that: (a) the SWAT model is an ideal tool to simulate rice development as well as hydrology; (b) there would be increases in rice yield ranged from +1.4 to +10.6% under climate projections of GFDL-ESM2M and IPSL-CM5A-LR but slight decreases ranged from −3.5 to −0.8% under that of HadGEM2-ES; (c) the yield and WFs of rice displayed clear differences in the catchment, with a characteristic that high in the south and low in the north, mainly due to the differences in climatic conditions, soil quality and fertilization amount; (d) there would be a decrease by 45.5% in blue WF with an increase by 88.1% in green WF, which could provide favorable conditions to enlarge irrigated areas and take technical measures for improving green water use efficiency of irrigation districts; (e) a clear rise in future grey WF would present enormous challenges for the protection of water resources and environmental pollution control in this catchment. So it should be to improved nutrient management strategies for the agricultural non-point source pollution control in irrigation districts, especially for the Hongchaojiang and Hepu irrigation districts.

1. Introduction

Over the past decades, the global climate has undergone significant changes with global warming [1]. Agriculture is the largest freshwater consumer, accounting for more than 70% of the world’s freshwater [2] and is extremely sensitive to climate change [3]. The change in the global climate will undoubtedly influence the future freshwater supply and food security by exerting pressure on agriculture and water resources, which presents a significant challenge to society [4]. It is necessary to research the impact of climate change on the utilization efficiency of agricultural water resources to mitigate the adverse impacts on agricultural production.
To assess water resource utilization during agricultural production processes, the water footprint (WF) concept has been proposed [5]. The calculation of a water footprint of a product can include the total amount of freshwater used along the supply chain of a product [6]. In the agricultural sector, WF of crop production has three components: (1) green WF: evaporation of water supplied from the rain during crop production; (2) blue WF: evaporation of the irrigation water supplied from groundwater or surface water sources; and (3) grey WF: volume of water required to dilute the pollution during the production process [7]. In recent years, various studies have focused mainly on the assessment of the water footprint in agricultural production on different scales (global, national, basin and district) from the perspective of production, climate change impact, or both [8,9,10].
Crop growth is a complex phenomenon that is closely related to water requirements and normally influenced by climate change. Rising temperature, elevated CO2 and changing precipitation patterns are important aspects of the changing climate with pronounced impacts on crops in general [11]. (1) Effect of CO2 increase on crop photosynthesis and respiration. Research has shown that the rates of crop photosynthesis and respiration increase along with the increment of the atmospheric CO2 concentration as well as crop yield. Unlike C4 crops (maize, sorghum, sugarcane), species of type C3 (wheat, rice, soybean) have lower rates of respiration, and the efficiency of photosynthesis and its metabolic pathway react positively to high levels of CO2 [8]. It had been reported that doubling the ambient CO2 levels stimulates the biomass yield of C3 plants by 40%, and data for C4 plants indicate a stimulation of 11% [12]; (2) Effect of temperature change. Temperature controls the rate of crop metabolic processes that ultimately influence transpiration and grain production [13]. And high temperatures may increase the length of the growing period. It has been revealed that wild oat grown at a high temperature of 23/19 °C (day/night) complete their development faster than those grown at a normal temperature 20/16 °C [14]; (3) Effect of precipitation change. Crops show a particular sensitivity to the water supply. Decreased precipitation may lead to water stress, reducing the rain-fed crop output and increasing irrigation requirements in irrigated crops. However, the crop response to the interaction of CO2, temperature and precipitation is complex. For instance, C3 crops usually benefit from rising CO2 levels physiologically, but rising temperatures can override the stimulating effects of CO2 on the photosynthesis of C3 crops [11]. Some studies have shown that high temperatures reduce or eliminate the high CO2 growth enhancement [14].
Rice is one of the largest water consumers, utilizing 25% of the global agricultural freshwater supply [7]. Approximately 75% of rice production is obtained from irrigated rice [15]. Large irrigation projects are often constructed to meet the water demand in rice production. Additionally, rice is the third largest food crop planted, feeding more than 3 billion people worldwide [16]. According to the latest statistics reported by the USDA (U.S. Department of Agriculture), the world’s rice-sown area expanded to more than 160 million hectares by 2016, and the rice harvest exceeded 484 million metric tons [17]. Rice production has been chiefly concentrated in East, South and Southeast Asia, with a proportional planting area of 88.4% and yields of 89.3% yields worldwide. China has a rice sown area of 30.16 million hectares with a total rice output of approximately 145 million metric tons in 2016, making it the largest rice producer in the world [17]. Approximately 52% of the global rice harvested areas are significantly influenced by climate, including China, India and Indonesia with approximately 75% of the global rice production [18].
The tools used for modeling yield and WF of crops are expected to have process descriptions defining the crop growth and the stresses (water, temperature, and nutrient) associated with crop growth. Various crop models have been used to assess climate change in response to crop production and water use efficiency but with mixed results, such as the APSIM [19] model, the CropSyst [20] model, the DSSAT [21] model and the EPIC [22] model. However, these field-scale models were usually only used at specific sites, but limited to study regional scientific questions. The Soil and Water Assessment Tool (SWAT) is based on the inclusion of a crop growth module (EPIC) into a physically-based semi-distributed hydrologic model, which has been widely used to assess the impact of prospective climate change on the water resources and agricultural production of various basins across the world [23,24].
There are many works regarding yield and WF of upland crops production focused on dry-farming areas [25,26,27], but few quantitatively studies have been done on the yield and WF of rice production in South China, our understanding on climate change in response to rice production and water resources utilization in this area is still insufficient. In this study, the SWAT model is used to assess the yield and water footprint of rice in the Nanliujiang Catchment. The Catchment is an important double-cropping rice production region and makes an important contribution to grain consumed within South China. Irrigation districts play an important role in rice production, and cropping therein requires large amounts of water for irrigation, which is potentially at stake under future water stress. Furthermore, the provincial agricultural department plans to continue enlarging irrigated areas of rice paddy in the next 15 years. Therefore, the question arises whether future climate change may be potentially harmful to rice production, agricultural irrigation water management as well as pollution control locally. In addition, one needs a tool that has good adaptability to simulate crop production including fertilization and irrigation, and assess crop yield and water use under future climate change. Consequently, the aims of this study were (1) to evaluate performance of the SWAT model; (2) to evaluate the future climate of the Nanliujiang catchment; (3) to assess rice yield under current and future climate conditions; and (4) to investigate quantitatively the WF of rice production under current and future climate conditions.

2. Materials and Methods

2.1. Study Area

The Nanliujiang catchment (21°35′54″–22°52′32″N, 109°00′03″–110°23′12″E) is located in the south of Guangxi province, China. It borders the Darong Mountain on the North and the North Bay on the South. The 9569 km2 catchment has a unique leaf-shaped physiographic characteristic with a drop of more than 1250 m, high-lying northeast, and low southwest (Figure 1). The Nanliujiang catchment belongs to the subtropics monsoon climate region, with an average annual temperature of 22 °C. The winter is mild with less rainfall, and the summer is hot and humid with abundant rainfall. The annual precipitation varies from 1596 mm in the northern mountains to 1804 mm on the southern coast. Having a proper soil and climate conditions, the catchment becomes one of the most important grain production bases in South China. The principal crops in this area are double-cropping rice. Early rice is usually planted in mid–March, until harvest in early July. Late rice is planted in mid–July and harvested around mid–November. The rainfall is heaving during the rainy season, about two–thirds of which falls during May–September; the maximum runoff usually occurs in July and August. However, the precipitation is less than the water requirement for rice growing in spring and autumn. Thus, irrigation is highly necessary for triggering rice growth and producing high yields. At present, the total cultivated land areas were expanded to 269.8 thousand hectares, approximately 70% of which were cultivated with rice (Figure 2a). There are 18 irrigating areas with a total area of 179.7 thousand hectares (Figure 2b) [28].

2.2. Description of the SWAT Model

The SWAT model is a physically semi-distributed hydrologic model, which has been widely used worldwide because of its powerful functionality for simulating water yield, sediment, water quality in basins [29,30], as well as predicting the impact of human activity and climate change on hydrology [31].
In the SWAT, the river basin is divided into subbasins, and each subbasin is further divided into several Hydrologic Response Units (HRUs). The hydrological component of SWAT calculates the daily water balance for each HRU and takes into account precipitation, evapotranspiration, soil moisture, surface runoff and percolation. The water balance equation is represented as
S W t = S W 0 + i = 1 T ( P d a y R s u r f E T W s e e p Q g w )
where SW is the soil water content; i is the time T (days) for the simulation period, Pday, Rsurf, ET, Wseep and Qgw are daily precipitation, surface runoff, evaporation, percolation and return flow (mm), respectively.
The plant growth module that is incorporated into SWAT is a biogeochemical-based EPIC model [32,33]. The daily heat unit accumulation is used to regulate plant growth. Plant growth and crop yield are constrained by temperature, water, atmospheric CO2 concentration and fertilizer application.
Plant biomass production is modeled by simulating leaf area development, radiation interception and conversion of intercepted light into dry biomass with the radiation-use efficiency of a specified plant. The equations are expressed as follows:
b i o = i = 1 d { R U E i 0.5 H d a y [ 1 exp ( k L A I ) ] }
where bio is the total dry biomass production in the growing season (kg·ha−1); RUEi is the radiation use efficiency on day i, kg·ha−1·(MJ·m−2)−1; Hday is the incident total solar (MJ·m−2); k is the light extinction coefficient; and LAI is the leaf area index.
The total leaf area index is calculated as follows:
L A I = { i = 1 d ( f r L A Im x , i f r L A Im x , i 1 ) L A I m x { 1 exp   [ 5   ( L A I i 1 L A I m x ) ] } f r P H U f r P H U , s e n L A I m x 1 f r P H U 1 f r P H U , s e n f r P H U > f r P H U , s e n
where frLAImx,i and frLAImx,i−1 are the fraction of the plant’s maximum leaf area index for day i and i − 1, respectively; LAImx is the maximum leaf area index of the plant; frPHU is the fraction of potential heat units accumulated for the plant on a given day in the growing season; frPHU,sen is the fraction of the growing season (PHU) at which senescence becomes the dominant growth process.
The crop yield is determined by the total crop biomass, and the harvest index and is calculated as follows:
Y i e l d = b i o g [ ( H I H I min ) γ γ + exp ( 6.13 0.883 γ ) + H I min ]
where Yield is the crop yield (kg·ha−1); biog is the aboveground biomass on the day of harvest (kg·ha−1); HI is the potential harvest index on the day of harvest; HImin is the minimum harvest index allowed for the plant; γ is the water deficiency factor.
The research shows that the change in the atmospheric CO2 concentrations has a large influence on the radiation use efficiency [34]. The formulas expressing the relationship between the radiation-use efficiency and CO2 concentrations have been incorporated into SWAT for climate change studies, and the equation is as follows:
R U E = 100 CO 2 CO 2 + exp ( r 1 r 2 CO 2 )
where CO2 is the concentration of carbon dioxide in the atmosphere (ppm); r1 and r2 are shape coefficients; and RUE is the radiation use efficiency, kg·ha−1⋅(MJ·m−2)−1.

2.3. Data and SWAT Model Setup

The data required in this paper include the following: (a) a DEM map provided by the Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences; (b) a land use map (Figure 2a) and (c) a soil map obtained from the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (d) daily meteorological data from 1995 to 2013 obtained from the China Meteorological Science Data Sharing Service Network, consisting of daily precipitation, maximum and minimum air temperature, relative humidity, sunshine duration and wind speed for 5 meteorological stations (Figure 1); (e) monthly pan evaporation data from 1995 to 2013 for 3 meteorological stations (Yulin, Bobai and Hepu) also obtained from the China Meteorological Science Data Sharing Service Network; (f) monthly measured streamflow data of three hydrometric stations (Hengjiang, Bobai and Changle) from 1995 to 2013; (g) the monthly outflow data of the 3 large reservoirs (Xiaojiang, Hongchaojiang and Wangshengjiang) during 1990–2013 provided by the Guangxi water conservancy and Electric Power Survey; (h) observed data of biomass and leaf area index for primary crops grown in the Hepu irrigation district from 2005 to 2007 provided by the Guangxi Academy of Agricultural Sciences, and rice yield data from 1999 to 2012 provided by the China Meteorological Science Data Sharing Service Network; (i) data for climate change scenarios from the output of three global climate models (Table 1) provided by the Inter-Sectoral Impact Model Inter-comparison Project (ISI-MIP), which had been linearly interpolated and verified by statistical bias the probability distribution [35]. CO2 was taken as a constant for each RCP (RCP2.6, CO2 = 443, RCP4.5, CO2 = 487 ppm, RCP8.5, CO2 = 540 ppm, against CO2 = 370 ppm currently) [8,24].
In plain agricultural regions with strong disturbance by human activity, the construction of ditch system projects profoundly alter the natural water flow paths and original sub-catchments [36,37]. Irrigation return water and runoff yield are mainly drained by the drainage canal system [38], making the irrigation district become an artificial subbasin [39,40]. Consequently, plain areas should not be divided into subbasins the watershed but the irrigation boundary controlled by drainage canal systems [41]. To show the real irrigation situation as closely as possible and reflect the differences between irrigation districts in yield and water use of rice, we used a combined method of automatic and manual delineation in this study. With the traditional watershed delineation method, the Nanliujiang catchment was divided into 26 subbasins according to the watershed (Figure 2c). Regarding the subbasins which contain both mountainous and plain, using manual delineation, 19 subbasins (1–2, 4–14, 16, 19–23) were redivided into 58 subbasins the boundaries of irrigation districts. Thus, a total number of 64 subbasins were identified for the catchment (Figure 2d), and 881 HRUs were defined therein. Depending on irrigation district location, irrigated agriculture was considered in subbasins 3, 6–7 and 24–59, and rain-fed cultivated land was considered in the others.
For each type of crop, the management practice was specified, such as planting, tillage, fertilization, irrigation and harvesting, among others. All operations were scheduled by the heat unit the cropping calendar. Fertilizers were applied with no effects of nitrogen and phosphorus stress through automatic fertilization application option in the SWAT model. The major sources of irrigation are reservoirs and river reaches (only irrigation district 4, 10 and 11 shown in Figure 2b). Irrigation water was extracted from the river reach or reservoir of a subbasin in response to the irrigation district location. Irrigation practices were accounted for in the simulations through automatic irrigation demand at the designated subbasins (it can be considered as sufficient irrigation without water stress), which aims to quantify the potential water requirement for rice growth under the influence of future climate change.

2.4. Crop Water Footprint

The water footprint (WF) is the total amount of fresh water needed to produce the goods and services related to a certain consumption pattern, which consists of green, blue, and grey WF [2]. The total WF of crop production can be calculated using the following equation:
W F = W F g r e e n + W F b l u e + W F g r e y
where WF is the total WF of crop production (m3·kg−1); and WFgreen, WFblue and WFgrey are the green, blue and grey WF (m3·kg−1), respectively.
WFgreen is calculated as follows:
W F g r e e n = W g r e e n / Y = 10 i = 1 T E T g r e e n / Y
where WFgreen is the green WF (m3·kg−1); Wgreen is the amount of green water use per unit area (m3·ha−1); Y is the crop yield per unit area (kg·ha−1); T is the length of the crop growing season (d); and ETgreen is the daily green water evapotranspiration (mm), which is the minor factor between the daily crop evapotranspiration ET and the daily effective precipitation Peff.
E T g r e e n = min ( E T , P e f f )
where ET is the daily evapotranspiration (mm) and Peff is the daily effective precipitation (mm).
Effective precipitation is the fraction of the total precipitation that is available to the crop and does not run off [42]. Based on this definition, Peff is calculated by using the SWAT model as follows [6]:
P e f f = P R O D P
where P is rainfall, RO is surface runoff and DP is deep percolation.
Wblue is calculated as follows [2]:
W F b l u e = W b l u e / Y = 10 i = 1 T E T b l u e / Y
where WFblue is the blue WF (m3·kg−1); Wblue is the amount of blue water use per unit area (m3·ha−1); and ETblue is the daily blue water evapotranspiration (mm). In irrigated systems, ETblue is the difference between ET and ETgreen, and ETblue = 0 in rain-fed systems [6,43].
E T b l u e = { E T E T g r e e n i n   i r r i g a t e d   s y s t e m s 0 i n   r a i n f e d   s y s t e m s
WFgrey is defined as the freshwater required to assess the load of pollutants during crop growth, and it expresses the degradative water use [7].
W F g r e y = a A R / ( c max c n a t ) / Y
where WFgrey is the grey WF; AR is the product of the fertilizer application rate in the field per unit area (kg·ha−1); and a is the leaching-runoff fraction. Nitrogen is the main pollutant considered for most water footprint studies, which is estimated due to its potential pollution effect via nitrate leaching into ground and surface waters, here a = 0.14 [43]; and cmax is the maximum acceptable concentration (10 mg·L−1) and cnat is the concentration in natural water (0.1 mg·L−1) [44].

2.5. Simulation of Future Crop Yield and Water Footprint

To quantify the impact of climate change on the crop yield and WF of rice in the Nanliujiang catchment, the bias-corrected output of climate variables (precipitation, temperature, and CO2) from GCMs were fed into the calibrated SWAT model, and thus rice yield and WF were calculated for future periods. Here, the reference period (2001–2010) was ten years, and therefore the future reference period of the same length (2041–2050) was used for improved comparability. It was assumed that there will be no change in land use for the project area in the future. Climate variables in the future, such as wind speed, relative humidity, and sunshine hours, were also assumed to be the same as the reference period.

3. Results and Discussion

3.1. Performance of the SWAT Model

3.1.1. LAI and Biomass

The crucial parameters that pertain to crop growth were calibrated on a daily time step to minimize the simulated and observed residuals for LAI and BIOM. The related parameters of early and late rice were displayed their initial and final values in Table 2. It can be seen that the values of BIO_E, FRGRW2, LAIMX2 and BLAI showed greater differences between early and late rice, which were caused mainly by rice varieties and climate variations in growing season. LAI and BIOM were adjusted observed measurements in the Hepu irrigation district from 2005–2007 (Figure 3). The performance indicators of calibration in this study were Nash–Sutcliffe (NS), root mean square error (PBIAS). Figure 3 and Table 3 display that both early and late rice are well depicted, with PBIAS values ranging from −8.61 to 2.13 and NS values from 0.83 to 0.97. Similarly, validation indicated an acceptable model performance with PBIAS values ranging from −7.43 to 8.03 and NS from 0.83 to 0.97. It can be seen that the SWAT model can be used to simulate the development of rice well.

3.1.2. Discharge and Evapotranspiration

The first years (1995–1997) used to warm-up the model, and the data from 1998 to 2006 were used for calibration and those from 2006 to 2013 for validation. The 12 hydrological parameters were selected for 3 hydrological stations (Hengjiang, Bobai and Changle), and their optimal values for SWAT calibration are shown in Table 4. The results showed that the accuracy of the simulated value for the monthly streamflow at the Nanliujiang catchment in the calibration and verification period affectively met the standard (Figure 4). PBIAS of the simulated streamflow at the three stations ranged from 2.8 to 11.7 and from 3.5 to 13.0, respectively, and NS of the simulated streamflow at the three stations ranged from 0.85 to 0.94 and from 0.73 to 0.82, respectively (Table 5).
A direct relationship has been between crop development and evapotranspiration [45], and thus the evapotranspiration is also necessary for calibration in this work. However, the data of the actual evapotranspiration is hard to be measured or recorded, and therefore the measured monthly pan evaporation data from 2000 to 2011 for 3 stations (Yulin, Bobai and Hepu) may be a good source of information. Here, the values of the actual evapotranspiration were estimated by using a simple conversion method. Xie et al., 2009 [46] examined the ratio between the actual evapotranspiration and pan evaporation in the Dongjiang catchment, which is near to the Nanliujiang catchment and shows similar climatic characteristics. The actual evapotranspiration was estimated by multiplying the measured pan evaporation with the ratio. The calibration results showed an acceptable performance of the model with PBIAS values from −4.2 to 10.4 and NS values from 0.70 to 0.85. Similarly, validation indicated an acceptable model performance with PBIAS values ranging from −8.42 to 13.2 and NS from 0.59 to 0.79 (Table 6).

3.1.3. Rice Yields

Root mean squared errors (RMSE) and mean relative error (MRE) were used to evaluate the SWAT performance. The parameters were calibrated measured yield data from two stations (Yulin and Hepu) for the period from 1999–2012. Table 7 showed acceptable performance by the model, with PBIAS values ranging from −13.13 to 7.90 and RMRE values from 0.24 to 0.83, and the MRE between the observed and simulated yield was less than 15%. All the results were considered acceptable for calibration of the model. Figure 4 showed the simulated and measured rice yields at Hepu station from 1999 to 2012.

3.2. Crop Yield and Water Footprint in the Reference Period

The rice yield in the Nanliujiang catchment was estimated using the validated SWAT model for the reference period (2001–2010), as shown in Figure 5a. The estimated total, green, blue and grey WF of rice in the 2000s calculated by Equations (6)–(12) is shown in Figure 5b–e respectively. Statistics (mean) of the yield and WF values are also shown in the figures. A combination of water use, soil and fertilization amount determined the spatial distribution characteristics of rice production in the catchment. Figure 5a shows that the yields of irrigated rice were higher than those of rain-fed rice. Since the rice in the rain-fed area was grown under rain-fed conditions during the growing period, there was a rapid decline in yield relative to that grown in irrigation districts under the same climate conditions. Rice yields varied from 4761 to 5264 kg·ha−1 in the irrigated area and from 3035 to 3378 kg·ha−1 in the rain-fed area. In the irrigated subbasins, the rice yields in the southern catchment were usually larger than that in northern catchment, and the rain-fed subbasins shared a similar characteristic. Hepu and Hongchaojiang irrigation districts (subbasin 24, 48–59) showed the maximum yield per acre in all irrigation districts, vary from 5174 to 5264 kg·ha−1. The Maolin irrigation district (subbasin 31) showed the minimum yield per acre at 4761 kg·ha−1.
The average total WF of rice production in the study area was 1064 m3·t−1 (Figure 5b), which was higher than the national average (971 m3·t−1) but below the global average (1450 m3·t−1) [47]. The largest share of the total WF corresponded to the green WF, accounting for 44.8%, followed by blue and grey WFs, which accounted for 29.1% and 26.1%, respectively (Figure 5d,e). The Hongchaojiang irrigation district (subbasin 56 and 59) showed the maximum WF in all irrigation districts of 1049 m3·t−1, with the green WF representing 179 m3·t−1, blue 591 m3·t−1, and grey 279 m3·t−1; the Longmen irrigation district (subbasin 27) showed the minimum total WF, consisting of a green WF of 200 m3·t−1, blue of 490 m3·t−1, and grey of 213 m3·t−1, resulting in a total WF of 903 m3·t−1. Subbasin 27 showed the minimum WF of 903 m3·t−1, and subbasin 9 showed the maximum WF of 1248 m3·t−1. The values of the WF of rain-fed subbasins with low rice production were higher than those of irrigated subbasins, confirming that a greater reduction in yield can result in large water footprints [6].

3.3. Future Climate Change

Since the south long-narrow catchment experiences a spatial and temporal variation in precipitation and temperature, the projection of future climate differs in different stations, which can be observed in the results provided by the different GCMs and RCPs. Future average precipitation and maximum and minimum temperature in the 2040s are shown in Table 8 during the growth season (March to November), under reference period (RE) and climate change scenarios. RCP8.5 of GFDL-ESM2M and RCP2.6 of HadGEM2-ES show a slight decrease in precipitation, whereas the others show an increase from 1390 mm currently to 1647 mm under RCP4.5 of GFDL-ESM2M at Yulin station, and from 1705 mm currently to 2043 mm under RCP2.6 of GFDL-ESM2M at Hepu station; RCP2.6 and 4.5 of GFDL-ESM2M show a decrease in temperature, whereas the others show an increase from 30.6 °C to 31.5 °C at a maximum temperature and from 22.8 °C to 23.8 °C at a minimum temperature under RCP8.5 of HadGEM2-ES at Yulin station, and from 29.8 °C currently to 32.3 °C at a maximum temperature and 23.5 °C currently to 25.6 °C at a minimum temperature under RCP8.5 of HadGEM2-ES at Hepu station.

3.4. Rice Yield and Water Footprint in the Future Climate

3.4.1. Rice Yield in the Future Climate

Rice is a C3 plant, and therefore an increase in CO2 can, in theory, increase its yield; however, yield is controlled by a combined effect of precipitation, temperature and CO2. A higher temperature can shorten the period of rice growth (according to the heat unit approach), and temperature stress can reduce the accumulation rate of dry materials, both of which may lead to low production. Figure 6 shows the potential rice yield in the 2040s under 9 climate change scenarios. On average, rice yield would increase by +3.0% in the future relative to the 2000s, ranging from −3.5% (HadGEM2-ES, RCP8.5) to +10.5% (GFDL-ESM2M, RCP4.5). Both GFDL-ESM2M and IPSL-CM5A-LR provide an increasing yield by different degrees under all RCPs, but HadGEM2-ES instead depict a decreasing yield under all RCPs (Figure 6). GFDL-ESM2M shows an increased yield up to +10.9% for irrigated rice and +10.3% for rain-fed rice under RCP4.5 in the 2040s relative to the 2000s, due to increasing precipitation and CO2. HadGEM2-ES provide the lowest rice yield due to the highest temperatures compared with the other models under all RCPs (Table 8), because rice yield will be largely impacted by temperature stress as well as a shortened growing period.
The spatial distributions of the average rice yields of the 3 GCMs under 3 RCPs in the 2040s were shown in Figure 7. It can be seen that the spatially clustered areas with high rice production are distributed in the northern (subbasin 3, 6, 7, 25–39), middle (subbasin 40–47) and southern (subbasin 24, 48–59) study area, which are the important irrigation districts and grain production base in South China. This finding can be attributed to strong factors such as a smooth topography, favorable soil quality and proper irrigation facilities. The average rice yield in the irrigation districts will be increased by 239 kg·ha−1 under RCP2.6, 380 kg·ha−1 under RCP4.5 and 199 kg·ha−1 under RCP8.5 relative to the 2000s, respectively. Subbasins with low rice yield are distributed in mountainous and hilly areas with rain-fed farming due to a shortage of water conservancy facilities. The mountainous subbasins with larger yields will be subbasin 14–16, 18–21 and 64, which profit from advantages in climate conditions. The average rice yield will be increased by 138 kg·ha−1 under RCP2.6, 223 kg·ha−1 under RCP4.5 and 105 kg·ha−1 under RCP8.5 relative to the 2000s, respectively. However, their rice yields will be no larger than 3400 kg·ha−1. These mountainous subbasins show a rugged topography with poor soils, mountainous cultivation could not only destroy mountain vegetation but also lead to soil erosion and environmental pollution. Therefore, reasonable adjustment and optimization of the agricultural structure is of great necessity in this area in the future. Mountain areas may be more adaptable to planting economic forests that can generate better economic returns, such as longan, litchi, mango and citrus.

3.4.2. Water Footprint in the Future Climate

Figure 8a shows that the average total WF of irrigated rice in the 2040s will decrease by different degrees under all scenarios. On average, the total WF will be reduced by 2.3% under RCP2.6, 4.9% under RCP4.5 and 5.4% under RCP8.5 relative to the 2000s, respectively. A projected rise in rice yield along with a decrease in blue water may be the main cause of this phenomenon. Due to the rise in precipitation during the growing season, rice has more green water to meet the increasing evapotranspiration demand triggered by climate change, leading to a substantial increase in green water demand and a reduction in blue water use. Hence, blue WF shows a significant decrease by 43.0% under RCP2.6, 48.3% under RCP4.5 and 45.2% under RCP8.5, respectively. Accordingly, green WF shows an increase from 207 m3·t−1 to 394 m3·t−1 under RCP2.6, 395 m3·t−1 under RCP4.5 and 382 m3·t−1 under RCP8.5, respectively. The amount of nitrogen fertilizer application will significantly increase with rice yield growth in the future, whereas the gray WF is likely to slightly increase by 4.0% under RCP2.6, 4.1% under RCP4.5 and 1.0% under RCP8.5, respectively, because rice yield improvements have helped to reduce obvious growth in the gray WF. Figure 8b shows that the mean total WF of rain-fed rice under 9 scenarios in the 2040s showed no obvious upward or downward trend. A slight decline in green WF and rise in gray WF may be the main causes. The green WF shows a small decrease by 4.0% under RCP2.6, 6.7% under RCP4.5 and 7.0% under RCP8.5, respectively. The amount of fertilizer applied to rain-fed rice will be important for growth along with rice yield growth in the 2040s, and the gray WF is likely to increase under all scenarios by 13.9% under RCP2.6, 17.1% under RCP4.5 and 10.2% under RCP8.5, respectively.
The spatial distributions of the average total, green, blue and gray WF of the 3 GCMs under 3 RCPs in the 2040s are shown in Figure 9. The WFs of rice production exhibit cleat differences among subbasins under 3 RCPs, mainly due to the differences in rice yields and water use levels. The rice yields and water use differ between subbasins and irrigation districts because of their variances in climatic conditions, soil quality and fertilizing amount caused by climatic conditions. Although the values of total WFs are different under the 3 RCPs as well as the green, blue and gray WFs, their distributions under different RCPs are similar (Figure 9a vs. Figure 9e vs. Figure 9i for total WF; Figure 9b vs. Figure 9f vs. Figure 9j for green WF; Figure 9c vs. Figure 9g vs. Figure 9k for blue WF; Figure 9d vs. Figure 9h vs. Figure 9l for gray WF). For instance, Figure 9a,e,i show that the subbasins in the northern, middle and southern catchment have a lower green, gray and total WF, representing the subbasins located in plains planted with irrigated rice. The irrigation districts located in the southern catchment usually have higher WFs than other irrigation districts due to the high precipitation and temperature. Hongchaojiang irrigation district (subbasin 56, 59) and Longmen irrigation district (subbasin 27) show the maximum and minimum total WF among all the irrigation districts, respectively. The spatial distribution of green WF is consistent with the distribution of precipitation. Due to abundant rainfall, the proportion of green WF in the Hepu (Subbasin 24, 48–55 and 57–58) and Hongchaojiang irrigation districts is relatively high. Farmers should take technical measures to increase green water use and reduce irrigation water use as more as possible. The southern irrigation districts also have the largest gray WF in the catchment, which increase by 10.8% under RCP2.6, 13.2% under RCP4.5 and 8.3% under RCP8.5, respectively. It will increase the pressure of water resource protection and water pollution control in the downstream area. So the fertilizers application should be regulate scientifically in the future so that grey WF can be kept at a proper level.

4. Conclusions

In this study, we used the SWAT model to quantify the yield and WFs of double cropping rice within the Nanliujiang catchment, China. We validated the simulations against observed leaf area index, biomass, runoff, evapotranspiration and rice yield in sequence. The impact of potential climate change on the yield and WF of rice was analyzed in the 2040s under 3 GCMs with 3 RCPs. The results showed that: (a) the SWAT model is an ideal tool to simulate rice development and WF under given and future climate conditions; (b) yield simulation showed significant variations under different future climate scenarios. It showed increased rice yields ranged from +1.4 to +10.6% under the GFDL-ESM2M and IPSL-CM5A-LR models but slightly decreased yields ranged from −3.5 to −0.8% under the HadGEM2-ES model; (c) both the yield and WF of rice exhibited obvious differences among subbasins, with characteristically high values in the south and low values in the north, mainly due to differences in climatic conditions, soil quality and fertilization amounts; (d) blue WF would significantly decrease by 45.5% and green WF would show a great increase by 88.1%, meaning green water would contribute the largest part in the WF of rice production in the catchment. These findings would be good news for local water resource development and utilization, and it would provide favorable conditions to expand the area under irrigation and take technical measures for improving green water use efficiency in irrigation districts; (e) the gray WF would increase under potential climate change, suggesting an increased requirement for fertilizer application in the future. The finding suggested that it was necessary to improved fertilization management strategies for the agricultural non-point source pollution control in irrigation districts, especially for the Hongchaojiang and Hepu irrigation districts. This work will also delineate some helpful information for the management and protection of water resources as well as irrigation planning.

Acknowledgments

This research was financially supported by the National Key Research and Development Program (during the 13th Five-year Plan) under Grant No. 2016YFA0601500, Ministry of Science and Technology, PRC, the Science and technology project of Guangxi Provincial Water Resources Bureau (201402), and the National Natural Science Foundation of China (No. 51679252). We are grateful to the Guangxi water conservancy and Electric Power Survey and Design Research Institute and China Meteorological Data Sharing Service System for their help in providing data.

Author Contributions

Weihua Xiao and Yong Zhao came up with the idea and designed the study; Mingzhi Yang performed the analyses and wrote the paper; Xudong Li helped with language editing and the WF calculation and analyses; Ya Huang and Fan Lu helped in collecting the data and partial analyses; Baodeng Hou and Baoqi Li contributed the rice yield analysis and figures production.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IPCC (Intergovernmental Panel on Climate Change). Climate Change 2014–Impacts, Adaptation and Vulnerability: Regional Aspects; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014. [Google Scholar]
  2. Hoekstra, A.Y.; Mekonnen, M.M. The water footprint of humanity. Proc. Natl. Acad. Sci. USA 2012, 109, 3232–3237. [Google Scholar] [CrossRef] [PubMed]
  3. Bocchiola, D.; Nana, E.; Soncini, A. Impact of climate change scenarios on crop yield and water footprint of maize in the Po valley of Italy. Agric. Water Manag. 2013, 116, 50–61. [Google Scholar] [CrossRef]
  4. Lobell, D.B.; Howden, S.M.; Smith, D.R.; Howden, S.M.; Smith, D.R.; Chhetri, N. A meta-analysis of crop yield under climate change and adaptation. Nat. Clim. Chang. 2014, 4, 287–291. [Google Scholar]
  5. Hoekstra, A.Y. Virtual Water Trade: Proceedings of the International Expert Meeting on Virtual Water Trade, Value of Water Research Report Series, No 12; UNESCO-IHE: Delft, The Netherlands, 2003. [Google Scholar]
  6. Dourte, D.R.; Fraisse, C.W.; Uryasev, O. Water Footprint on AgroClimate: A dynamic, web-based tool for comparing agricultural systems. Agric. Syst. 2014, 125, 33–41. [Google Scholar] [CrossRef]
  7. Chapagain, A.K.; Hoekstra, A.Y. The blue, green and grey water footprint of rice from production and consumption perspectives. Ecol. Econ. 2011, 70, 749–758. [Google Scholar] [CrossRef]
  8. Bocchiola, D. Impact of potential climate change on crop yield and water footprint of rice in the Po valley of Italy. Agric. Syst. 2015, 139, 223–237. [Google Scholar] [CrossRef]
  9. Deng, G.; Ma, Y.; Li, X. Regional water footprint evaluation and trend analysis of China—Based on interregional input–output model. J. Clean. Prod. 2016, 112, 4674–4682. [Google Scholar] [CrossRef]
  10. Shrestha, S.; Chapagain, R.; Babel, M.S. Quantifying the impact of climate change on crop yield and water footprint of rice in the Nam Oon Irrigation Project, Thailand. Sci. Total Environ. 2017, 599–600, 689–699. [Google Scholar] [CrossRef] [PubMed]
  11. Jinger, D.; Kaur, R.; Kaur, N.; Dass, A. Weed Dynamics under Changing Climatic Scenario: A Review. Int. J. Curr. Microbiol. Appl. Sci. 2017, 6, 2376–2388. [Google Scholar]
  12. Kimball, B.A. Carbon dioxide and agricultural yield: An assemblage and analysis of prior observations. Agron. J. 1983, 75, 779–788. [Google Scholar] [CrossRef]
  13. Teixeira, E.I.; Fischer, G.; Van, V.H.; Walter, C.; Ewert, F. Global hot-spots of heat stress on agricultural crops due to climate change. Agric. For. Meteorol. 2013, 170, 206–215. [Google Scholar] [CrossRef]
  14. Hanson, M.R.; Lin, M.T.; Carmo-Silva, A.E.; Parry, M.A. Towards engineering carboxysomes into C3 plants. Plant J. 2016, 87, 38–50. [Google Scholar] [CrossRef] [PubMed]
  15. Thanawong, K.; Perret, S.R.; Basset-Mens, C. Eco-efficiency of paddy rice production in Northeastern Thailand: Acomparison of rain-fed and irrigated cropping systems. J. Clean. Prod. 2014, 73, 204–217. [Google Scholar] [CrossRef]
  16. Bazargan, A.; Tan, J.; Hui, C.W.; McKay, G. Utilization of rice husks for the production of oil sorbent materials. Cellulose 2014, 21, 1679–1688. [Google Scholar] [CrossRef]
  17. United States of Agricultural Department (USDA): Grain: World Markets and Trade, 2017. Available online: http://usda.mannlib.cornell.edu/usda/fas/grain-market//2010s/2016/grain-market-12-09-2016.pdf (accessed on 10 November 2017).
  18. Ray, D.K.; Gerber, J.S.; Macdonald, G.K.; West, P.C. Climate variation explains a third of global crop yield variability. Nat. Commun. 2016, 6, 5989. [Google Scholar] [CrossRef] [PubMed]
  19. Keating, B.A.; Carberry, P.S.; Hammer, G.L.; Probert, M.E.; Robertson, M.J.; Holzworth, D.; Huth, N.I.; Hargreaves, J.N.G.; Meinke, H.; Hochman, Z.; et al. An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 2003, 18, 267–288. [Google Scholar] [CrossRef]
  20. Stöckle, C.O.; Donatelli, M.; Nelson, R. CropSyst, a cropping systems simulation model. Eur. J. Agron. 2003, 18, 289–307. [Google Scholar] [CrossRef]
  21. Jones, J.W.; Hoogenboom, G.; Porter, C.H.; Boote, K.J.; Batchelor, W.D.; Hunt, L.A.; Wilkens, P.W.; Singh, U.; Gijsman, A.J.; Ritchie, J.T. The DSSAT cropping system model. Eur. J. Agron. 2003, 18, 235–265. [Google Scholar] [CrossRef]
  22. Williams, J.R.; Jones, C.A.; Dyke, P.T. A modelling approach to determining the relationship between erosion and soil productivity. Trans. ASAE 1984, 27, 129–144. [Google Scholar] [CrossRef]
  23. Ahmadzadeh, H.; Morid, S.; Delavar, M.; Srinivasan, R. Using the SWAT model to assess the impacts of changing irrigation from surface to pressurized systems on water productivity and water saving in the Zarrineh Rud catchment. Agric. Water Manag. 2015, 175, 15–28. [Google Scholar] [CrossRef]
  24. Palazzoli, I.; Maskey, S.; Uhlenbrook, S.; Nana, E.; Bocchiola, D. Impact of prospective climate change on water resources and crop yields in the Indrawati basin, Nepal. Agric. Syst. 2015, 133, 143–157. [Google Scholar] [CrossRef]
  25. Sun, S.; Wu, P.; Wang, Y.; Zhao, X.; Liu, J.; Zhang, X. The impacts of interannual climate variability and agricultural inputs on water footprint of crop production in an irrigation district of China. Sci. Total Environ. 2013, 444, 498. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, X.; Li, X.; Fischer, G.; Sun, L.; Tan, M.; Xin, L.; Liang, Z. Impact of the changing area sown to winter wheat on crop water footprint in the North China Plain. Ecol. Indic. 2015, 57, 100–109. [Google Scholar] [CrossRef]
  27. Zhuo, L.; Mekonnen, M.M.; Hoekstra, A.Y. Sensitivity and uncertainty in crop water footprint accounting: A case study for the Yellow River basin. Hydrol. Earth Syst. Sci. 2017, 11, 2219–2234. [Google Scholar]
  28. Guangxi Statistical Bureau. Guangxi Statistical Yearbook; Guangxi Statistical Bureau: Nanning, China, 2015. [Google Scholar]
  29. Vaghefi, S.A.; Mousavi, S.J.; Abbaspour, K.C.; Srinivasan, R.; Yang, H. Analyses of the impact of climate change on water resources components, drought and wheat yield in semiarid regions: Karkheh River Basin in Iran. Hydrol. Process. 2014, 28, 2018–2032. [Google Scholar] [CrossRef]
  30. Vaghefi, S.A.; Mousavi, S.J.; Abbaspour, K.C.; Srinivasan, R.; Arnold, J.R. Integration of hydrologic and water allocation models in basin-scale water resources management considering crop pattern and climate change: Karkheh River Basin in Iran. Reg. Environ. Chang. 2015, 15, 475–484. [Google Scholar] [CrossRef]
  31. Zhang, A.; Zhang, C.; Fu, G.; Wang, B.; Bao, Z.; Zheng, H. Assessments of impacts of climate change and human activities on runoff with SWAT for the Huifa river basin, northeast China. Water Resour. Manag. 2012, 26, 2199–2217. [Google Scholar] [CrossRef]
  32. Balkovič, J.; van der Velde, M.; Schmid, E.; Skalský, R.; Khabarov, N.; Obersteiner, M.; Stürmerbe, B.; Xionga, W. Pan-European crop modelling with EPIC: Implementation, up-scaling and regional crop yield validation. Agric. Syst. 2013, 120, 61–75. [Google Scholar] [CrossRef]
  33. Xiong, W.; Skalský, R.; Porter, C.H.; Balkovič, J.; Jones, J.W.; Yang, D. Calibratio-induced uncertainty of the EPIC model to estimate climate change impact on global maize yield. J. Adv. Model. Earth Syst. 2016, 8, 1358–1375. [Google Scholar] [CrossRef]
  34. Stockle, C.O.; Williams, J.R.; Rosenburg, N.J.; Jones, C.A. A method for estimating the direct and climatic effects of rising atmospheric carbon dioxide on growth and yield of crops: Part 1—Modification of the EPIC model for climate change analysis. Agric. Syst. 1992, 38, 225–238. [Google Scholar] [CrossRef]
  35. Hagemann, S.; Chen, C.; Haerter, J.O.; Heinke, J.; Gerten, D.; Piani, C. Impact of a statistical bias correction on the projected hydrological changes obtained from three GCMs and two hydrology models. J. Hydrometeorol. 2011, 12, 556–578. [Google Scholar] [CrossRef]
  36. Loon, A.F.V.; Gleeson, T.; Clark, J.; Van Dijk, A.I.; Stahl, K.; Hannaford, J.; di Baldassarre, G.; Teuling, A.J.; Tallaksen, L.M.; Uijlenhoet, R.; et al. Drought in the Anthropocene. Nat. Geosci. 2016, 9, 89–91. [Google Scholar] [CrossRef]
  37. Loon, A.F.V.; Stahl, K.; Baldassarre, G.D.; Clark, J.; Rangecroft, S.; Wanders, N.; Gleeson, T.; van Dijk, A.I.J.M.; Tallaksen, L.M.; Hannaford, J.; et al. Drought in a human-modified world: Reframing drought definitions, understanding, and analysis approaches. Hydrol. Earth Syst. Sci. 2016, 20, 3631–3650. [Google Scholar] [CrossRef]
  38. Zheng, J.; Li, G.Y.; Han, Z.Z.; Meng, G.X. Application of modified SWAT model in plain irrigation district. J. Hydraul. Eng. 2011, 42, 88–97. (In Chinese) [Google Scholar]
  39. Wu, Y.; Li, C.; Zhang, C.; Shi, X.; Zhao, S.; Lin, T. A watershed delineation method for mountains, plains complex landform area ArcGIS and SWAT. Arid Land Geogr. 2016, 39, 413–419. (In Chinese) [Google Scholar]
  40. Awan, U.K.; Liaqat, U.W.; Choi, M.; Ismaeel, A. A SWAT modeling approach to assess the impact of climate change on consumptive water use in Lower Chenab Canal area of Indus basin. Hydrol. Res. 2016, 47, 1025–1037. [Google Scholar] [CrossRef]
  41. Zhao, Y.; Zhang, J.; Pei, Y. Distributed water cycle simulation for plain area in Ningxia autonomous region. J. Hydraul. Eng. 2007, 38, 498–505. (In Chinese) [Google Scholar]
  42. Döll, P. Impact of Climate Change and Variability on Irrigation Requirements: A Global Perspective. Clim. Chang. 2002, 54, 269–293. [Google Scholar] [CrossRef]
  43. Mekonnen, M.M.; Hoekstra, A.Y. The green, blue and grey water footprint of crops and derived crop products. Hydrol. Earth Syst. Sci. 2011, 15, 1577–1600. [Google Scholar] [CrossRef] [Green Version]
  44. Gil, R.; Bojacá, C.R.; Schrevens, E. Uncertainty of the Agricultural Grey Water Footprint Based on High Resolution Primary Data. Water Resour. Manag. 2017, 31, 1–12. [Google Scholar] [CrossRef]
  45. Rouholahnejad, E.; Abbaspour, K.C.; Srinivasan, R.; Bacu, V.; Lehmann, A. Water resources of the Black Sea Basin at high spatial and temporal resolution. Water Resour. Res. 2015, 50, 5866–5885. [Google Scholar] [CrossRef]
  46. Xie, P.; Chen, X.; Wang, Z.; Xie, Y. Comparison of Actual Evapotranspiration and Pan Evaporation. Acta Geogr. Sin. 2009, 64, 270–277. (In Chinese) [Google Scholar]
  47. Chapagain, A.K.; Hoekstra, A.Y. The Green, Blue and Grey Water Footprint of Rice from Both a Production and Consumption Perspective; Value of Water Research Report Series No. 40; UNESCO-IHE Institute for Water Education: Delft, The Netherlands, 2010. [Google Scholar]
Figure 1. Location of the Nanliujiang catchment in China.
Figure 1. Location of the Nanliujiang catchment in China.
Sustainability 10 00242 g001
Figure 2. (a) Land use; (b) 18 irrigating areas and distribution of the main irrigation and drainage canal; (c) auto-delineated and (d) hand-delineated subbasins in the Nanliujiang catchment.
Figure 2. (a) Land use; (b) 18 irrigating areas and distribution of the main irrigation and drainage canal; (c) auto-delineated and (d) hand-delineated subbasins in the Nanliujiang catchment.
Sustainability 10 00242 g002
Figure 3. Comparison of simulated vs. observed BIOM and LAI for rice grown in the Hepu irrigation district.
Figure 3. Comparison of simulated vs. observed BIOM and LAI for rice grown in the Hepu irrigation district.
Sustainability 10 00242 g003
Figure 4. Comparison of simulated vs. observed rice yield in Hepu station.
Figure 4. Comparison of simulated vs. observed rice yield in Hepu station.
Sustainability 10 00242 g004
Figure 5. Spatial variations in the (a) yield and (b) total; (c) green; (d) blue and (e) grey WFs of rice from 2001–2010.
Figure 5. Spatial variations in the (a) yield and (b) total; (c) green; (d) blue and (e) grey WFs of rice from 2001–2010.
Sustainability 10 00242 g005
Figure 6. Average rice yield with RE and climate change scenarios (3 GCMs with 3 RCPs) from 2041–2050.
Figure 6. Average rice yield with RE and climate change scenarios (3 GCMs with 3 RCPs) from 2041–2050.
Sustainability 10 00242 g006
Figure 7. Spatial variations in rice yields under RCP2.6, 4.5 and 8.5 in the Nanliujiang catchment from 2041–2050.
Figure 7. Spatial variations in rice yields under RCP2.6, 4.5 and 8.5 in the Nanliujiang catchment from 2041–2050.
Sustainability 10 00242 g007
Figure 8. Average WFs of (a) irrigated and (b) rain-fed rice with RE and climate change scenarios (3 GCMs with 3 RCPs) from 2041–2050.
Figure 8. Average WFs of (a) irrigated and (b) rain-fed rice with RE and climate change scenarios (3 GCMs with 3 RCPs) from 2041–2050.
Sustainability 10 00242 g008
Figure 9. Spatial variations of (a,e,i) total; (b,f,j) green; (c,g,k) blue and (d,h,l) gray WFs under RCP2.6, 4.5 and 8.5 from 2041–2050.
Figure 9. Spatial variations of (a,e,i) total; (b,f,j) green; (c,g,k) blue and (d,h,l) gray WFs under RCP2.6, 4.5 and 8.5 from 2041–2050.
Sustainability 10 00242 g009
Table 1. Descriptions of the 3 regional climate change scenarios used in this study.
Table 1. Descriptions of the 3 regional climate change scenarios used in this study.
ContentDescription
Model nameGFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR
ScenarioRCP2.6, RCP4.5, RCP8.5
Spatial resolution0.5° × 0.5°
Meteorological variablesPrecipitation, maximum and minimum air temperature
Temporal scaleDaily data from 1 January 1960 to 31 December 2050
Table 2. Results from the calibration of SWAT rice growth parameters.
Table 2. Results from the calibration of SWAT rice growth parameters.
Parameter NameParameter DefinitionEarly RiceLate Rice
Initial ValueFinal ValueInitial ValueFinal Value
BIO_EBiomass/energy ratio2224.872228.71
HVSTIHarvest index0.50.540.50.55
FRGRW1Fraction of the plant growing season corresponding to the first point on the optimal leaf area development curve (OLADC)0.30.310.30.34
FRGRW2Fraction of the plant growing season corresponding to the second point on the OLADC0.70.510.70.38
LAIMX1Fraction of the maximum leaf area index corresponding to the first point on the OLADC0.010.280.010.34
LAIMX2Fraction of the maximum leaf area index corresponding to the second point on the OLADC0.950.990.950.58
BLAIMaximum leaf area index56.3057.28
DLAIFraction of growing season when leaf area starts declining0.80.790.80.76
EXT_COFFLight extinction coefficient0.350.570.350.56
Table 3. Model performance for the simulation of rice in the Nanliujiang catchment. PBIAS—Percent bias; NS—Nash-Sutcliffe coefficient.
Table 3. Model performance for the simulation of rice in the Nanliujiang catchment. PBIAS—Percent bias; NS—Nash-Sutcliffe coefficient.
Crop TypeVariable NameCalibration (2005–2006)Validation (2007)
PBIASNSPBIASNS
Early RiceBIOM−5.620.94−7.430.97
LAI1.830.83−3.820.87
Late RiceBIOM2.130.978.030.95
LAI−8.610.957.720.83
Table 4. Results from the calibration of SWAT hydrological parameters.
Table 4. Results from the calibration of SWAT hydrological parameters.
Parameter NameParameter DefinitionInitial RangeCalibrated Value
HengjiangBobaiChangle
CN2SCS runoff curve numberr (−0.5, 0.5)0.10−0.27−0.24
ESCOSoil evaporation compensation factorv (0, 1)0.210.240.31
EPCOPlant uptake compensation factorv (0, 1)0.730.650.86
GWQMNThreshold depth of water in the shallow aquifer required for return flow to occurv (0, 5000)1615.41248.51192.5
GW_REVAPGroundwater “revap” coefficientv (0.02, 2)0.480.301.67
REVAPMNThreshold depth of water in the shallow aquifer for “revap” to occurv (0, 500)105.67334.0164.39
SOL_AWCAvailable water capacity of the soil layerr (−0.5, 0.5)0.280.270.28
SOL_KSaturated hydraulic conductivityr (−0.5, 0.5)0.30−0.140.17
ALPHA_BFBase flow alpha factorv (0, 1)0.480.900.55
CANMXMaximum canopy storager (−1, 1)0.890.690.12
GW_DELAYGroundwater delayv (0, 500)35.7821.2926.03
SLSUBBSNAverage slope lengthr (−0.5, 0.5)−0.26−0.10−0.25
Notes: r indicates the existing parameter value is multiplied by (1 + a given value), v indicates the default parameter is replaced by a given value.
Table 5. Model performance for the simulation of monthly discharge in the Nanliujiang catchment.
Table 5. Model performance for the simulation of monthly discharge in the Nanliujiang catchment.
Station NameCalibration (1998–2005)Validation (2006–2013)
PBIASNSPBIASNS
Hengjiang11.720.8513.020.73
Bobai3.220.903.540.87
Changle2.810.944.330.82
Table 6. Model performance for the simulation of monthly evapotranspiration in the Nanliujiang catchment.
Table 6. Model performance for the simulation of monthly evapotranspiration in the Nanliujiang catchment.
Station NameCalibration (2000–2005)Validation (2006–2011)
PBIASNSPBIASNS
Yulin−4.210.70−6.530.59
Bobai10.40.7413.20.61
Hepu−7.130.85−8.420.79
Table 7. Model performance for the simulation of rice yield.
Table 7. Model performance for the simulation of rice yield.
Station NameCrop TypeCalibration (1999–2006)Validation (2007–2012)
PBIASRMSE/tonMRE/%PBIASRMSE/tonMRE/%
YulinEarly Rice−13.130.8313.11−8.550.596.11
Late Rice7.900.6311.421.490.628.98
HepuEarly Rice0.310.247.154.580.298.61
Late Rice−6.630.566.832.010.495.02
Table 8. Average values of climate variables during the growth season for the 2000s and 2040s under 9 climate change scenarios.
Table 8. Average values of climate variables during the growth season for the 2000s and 2040s under 9 climate change scenarios.
PeriodScenariosPrecipitation (mm)Temperature (°C)
MaximumMinimum
YulinHepuYulinHepuYulinHepu
2000sRE1390170530.629.822.823.5
2040sGFDL-ESM2MRCP2.61553204329.630.522.324.0
RCP4.51647195329.430.422.524.2
RCP8.51377165530.631.623.124.9
HadGEM2-ESRCP2.61364160531.332.223.325.3
RCP4.51423178531.232.023.525.4
RCP8.51450177131.532.323.825.6
IPSL-CM5A-LRRCP2.61420170630.231.222.924.6
RCP4.51451187130.431.323.225.0
RCP8.51413173131.032.023.725.5

Share and Cite

MDPI and ACS Style

Yang, M.; Xiao, W.; Zhao, Y.; Li, X.; Huang, Y.; Lu, F.; Hou, B.; Li, B. Assessment of Potential Climate Change Effects on the Rice Yield and Water Footprint in the Nanliujiang Catchment, China. Sustainability 2018, 10, 242. https://0-doi-org.brum.beds.ac.uk/10.3390/su10020242

AMA Style

Yang M, Xiao W, Zhao Y, Li X, Huang Y, Lu F, Hou B, Li B. Assessment of Potential Climate Change Effects on the Rice Yield and Water Footprint in the Nanliujiang Catchment, China. Sustainability. 2018; 10(2):242. https://0-doi-org.brum.beds.ac.uk/10.3390/su10020242

Chicago/Turabian Style

Yang, Mingzhi, Weihua Xiao, Yong Zhao, Xudong Li, Ya Huang, Fan Lu, Baodeng Hou, and Baoqi Li. 2018. "Assessment of Potential Climate Change Effects on the Rice Yield and Water Footprint in the Nanliujiang Catchment, China" Sustainability 10, no. 2: 242. https://0-doi-org.brum.beds.ac.uk/10.3390/su10020242

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop