Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Changes in the Potential Multiple Cropping System in Response to Climate Change in China from 1960–2010</title>
<author><name sortKey="Liu, Luo" sort="Liu, Luo" uniqKey="Liu L" first="Luo" last="Liu">Luo Liu</name>
<affiliation><nlm:aff id="aff1"><addr-line>State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2"><addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><addr-line>University of Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Xu, Xinliang" sort="Xu, Xinliang" uniqKey="Xu X" first="Xinliang" last="Xu">Xinliang Xu</name>
<affiliation><nlm:aff id="aff2"><addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zhuang, Dafang" sort="Zhuang, Dafang" uniqKey="Zhuang D" first="Dafang" last="Zhuang">Dafang Zhuang</name>
<affiliation><nlm:aff id="aff2"><addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Chen, Xi" sort="Chen, Xi" uniqKey="Chen X" first="Xi" last="Chen">Xi Chen</name>
<affiliation><nlm:aff id="aff1"><addr-line>State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Li, Shuang" sort="Li, Shuang" uniqKey="Li S" first="Shuang" last="Li">Shuang Li</name>
<affiliation><nlm:aff id="aff2"><addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><addr-line>University of Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">24312511</idno>
<idno type="pmc">3849037</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3849037</idno>
<idno type="RBID">PMC:3849037</idno>
<idno type="doi">10.1371/journal.pone.0080990</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">000177</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000177</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Changes in the Potential Multiple Cropping System in Response to Climate Change in China from 1960–2010</title>
<author><name sortKey="Liu, Luo" sort="Liu, Luo" uniqKey="Liu L" first="Luo" last="Liu">Luo Liu</name>
<affiliation><nlm:aff id="aff1"><addr-line>State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2"><addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><addr-line>University of Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Xu, Xinliang" sort="Xu, Xinliang" uniqKey="Xu X" first="Xinliang" last="Xu">Xinliang Xu</name>
<affiliation><nlm:aff id="aff2"><addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zhuang, Dafang" sort="Zhuang, Dafang" uniqKey="Zhuang D" first="Dafang" last="Zhuang">Dafang Zhuang</name>
<affiliation><nlm:aff id="aff2"><addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Chen, Xi" sort="Chen, Xi" uniqKey="Chen X" first="Xi" last="Chen">Xi Chen</name>
<affiliation><nlm:aff id="aff1"><addr-line>State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Li, Shuang" sort="Li, Shuang" uniqKey="Li S" first="Shuang" last="Li">Shuang Li</name>
<affiliation><nlm:aff id="aff2"><addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><addr-line>University of Chinese Academy of Sciences, Beijing, China</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint><date when="2013">2013</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>The multiple cropping practice is essential to agriculture because it has been shown to significantly increase the grain yield and promote agricultural economic development. In this study, potential multiple cropping systems in China are calculated based on meteorological observation data by using the Agricultural Ecology Zone (AEZ) model. Following this, the changes in the potential cropping systems in response to climate change between the 1960s and the 2010s were subsequently analyzed. The results indicate that the changes of potential multiple cropping systems show tremendous heterogeneity in respect to the spatial pattern in China. A key finding is that the magnitude of change of the potential cropping systems showed a pattern of increase both from northern China to southern China and from western China to eastern China. Furthermore, the area found to be suitable only for single cropping decreased, while the area suitable for triple cropping increased significantly from the 1960s to the 2000s. During the studied period, the potential multiple cropping index (PMCI) gap between rain-fed and irrigated scenarios increased from 18% to 24%, which indicated noticeable growth of water supply limitations under the rain-fed scenario. The most significant finding of this research was that from the 1960s to the 2000s climate change had led to a significant increase of PMCI by 13% under irrigated scenario and 7% under rain-fed scenario across the whole of China. Furthermore, the growth of the annual mean temperature is identified as the main reason underlying the increase of PMCI. It has also been noticed that across China the changes of potential multiple cropping systems under climate change were different from region to region.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Ziqiang, F" uniqKey="Ziqiang F">F Ziqiang</name>
</author>
<author><name sortKey="Yunlong, C" uniqKey="Yunlong C">C Yunlong</name>
</author>
<author><name sortKey="Youxiao, Y" uniqKey="Youxiao Y">Y Youxiao</name>
</author>
<author><name sortKey="Erfu, D" uniqKey="Erfu D">D Erfu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hong, Y" uniqKey="Hong Y">Y Hong</name>
</author>
<author><name sortKey="Xiubin, Y" uniqKey="Xiubin Y">Y Xiubin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Von, B" uniqKey="Von B">B Von</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Guo, B" uniqKey="Guo B">B Guo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jinshen, Y" uniqKey="Jinshen Y">Y Jinshen</name>
</author>
<author><name sortKey="Limin, S" uniqKey="Limin S">S Limin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shulan, J" uniqKey="Shulan J">J Shulan</name>
</author>
<author><name sortKey="Lichun, H" uniqKey="Lichun H">H Lichun</name>
</author>
<author><name sortKey="Lei, X" uniqKey="Lei X">X Lei</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shulan, J" uniqKey="Shulan J">J Shulan</name>
</author>
<author><name sortKey="Caiqiu, X" uniqKey="Caiqiu X">X Caiqiu</name>
</author>
<author><name sortKey="Huahua, P" uniqKey="Huahua P">P Huahua</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jing, L" uniqKey="Jing L">L Jing</name>
</author>
<author><name sortKey="Zhiyuan, R" uniqKey="Zhiyuan R">R Zhiyuan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dailiang, P" uniqKey="Dailiang P">P Dailiang</name>
</author>
<author><name sortKey="Jingfeng, H" uniqKey="Jingfeng H">H Jingfeng</name>
</author>
<author><name sortKey="Huimin, J" uniqKey="Huimin J">J Huimin</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Huimin, Y" uniqKey="Huimin Y">Y Huimin</name>
</author>
<author><name sortKey="Jiyuan, L" uniqKey="Jiyuan L">L Jiyuan</name>
</author>
<author><name sortKey="Mingkui, C" uniqKey="Mingkui C">C Mingkui</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hongping, D" uniqKey="Hongping D">D Hongping</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shumin, L" uniqKey="Shumin L">L Shumin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mao, L" uniqKey="Mao L">L Mao</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Xin, J" uniqKey="Xin J">J Xin</name>
</author>
<author><name sortKey="Yua, Zr" uniqKey="Yua Z">ZR Yua</name>
</author>
<author><name sortKey="Van, Leeuwen" uniqKey="Van L">Leeuwen Van</name>
</author>
<author><name sortKey="Driessen, Pm" uniqKey="Driessen P">PM Driessen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huimin, Y" uniqKey="Huimin Y">Y Huimin</name>
</author>
<author><name sortKey="Yuling, F" uniqKey="Yuling F">F Yuling</name>
</author>
<author><name sortKey="Xiangming, X" uniqKey="Xiangming X">X Xiangming</name>
</author>
<author><name sortKey="Heqing, H" uniqKey="Heqing H">H Heqing</name>
</author>
<author><name sortKey="Honglin, H" uniqKey="Honglin H">H Honglin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jinlong, F" uniqKey="Jinlong F">F Jinlong</name>
</author>
<author><name sortKey="Bingfang, W" uniqKey="Bingfang W">W Bingfang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Xiaolin, Z" uniqKey="Xiaolin Z">Z Xiaolin</name>
</author>
<author><name sortKey="Qiang, L" uniqKey="Qiang L">L Qiang</name>
</author>
<author><name sortKey="Miaogen, S" uniqKey="Miaogen S">S Miaogen</name>
</author>
<author><name sortKey="Jin, C" uniqKey="Jin C">C Jin</name>
</author>
<author><name sortKey="Jin, W" uniqKey="Jin W">W Jin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ashton, S" uniqKey="Ashton S">S Ashton</name>
</author>
<author><name sortKey="Joseph, M" uniqKey="Joseph M">M Joseph</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zandbergen, P" uniqKey="Zandbergen P">P Zandbergen</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Jiyuan, L" uniqKey="Jiyuan L">L Jiyuan</name>
</author>
<author><name sortKey="Mingliang, L" uniqKey="Mingliang L">L Mingliang</name>
</author>
<author><name sortKey="Dafang, Z" uniqKey="Dafang Z">Z Dafang</name>
</author>
<author><name sortKey="Zengxiang, Z" uniqKey="Zengxiang Z">Z Zengxiang</name>
</author>
<author><name sortKey="Xiangzheng, D" uniqKey="Xiangzheng D">D Xiangzheng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jiyuan, L" uniqKey="Jiyuan L">L Jiyuan</name>
</author>
<author><name sortKey="Mingliang, L" uniqKey="Mingliang L">L Mingliang</name>
</author>
<author><name sortKey="Hanqin, T" uniqKey="Hanqin T">T Hanqin</name>
</author>
<author><name sortKey="Dafang, Z" uniqKey="Dafang Z">Z Dafang</name>
</author>
<author><name sortKey="Zengxiang, Z" uniqKey="Zengxiang Z">Z Zengxiang</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Jiyuan, L" uniqKey="Jiyuan L">L Jiyuan</name>
</author>
<author><name sortKey="Zengxiang, Z" uniqKey="Zengxiang Z">Z Zengxiang</name>
</author>
<author><name sortKey="Xinliang, X" uniqKey="Xinliang X">X Xinliang</name>
</author>
<author><name sortKey="Wenhui, K" uniqKey="Wenhui K">K Wenhui</name>
</author>
<author><name sortKey="Wancun, Z" uniqKey="Wancun Z">Z Wancun</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jiyuan, L" uniqKey="Jiyuan L">L Jiyuan</name>
</author>
<author><name sortKey="Zhang, Q" uniqKey="Zhang Q">Q Zhang</name>
</author>
<author><name sortKey="Hu, Y" uniqKey="Hu Y">Y Hu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutchinson, Mf" uniqKey="Hutchinson M">MF Hutchinson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutchinson, Mf" uniqKey="Hutchinson M">MF Hutchinson</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Monteith, Jl" uniqKey="Monteith J">JL Monteith</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jinlong, F" uniqKey="Jinlong F">F Jinlong</name>
</author>
<author><name sortKey="Bingfang, W" uniqKey="Bingfang W">W Bingfang</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group><journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">24312511</article-id>
<article-id pub-id-type="pmc">3849037</article-id>
<article-id pub-id-type="publisher-id">PONE-D-13-06797</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0080990</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Changes in the Potential Multiple Cropping System in Response to Climate Change in China from 1960–2010</article-title>
<alt-title alt-title-type="running-head">Multiple Cropping System to Climate Changes</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Liu</surname>
<given-names>Luo</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Xu</surname>
<given-names>Xinliang</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Zhuang</surname>
<given-names>Dafang</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Chen</surname>
<given-names>Xi</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Li</surname>
<given-names>Shuang</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>University of Chinese Academy of Sciences, Beijing, China</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Yan</surname>
<given-names>Wengui</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>National Rice Research Center, United States of America</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>xuxl@lreis.ac.cn</email>
</corresp>
<fn fn-type="conflict"><p><bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p>Conceived and designed the experiments: XX DZ XC. Performed the experiments: SL. Analyzed the data: XX. Wrote the paper: LL.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><year>2013</year>
</pub-date>
<pub-date pub-type="epub"><day>3</day>
<month>12</month>
<year>2013</year>
</pub-date>
<volume>8</volume>
<issue>12</issue>
<elocation-id>e80990</elocation-id>
<history><date date-type="received"><day>1</day>
<month>2</month>
<year>2013</year>
</date>
<date date-type="accepted"><day>15</day>
<month>10</month>
<year>2013</year>
</date>
</history>
<permissions><copyright-year>2013</copyright-year>
<copyright-holder>Liu et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/3.0/"><license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/3.0/">Creative Commons Attribution License</ext-link>
, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<abstract><p>The multiple cropping practice is essential to agriculture because it has been shown to significantly increase the grain yield and promote agricultural economic development. In this study, potential multiple cropping systems in China are calculated based on meteorological observation data by using the Agricultural Ecology Zone (AEZ) model. Following this, the changes in the potential cropping systems in response to climate change between the 1960s and the 2010s were subsequently analyzed. The results indicate that the changes of potential multiple cropping systems show tremendous heterogeneity in respect to the spatial pattern in China. A key finding is that the magnitude of change of the potential cropping systems showed a pattern of increase both from northern China to southern China and from western China to eastern China. Furthermore, the area found to be suitable only for single cropping decreased, while the area suitable for triple cropping increased significantly from the 1960s to the 2000s. During the studied period, the potential multiple cropping index (PMCI) gap between rain-fed and irrigated scenarios increased from 18% to 24%, which indicated noticeable growth of water supply limitations under the rain-fed scenario. The most significant finding of this research was that from the 1960s to the 2000s climate change had led to a significant increase of PMCI by 13% under irrigated scenario and 7% under rain-fed scenario across the whole of China. Furthermore, the growth of the annual mean temperature is identified as the main reason underlying the increase of PMCI. It has also been noticed that across China the changes of potential multiple cropping systems under climate change were different from region to region.</p>
</abstract>
<funding-group><funding-statement>This research was supported and funded by National Program on Key Basic Research Project (973 Program) (Grant No. 2010CB950901), Project of Remote sensing investigation and assessment for Ecological environment changes in China(No. STSN-14-00) and project of CAS action-plan for West Development (No. KZCX2-XB3-08-01) The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts><page-count count="12"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>Multiple cropping is a general cropping practice which involves growing two or more crops on the same field in a given year. It is one of the most effective ways to increase grain yield and to promote agricultural economic development. There is a common belief that food production and food security are continuously threatened by the increasing price of crude oil, extreme weather, and the diversion of crop land for bio-fuel production. Therefore, food security has long been an issue which has attracted global attention under a situation of urbanization by which the total area of active farmland continues to decrease. Insufficient attention has been paid to processes that increase crop yield, leading to an increase in the multiple cropping index (MCI; ratio of sown area to cropland area) <xref ref-type="bibr" rid="pone.0080990-Ziqiang1">[1]</xref>
–<xref ref-type="bibr" rid="pone.0080990-Hong1">[2]</xref>
. Due to environmental limitations, increased food production is unlikely to be achieved from farmland expansion in many regions <xref ref-type="bibr" rid="pone.0080990-Von1">[3]</xref>
. In this situation, multiple cropping is an effective method for improving crop productivity as it is an efficient use of natural resources (i.e. land, light, heat and so on) and human resources. This practice also helps to alleviate competition for land use between food production and economic crops.</p>
<p>China is one of the largest countries in the world and has succeeded in achieving the highest MCI. Nearly half of the cultivated land in China is subject to multiple cropping practices. This has resulted in a 25% increase in crop yields in the past few decades, which has enabled China to feed 22% of the world's population using only 7% of the world's cultivated land <xref ref-type="bibr" rid="pone.0080990-Guo1">[4]</xref>
. Thus, an increase in the use of multiple cropping practices has assisted China in successfully addressing the issue of food security. Two factors which further underline the importance of multiple cropping practices are the global trends of climate change and urbanization. Both of these factors have caused cropping systems to undergo significant change during the last several decades <xref ref-type="bibr" rid="pone.0080990-Jinshen1">[5]</xref>
–<xref ref-type="bibr" rid="pone.0080990-Zhiguo1">[10]</xref>
. In the case of China, this has resulted in a decrease in cropland area by 15% <xref ref-type="bibr" rid="pone.0080990-Huimin1">[11]</xref>
(primarily distributed in developed areas).</p>
<p>The statistical method used for acquiring information on cropping system changes was -based on administrative boundaries from the government, which ignored the spatial heterogeneity within the same administrative area and thus failed to represent the spatial characteristics of the multiple cropping systems <xref ref-type="bibr" rid="pone.0080990-Hongping1">[12]</xref>
–<xref ref-type="bibr" rid="pone.0080990-Mao1">[14]</xref>
. These systems can be measured according to the phenological and ecological characteristics of crop growth. This is done by identifying the periods of cropping cycles from time series data collected by remote sensing systems including AVHRR/NDVI, SPOT VEGETATION and MODIS <xref ref-type="bibr" rid="pone.0080990-Xin1">[15]</xref>
–<xref ref-type="bibr" rid="pone.0080990-Xiaolin1">[18]</xref>
. This method is restrained by the temporal resolution of remote sensing data, and therefore cannot capture the changes in the multiple cropping systems over long periods of time. However it is suitable for food security assessment and agricultural development planning, which aims to quantify long term sequenced MCI changes and investigate its response to global climate change.</p>
<p>In this study, the potential multiple cropping systems in China were calculated using the method of AEZ based on meteorological observation data. The spatio-temporal change patterns between 1960 and -2010 were also evaluated, via the spatial analysis techniques of the Geographic Information System (GIS). Finally, the characteristic variations of the potential cropping systems were analyzed, and their responses to climate change in different regions were illustrated.</p>
</sec>
<sec id="s2"><title>Data and Methods</title>
<sec id="s2a"><title>Data sources</title>
<p>The input data for this study included terrain elevation, soil, land use and meteorological data (including the monthly maximum air temperature, minimum air temperature, precipitation, relative humidity, wind speed and sunshine hours).</p>
<sec id="s2a1"><title>Terrain Elevation Data</title>
<p>The terrain elevation dataset derived from the Shuttle Radar Topography Mission C-band (SRTM) was the first publicly available near-global, high resolution raster Digital Elevation Model (DEM) <xref ref-type="bibr" rid="pone.0080990-Ashton1">[19]</xref>
. The SRTM data was distributed with a 90 m spatial resolution around the Earth, reduced from the original 30 m resolution via averaging and sub-sampling. Numerous worldwide applied studies have used SRTM data for environmental analysis <xref ref-type="bibr" rid="pone.0080990-Zandbergen1">[20]</xref>
–<xref ref-type="bibr" rid="pone.0080990-Farr1">[21]</xref>
. The SRTM Version 2 data (by naming convention SRTM3 for 3 arc sec data) was sourced from the Jet Propulsion Laboratory (<ext-link ext-link-type="uri" xlink:href="http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp">http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp</ext-link>
).</p>
</sec>
<sec id="s2a2"><title>Soil Data</title>
<p>Soil quality was determined by several parameters, including soil type, effective soil depth and soil water holding capacity. A nation-wide soil dataset at the scale of 1∶1,000,000 was provided by the Data Center for Resources and Environmental Sciences at the Chinese Academy of Sciences (RESDC) (<ext-link ext-link-type="uri" xlink:href="http://www.resdc.cn/first.asp">http://www.resdc.cn/first.asp</ext-link>
). Soil data was used to calculate the soil-water balance, which determined the potential and actual evapotranspiration for a reference crop and the length of its growing period (LGP, days).</p>
</sec>
<sec id="s2a3"><title>Land Use Data</title>
<p>A land-use database with a mapping scale of 1∶1,000,000 was developed by CAS (The Chinese Academy of Sciences). This database included five time periods: the late 1980s, the mid-1990s, the late 1990s, the mid-2000s and 2008. The primary data source for the land use database was Landsat MSS/TM/ETM CCD digital images. CBERS (the China-Brazil Earth Resources Satellite) and HJ-1 (small satellite constellation for environment and disaster monitoring) images were also used as a supplement for the areas uncovered by Landsat. The land-use data were classified into 25 categories, which were subsequently grouped into six classes: cropland, woodland, grassland, water body, built-up area and unused land. Detailed information on this database can be found in previous papers <xref ref-type="bibr" rid="pone.0080990-Jiyuan1">[22]</xref>
–<xref ref-type="bibr" rid="pone.0080990-Jiyuan5">[26]</xref>
. In this study, farmland data were extracted from 2008 land use data.</p>
</sec>
<sec id="s2a4"><title>Meteorological Data</title>
<p>Meteorological data, which included the monthly maximum air temperature, minimum air temperature, precipitation, relative humidity, wind speed at 10 m height and sunshine hours from 1960 to 2010, were obtained from the national agro-meteorological stations of China maintained by Chinese Meteorological Administration (CMA) (<ext-link ext-link-type="uri" xlink:href="http://cdc.cma.gov.cn">http://cdc.cma.gov.cn</ext-link>
). Because of the diverse terrain across China, the impact of topographic conditions on the interpolation of the meteorological data was also considered. The ANUSPLIN software <xref ref-type="bibr" rid="pone.0080990-Hutchinson1">[27]</xref>
–<xref ref-type="bibr" rid="pone.0080990-Hutchinson2">[28]</xref>
, which was designed for spatial interpolation of climate data, was used to interpolate the meteorological data with the DEM mentioned in Section 2.1.1. These data measured monthly for the above six key factors of plant growth were then interpolated by ANUSPLIN to individual 10 km resolutions based on the digital terrain model of China (excluding Taiwan).</p>
<p>Before further processing, all source data were resampled into a raster dataset with a 10 km spatial resolution, and the data were transformed into the same coordinate system (Krasovsky_1940_Albers projection system, with central_meridian 105°E, double standard parallel of 27°N and 45°N, and D_Krasovsky_1940 datum).</p>
</sec>
</sec>
<sec id="s2b"><title>Method</title>
<p>The calculation method for the potential multiple cropping systems were derived from the Global Agro-Ecological Zones Model (AEZ) which was developed in the 1970s, and updated in 2010 by the Food and Agriculture Organization (FAO) and the International Institute for Applied Systems Analysis (IIASA) <xref ref-type="bibr" rid="pone.0080990-FAO1">[29]</xref>
–<xref ref-type="bibr" rid="pone.0080990-Fischer2">[33]</xref>
. The crop growth process was primarily affected by the local climate conditions during the growth period. The potential multiple cropping index (PMCI) is directly determined by the light-temperature-water condition. <xref ref-type="fig" rid="pone-0080990-g001">Figure 1</xref>
shows the technological framework designed for the calculation of the multiple cropping systems.</p>
<fig id="pone-0080990-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g001</object-id>
<label>Figure 1</label>
<caption><title>Technological framework designed for potential estimation of multiple cropping systems.</title>
</caption>
<graphic xlink:href="pone.0080990.g001"></graphic>
</fig>
<sec id="s2b1"><title>Multiple cropping zones</title>
<p>To assess the potential multiple cropping systems, a number of multiple cropping zones were defined by matching both the thermal and water requirements with the time required for crop growth in different latitudinal thermal climates. The following parameters were used to define the cropping zones (<xref ref-type="table" rid="pone-0080990-t001">Table 1</xref>
, <xref ref-type="table" rid="pone-0080990-t002">Table 2</xref>
): length of the growing period (LGP), number of days with mean daily temperatures above 5°C and 10°C (LGP<sub>t = 5</sub>
and LGP<sub>t = 10</sub>
), accumulated temperature on days with a mean daily temperature ≥0°C and ≥10°C (TS<sub>t = 0</sub>
and TS<sub>t = 10</sub>
) and accumulated temperature during the growing period with a mean daily temperature ≥5°C and 10°C (TS-G<sub>t = 5</sub>
and TS-G<sub>t = 10</sub>
).</p>
<table-wrap id="pone-0080990-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.t001</object-id>
<label>Table 1</label>
<caption><title>Delineation of multiple cropping zones in tropical areas.</title>
</caption>
<alternatives><graphic id="pone-0080990-t001-1" xlink:href="pone.0080990.t001"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Zone</td>
<td align="left" rowspan="1" colspan="1">LGP</td>
<td align="left" rowspan="1" colspan="1">LGP<sub>t = 5</sub>
</td>
<td align="left" rowspan="1" colspan="1">LGP<sub>t = 10</sub>
</td>
<td align="left" rowspan="1" colspan="1">TS<sub>t = 0</sub>
</td>
<td align="left" rowspan="1" colspan="1">TS<sub>t = 10</sub>
</td>
<td align="left" rowspan="1" colspan="1">TS-G<sub>t = 5</sub>
</td>
<td align="left" rowspan="1" colspan="1">TS-G<sub>t = 10</sub>
</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">A</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">B</td>
<td align="left" rowspan="1" colspan="1">≥45</td>
<td align="left" rowspan="1" colspan="1">≥120</td>
<td align="left" rowspan="1" colspan="1">≥90</td>
<td align="left" rowspan="1" colspan="1">≥1600</td>
<td align="left" rowspan="1" colspan="1">≥1000</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">C</td>
<td align="left" rowspan="1" colspan="1">≥200</td>
<td align="left" rowspan="1" colspan="1">≥200</td>
<td align="left" rowspan="1" colspan="1">≥120</td>
<td align="left" rowspan="1" colspan="1">≥6400</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">≥3200</td>
<td align="left" rowspan="1" colspan="1">≥2700</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">D</td>
<td align="left" rowspan="1" colspan="1">≥240</td>
<td align="left" rowspan="1" colspan="1">≥240</td>
<td align="left" rowspan="1" colspan="1">≥165</td>
<td align="left" rowspan="1" colspan="1">≥6400</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">≥4000</td>
<td align="left" rowspan="1" colspan="1">≥3200</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">E</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">F</td>
<td align="left" rowspan="1" colspan="1">≥300</td>
<td align="left" rowspan="1" colspan="1">≥300</td>
<td align="left" rowspan="1" colspan="1">≥240</td>
<td align="left" rowspan="1" colspan="1">≥7200</td>
<td align="left" rowspan="1" colspan="1">≥7000</td>
<td align="left" rowspan="1" colspan="1">≥5100</td>
<td align="left" rowspan="1" colspan="1">≥4800</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">G</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">H</td>
<td align="left" rowspan="1" colspan="1">≥360</td>
<td align="left" rowspan="1" colspan="1">≥360</td>
<td align="left" rowspan="1" colspan="1">≥360</td>
<td align="left" rowspan="1" colspan="1">≥7200</td>
<td align="left" rowspan="1" colspan="1">≥7000</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<table-wrap id="pone-0080990-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.t002</object-id>
<label>Table 2</label>
<caption><title>Delineation of multiple cropping zones in non-tropical areas.</title>
</caption>
<alternatives><graphic id="pone-0080990-t002-2" xlink:href="pone.0080990.t002"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Zone</td>
<td align="left" rowspan="1" colspan="1">LGP</td>
<td align="left" rowspan="1" colspan="1">LGP<sub>t = 5</sub>
</td>
<td align="left" rowspan="1" colspan="1">LGP<sub>t = 10</sub>
</td>
<td align="left" rowspan="1" colspan="1">TS<sub>t = 0</sub>
</td>
<td align="left" rowspan="1" colspan="1">TS<sub>t = 10</sub>
</td>
<td align="left" rowspan="1" colspan="1">TS-G<sub>t = 5</sub>
</td>
<td align="left" rowspan="1" colspan="1">TS-G<sub>t = 10</sub>
</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">A</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">B</td>
<td align="left" rowspan="1" colspan="1">≥45</td>
<td align="left" rowspan="1" colspan="1">≥120</td>
<td align="left" rowspan="1" colspan="1">≥90</td>
<td align="left" rowspan="1" colspan="1">≥1600</td>
<td align="left" rowspan="1" colspan="1">≥1000</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">C</td>
<td align="left" rowspan="1" colspan="1">≥180</td>
<td align="left" rowspan="1" colspan="1">≥200</td>
<td align="left" rowspan="1" colspan="1">≥120</td>
<td align="left" rowspan="1" colspan="1">≥3600</td>
<td align="left" rowspan="1" colspan="1">≥3000</td>
<td align="left" rowspan="1" colspan="1">≥3200</td>
<td align="left" rowspan="1" colspan="1">≥2900</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">D</td>
<td align="left" rowspan="1" colspan="1">≥210</td>
<td align="left" rowspan="1" colspan="1">≥240</td>
<td align="left" rowspan="1" colspan="1">≥165</td>
<td align="left" rowspan="1" colspan="1">≥4500</td>
<td align="left" rowspan="1" colspan="1">≥3600</td>
<td align="left" rowspan="1" colspan="1">≥4000</td>
<td align="left" rowspan="1" colspan="1">≥3200</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">E</td>
<td align="left" rowspan="1" colspan="1">≥240</td>
<td align="left" rowspan="1" colspan="1">≥270</td>
<td align="left" rowspan="1" colspan="1">≥180</td>
<td align="left" rowspan="1" colspan="1">≥4800</td>
<td align="left" rowspan="1" colspan="1">≥4500</td>
<td align="left" rowspan="1" colspan="1">≥4300</td>
<td align="left" rowspan="1" colspan="1">≥4000</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">F</td>
<td align="left" rowspan="1" colspan="1">≥300</td>
<td align="left" rowspan="1" colspan="1">≥300</td>
<td align="left" rowspan="1" colspan="1">≥240</td>
<td align="left" rowspan="1" colspan="1">≥5400</td>
<td align="left" rowspan="1" colspan="1">≥5100</td>
<td align="left" rowspan="1" colspan="1">≥5100</td>
<td align="left" rowspan="1" colspan="1">≥4800</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">G</td>
<td align="left" rowspan="1" colspan="1">≥330</td>
<td align="left" rowspan="1" colspan="1">≥330</td>
<td align="left" rowspan="1" colspan="1">≥270</td>
<td align="left" rowspan="1" colspan="1">≥5700</td>
<td align="left" rowspan="1" colspan="1">≥5500</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">H</td>
<td align="left" rowspan="1" colspan="1">≥360</td>
<td align="left" rowspan="1" colspan="1">≥360</td>
<td align="left" rowspan="1" colspan="1">≥330</td>
<td align="left" rowspan="1" colspan="1">≥7200</td>
<td align="left" rowspan="1" colspan="1">≥7000</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>To reveal the impact of water conditions on the multiple cropping systems, two scenarios (irrigated and rain-fed) were considered for the calculation of the potential multiple cropping systems. The rain-fed scenarios <xref ref-type="bibr" rid="pone.0080990-FAO1">[29]</xref>
were obtained from the light-temperature-water condition, while the irrigated scenarios <xref ref-type="bibr" rid="pone.0080990-FAO2">[30]</xref>
were calculated using only the light-temperature conditions, which assumed sufficient water for crop growth. The two scenarios were consistent with two existing agricultural management methods in China (naturally rain-fed farmland and irrigated farmland). The water-limited condition referred to the LGP and TS-G mentioned previously with an additional requirement that indicated moisture was insufficient (ETa <0.4 ET<sub>0</sub>
) <xref ref-type="bibr" rid="pone.0080990-FAO3">[31]</xref>
. According to the above criteria (<xref ref-type="table" rid="pone-0080990-t001">Table 1</xref>
, <xref ref-type="table" rid="pone-0080990-t002">Table 2</xref>
), the following eight zones were classified and mapped:</p>
<list list-type="order"><list-item><p>A. <italic>Zone of no cropping</italic>
(too cold or too dry for rain-fed crops)</p>
</list-item>
<list-item><p>B. <italic>Zone of single cropping</italic>
</p>
</list-item>
<list-item><p>C. <italic>Zone of limited double cropping</italic>
(relay cropping; single wetland rice cropping may be possible)</p>
</list-item>
<list-item><p>D. <italic>Zone of double cropping</italic>
(sequential cropping; double cropping with wetland rice cropping not possible)</p>
</list-item>
<list-item><p>E. <italic>Zone of double cropping</italic>
(sequential cropping; wetland rice cropping possible)</p>
</list-item>
<list-item><p>F. <italic>Zone of limited triple cropping</italic>
(partial relay cropping; no third cropping possible in case of two wetland rice crops)</p>
</list-item>
<list-item><p>G. <italic>Zone of triple cropping</italic>
(sequential cropping of three short-cycle crops; two wetland rice crops possible)</p>
</list-item>
<list-item><p>H. <italic>Zone of triple rice cropping</italic>
(sequential cropping of three wetland rice crops possible)</p>
</list-item>
</list>
</sec>
<sec id="s2b2"><title>Latitudinal thermal classification</title>
<p>The latitudinal thermal classification was the basis for the calculation of potential multiple cropping systems (<xref ref-type="table" rid="pone-0080990-t003">Table 3</xref>
). The delineation of a latitudinal thermal regime was based on the total affected accumulation temperature above 10°C (T<sub>sum10</sub>
), which accumulates the daily average temperatures for days in which the daily average temperature is above the appropriate threshold temperature of 10°C <xref ref-type="bibr" rid="pone.0080990-Bingwei1">[34]</xref>
.</p>
<table-wrap id="pone-0080990-t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.t003</object-id>
<label>Table 3</label>
<caption><title>Classification of latitudinal thermal climates.</title>
</caption>
<alternatives><graphic id="pone-0080990-t003-3" xlink:href="pone.0080990.t003"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Latitudinal thermal regime</td>
<td align="left" rowspan="1" colspan="1">Condition</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">Tropics</td>
<td align="left" rowspan="1" colspan="1">Tsum<sub>10</sub>
>8000</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Subtropics</td>
<td align="left" rowspan="1" colspan="1">4500< Tsum<sub>10</sub>
<8000</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Temperate</td>
<td align="left" rowspan="1" colspan="1">1500< Tsum<sub>10</sub>
<4500</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Boreal</td>
<td align="left" rowspan="1" colspan="1">0< Tsum<sub>10</sub>
<1500</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Arctic</td>
<td align="left" rowspan="1" colspan="1">Tsum<sub>10</sub>
= 0</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</sec>
<sec id="s2b3"><title>Evapotranspiration</title>
<p>The reference evapotranspiration (ET<sub>0</sub>
,mm/day) and actual evapotranspiration (ET<sub>a</sub>
, mm/day) values were calculated using the soil-water balance calculation procedures following the methodologies outlined in “CROPWAT” <xref ref-type="bibr" rid="pone.0080990-FAO3">[31]</xref>
and “Crop Evapotranspiration” <xref ref-type="bibr" rid="pone.0080990-FAO4">[35]</xref>
for crop water requirements and growing period suitability.</p>
<p>The reference evapotranspiration (ET<sub>0</sub>
, mm/day) represents the evapotranspiration from a defined reference surface which is similar to an extensive surface of green and well-watered grass of uniform height (0.12 meters) that is actively growing and completely shades the ground. The ET<sub>0</sub>
was calculated according to the Penman-Monteith equation <xref ref-type="bibr" rid="pone.0080990-FAO3">[31]</xref>
, <xref ref-type="bibr" rid="pone.0080990-Monteith1">[36]</xref>
, <xref ref-type="bibr" rid="pone.0080990-Monteith2">[37]</xref>
.</p>
<p>The calculation of ET<sub>m</sub>
for a ‘reference crop’ is based on the assumption that sufficient water is available for uptake in the rooting zone. The value of ET<sub>m</sub>
is related to ET<sub>0</sub>
by applying the crop coefficients for the water requirements (<italic>kc</italic>
, fractional) (<xref ref-type="disp-formula" rid="pone.0080990.e001">Equation 1</xref>
). The <italic>kc</italic>
factors are determined by phenological development and leaf area <xref ref-type="bibr" rid="pone.0080990-FAO4">[35]</xref>
. <disp-formula id="pone.0080990.e001"><graphic xlink:href="pone.0080990.e001"></graphic>
<label>(1)</label>
</disp-formula>
</p>
<p>The actual uptake of water for the ‘reference’ crop was characterized by the actual evapotranspiration (ET<sub>a</sub>
, mm/day). The calculation of ET<sub>a</sub>
differentiated between two possible cases depending on the availability of water for plant extraction: (i) adequate soil water availability (ET<sub>a</sub>
= ET<sub>m</sub>
) and (ii) limited soil-water availability (ET<sub>a</sub>
<sub>m</sub>
) <xref ref-type="bibr" rid="pone.0080990-FAO4">[35]</xref>
.</p>
<p>Under adequate soil-water conditions, the value of ET<sub>a</sub>
was assumed to be equal to ET<sub>m</sub>
as long as the water balance (Wb, mm) was greater than or equal to the “readily” available soil-water (W<sub>read</sub>
, mm). This requirement characterizes a situation in which crops are able to “easily” extract sufficient water and therefore no water stress occurs (<xref ref-type="disp-formula" rid="pone.0080990.e002">Equation 2</xref>
).<disp-formula id="pone.0080990.e002"><graphic xlink:href="pone.0080990.e002"></graphic>
<label>(2)</label>
</disp-formula>
when <inline-formula><inline-graphic xlink:href="pone.0080990.e003.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pone.0080990.e004.jpg"></inline-graphic>
</inline-formula>
</p>
<p>For limited water conditions, the ET<sub>a</sub>
can be calculated as the product of ET<sub>m</sub>
and the variable ρ (<xref ref-type="disp-formula" rid="pone.0080990.e005">Equation 3</xref>
).<disp-formula id="pone.0080990.e005"><graphic xlink:href="pone.0080990.e005"></graphic>
<label>(3)</label>
</disp-formula>
</p>
<p>The <italic>pre</italic>
represents the precipitation, and ρ is the quotient of the current water balance (W<italic>b</italic>
, mm) and the readily available soil water (W<italic><sub>read</sub>
</italic>
, mm). <disp-formula id="pone.0080990.e006"><graphic xlink:href="pone.0080990.e006"></graphic>
<label>(4)</label>
</disp-formula>
</p>
</sec>
<sec id="s2b4"><title>Soil-water balance</title>
<p>The volume of water available for plant uptake is calculated using a daily soil-water balance (W<italic>b</italic>
, mm). The W<italic>b</italic>
accounts for the accumulated daily water inflow from precipitation (pre, mm/day) or snowmelt (snm, mm/day) and the outflow from actual evapotranspiration (ETa, mm/day) and excess water loss due to runoff and deep percolation.<disp-formula id="pone.0080990.e007"><graphic xlink:href="pone.0080990.e007"></graphic>
<label>(5)</label>
</disp-formula>
</p>
<p>where j is the day of the year, and W<sub>max</sub>
(W<sub>max</sub>
, mm) is the maximum soil water storage capacity.</p>
</sec>
</sec>
</sec>
<sec id="s3"><title>Results and Analysis</title>
<sec id="s3a"><title>Comparison and Verification</title>
<p>To verify the accuracy of the calculated results, the calculated potential multiple cropping system in 2000 was compared with the actual result from remote sensing monitoring based on MODIS data <xref ref-type="bibr" rid="pone.0080990-Huimin1">[11]</xref>
. Three types of cropping systems were used in the actual situation. The potential multiple cropping system was also classified into three corresponding types: single cropping per year, double cropping per year, and triple cropping per year. <xref ref-type="fig" rid="pone-0080990-g002">Figure 2</xref>
demonstrates that the potential multiple cropping system was consistent with the actual system overall. For example, regions that use triple cropping (the area of 93,445 square kilometers) exist in Sichuan Basin, the Middle-Lower Yangtze Plain, the south of the Yangtze River and Hainan Island are completely simulated as triple cropping in the potential cropping system. Similarly, the single cropping in the potential cropping system (the area of 870,093 square kilometers) duplicated the actual single cropping regions. However, as shown in <xref ref-type="table" rid="pone-0080990-t004">Table 4</xref>
, the simulated MCI was larger than the actual MCI, which revealed the difference between the calculated potential multiple cropping systems and the actual systems. Farmers tended to select their multiple cropping practices following the traditional pattern which was based on previous climate conditions rather than the optimal multiple cropping systems (which are consistent with the actual local climate conditions). Moreover, many farmers preferred to grow plants that yielded greater economic profits, and the lack of young and middle-aged rural laborers who prefer to work in cities further reduced the MCI.</p>
<fig id="pone-0080990-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g002</object-id>
<label>Figure 2</label>
<caption><title>Comparison between potential and actual multiple cropping systems.</title>
</caption>
<graphic xlink:href="pone.0080990.g002"></graphic>
</fig>
<table-wrap id="pone-0080990-t004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.t004</object-id>
<label>Table 4</label>
<caption><title>Comparison table between potential and actual multiple cropping systems (unit: km<sup>2</sup>
).</title>
</caption>
<alternatives><graphic id="pone-0080990-t004-4" xlink:href="pone.0080990.t004"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">actual</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">single</td>
<td align="left" rowspan="1" colspan="1">double</td>
<td align="left" rowspan="1" colspan="1">triple</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">single</td>
<td align="left" rowspan="1" colspan="1">870,093</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">potential</td>
<td align="left" rowspan="1" colspan="1">double</td>
<td align="left" rowspan="1" colspan="1">981,026</td>
<td align="left" rowspan="1" colspan="1">301,043</td>
<td align="left" rowspan="1" colspan="1">0</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">triple</td>
<td align="left" rowspan="1" colspan="1">629,786</td>
<td align="left" rowspan="1" colspan="1">674,971</td>
<td align="left" rowspan="1" colspan="1">93,445</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</sec>
<sec id="s3b"><title>Spatial patterns of the potential multiple cropping system in China</title>
<p>Based on the method in Section 2.2, the potential multiple cropping systems were calculated using the average climatic conditions at intervals of 10 years from the 1960s to the 2000s. <xref ref-type="fig" rid="pone-0080990-g003">Figure 3</xref>
and <xref ref-type="fig" rid="pone-0080990-g004">Figure 4</xref>
show the distribution maps of the multiple cropping systems under the rain-fed and irrigated scenarios during each 10-year period.</p>
<fig id="pone-0080990-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g003</object-id>
<label>Figure 3</label>
<caption><title>Distribution of multiple cropping systems under the rain-fed scenario.</title>
</caption>
<graphic xlink:href="pone.0080990.g003"></graphic>
</fig>
<fig id="pone-0080990-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g004</object-id>
<label>Figure 4</label>
<caption><title>Distribution of multiple cropping systems under the irrigated scenario.</title>
</caption>
<graphic xlink:href="pone.0080990.g004"></graphic>
</fig>
<p>In the 2000s, the complexity of the potential cropping system pattern and the spatial pattern of the potential multiple cropping systems increased, as manifested from north to south and from west to east (<xref ref-type="fig" rid="pone-0080990-g003">Figure 3</xref>
, <xref ref-type="fig" rid="pone-0080990-g004">Figure 4</xref>
). Due to the relatively high temperature and abundant rainfall, the PMCI in the southeastern areas was high. The PMCI declined with increasing distance from the coastal line. In the northwestern portion of the country, potential multiple cropping was constrained by low temperatures and rainfall.</p>
<p>Under the irrigated scenario, the triple cropping zones (including limited triple cropping, triple cropping with rice and triple rice cropping) were primarily distributed over South China and the middle and lower reaches of the Yangtze River, which cover 36.39% of the total land area of China (<xref ref-type="fig" rid="pone-0080990-g005">Figure 5</xref>
). The double cropping zones (including limited double cropping, double cropping and double cropping with rice) were primarily spread across the Loess Plateau, the Huang-Huai-Hai plain and most of the Gansu-Xinjiang regions, which occupied 32.49% of the total land area of China. The single cropping zones were primarily distributed in the north of Inner Mongolia and along the Great Wall region (denoted as f) and a small portion of the Qinghai-Tibet Plateau, which occupy 31.11% of the total land area of China.</p>
<fig id="pone-0080990-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g005</object-id>
<label>Figure 5</label>
<caption><title>Area ratio difference of the multiple cropping systems under two scenarios in the 2000s.</title>
</caption>
<graphic xlink:href="pone.0080990.g005"></graphic>
</fig>
<p>Under the rain-fed scenario, all multiple cropping zones were shifted to the south and the east, and the areas for the double and triple cropping zones were decreased by 7.84% and 8.95% respectively (<xref ref-type="fig" rid="pone-0080990-g005">Figure 5</xref>
). The limited rainfall led to a decrease of the PMCI from 205% to 181%. The data showed that the irrigated condition was crucial for multiple cropping in China. Therefore, improving the construction of water facilities for farmland will benefit food productivity in China.</p>
</sec>
<sec id="s3c"><title>Changes in the potential multiple cropping systems</title>
<p>The potential multiple cropping systems have undergone significant changes between the 1960s -and the 2000s. The fundamental characteristics of the changes in the potential multiple cropping systems were that the single cropping area decreased and the triple cropping area increased. Moreover, the PMCI showed an increasing trend, with the changes primarily occurring after the 1980s.</p>
<p>Under rain-fed scenarios from the 1960s to the 2000s, the single cropping area decreased by 3.23% (approximately 173,900 km<sup>2</sup>
), whereas the triple cropping area increased by 3.21% (approximately 173,000 km<sup>2</sup>
) and the PMCI increased by 7% (from 174% to 181%). No significant changes occurred from the 1960s to the 1980s, but shifts occurred in all areas and the PMCI grew rapidly from the 1980s to the 2000s. The major conversion from double cropping to triple cropping and from single cropping to double cropping took place in the time period from the 1980s to the 1990s and from the 1990s to the 2000s respectively. The size of the conversion areas was approximately equal during the two decades (<xref ref-type="fig" rid="pone-0080990-g006">Figure 6</xref>
).</p>
<fig id="pone-0080990-g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g006</object-id>
<label>Figure 6</label>
<caption><title>Percentage changes in the multiple cropping area and PMCI under two scenarios from the 1960s to the 2000s.</title>
</caption>
<graphic xlink:href="pone.0080990.g006"></graphic>
</fig>
<p>Under the irrigated scenarios, the single cropping area declined by 7.28% (approximately 391,600 km<sup>2</sup>
), the triple cropping areas increased by 6.38% (approximately 343,300 km<sup>2</sup>
), and the PMCI increased by 13% (from 192% to 205%). Major changes occurred from the 1980s to the 2000s, and the PMCI grew much more rapidly. During this period, the increase in the triple cropping area was nearly equivalent to the decrease in single cropping area, whereas the double cropping area remained unchanged as a whole(<xref ref-type="fig" rid="pone-0080990-g006">Figure 6</xref>
).</p>
</sec>
<sec id="s3d"><title>Response of potential multiple cropping systems to climate change</title>
<p>As presented in <xref ref-type="fig" rid="pone-0080990-g007">figure 7(a, b, c)</xref>
, the climate change from the 1960s to the 2000s showed that the annual mean temperature, precipitation and radiation data take on different spatial change patterns. The annual mean temperature generally increased across the major cultivation regions. The line chart from the 1960s to the 2000s presented in <xref ref-type="fig" rid="pone-0080990-g008">Figure 8</xref>
reflects an obvious upward trend in the temperature change. More specifically, the temperature increased by more than 2°C in parts of the northeast, the north and the northwest of China (<xref ref-type="fig" rid="pone-0080990-g007">Figure 7a</xref>
). The rate of temperature change after the 1980s (0.41°C/decade) increased more dramatically than prior to the 1980s (0.11°C/decade).</p>
<fig id="pone-0080990-g007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g007</object-id>
<label>Figure 7</label>
<caption><title>Spatial distribution pattern of climate changes from the 1960s to the 2000s in China.</title>
</caption>
<graphic xlink:href="pone.0080990.g007"></graphic>
</fig>
<fig id="pone-0080990-g008" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g008</object-id>
<label>Figure 8</label>
<caption><title>Climate changes from the 1960s to the 2000s in China.</title>
</caption>
<graphic xlink:href="pone.0080990.g008"></graphic>
</fig>
<p>Precipitation was reduced in most regions during the entire period, especially in Southwest China, the Loess Plateau, Inner Mongolia and along the Great Wall region where the precipitation decreased by more than 10%. The areas where precipitation increased were primarily distributed in the northwest of China and the middle and lower reaches of the Yangtze River Region (<xref ref-type="fig" rid="pone-0080990-g007">Figure 7b</xref>
). As an illustration, the precipitation in the Gansu-Xinjiang and Southern China region increased by 13.7% (approximately 21 mm) and 2% (approximately 31 mm), respectively. In contrast, the temporal change process showed that precipitation did not change significantly before the 1990s but took on an obvious decreasing trend (31 mm) from the 1990s to the 2000s.</p>
<p>Solar radiation decreased in large areas across Huang-Huai-Hai plain and most of southern China but increased in parts of Northeastern China, the Loess Plateau and the western area of Southern China (<xref ref-type="fig" rid="pone-0080990-g007">Figure 7c</xref>
).The amount of radiation dropped at an average rate of 182 MJ/m<sup>2</sup>
·decade from the 1960s to the 1980s and at an average rate of approximately 22 MJ/m<sup>2</sup>
·decade from the 1980s to the 2010s.</p>
<p>The obvious spatial heterogeneity of the climate changes from the 1960s to the 2000s has an important influence on the potential multiple cropping systems in China. Temperature and precipitation were two leading factors that affected the changes of the PMCI, and <xref ref-type="fig" rid="pone-0080990-g008">Figure 8</xref>
and <xref ref-type="fig" rid="pone-0080990-g009">Figure 9</xref>
show that climate change increased the PMCI under the irrigated and rain-fed scenarios by 13% and 7% respectively in China. The increase in the PMCI was mostly due to the growth of the annual mean temperature. The gap of PMCI between rain-fed scenario and irrigated scenario was 18% in the 1960s, while this gap had increased to 24% in the 2000s.The reason for the increased gap was that the crops could not get enough water due to the reduction in rainfall, which limited the increase of the PMCI under rain-fed conditions. Consequently, in China, moisture was one of the most significant factors for the PMCI. The moisture limit will be the main factor restraining the PMCI under rain-fed scenarios in the present and perhaps also in the future.</p>
<fig id="pone-0080990-g009" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080990.g009</object-id>
<label>Figure 9</label>
<caption><title>The relationship between the PMCI and climate changes in China from the 1960s to the 2000s.</title>
</caption>
<graphic xlink:href="pone.0080990.g009"></graphic>
</fig>
<p>In different regions of China, the responses of the potential multiple cropping systems to climate change took on different characteristics. In Northeast China (a), the PMCI continuously increased in two scenarios with rising annual mean temperature and declining radiation throughout the entire period. The gap in the PMCI under the two scenarios remained nearly the same, even when the precipitation sharply dropped by more than 50 mm from the 1980s to the 2000s in the rain-fed scenario. This result indicated that heat is a major factor in limiting the potential multiple cropping systems, whereas precipitation had a lesser influence.</p>
<p>In the Gansu-Xinjiang region (b), the potential multiple cropping systems under the rain-fed scenario was always one crop per year, whereas the PMCI under the irrigated condition rose by 0.178. This result shows that the temperature increase had a positive impact on the PMCI under the irrigated condition. The largest gap of the PMCI under the two scenarios was approximately 0.7. However, a small increase in precipitation did not cause the increase in the PMCI under the rain-fed condition due to the extreme deficiency of rainfall.</p>
<p>In Southern China (c), the PMCI under the irrigated scenario was close to 3, which is primarily due to the effect of abundant heat on crop growth. The PMCI under the rain-fed scenario took on an obvious decreasing trend, and the PMCI gap under the two scenarios fluctuated between 0.26 and 0.5. This result was primarily due to the highly uneven rainfall.</p>
<p>In Huang-Huai-Hai plain (d), the annual mean temperature increased by 1.5°C from the 1960s to the 2000s, which led to an increase in the PMCI of 0.103 under the irrigated scenario. Under the rain-fed scenarios, the PMCI increased by 0.266, and displayed a fluctuation that first increased, subsequently decreased, and finally increased again. The change in the PMCI was primarily due to the sharp decrease in radiation by more than 500 MJ/m<sup>2</sup>
·year and the increase in precipitation during certain periods. The gap between the two scenarios was approximately 0.35, which indicates that water restricts crop growth in this region to a certain extent.</p>
<p>In the Loess Plateau (e), Inner Mongolia and along the Great Wall region (f), the annual mean temperature and PMCI under irrigated conditions showed a simultaneous growth trend, whereas the PMCI under the rain-fed condition changed little, resulting in larger gaps. A potential reason for this finding was that the declining precipitation caused a moisture shortage for the crops.</p>
<p>The PMCI of the Qinghai-Tibet Plateau (h) under the two scenarios were both lower than those of other regions, which rose slightly with increasing temperatures, and the gaps were rarely affected by a reduction in precipitation. The main reason for this appeared to be that the temperature was too low (below 5°C) for normal crop growth.</p>
<p>In Southwest China (i) and the middle and lower reaches of the Yangtze River (j), the annual mean temperature and PMCI under the two scenarios showed a concurrent growth trend. Although the precipitation fluctuated, the PMCI gap under the two scenarios was almost unchanged. This result shows that these areas received sufficient rain for plant growth under the prevailing heat condition, which implies that these areas are the most suitable zones for plant growth in China.</p>
</sec>
</sec>
<sec id="s4"><title>Conclusions and Discussions</title>
<p>The multiple cropping systems demonstrate the efficient use of water, soil, light energy, and other natural resources <xref ref-type="bibr" rid="pone.0080990-Jinlong2">[38]</xref>
, all of which are important for agricultural production and food security. Wherever the natural conditions allow, it is necessary to adopt high intensity cultivation to secure efficient food supply in China <xref ref-type="bibr" rid="pone.0080990-Hong1">[2]</xref>
.</p>
<p>This study focused on the calculation and analysis of potential multiple cropping systems in China based on multi-source data. We also analyzed the changes in the potential multiple cropping systems in response to climate change in China from the 1960s to the 2010s. Conclusions from this work are as follows <xref ref-type="bibr" rid="pone.0080990-Li1">[39]</xref>
–<xref ref-type="bibr" rid="pone.0080990-ZhaoD1">[42]</xref>
:</p>
<p>(1) The spatial pattern of potential multiple cropping systems in China displays tremendous heterogeneity, and an increasing complexity from northern China to southern China, and from western China to eastern China. The decrease in the single cropping area, and the increase in the triple cropping area, are the fundamental characteristics of the potential multiple cropping systems change from the 1960s to the 2000s. The rate of increase of PMCI after the 1980s was far larger than that before the 1980s. During the studied period, the PMCI gap between rain-fed and irrigated scenarios increased from 18% to 24%, which indicated noticeable growth of water supply limitations under the rain-fed scenario. Irrigation conditions were important for multiple cropping practices, underlining the importance of the construction of water facilities would be beneficial for food production in China.</p>
<p>(2) The obvious spatial heterogeneity of climate change from the 1960s to the 2000s had a significant influence on the potential multiple cropping systems in China. The climate change caused the PMCI under the irrigated and rain-fed scenarios to increase by 13% and 7% respectively across the whole of China. Furthermore, the growth of the annual mean temperature was identified as the main reason for the increase of the PMCI.</p>
<p>(3) In the different regions of China, the response of the potential multiple cropping systems to climate change took on different characteristics. In Southern China, the PMCI under the rain-fed scenario took on an obvious decreasing trend because of highly uneven rainfall, although precipitation was abundant. In the Northeast and Qinghai-Tibet Plateau, the heat severely restricted crop growth, and in Gansu-Xinjiang region, Loess Plateau, Inner Mongolia and along the Great Wall region, precipitation was the major factor that limited the potential multiple cropping systems. In Southwest China and the middle and lower reaches of the Yangtze River, the PMCI showed a simultaneous growth trend under both the irrigation and rain-fed scenarios, although precipitation fluctuated.</p>
<p>This study calculated the potential multiple cropping systems based on meteorological data, and subsequently analyzed the relationships between the potential multiple cropping systems and climate change, which assist in dealing with the issue of food security in China. However, certain limitations exist in this research. For instance, (1) due to the lack of irrigation data, we were not able to simulate the actual multiple cropping systems according to real irrigation conditions, instead we simulated two different scenarios (rain-fed and irrigated scenarios); (2) it was assumed that the irrigated scenarios would supply sufficient water for crop growth, but in practice the availability of water under the irrigation scenario still poses limitations to crop growth; (3) extreme weather conditions (i.e., freezing or high temperatures, rainstorms and heavy snow) may have an extreme effect on the multiple cropping systems, but these factors were not taken into consideration in our research.</p>
</sec>
</body>
<back><ref-list><title>References</title>
<ref id="pone.0080990-Ziqiang1"><label>1</label>
<mixed-citation publication-type="journal"><name><surname>Ziqiang</surname>
<given-names>F</given-names>
</name>
, <name><surname>Yunlong</surname>
<given-names>C</given-names>
</name>
, <name><surname>Youxiao</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Erfu</surname>
<given-names>D</given-names>
</name>
(<year>2001</year>
) <article-title>Research on the relationship of cultivated land change and food security in China</article-title>
. <source>Journal of Natural Resources</source>
<volume>16</volume>
: <fpage>313</fpage>
–<lpage>319</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Hong1"><label>2</label>
<mixed-citation publication-type="journal"><name><surname>Hong</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Xiubin</surname>
<given-names>Y</given-names>
</name>
(<year>2000</year>
) <article-title>Cultivated land and food supply in China</article-title>
. <source>Land Use Policy</source>
<volume>17</volume>
: <fpage>73</fpage>
–<lpage>88</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Von1"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Von</surname>
<given-names>B</given-names>
</name>
(<year>2007</year>
) <article-title>Rising world food prices: impact on the poor</article-title>
. <source>The International Journal for Rural Development (Rural 21)</source>
<volume>42</volume>
: <fpage>19</fpage>
–<lpage>21</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Guo1"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Guo</surname>
<given-names>B</given-names>
</name>
(<year>1997</year>
) <article-title>The analysis for the change of multiple cropping index in China</article-title>
. <source>Economic Geography</source>
<volume>17</volume>
: <fpage>8</fpage>
–<lpage>13</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Jinshen1"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Jinshen</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Limin</surname>
<given-names>S</given-names>
</name>
(<year>2000</year>
) <article-title>The transition of cropping system in Hebei Province</article-title>
. <source>Acta Agriculture Boreali-sinica</source>
<volume>4</volume>
: <fpage>126</fpage>
–<lpage>130</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Shulan1"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Shulan</surname>
<given-names>J</given-names>
</name>
, <name><surname>Lichun</surname>
<given-names>H</given-names>
</name>
, <name><surname>Lei</surname>
<given-names>X</given-names>
</name>
(<year>2011a</year>
) <article-title>The Relationship of Multiple Cropping Index of Arableland Change and National Food Security in the Middle and Lower Reaches of Yangtze River</article-title>
. <source>Chinese Agricultural Science Bulletin</source>
<volume>27</volume>
: <fpage>208</fpage>
–<lpage>212</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Shulan2"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Shulan</surname>
<given-names>J</given-names>
</name>
, <name><surname>Caiqiu</surname>
<given-names>X</given-names>
</name>
, <name><surname>Huahua</surname>
<given-names>P</given-names>
</name>
(<year>2011b</year>
) <article-title>Change Characteristics and Potential Analysis for Multiple Cropping Index of Cropland in Major Grain Producing Areas in China</article-title>
. <source>Guizhou Agricultural Sciences</source>
<volume>39</volume>
: <fpage>201</fpage>
–<lpage>204</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Jing1"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>Jing</surname>
<given-names>L</given-names>
</name>
, <name><surname>Zhiyuan</surname>
<given-names>R</given-names>
</name>
(<year>2011</year>
) <article-title>The Monitoring for cropping index of arable land in northwest region using SPOTNDVI-A case of Shaanxi Province</article-title>
. <source>Journal of Arid Land Resources and Environment</source>
<volume>25</volume>
: <fpage>86</fpage>
–<lpage>91</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Dailiang1"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Dailiang</surname>
<given-names>P</given-names>
</name>
, <name><surname>Jingfeng</surname>
<given-names>H</given-names>
</name>
, <name><surname>Huimin</surname>
<given-names>J</given-names>
</name>
(<year>2007</year>
) <article-title>Monitoring the Sequential Cropping Index of Arable Land in Zhejiang Province of China Using MODIS-NDVI</article-title>
. <source>Agricultural Sciences in China</source>
<volume>6</volume>
: <fpage>208</fpage>
–<lpage>213</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Zhiguo1"><label>10</label>
<mixed-citation publication-type="other">Zhiguo Z (2011) Spatial-temporal Characteristics of Multiple Cropping Index and Relationship between Multiple Cropping Index and Grain Yield in Henan Province. Hubei Agricultural Sciences: 3653–3656.</mixed-citation>
</ref>
<ref id="pone.0080990-Huimin1"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Huimin</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Jiyuan</surname>
<given-names>L</given-names>
</name>
, <name><surname>Mingkui</surname>
<given-names>C</given-names>
</name>
(<year>2005</year>
) <article-title>Remotely Sensed Multiple Cropping Index Variations in China during 1981-2000</article-title>
. <source>ACTA Geographica SINICA</source>
<volume>60</volume>
: <fpage>559</fpage>
–<lpage>566</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Hongping1"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Hongping</surname>
<given-names>D</given-names>
</name>
(<year>2001</year>
) <article-title>Analysis of the development of Hunan farming system from 1949 to 1998</article-title>
. <source>Cropping and Planting</source>
<volume>3</volume>
: <fpage>1</fpage>
–<lpage>4</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Shumin1"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Shumin</surname>
<given-names>L</given-names>
</name>
(<year>2007</year>
) <article-title>Research on the potential of the multiple cropping index and the behavior of multiple cropping of China</article-title>
. <source>Problems of Agricultural Economy</source>
<volume>5</volume>
: <fpage>85</fpage>
–<lpage>90</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Mao1"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Mao</surname>
<given-names>L</given-names>
</name>
(<year>2002</year>
) <article-title>Gray incidence analysis between cultivated area and grain yield and its affecting factors of Henan Province</article-title>
. <source>Progress in Geography</source>
<volume>21</volume>
: <fpage>163</fpage>
–<lpage>172</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Xin1"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Xin</surname>
<given-names>J</given-names>
</name>
, <name><surname>Yua</surname>
<given-names>ZR</given-names>
</name>
, <name><surname>Van</surname>
<given-names>Leeuwen</given-names>
</name>
, <name><surname>Driessen</surname>
<given-names>PM</given-names>
</name>
(<year>2002</year>
) <article-title>Mapping crop key phonological stages in the North China Plain using NOAA time series images. Int. J.Appl. Earth Obs</article-title>
. <source>Geoinf</source>
<volume>4</volume>
: <fpage>109</fpage>
–<lpage>117</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Huimin2"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Huimin</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Yuling</surname>
<given-names>F</given-names>
</name>
, <name><surname>Xiangming</surname>
<given-names>X</given-names>
</name>
, <name><surname>Heqing</surname>
<given-names>H</given-names>
</name>
, <name><surname>Honglin</surname>
<given-names>H</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Modeling gross primary productivity for winter wheat–maize double cropping system using MODIS time series and CO<sub>2</sub>
eddy flux tower data</article-title>
. <source>Agriculture Ecosystems and Environment</source>
<volume>129</volume>
: <fpage>391</fpage>
–<lpage>400</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Jinlong1"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Jinlong</surname>
<given-names>F</given-names>
</name>
, <name><surname>Bingfang</surname>
<given-names>W</given-names>
</name>
(<year>2004a</year>
) <article-title>A study on cropping index potential based on GIS</article-title>
. <source>Journal of Remote Sensing</source>
<volume>8</volume>
: <fpage>637</fpage>
–<lpage>644</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Xiaolin1"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Xiaolin</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Qiang</surname>
<given-names>L</given-names>
</name>
, <name><surname>Miaogen</surname>
<given-names>S</given-names>
</name>
, <name><surname>Jin</surname>
<given-names>C</given-names>
</name>
, <name><surname>Jin</surname>
<given-names>W</given-names>
</name>
(<year>2008</year>
) <article-title>A methodology for multiple cropping index extraction based on NDVI time series</article-title>
. <source>Journal of Natural Resources</source>
<volume>23</volume>
: <fpage>534</fpage>
–<lpage>544</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Ashton1"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Ashton</surname>
<given-names>S</given-names>
</name>
, <name><surname>Joseph</surname>
<given-names>M</given-names>
</name>
(<year>2011</year>
) <article-title>Spatial structure and landscape associations of SRTM error</article-title>
. <source>Remote Sensing of Environment</source>
<volume>115</volume>
: <fpage>1576</fpage>
–<lpage>1587</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Zandbergen1"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Zandbergen</surname>
<given-names>P</given-names>
</name>
(<year>2008</year>
) <article-title>Applications of Shuttle Radar Topography Mission Elevation Data</article-title>
. <source>Geography Compass</source>
<volume>2</volume>
: <fpage>1404</fpage>
–<lpage>1431</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Farr1"><label>21</label>
<mixed-citation publication-type="journal">Farr TG, Rosen PA, Caro E, Crippen R, Duren R, <etal>et al</etal>
. (2007) The Shuttle Radar Topography mission. Reviews of Geophysics <volume>45</volume>
.</mixed-citation>
</ref>
<ref id="pone.0080990-Jiyuan1"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Jiyuan</surname>
<given-names>L</given-names>
</name>
, <name><surname>Mingliang</surname>
<given-names>L</given-names>
</name>
, <name><surname>Dafang</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Zengxiang</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Xiangzheng</surname>
<given-names>D</given-names>
</name>
, <etal>et al</etal>
(<year>2003</year>
) <article-title>Study on spatial pattern of land-use change in China during 1995–2000</article-title>
. <source>Science in China (Series D)</source>
<volume>46</volume>
: <fpage>373</fpage>
–<lpage>384</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Jiyuan2"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>Jiyuan</surname>
<given-names>L</given-names>
</name>
, <name><surname>Mingliang</surname>
<given-names>L</given-names>
</name>
, <name><surname>Hanqin</surname>
<given-names>T</given-names>
</name>
, <name><surname>Dafang</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Zengxiang</surname>
<given-names>Z</given-names>
</name>
, <etal>et al</etal>
(<year>2005a</year>
) <article-title>Spatial and temporal patterns of China's cropland during 1990–2000</article-title>
. <source>Remote Sensing of Environment</source>
<volume>98</volume>
: <fpage>442</fpage>
–<lpage>456</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Jiyuan3"><label>24</label>
<mixed-citation publication-type="journal">Jiyuan L, Tian H, Mingliang L, Dafang Z, Jerry M, <etal>et al</etal>
. (2005b) China's changing landscape during the 1990s: Large-scale land transformations estimated with satellite data. Geophysical Research Letters <volume>32</volume>
.</mixed-citation>
</ref>
<ref id="pone.0080990-Jiyuan4"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Jiyuan</surname>
<given-names>L</given-names>
</name>
, <name><surname>Zengxiang</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Xinliang</surname>
<given-names>X</given-names>
</name>
, <name><surname>Wenhui</surname>
<given-names>K</given-names>
</name>
, <name><surname>Wancun</surname>
<given-names>Z</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Spatial patterns and driving forces of Land use change in China during the early 21st century</article-title>
. <source>Journal of Geographical Sciences</source>
<volume>20</volume>
: <fpage>483</fpage>
–<lpage>494</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Jiyuan5"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Jiyuan</surname>
<given-names>L</given-names>
</name>
, <name><surname>Zhang</surname>
<given-names>Q</given-names>
</name>
, <name><surname>Hu</surname>
<given-names>Y</given-names>
</name>
(<year>2012</year>
) <article-title>Regional differences of China's urban expansion from late 20th to early 21st century based on remote sensing information</article-title>
. <source>Chinese Geographical Science</source>
<volume>22</volume>
: <fpage>1</fpage>
–<lpage>14</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Hutchinson1"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Hutchinson</surname>
<given-names>MF</given-names>
</name>
(<year>1998a</year>
) <article-title>Interpolation of rainfall data with thin plate smoothing splines I. Two dimensional smoothing of data with short range correlation</article-title>
. <source>Geographic Information Decision Analysis</source>
<volume>2</volume>
: <fpage>153</fpage>
–<lpage>167</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Hutchinson2"><label>28</label>
<mixed-citation publication-type="journal"><name><surname>Hutchinson</surname>
<given-names>MF</given-names>
</name>
(<year>1998b</year>
) <article-title>Interpolation of rainfall data with thin plate smoothing splines II. Analysis of topographic dependence</article-title>
. <source>Geographic Information Decision Analysis</source>
<volume>2</volume>
: <fpage>168</fpage>
–<lpage>185</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-FAO1"><label>29</label>
<mixed-citation publication-type="other">FAO (1984) Guidelines: Land Evaluation for Rain-fed Agriculture. FAO Soils Bulletin. 52.</mixed-citation>
</ref>
<ref id="pone.0080990-FAO2"><label>30</label>
<mixed-citation publication-type="other">FAO (1985) Guidelines: Land Evaluation for Irrigated Agriculture. FAO Soils Bulletin. 55.</mixed-citation>
</ref>
<ref id="pone.0080990-FAO3"><label>31</label>
<mixed-citation publication-type="other">FAO (1992) CROPWAT: A Computer Program for Irrigation Planning and Management. FAO Irrigation and Drainage Paper no 46. Land and Water Development Division, Rome, Italy.</mixed-citation>
</ref>
<ref id="pone.0080990-Fischer1"><label>32</label>
<mixed-citation publication-type="other">Fischer G, M Shah, Van Velthuizen H, Nachtergaele O. (2002a) Global Agro-Ecological Assessment for Agriculture in the 21st Century: Methodology and Results. IIASA RR-02-02, IIASA, Laxenburg, Austria.</mixed-citation>
</ref>
<ref id="pone.0080990-Fischer2"><label>33</label>
<mixed-citation publication-type="other">Fischer G, M Shah, Van Velthuizen H. (2002b) Climate Change and Agricultural Vulnerability. Special Report as contribution to the World Summit on Sustainable Development, Johannesburg 2002. International Institute for Applied Systems Analysis, Laxenburg, Austria: 150–152.</mixed-citation>
</ref>
<ref id="pone.0080990-Bingwei1"><label>34</label>
<mixed-citation publication-type="journal">Bingwei H (1958) The preliminary draft of comprehensive natural divisions in China. Acta Geographica Sinica <volume>4</volume>
.</mixed-citation>
</ref>
<ref id="pone.0080990-FAO4"><label>35</label>
<mixed-citation publication-type="other">FAO (1998) Crop Evapotranspiration. FAO Irrigation and Drainage Paper no. 56 Rome, Italy.</mixed-citation>
</ref>
<ref id="pone.0080990-Monteith1"><label>36</label>
<mixed-citation publication-type="other">Monteith JL (1965) Evapotranspiration and the environment. In The State and Movement of Water in Living Organisms: 205–234.</mixed-citation>
</ref>
<ref id="pone.0080990-Monteith2"><label>37</label>
<mixed-citation publication-type="journal"><name><surname>Monteith</surname>
<given-names>JL</given-names>
</name>
(<year>1981</year>
) <article-title>Evapotranspiration and surface temperature</article-title>
. <source>Quarterly Journal Royal Meteorological Society</source>
<volume>107</volume>
: <fpage>1</fpage>
–<lpage>27</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Jinlong2"><label>38</label>
<mixed-citation publication-type="journal"><name><surname>Jinlong</surname>
<given-names>F</given-names>
</name>
, <name><surname>Bingfang</surname>
<given-names>W</given-names>
</name>
(<year>2004b</year>
) <article-title>A Methodology for Retrieving Cropping Index from NDVI Profile</article-title>
. <source>Journal of Remote Sensing</source>
<volume>8</volume>
: <fpage>628</fpage>
–<lpage>636</lpage>
</mixed-citation>
</ref>
<ref id="pone.0080990-Li1"><label>39</label>
<mixed-citation publication-type="journal">Li P, Feng Z, Jiang L, Yujie L, Xiangming X (2012). Changes in rice cropping systems in the Poyang Lake Region, China during 2004–2010. Journal of Geographical Sciences, <volume>22(4)</volume>
: , 653–668.</mixed-citation>
</ref>
<ref id="pone.0080990-HoqueM1"><label>40</label>
<mixed-citation publication-type="other">Hoque M Z. (1984). Cropping systems in Asia: on-farm research and management. Irri.</mixed-citation>
</ref>
<ref id="pone.0080990-TianXiang1"><label>41</label>
<mixed-citation publication-type="other">TianXiang Y, Na Z, Douglas R, ChenLiang W, ZeMeng F, <etal>et al</etal>
. (2013). Climate change trend in China, with improved accuracy. Climatic Change, 1–15.</mixed-citation>
</ref>
<ref id="pone.0080990-ZhaoD1"><label>42</label>
<mixed-citation publication-type="journal">Zhao D S, Wu S H, Yin Y H. (2013). Responses of Terrestrial Ecosystems' Net Primary Productivity to Future Regional Climate Change in China. PloS one, <volume>8(4)</volume>
: , e60849.</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>