1
Supplementary Material
2
Agronomy Journal
3
Winter wheat yield gaps and patterns in China
4
Shuang Sun1,2, Xiaoguang Yang1*, Xiaomao Lin2*, Gretchen F. Sassenrath2, and
5
Kenan Li3
6 7
1
8 9
2
10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
College of Resources and Environmental Sciences, China Agricultural University, No.2 Yuanmingyuan West Rd., Haidian District, Beijing 100193, China Department of Agronomy, Kansas State University, 2108 Throckmorton Plant Sciences Center, Manhattan, KS 66506, USA 3
College of Air Traffic Management, Civil Aviation University of China, No.2898 Jinbei Rd., Dongli District, Tianjin 300300, China
*
Corresponding author: Xiaoguang Yang College of Resources and Environmental Sciences China Agricultural University, Beijing 100193, China Tel: +86-10-6273 3939, Fax: +86-10-6273 3939 Email:
[email protected] Xiaomao Lin Department of Agronomy Kansas State University, Kansas 66506, USA Tel: +1 7855326816, Fax: +1 7855326094 Email:
[email protected]
25 26 27
This file includes:
28
Classification of yield trend pattern and data analysis
29
Example of using more than one cultivar in simulation
30
How to use the weighting factors
31
Supplemental Tables S1–S3
32
Supplemental Figure S1
33
References
1
34
Classification of yield trend pattern and data analysis
35
The wheat yield trends were classified for each county according to the global crop
36
yield trends identified by Ray et al. (2012): yields never improved, yields stagnated,
37
yields collapsed, and yields increasing. Regions where no significant yield improvement
38
occurred over the 30-year study period were classified as ‘yields never improved’.
39
Regions in which previous yield improvements have slowed or ceased were classified as
40
‘yields stagnated’. Regions where yields decreased since the 1980s or initially increased
41
and then declined to the 1980s level were classified as ‘yields collapsed’. If yields were
42
still increasing, the regions were classified as ‘yields increasing’ (Ray et al. 2012). For
43
each yield time series for each county, four linear regression models, an intercept-only
44
model (Eq. 1), a linear model (Eq. 2), a quadratic model (Eq. 3), and a cubic model (Eq.
45
4), were fitted for each county:
46 47 48 49
Yield = k
(1)
Yield = at 2 + bt + k
(3)
Yield = at + k
(2)
Yield = at 3 + bt 2 + ct + k
(4)
50
where Yield predictant represents the actual county-average yield in kg ha-1 and t
51
represents year; a, b, and c are the regression coefficients and k is the intercept. The best-
52
fit model for yield data was chosen based on the Akaike Information Criterion (AIC) for
53
each model (Akaike 1974) (Eq. 5); the model with the minimum AIC was chosen as the
54
best for the county evaluated:
55
RSS
AIC=nlog(
n
)+2
(5)
2
56
where RSS denotes residual sum of squares for a sample size of n and
57
parameters in regression.
is the number of
58
The accuracy and significance of model fit was then determined by the F-test with p-
59
value