Next Article in Journal
Unsupervised Saliency Model with Color Markov Chain for Oil Tank Detection
Next Article in Special Issue
Long-Term Monitoring of Cropland Change near Dongting Lake, China, Using the LandTrendr Algorithm with Landsat Imagery
Previous Article in Journal
Automatic Data Selection and Boresight Adjustment of LiDAR Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved CASA Model for Estimating Winter Wheat Yield from Remote Sensing Images

1
Key Laboratory of Quantitative Remote Sensing in Agriculture of Ministry of Agriculture, Beijing Research Center for Information Technology in Agriculture, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China
2
National Engineering Research Center for Agro-Ecological Big Data Analysis & Application, Anhui University, Hefei 230601, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(9), 1088; https://doi.org/10.3390/rs11091088
Submission received: 4 April 2019 / Revised: 27 April 2019 / Accepted: 30 April 2019 / Published: 7 May 2019

Abstract

:
The accurate and timely monitoring and evaluation of the regional grain crop yield is more significant for formulating import and export plans of agricultural products, regulating grain markets and adjusting the planting structure. In this study, an improved Carnegie–Ames–Stanford approach (CASA) model was coupled with time-series satellite remote sensing images to estimate winter wheat yield. Firstly, in 2009 the entire growing season of winter wheat in the two districts of Tongzhou and Shunyi of Beijing was divided into 54 stages at five-day intervals. Net Primary Production (NPP) of winter wheat was estimated by the improved CASA model with HJ-1A/B satellite images from 39 transits. For the 15 stages without HJ-1A/B transit, MOD17A2H data products were interpolated to obtain the spatial distribution of winter wheat NPP at 5-day intervals over the entire growing season of winter wheat. Then, an NPP-yield conversion model was utilized to estimate winter wheat yield in the study area. Finally, the accuracy of the method to estimate winter wheat yield with remote sensing images was verified by comparing its results to the ground-measured yield. The results showed that the estimated yield of winter wheat based on remote sensing images is consistent with the ground-measured yield, with R2 of 0.56, RMSE of 1.22 t ha−1, and an average relative error of −6.01%. Based on time-series satellite remote sensing images, the improved CASA model can be used to estimate the NPP and thereby the yield of regional winter wheat. This approach satisfies the accuracy requirements for estimating regional winter wheat yield and thus may be used in actual applications. It also provides a technical reference for estimating large-scale crop yield.

1. Introduction

Food is the basis of human survival and development. However, in recent years, with the increase in population and human activity, the global ecosystem has been threatened. Air pollution and the loss of cultivated land have had significant impacts on grain crops, creating a major worldwide challenge in the form of food shortages and food security [1,2,3]. Accurate and timely monitoring and evaluation of regional grain yield is vital to the scientific formulation of agricultural import and export plans, the regulation of the grain market, and the guidance and adjustment of planting structures [4,5,6,7]. At present, with the rapid development of remote sensing technology, the traditional time-consuming methods for estimating grain yield, such as statistical surveying, agronomy forecasting, and agrometeorology, have gradually integrated the use of remote sensing images to more rapidly estimate grain yield and other factors over large areas [8,9,10,11].
In recent years, research has progressed in the field of estimating crop yield based on remote sensing images [12,13,14,15,16]. In particular, research has focused on methods to calculate the net primary production (NPP) of crops based on remote sensing data to estimate the mass of crop dry matter quality and the associated crop yield [17]. The conventional models used to estimate models crop yield from remote sensing data include statistical, crop-growth, and parameter models [18]. Examples of the statistical model include the work of Xu et al. [19], who used Landsat5 TM remote sensing images as input into an iteration algorithm of weighted optimal combination (WOC) to assign the optimal weight to barley yield over multiple time periods. Lewis et al. [20] then established a more accurate statistical model to estimate yield by using the maximum Normalized Difference Vegetation Index (NDVI) during the crop-growing season as the yield-evaluation index. Kalubarme et al. [21] used a time-series NDVI curve to extract numerous characteristic parameters and constructed a model based on multiple linear regression analysis of yield to estimate the winter wheat yield. However, these statistical models for estimating yield are not based on crop physiology, which causes problems such as low stability and poor spatial and temporal expansion. In contrast, the high-precision crop-growth models: Crop Environment Resource Synthesis (CERES) [22], Soil Water Atmosphere Plant (SWAP) [23], and World Food Study (WOFOST) [24,25] take full account of the physiological process of crop formation. Unfortunately, these models are too complex and involve a large number of input parameters, so they are only suitable for small areas [8,26]. Parameter model combines the characteristics of empirical and mechanistic models and therefore offers the advantages of both types of models, such as simple structure, less complicated parameters, easy access, and strong applicability. It not only meets the requirements of national and even global large-scale NPP estimates, but can also be used on regional scales. It is a good model for estimating yield by simulating crop NPP [27,28,29].
A typical parameter model is the Carnegie–Ames–Stanford approach (CASA) model, which was originally developed from the model proposed by Monteith in 1972. The CASA model estimates the NPP by focusing on the driving role of the absorbed photosynthetically active radiation (APAR) and the light use efficiency ( ε ) absorbed by vegetation [30]. Since the development of the CASA model, it has been widely used to estimate the NPP [31,32,33,34]. However, the original CASA model is established based on estimates of the NPP of various types of vegetation in North America, and numerous deficiencies remain in the calculation and processing of some of its parameters. For example, the original CASA model used the maximum light use efficiency ( ε * ) = 0.389 g C MJ−1 as the maximum light use efficiency and did not distinguish between vegetation types. For the fraction of absorbed photosynthetically active radiation (fPAR), the original CASA model only established a linear relationship with the Ratio Vegetation Index (SR), which did not properly represent the internal relationship between fPAR and vegetation. In addition, the original CASA model used the soil-moisture model to calculate the water stress factor which was based on the light use efficiency, but the parameters of the soil-moisture model were difficult to obtain, and complex to calculate—making this approach unsuitable for application over a small region. Thus, with the change of environment and climatic conditions and plant characteristics, the original CASA model is not suitable to estimate the NPP of vegetation. It is vital to improve the original CASA model based on actual production and applications, which may be done by incorporating high-quality remote sensing data.
Given this motivation, the major objectives of this study are (1) to optimize the important parameters of the original CASA model in order to fit the growth of winter wheat in the study area; (2) to obtain the high temporal and spatial distribution of winter wheat NPP at five-day cycle intervals in the growing season based on the improved CASA model; (3) to analyze the distribution of winter wheat yield in 2009 combined with the NPP-yield conversion model. Therefore, two multispectral data: HJ-1A/B remote sensing images and the MOD17A2H data products, as data sources in this study. In addition, the accuracy of this estimation was verified by the measured yield data.

2. Study Area and Data

2.1. Study Area

The winter wheat experiment was carried out in the Tongzhou and Shunyi districts of Beijing from 2008 to 2009, as indicated in Figure 1.
Tongzhou District is located southeast of Beijing. It is 36.5 km wide from east to west and 48 km long from north to south, for an area of 906 km2. The highest altitude is 27.6 m and the lowest altitude is 8.2 m. It has a typical continental monsoon climate with distinct seasons; with rain and heat in the same season, a large temperature difference between day and night, an annual average temperature of 10–12 °C, an annual average precipitation of 620 mm, and an average frost-free period of about 180 days. Winter wheat was sown over an area of 156 km2.
Shunyi District is located northeast of Beijing. It is 45 km from east to west and 30 km from north to south, for a total area of 1021 km2. The plain area accounts for 95.7%. The terrain is high in the north and low in the south. The highest point in the northern mountains is 637 m above sea level. It has a temperate continental semi-humid monsoon climate with distinct seasons; with cold and dry in winter, hot and rainy in summer, an annual average temperature of 11.5 °C, an annual average precipitation of 610 mm, and an average frost-free of about 195 days. The average elevation is 35 m and the planting area covers 139 km2.
The main grain crops in the study area are winter wheat and summer maize. Winter wheat is the main crop of summer grain in both Tongzhou District and Shunyi District, and its sown area and harvested yield are in the forefront in Beijing.

2.2. Data and Processing

2.2.1. Remote Sensing Data

The remote sensing data used in this study came from two sources: HJ-1A/B satellite images with a spatial resolution of 30 m and MOD17A2H data with a spatial resolution of 500 m from 2008 to 2009.
(1) HJ-1A/B Data Product
China’s independently developed environment and disaster monitoring and forecasting small satellite constellations A and B, called HJ-1A/B, are equipped with CCD sensors with a spatial resolution of 30 m, radiometric resolution of 8 bit, spectral resolution (visible and near-infrared bands: 430–520 nm, 520–600 nm, 630–690 nm, and 760–900 nm). When satellites A and B are used in parallel, the return visit period can attain two days [35,36]. The level-2 products of the HJ-1A/B satellite used in this work are taken from the China Centre For Resources Satellite Data and Application website (http://www.cresda.com) [37]. From the beginning of October 2008 to the end of June 2009, high-quality images were downloaded following the principle of no cloud or less cloud in the study area. High-quality remote sensing image data for 39 scenes were screened and reserved for atmospheric radiation correction, geometric correction, image splicing, clipping, and other image preprocessing operations.
According to the growth characteristics of winter wheat, the winter wheat growing season was divided into 54 stages at 5-day cycle intervals, 39 of which could be applied to the 39 scenes of the preprocessed HJ-1A/B images and used as input data for the CASA model to extract the NDVI of internal stages. In addition, Savitzky–Golay filter was applied to eliminate image noise for the preprocessed images [38,39,40].
Then, a supervised classification method of maximum likelihood combined with decision tree classification was used to complete the extraction of the winter wheat planting area based on processed HJ-1A/B images. The wheat extraction results are shown in Figure 2.
(2) MOD17A2H Data Product
MOD17A2 is the primary data product of the Moderate Resolution Imaging Spectrometer (MODIS) with a spectrum range of 400–14400 nm, 36 spectral bands, and a radiometric resolution of 12-bit. It offers the first periodic, near-real-time, repetitive monitoring of global vegetation, and provides Gross Primary Production (GPP) synthesized data products with a resolution of 1 km and 8-day [41,42]. MOD17A2H is an upgrade of MOD17A2, which offers level-4 products with a spatial resolution of 500 m and temporal resolution of 8-day.
For the 15 stages of the study not covered by the HJ-1A/B data with medium spatial resolution, MODIS land standard MOD17A2H data products were used for interpolation in 2008 and 2009, which were downloaded from the NASA website (https://ladsweb.modaps.eosdis.nasa.gov/) [43]. Its data preprocessing mainly included projection conversion, data mosaic, clipping, and so on. The MOD17A2H data products of spatial resolution were changed from 500 m to 30 m by resampling in order to extract the winter wheat images of single vegetation. In addition, the MOD17A2H value represented the daily cumulative value of GPP in 8-day, which was multiplied by 5 8 can be considered as the cumulative value of GPP for five days.
Gross primary production (GPP) is the total amount of organic matter produced by photosynthesis of green plants, which is the sum of NPP and autotrophic respiration [44]. Since the NPP model was used to estimate the winter wheat yield in this work, the GPP generated by winter wheat photosynthesis was transformed into the NPP. According to numerous studies, the proportion in the GPP of organic mass consumed by autotrophic respiration during crop growth can be regarded as a fixed [45], with a constant of proportionality of 0.42 [46,47]. Therefore, the value 0.58 can be used as the conversion coefficient between GPP and NPP. Processing and calculating the selected MOD17A2H data products gave the NPP of winter wheat with a time resolution of 5-day, which was used in our model to estimate the yield of winter wheat.

2.2.2. Meteorological Data

The average daily temperature, hours of sunshine, and daily precipitation used herein were provided by the National Meteorological Information Center (http://data.cma.cn/) [48]. The processing consisted of the following sections:
(1)
The meteorological data from five meteorological stations around the study area was selected, and then it was arranged with 5-day intervals to obtain the 5-day-average temperature, total hours of sunshine, and total precipitation.
(2)
The spatial coordinate information based on the longitude and latitude of each station was assigned to meteorological data.
(3)
The meteorological data with a spatial resolution of 30 m was produced by spatial interpolation of the date from the five meteorological stations.

2.2.3. Measured Yield Data

A total of 29 samples were selected from the study area, 15 from the Tongzhou District and 14 from the Shunyi District. The location of the sample area was determined by considering the representative growth status and regional distribution of winter wheat, and the size of each sample area was not less than 100 m × 100 m. The area of each sampling point was 1 m × 1 m. Five-point sampling was used for sampling. After drying and weighing, the average weight of winter wheat from the five sampling points was used as the yield data for the given sample area.

3. Study Methods

3.1. Construction of Improved CASA Model

Net primary production (NPP) refers to the total amount of organic matter accumulated by plants in the primary production stage per unit time and unit area. It is an essential indicator of plant-growth status [34]. The improved CASA model was used in this study to simulate the winter wheat NPP [49], and then the yield of the study area was estimated. After years of development, the parameters of the CASA model can be obtained by the inversion of remote sensing data and processing of meteorological data. This is done by using
NPP ( x , t ) = APAR ( x , t ) × ε ( x , t ) ,
where t is time, x is the pixel number, NPP(x, t) [g C m−2] is the net primary production of pixel x at time t, APAR(x, t) [MJ m−2] is the absorbed photosynthetically active radiation of pixel x at time t, and ε ( x , t ) [g C MJ−1] is the light use efficiency for pixel x at time t.

3.2. Determination of Absorbed Photosynthetically Active Radiation

The absorbed photosynthetically active radiation by actual vegetation (APAR) is mainly related to two factors: the total solar radiation SOL and the fraction of absorption by vegetation to APAR. This is calculated as follows:
APAR ( x , t ) = 1 2 × S O L ( x , t ) × fPAR ( x , t ) ,
where SOL(x, t) [MJ m−2] is the total solar radiation of pixel x at time t, fPAR(x, t) is the fraction of absorbed photosynthetically active radiation, and the factor 1 2 is the ratio of effective solar radiation that can be used by vegetation to the total solar radiation.

3.2.1. Determination of Total Solar Radiation

The direct calculation of the total solar radiation was not included in this study because the small number of meteorological stations that provide total solar radiation. Therefore, the daily solar radiation was calculated by using the hours of sunshine, and the total solar radiation at 5-day intervals is incorporated separately by using
S O L ( x , t ) = ( a + b n N ) R a ,
where SOL(x, t) is the same as for Equation (2), and a and b are empirical coefficients that are determined from the measured solar radiation. In this study, the method of Zuo et al. [50] was used to calculate a and b. For the Beijing area, a = 0.22 and b = 0.72. The ratio n N , where n is the hours of sunshine and N = 24 h, gives the fraction of sunshine time for one day. Finally, R a [MJ m−2] is the solar intensity at ground level.
The quantity R a can be calculated from technical documents published by Food and Agriculture Organization of the United Nations (FAO) [51], which is related to the distance d between the Sun and Earth, the latitude φ , the solar sunset angle ω 0 , the solar constant S 0 = 0.082 MJ m−2 min, the solar declination delta δ , and the number J of days. R a is calculated by using
R a = 24 60 π × S 0 d ( ω 0 s i n φ s i n δ + c o s φ c o s δ s i n ω 0 ) ,
In Equation (4), the parameters are given by
d = 1 + 0.033 cos ( 2 π 365 J ) ,
ω 0 = a r c c o s ( t a n φ t a n δ ) ,
δ = 0.409 sin ( 2 π 365 J 1.39 ) ,

3.2.2. Improved Calculation of fPAR

For winter wheat, fPAR depends on the growing season and the growing period. References [18,52,53,54] showed the NDVI correctly represents the vegetation type and growth situation, and it is linear in fPAR. In the present study, the maximum and minimum NDVI and fPAR can be determined by HJ-1A/B images to calculate fPAR(x, t) as follows:
fPAR ( x , t ) = [ NDVI ( x , t ) NDVI min ] ( fPAR max fPAR min ) ( NDVI max NDVI min ) + fPAR min ,
In addition, Field et al. [55] found that fPAR is linear in SR:
fPAR ( x , t ) = [ SR ( x , t ) SR min ] ( fPAR max fPAR min ) ( SR max SR min ) + fPAR min ,
SR = 1 + NDVI ( x , t ) 1 NDVI ( x , t ) ,
where fPAR max = 0.950 and fPAR min = 0.001 are the constant maximum and minimum values of fPAR.
The constants NDVI max and NDVI min are the maximum and minimum NDVI for winter wheat. The determination of NDVI maximum and minimum has always been controversial. After making a comparative analysis, the 95% lower-side value and 5% lower-side value of the winter wheat NDVI maximum and minimum probability distribution for NDVI max and NDVI min was used in this study, respectively [56]. However, the present study covers the entire growing season of winter wheat, which translates into 39 scenes of HJ-1A/B data and 15 scenes of MOD17A2H data, with more processed images than usual for such studies. Thus, the available HJ-1A/B images for each month with the NDVIs for winter wheat were applied, which may differ greatly between months during the growing season of winter wheat, to calculate the monthly maximum and minimum NDVI by the calculation method mentioned above. SR max and SR min are determined from NDVI max and NDVI min , respectively.
With the improved CASA model proposed herein, both vegetation indices NDVI and SR can be used to estimate fPAR. According to previous research, the calculation of fPAR based on NDVI is often greater than the NVDI based on measurements, and the calculation of fPAR based on SR is often lower than the NVDI based on measurements. Therefore, the same weight α = 0.5 was given to the two calculation methods:
fPAR ( x , t ) = α fPAR NDVI + ( 1 α ) fPAR SR ,

3.3. Improved Estimation of Light Use Efficiency

Ideally, plant vegetation has maximum light use efficiency but, in reality, the light use efficiency is affected by many external environmental factors: mainly temperature and precipitation. The light use efficiency thus was calculated by using the model established by Potter and Field [49,55] combined with the regional actual evapotranspiration model proposed by Zhou and Zhang in 1995 [57]. Specifically, the ε is calculated as follows:
ε ( x , t ) = T ε 1 ( x , t ) × T ε 2 ( x , t ) × W ε ( x , t ) × ε * ,
T ε 1 ( x , t ) = 0.8 + 0.002 T o p t ( x , t ) 0.0005 T 2 o p t ( x , t ) ,
T ε 2 ( x , t ) = 1.184 { 1 + e x p [ 0.2 ( T o p t ( x , t ) 10 T ( x , t ) ) ] } × { 1 + e x p [ 0.3 ( T o p t ( x , t ) 10 + T ( x , t ) ) ] } ,
W ε ( x , t ) = ( 1 2 ) E E T ( x , t ) P E T ( x , t ) + 1 2 ,
where T ε 1 ( x , t ) and T ε 2 ( x , t ) are the temperature stress factors of the light use efficiency for pixel x at time t, W ε ( x , t ) is the water stress factor of the light use efficiency for pixel x at time t, ε * [g C MJ−1] is the maximum light use efficiency under ideal conditions, T ( x , t ) is the average temperature for pixel x at time t, T opt ( x , t ) is the appropriate temperature for plant growth for pixel x at time t; E E T ( x , t ) [mm] is the actual evapotranspiration and based on using Equation (16), which in turn depends on the Penman Equation (17) [44]. P E T ( x , t ) [mm] is the potential evapotranspiration, which is calculated as follows:
E E T ( x , t ) = { P ( x , t ) × R n ( x , t ) × [ P 2 ( x , t ) + R 2 n ( x , t ) + P ( x , t ) × R n ( x , t ) ] } { [ P ( x , t ) + R n ( x , t ) ] × [ P 2 ( x , t ) + R 2 n ( x , t ) ] }
R n ( x , t ) = [ E p 0 ( x , t ) × P ( x , t ) ] 1 2 × { 0.369 + 0.598 × [ E p 0 ( x , t ) P ( x , t ) ] 1 2 }
P E T ( x , t ) = E p o ( x , t ) + E E T ( x , t ) 2
where, P ( x , t ) [mm] is the precipitation for pixel x at time t, R n ( x , t ) [MJ m−2] is the net radiation for pixel x at time t, E p 0 ( x , t ) is the local potential evapotranspiration, which is related to the temperature for pixel x at time t, and Reference [58] is referenced to calculate E p 0 ( x , t ) based on vegetation-climate method established by Thornthwaite.
An important parameter of the CASA model is ε * , the ideal-condition maximum light use efficiency, which can be assigned a fixed value when processing data corresponding to a small area [55,59]. ε * directly affects the value calculated for the light use efficiency. Without distinguishing vegetation types and geographical types, Field et al. [55] used 0.389 g C MJ−1 as the global ε * of vegetation. Hunt [60] assigned an upper limit of 3.5 g C MJ−1 to the ε , and Zhu et al. [56] estimated the ε * for various typical types of vegetation in China based on NPP estimated by remote sensing models. Goetz et al. [61,62] claimed that the ε * ranges from 0.42 to 3.8 g C MJ−1, whereas, other results suggested that the ε * of wheat ranges from 1.46 to 2.93 g C MJ−1 [63,64]. Different types of vegetation have different ecological and physiological structures, so their ε * should differ. Since the present study focuses on winter wheat, the environmental conditions of the study area along with the results of previous research are considered to fix the ε * of winter wheat at ranges from 0.42 to 2.93 g C MJ−1. Then, the average of the range was calculated and rounded to get ε * of 1.7 g C MJ−1, which had a good result for the improved CASA model in this study.

3.4. NPP-Yield Conversion Model

Winter wheat yield is the grain harvest of winter wheat, which often refers to the economic yield of winter wheat. Along with Ren et al. [65], the method was used to estimate winter wheat yield, which is
Yield = α NPP × p × HI 1 ω × 10 2
where Yield is the winter wheat yield [t ha−1], α is the carbon-conversion coefficient (the carbon content of winter wheat is about 45%, so the α = 2.22 [66]), NPP is the cumulative net primary production of organic matter over the entire growing season of winter wheat, p = 0.9 is the distribution coefficient of aboveground parts, ω = 12.5% is the water-content coefficient of wheat grain during the storage period after harvest [67], and HI = 0.45 is the harvest index (this value is based on the growth characteristics of winter wheat in the Tongzhou and Shunyi districts of Beijing, as determined by previous research), 10 2 is the conversion factor between the unit g m−2 and t ha−1.

4. Results and Analysis

With the gradual increase of the hours of sunshine, temperature and precipitation; the winter wheat enters a stage of rapid growth, especially from March to May, which is a critical period for accumulating dry matter quality. Although this study analyzes the distribution of winter wheat yield over the entire growing season, the parameters fPAR, the light use efficiency, and the NPP from March to May are focused on to parametrize the improved CASA model.

4.1. Fraction of Absorbed Photosynthetically Active Radiation

The monthly NDVI and SR maximum and minimum values need to be determined before calculating fPAR based on Equation (11), which are listed in Table 1.
Figure 3 shows the results for fPAR from March to May 2009. The fPAR results for March, which average 0.47, were less than those for April and May, which fell mostly between 0.32 and 0.57. The fPAR for winter wheat in April continued to increase due to the work of photosynthesis and generally exceeds 0.40. The process of photosynthesis of winter wheat gradually weakened as May approaches, so the fPAR values concentrated in the range 0.40–0.60. The greatest average fPAR of 0.59 occurred in April 2009, whereas the average fPAR was 0.53 in May of the same year. These results reflected the mild weather during April, in which hydrothermal and sunlight conditions were suitable for winter wheat growth, resulting in increased photosynthesis of winter wheat and the subsequent increase in growth. In addition, the April and May fPAR results for Tongzhou and Shunyi districts did not differ significantly, which meant that the geographical environment of these two districts was suitable for winter wheat growth.

4.2. Light Use Efficiency

In 2009, the light use efficiency for winter wheat in the Tongzhou and Shunyi districts of Beijing remained relatively constant over each given month, although slight differences appeared due to differences in climate and environment. As shown in Figure 4, the ε in March ranged from 0.67 to 0.71 g C MJ−1, and the average monthly light use efficiency was 0.69 g C MJ−1. In April, the ε increased to reach 0.80–1.0 g C MJ−1, and the average monthly ε was 0.92 g C MJ−1. Environmental conditions such as sunshine, temperature, and precipitation in May provided improved growing conditions for winter wheat, so the ε in May ranged from 1.06 to 1.15 g C MJ−1, with an average of 1.11 g C MJ−1. Figure 4 also shows that the ε in March and April was approximately the same in Tongzhou District as in Shunyi District. In May, however, the ε in Shunyi District was significantly greater than in Tongzhou District. This result reflected the better environmental conditions for winter wheat growth, such as temperature and precipitation, in Shunyi District from March to May, which allowed superior growth and development of winter wheat.

4.3. Net Primary Production

Figure 5 shows the spatial distribution of NPP in the study area from March to May, 2009. The NPP generally increased from March to May. In March, winter wheat gradually enters the greening stage and the erecting stage. The results for NPP in March were all less than 100 g C m−2, and the average NPP was 25 g C m−2. In April, most of the winter wheat was in the jointing stage, booting stage, and heading stage and the growth rate continued to accelerate. The NPP was mainly concentrated in the range 50–150 g C m−2, with an average of 88 g C m−2. In May, the average NPP reached a maximum of 128 g C m−2 and mostly remained above 100 g C m−2. The winter wheat growth increased continuously, and the dry matter quality accumulated rapidly.
The Tongzhou and Shunyi districts in Beijing have typical continental monsoon climates. The seasonal climate directly affects the formation of winter wheat NPP. Figure 6 shows the mean monthly NPP of winter wheat in the study area as a function of time over the entire growing season. From October to December 2008, the winter wheat NPP in the seedling stage decreased monotonically, primarily because, as the study area gradually entered the winter season, climatic factors such as solar radiation, temperature, and precipitation reached their low point. Thus, winter wheat growth was limited; the growth rate was slow, and the rate of NPP formation was relatively low. From January to March 2009, with the improvement of climatic conditions, the winter wheat gradually passed through the over-wintering stage and entered the stage of rapid growth and development, so the NPP began to increase. From April to June 2009, the winter wheat leaf area continued to increase, photosynthesis, respiration, and other growth activities continued to increase, and dry matter quality accumulated rapidly. Finally, in May, the winter wheat growth and NPP peaked.
Figure 7 shows the spatial distribution of the total winter wheat NPP for the entire growing season, which was obtained by summing the NPP of each month. The overall distribution of winter wheat NPP in Tongzhou District was similar to that in Shunyi District and was concentrated in the range 300–700 g C m−2, with an average NPP of 464 g C m−2.

4.4. Estimation of Winter Wheat Yield

In this study, the winter wheat NPP was used to estimate the winter wheat economic yield. First, the improved CASA model was applied to estimate the winter wheat NPP in the Tongzhou and Shunyi districts of Beijing for the 2008–2009 growing season. Next, the NPP-yield conversion model [Equation (19)] was used to convert the winter wheat NPP into the winter wheat economic yield.
Figure 8 shows the spatial distribution of the estimated winter wheat yield in the study area, which ranged from 3 to 8 t ha−1 over more than 90% of the total study area. The output of the entire study area was about 50% for the range 3–5 t ha−1 and 35% for the range 5–7 t ha−1. The results show that, during the 2008–2009 growing season, the Tongzhou and Shunyi districts of Beijing had similar winter wheat growth trends, however, the winter wheat yield in the central part of the Tongzhou District and in the western part of the Shunyi District was higher than that in the northern part of the Tongzhou District and in the western part of the Shunyi District. Overall, the growth of winter wheat in the study area was ideal.

4.5. Verification of Estimated Yield

A total of 29 sampling points were selected to verify that the improved CASA model and the NPP-yield conversion model provided accurate estimates of winter wheat yield. Each pixel for the sampling points was averaged with the eight adjacent pixels to obtain the estimated yield for the given sampling point. Figure 9 shows the estimated yield as a function of the measured yield.
The determination coefficient R2 = 0.56 and RMSE was 1.22 t ha−1. This analysis gave an average absolute error of −0.57 t ha−1 and an average relative error of −6.01%, which satisfied the accuracy requirements for estimating regional winter wheat yield. Thus, the proposed method may have good application potential.

5. Discussion

Accurate and timely monitoring and evaluation of regional grain crop yield is a research priority of numerous countries, and remote sensing technology is widely used in a wide range of techniques to estimate crop yield [68,69,70]. MODIS remote sensing images are often used as a data source in many studies. These images may be downloaded for free and offer multiple data bands and a fast update frequency [71,72]. However, the choice of spatial resolution for MODIS images is only 250, 500, and 1000 m. Thus, single MODIS remote sensing images with high temporal resolution and low spatial resolution are often applicable for large areas. The Chinese domestic satellite HJ-1A/B, however, offers a relatively high spatial resolution of 30 m, but the return-time period is inferior to that of MODIS. Therefore, to study the Tongzhou and Shunyi districts of Beijing, the improved CASA model was developed to estimate winter wheat NPP based on HJ-1A/B images with high temporal and medium spatial resolution. At the same time, MOD17A2H images were incorporated to fully exploit the advantages of each remote sensing data source and obtain the spatial distribution of winter wheat NPP with a high spatial and temporal resolution (5-day intervals over the entire winter wheat growing season). The results provide support for monitoring and evaluation of winter wheat yield.
Thus, by using the improved CASA model and the NPP-yield conversion model, the winter wheat yield in the study area was estimated based on remote sensing images. However, the parameters used in the two models were often affected by factors such as vegetation type, geographical conditions, climate, and environment. Therefore, based on previous studies, some of the model parameters were optimized based on the actual environmental conditions in the study area during the winter wheat growing season. The most influential parameters for the improved CASA model were fPAR and the light use efficiency ( ε ).
The original CASA model only establishes a linear relationship between fPAR and SR. In this study, the two vegetation indices NDVI and SR are used to calculate fPAR NDVI and fPAR SR , respectively, with each weighted the same, to estimate fPAR. In addition, although most studies calculate the maximum and minimum NDVI by using the 95% lower side value ( NDVI max ) and the 5% lower side value ( NDVI min ) of the NDVI maximum probability distribution extracted from remote sensing images [56], respectively a single constant was assigned to the NDVI maximum and minimum monthly because of the long growing season of winter wheat and the large number of remote sensing images involved in the calculations. In view of the actual situation of winter wheat over the entire growing season, this approach was justified because it reduced the parameter error and improves the estimation accuracy.
The original CASA model used a soil-moisture model to calculate the water stress factor that appeared in the light use efficiency ( ε ) . The maximum light use efficiency ( ε * ) was assigned a value of 0.389 g C MJ−1. In addition, many studies have been done to optimize the parameters of the model Reference [73], which showed that remote sensing and optimization techniques can be used to improve the precision of crop yield estimation, where the water stress factor depends on energy balance equation and   ε * is equal to 0.29 g C/MJ for C3 crops. However, both the structure of the soil-moisture model and energy balance equation are complex and the parameters are not easy to determine, so ε in this study was calculated by using the actual regional evapotranspiration model proposed by Zhou et al. [57], where temperature and precipitation can be used as input data to calculate parameters: the net radiation ( R n ) and the local potential evapotranspiration ( E p 0 ).
The value of ε * has a great influence on the simulation of NPP, and its size is controversial [55]. Therefore, designing experiment is the good way to determine the value of ε * . Li et al. [74] applied Particle Swarm Optimizer (PSO) to search for the optimization of LUE retrievals through the CASA model combined with time-series NDVI and NPP ground measurements, which not only proved the significant difference of among various vegetation types, but also determined the specific value of ε * in grassland. However, the objective of this study was to estimate NPP of winter wheat based on the improved CASA model, so no specific experiments were designed to determine the value of ε * . Upon combining this with the growth characteristics of winter wheat in the study area, the range of maximum light use efficiency ( ε * ) was determined. Then, the average of the range was calculated and rounded to get ε * of 1.7 g C MJ−1.
This study develops the improved CASA model and combines it with the NPP to estimate winter wheat yield from time-series satellite remote sensing images. The comparison of the estimated results with the ground-measured-yield data shows that the estimation accuracy satisfies the requirements for estimating regional winter wheat. Compared with the study by Ren et al. [65], both methods have good results, but this model has obvious advantages in particular, widely used remote sensing images, NPP distribution with high temporal and spatial resolution, and simpler input parameters. It is still an important objective to improve the accuracy of the estimated yield of crop, which can be considered from two aspects. One is the optimization of data sources. The estimation of crop yield depends partly on the quality and quantity of the input data. In this study, the mature integrated GPP product (MOD17A2H) was applied, which existed objectively and was widely used. Reference [73] creates an optimization model to compensate for missing data, which can improve the precision of crop yield estimation combined with Monteith model (CASA). Although the study object is not winter wheat, it is still worth paying attention to this optimization method. The other is the essential parameters in the model for estimating winter wheat yield based on remote sensing images: the ε * and the HI, which, based on previous research, were assigned constant values in the present study. However, these two parameters are affected by other factors, so variable processing must be done in addition to relevant experimentation, which will be the subject of future research.
A total of 29 ground-measured winter wheat yield data were used in this study and winter wheat yield were concentrated in the relatively small range of 3–8 t ha−1. An ideal result was obtained by using 29 sampling points in this study area, but if the study area scale is continuously expanded, further research is required to verify the accuracy and applicability of the proposed model.

6. Conclusions

This study combined two kinds of satellite data: HJ-1A/B remote sensing images and MOD17A2H data products as input for an improved CASA model to estimate the net primary production (NPP) during the entire 2008–2009 winter wheat growing season. The NPP-yield conversion model was then used to estimate winter wheat yield on a regional scale. The results give a determination coefficient R2 = 0.56 between the measured and estimated winter wheat yield, RMSE of 1.22 t ha−1, and an average relative error of −6.01%. Thus, the improved CASA model combined with the NPP-yield conversion model satisfied the requirements for regional-scale estimates of winter wheat yield from remote sensing images. This approach may thus be applied in the field and can also provide a benchmark for a wide range of models for estimating crop yield.

Author Contributions

Y.W.; processed and analyzed the data and drafted the manuscript. X.X.; guided the ideal, designed the experiment, advised on data analysis and revised the manuscript. G.Y.; L.H.; L.F.; P.W. and G.C. were involved in the experiments, ground data collection, and manuscript revision. All authors read and approved the final version.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 41571416) and the National Key Research and Development Program (Grant No. 2017YFD0201501).

Acknowledgments

We appreciate Hong Chang and Weiguo Li for their help during field data collection. The funders had no role in choosing the study design or in the collection, analysis and interpretation of the data, in the writing of the report, or in the decision to submit the article for publication.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Godfray, H.C.J.; Beddington, J.R.; Crute, I.R.; Haddad, L.; Lawrence, D.; Muir, J.F.; Pretty, J.; Robinson, S.; Thomas, S.M.; Toulmin, C. Food Security: The Challenge of Feeding 9 Billion People. Science 2010, 327, 812–818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Atzberger, C. Advances in Remote Sensing of Agriculture: Context Description, Existing Operational Monitoring Systems and Major Information Needs. Remote Sens. 2013, 5, 949–981. [Google Scholar] [CrossRef] [Green Version]
  3. Tao, F.; Yokozawa, M.; Zhang, Z.; Xu, Y.; Hayashi, Y. Remote sensing of crop production in China by production efficiency models: Models comparisons, estimates and uncertainties. Ecol. Model. 2005, 183, 385–396. [Google Scholar] [CrossRef]
  4. Rae, A.; Pardey, P. Global food security-introduction. Aust. J. Agric. Resour. Econ. 2014, 58, 499–503. [Google Scholar] [CrossRef]
  5. Lipper, L.; Thornton, P.; Campbell, B.M.; Baedeker, T.; Braimoh, A.; Bwalya, M.; Caron, P.; Cattaneo, A.; Garrity, D.; Henry, K.; et al. Climate-smart agriculture for food security. Nat. Clim. Chang. 2014, 4, 1068–1072. [Google Scholar] [CrossRef]
  6. Zhao, Y.; Chen, S.; Shen, S. Assimilating remote sensing information with crop model using Ensemble Kalman Filter for improving LAI monitoring and yield estimation. Ecol. Model. 2013, 270, 30–42. [Google Scholar] [CrossRef]
  7. Li, Y.; Zhou, Q.; Zhou, J.; Zhang, G.; Chen, C.; Wang, J. Assimilating remote sensing information into a coupled hydrology-crop growth model to estimate regional maize yield in arid regions. Ecol. Model. 2014, 291, 15–27. [Google Scholar]
  8. Huang, J.; Tian, L.; Liang, S.; Ma, H.; Becker-Reshef, I.; Huang, Y.; Su, W.; Zhang, X.; Zhu, D.; Wu, W. Improving winter wheat yield estimation by assimilation of the leaf area index from Landsat TM and MODIS data into the WOFOST model. Agric. For. Meteorol. 2015, 204, 106–121. [Google Scholar] [CrossRef]
  9. Singh, R.; Semwal, D.P.; Rai, A.; Chhikara, R.S. Small area estimation of crop yield using remote sensing satellite data. Int. J. Remote Sens. 2002, 23, 49–56. [Google Scholar] [CrossRef]
  10. Fang, H.; Liang, S.; Hoogenboom, G. Integration of MODIS LAI and vegetation index products with the CSM–CERES–Maize model for corn yield estimation. Int. J. Remote Sens. 2011, 32, 1039–1065. [Google Scholar] [CrossRef] [Green Version]
  11. Doraiswamy, P.C.; Moulin, S.; Cook, P.W.; Stern, A. Crop Yield Assessment from Remote Sensing. Photogramm. Eng. Remote Sens. 2003, 69, 665–674. [Google Scholar] [CrossRef] [Green Version]
  12. Lobell, D.B. The use of satellite data for crop yield gap analysis. Field Crop. Res. 2013, 143, 56–64. [Google Scholar] [CrossRef] [Green Version]
  13. Prasad, A.K.; Chai, L.; Singh, R.P.; Kafatos, M. Crop yield estimation model for Iowa using remote sensing and surface parameters. Int. J. Appl. Earth Obs. Geoinf. 2006, 8, 26–33. [Google Scholar] [CrossRef]
  14. Nemani, R.R.; Keeling, C.D.; Hashimoto, H.; Jolly, W.M.; Piper, S.C.; Tucker, C.J.; Myneni, R.B.; Running, S.W. Climate-Driven Increases in Global Terrestrial Net Primary Production from 1982 to 1999. Science 2003, 300, 1560–1563. [Google Scholar] [CrossRef] [Green Version]
  15. Thorp, K.R.; Wang, G.; West, A.L.; Moran, M.S.; Bronson, K.F.; White, J.W.; Mon, J. Estimating crop biophysical properties from remote sensing data by inverting linked radiative transfer and ecophysiological models. Remote. Sens. Environ. 2012, 124, 224–233. [Google Scholar] [CrossRef] [Green Version]
  16. Mirschel, W.; Wieland, R.; Wenkel, K.-O.; Nendel, C.; Guddat, C. YIELDSTAT—A spatial yield model for agricultural crops. Eur. J. Agron. 2014, 52, 33–46. [Google Scholar] [CrossRef]
  17. Xu, X.G.; Wu, B.F.; Meng, J.H.; Li, Q.Z.; Huang, W.J.; Liu, L.Y.; Wang, J.H. Research advances in crop yield estimation models based on remote sensing. Trans. CSAE 2008, 24, 290–298, (In Chinese with English abstract). [Google Scholar]
  18. Ruimy, A.; Saugier, B. Methodology for the estimation of terrestrial net primary production from remotely sensed data. J. Geophys. Res. 1994, 97, 18515–18521. [Google Scholar] [CrossRef]
  19. Xu, X.G.; Wang, J.H.; Huang, W.J.; Li, C.J.; Yang, X.D.; Gu, X.H. Estimation of crop yield based on weight optimization combination and multi-temporal remote sensing data. Trans. CSAE 2009, 25, 137–142, (In Chinese with English abstract). [Google Scholar]
  20. Lewis, J.E.; Rowland, J.; Nadeau, A. Estimating maize production in Kenya using NDVI: Some statistical considerations. Int. J. Remote Sens. 1998, 19, 2609–2617. [Google Scholar] [CrossRef]
  21. Kalubarme, M.H.; Potdar, M.B.; Manjunath, K.R.; Mahey, R.K.; Siddhu, S.S. Growth profile based crop yield models: A case study of large area wheat yield modelling and its extendibility using atmospheric corrected NOAA AVHRR data. Int. J. Remote Sens. 2003, 24, 2037–2054. [Google Scholar] [CrossRef]
  22. Jones, J.W.; Hoogenboom, G.; Porter, C.H.; Boote, K.J.; Batchelor, W.D.; Hunt, L.A.; Wilkens, P.W.; Singh, U.; Gijsman, A.J.; Ritchie, J.T. The DSSAT cropping system model. Eur. J. Agron. 2003, 18, 235–265. [Google Scholar] [CrossRef]
  23. Huang, J.; Ma, H.; Su, W.; Zhang, X.; Huang, Y.; Fan, J.; Wu, W. Jointly assimilating MODIS LAI and ET products into the SWAP model to estimate winter wheat yield. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 4060–4071. [Google Scholar] [CrossRef]
  24. Huang, J.; Fernando, S.; Huang, Y.; Ma, H.; Li, X.; Liang, S.; Zhang, X.; Fan, J.; Wu, W. Assimilating a synthetic Kalman filter leaf area index series into the WOFOST model to estimate regional winter wheat yield. Agric. For. Meteorol. 2016, 216, 188–202. [Google Scholar] [CrossRef]
  25. Huang, J.; Ma, H.; Sedano, F.; Lewis, P.; Liang, S.; Wu, Q.; Su, W.; Zhang, X.; Zhu, D. Evaluation of regional estimates of winter wheat yield by assimilating three remotely sensed reflectance datasets into the coupled WOFOST–PROSAIL model. Eur. J. Agron. 2019, 102, 1–13. [Google Scholar] [CrossRef]
  26. Dorigo, W.A.; Zurita-Milla, R.; De Wit, A.J.W.; Brazile, J.; Singh, R.; Schaepman, M.E. A review on reflective remote sensing and data assimilation techniques for enhanced agroecosystem modeling. Int. J. Appl. Earth Obs. Geoinf. 2007, 9, 165–193. [Google Scholar] [CrossRef]
  27. Bastiaanssen, W.G.M.; Ali, S. A new crop yield forecasting model based on satellite measurements applied across the Indus Basin, Pakistan. Agric. Ecosyst. Environ. 2003, 94, 321–340. [Google Scholar] [CrossRef]
  28. Lobell, D.B.; Asner, G.P.; Ortiz-Monasterio, J.I.; Benning, T.L. Remote sensing of regional crop production in the Yaqui Valley, Mexico: Estimates and uecertainties. Agric. Ecosyst. Environ. 2003, 94, 205–220. [Google Scholar] [CrossRef]
  29. Wang, P.; Sun, R.; Zhang, J.; Zhou, Y.; Xie, D.; Zhu, Q. Yield estimation of winter wheat in the North China Plain using the remote-sensing–photosynthesis–yield estimation for crops (RS–P–YEC) model. Int. J. Remote Sens. 2011, 32, 6335–6348. [Google Scholar] [CrossRef]
  30. Monteith, J.L. Climate and the efficiency of crop production in Britain. Philos. Trans. R. Soc. B 1977, 28l, 277–294. [Google Scholar] [CrossRef]
  31. Nayak, R.K.; Patel, N.R.; Dadhwal, V.K. Estimation and analysis of terrestrial net primary productivity over India by remote-sensing-driven terrestrial biosphere model. Environ. Monit. Assess. 2010, 170, 195–213. [Google Scholar] [CrossRef] [PubMed]
  32. Xing, X.; Xu, X.; Zhang, X.; Zhou, C.; Song, M.; Shao, B.; Ouyang, H. Simulating net primary production of grasslands in northeastern Asia using MODIS data from 2000 to 2005. J. Geogr. Sci. 2010, 20, 193–204. [Google Scholar] [CrossRef]
  33. Hicke, J.A.; Asner, G.P.; Randerson, J.T.; Tucker, C.J. Trends in North American net primary productivity derived from satellite obser-vations, 1982–1998. Glob. Biogeochem. Cycles 2002, 16, 1019–1040. [Google Scholar] [CrossRef]
  34. Field, C.B.; Behrenfeld, M.J.; Randerson, J.T.; Falkowski, P. Primary Production of the Biosphere: Integrating Terrestrial and Oceanic Components. Science 1998, 281, 237–240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wang, Q. Technical system design and construction of China’s HJ-1 satellites. Int. J. Digit. Earth 2012, 5, 202–216. [Google Scholar] [CrossRef]
  36. Yu, B.; Shang, S. Multi-Year Mapping of Maize and Sunflower in Hetao Irrigation District of China with High Spatial and Temporal Resolution Vegetation Index Series. Remote. Sens. 2017, 9, 855. [Google Scholar] [Green Version]
  37. China Centre for Resources Satellite Data and Application (CRESDA). Available online: http://www.cresda.com/ (accessed on 15 November 2018).
  38. Savitzky, A.; Golay, M.J.E. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  39. Viovy, N.; Arino, O.; Belward, A.S. The Best Index Slope Extraction (BISE): A method for reducing noise in NDVI time-series. Int. J. Remote Sens. 1992, 13, 1585–1590. [Google Scholar] [CrossRef]
  40. Chen, J.; Jönsson, P.; Tamura, M.; Gu, Z.; Matsushita, B.; Eklundh, L. A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter. Remote Sens. Environ. 2004, 91, 332–344. [Google Scholar] [CrossRef]
  41. Running, S.W.; Nemani, R.R.; Heinsch, F.A.; Zhao, M.; Reeves, M.; Hashimoto, H. A Continuous Satellite-Derived Measure of Global Terrestrial Primary Production. BioScience 2004, 54, 547–560. [Google Scholar] [CrossRef]
  42. Running, S.W.; Zhao, M. Daily GPP and Annual NPP (MOD17A2/A3) Products NASA Earth Observing System MODIS Land Algorithm; version 3.0 for Collection 6; NASA: Greenbelt, MD, USA, 2015.
  43. Level-1 and Atmosphere Archive & Distribution System (LAADS) Distributed Active Archive Center (DAAC). Available online: https://ladsweb.modaps.eosdis.nasa.gov/ (accessed on 20 November 2018).
  44. Lieth, H.; Whittaker, R.H. Primary Productivity of the Biosphere; Springer: Berlin, Germany, 1975; pp. 237–263. [Google Scholar]
  45. Albrizio, R.; Steduto, P. Photosynthesis, respiration and conservative carbon use efficiency of four field grown crops. Agric. For. Meteorol. 2003, 116, 19–36. [Google Scholar] [CrossRef]
  46. Gifford, R.M. Whole plant respiration and photosynthesis of wheat under increased CO2 concentration and temperature: Long-term vs. short-term distinctions for modelling. GCB Biol. 1995, 1, 385–396. [Google Scholar] [CrossRef]
  47. Cheng, W.X.; Sims, D.A.; Luo, Y.Q.; Coleman, J. Photosynthesis, respiration, and net primary production of sunflower stands in ambient and elevated atmospheric CO2 concentrations. GCB Biol. 2000, 6, 931–941. [Google Scholar]
  48. National Meteorological Information Center (NMC). Available online: http://data.cma.cn/ (accessed on 11 December 2018).
  49. Potter, C.S.; Randerson, J.T.; Field, C.B.; Matson, P.A.; Vitousek, P.M.; Mooney, H.A.; Klooster, S.A. Terrestrial ecosystem production: A process model based on global satellite and surface data. Glob. Biogeochem. Cycles 1993, 7, 811–841. [Google Scholar] [CrossRef]
  50. Zuo, D.K.; Wang, Y.X.; Chen, J.S. Characteristics of the distribution of total radiation in China. Acta Meteorol. Sin. 1963, 33, 78–96, (In Chinese with English abstract). [Google Scholar]
  51. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements; Irrigation and Drainage Paper No. 56; FAO: Rome, Italy, 1998; p. 300. [Google Scholar]
  52. Fensholt, R.; Sandholt, I.; Rasmussen, M.S. Evaluation of Modis LAI and the relation between fAPAR and NDVI in a semi-arid environment using situ measurements. Remote Sens. Environ. 2004, 91, 490–507. [Google Scholar] [CrossRef]
  53. Myneni, R.B.; Hoffman, S.; Knyazikhin, Y.; Privette, J.L.; Glassy, J.; Tian, Y.; Wang, Y.; Song, X.; Zhang, Y.; Smith, G.R.; et al. Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data. Remote Sens. Environ. 2002, 83, 214–231. [Google Scholar] [CrossRef] [Green Version]
  54. Myneni, R.B.; Williams, D.L. On the relationship between fAPAR and NDVI. Remote Sens. Environ. 1994, 49, 200–211. [Google Scholar] [CrossRef]
  55. Field, C.B.; Randerson, J.T.; Malmstrom, C.M. Global net primary production: Combining ecology and remote sensing. Remote. Sens. Environ. 1995, 51, 74–88. [Google Scholar] [CrossRef] [Green Version]
  56. Zhu, W.Q.; Pan, Y.Z.; He, H.; Yu, D.Y.; Fu, H.D. China’s largest light utilization simulation of typical vegetation. Chin. Sci. Bull. 2006, 51, 700–706, (In Chinese with English abstract). [Google Scholar] [CrossRef]
  57. Zhou, G.S.; Zhang, X.S. A natural vegetation NPP Model. J. Plant Ecol. 1995, 19, 193–200, (In Chinese with English abstract). [Google Scholar]
  58. Thornthwaite, C.W. An Approach toward a Rational Classification of Climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  59. Russell, G.; Marshall, B.; Jarvis, P.G. Absorption of Radiation by Canopies and Stand Growth; Elsevier: Cambridge, UK, 1989; Volume 10, pp. 21–40. [Google Scholar]
  60. Hunt, E.R., Jr. Relationship between woody biomass and PAR conversion efficiency for estimating net primary production from NDVI. Int. J. Remote Sens. 1994, 15, 1725–1729. [Google Scholar]
  61. Goetz, S.J.; Prince, S.D. Modelling Terrestrial Carbon Exchange and Storage: Evidence and Implications of Functional Convergence in Light-use Efficiency. Adv. Ecol. Res. 1999, 28, 57–92. [Google Scholar]
  62. Goetz, S.J.; Prince, S.D.; Goward, S.N.; Thawley, M.M.; Small, J. Satellite remote sensing of primary production: An improved production efficiency modeling approach. Ecol. Model. 1999, 122, 239–255. [Google Scholar] [CrossRef]
  63. Gregory, P.J.; Tennant, D.; Belford, R.K. Root and shoot growth, and water and light use efficiency of barley and wheat crops grown on a shallow duplex soil in a mediterranean-type environment. Aust. J. Agric. Res. 1992, 43, 555–573. [Google Scholar] [CrossRef]
  64. Yunusa, I.A.M.; Siddique, K.H.M.; Belford, R.K.; Karimi, M.M. Effect of canopy structure on efficiency of radiation interception and use in spring wheat cultivars during the pre-anthesis period in a mediterranean-type environment. Field Crop. Res. 1993, 35, 113–122. [Google Scholar] [CrossRef]
  65. Ren, J.Q.; Chen, Z.X.; Tang, H.J.; Shi, R.X. Regional yield estimation for winter wheat based on net primary production model. Trans. CSAE 2006, 22, 111–116, (In Chinese with English abstract). [Google Scholar]
  66. Schlesinger, W.H. Biogeochemistry: An Analysis of Global Change. Science 1991, 253, 686–687. [Google Scholar]
  67. Lobell, D.B.; Hicke, J.A.; Asner, G.P.; Field, C.B.; Tucker, C.J.; Los, S.O. Satellite estimates of productivity and light use efficiency in the United States agriculture, 1982–1998. GCB Biol. 2002, 8, 722–735. [Google Scholar]
  68. Doraiswamy, P.C.; Sinclair, T.R.; Hollinger, S.; Akhmedov, B.; Stern, A.; Prueger, J. Application of MODIS derived parameters for regional crop yield assessment. Remote. Sens. Environ. 2005, 97, 192–202. [Google Scholar] [CrossRef]
  69. Launay, M.; Guérif, M. Assimilating remote sensing data into a crop model to improve predictive performance for spatial applications. Agric. Ecosyst. Environ. 2005, 111, 321–339. [Google Scholar] [CrossRef]
  70. Ren, J.Q.; Chen, Z.X.; Zhou, Q.B.; Tang, H.J. Regional yield estimation for winter wheat with MODIS-NDVI data in Shandong, China. Int. J. Appl. Earth Obs. Geoinf. 2008, 10, 403–413. [Google Scholar] [CrossRef]
  71. Wardlow, B.; Egbert, S.; Kastens, J. Analysis of time-series MODIS 250 m vegetation index data for crop classification in the U.S. Central Great Plains. Remote Sens. Environ. 2007, 108, 290–310. [Google Scholar] [CrossRef] [Green Version]
  72. Fritz, S.; Massart, M.; Savin, I.; Gallego, J.; Rembold, F. The use of MODIS data to derive acreage estimations for larger fields: A case study in the south-western Rostov region of Russia. Int. J. Appl. Earth Obs. Geoinf. 2008, 10, 453–466. [Google Scholar] [CrossRef]
  73. Awad, M.M. Toward Precision in Crop Yield Estimation Using Remote Sensing and Optimization Techniques. Agriculture 2019, 9, 54. [Google Scholar] [CrossRef]
  74. Li, A.; Bian, J.; Lei, G.; Huang, C. Estimating the Maximal Light Use Efficiency for Different Vegetation through the CASA Model Combined with Time-Series Remote Sensing Data and Ground Measurements. Remote. Sens. 2012, 4, 3857–3876. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Study area: (a) Beijing in China; (b) Tongzhou and Shunyi districts in Beijing. Green dots are sampling points.
Figure 1. Study area: (a) Beijing in China; (b) Tongzhou and Shunyi districts in Beijing. Green dots are sampling points.
Remotesensing 11 01088 g001
Figure 2. Winter wheat planting area in Tongzhou and Shunyi districts of Beijing in 2009.
Figure 2. Winter wheat planting area in Tongzhou and Shunyi districts of Beijing in 2009.
Remotesensing 11 01088 g002
Figure 3. Spatial distribution of winter wheat fPAR from March to May, 2009: (a) March; (b) April; (c) May.
Figure 3. Spatial distribution of winter wheat fPAR from March to May, 2009: (a) March; (b) April; (c) May.
Remotesensing 11 01088 g003
Figure 4. Spatial distribution of light use efficiency of winter wheat from March to May 2009: (a) March; (b) April; (c) May.
Figure 4. Spatial distribution of light use efficiency of winter wheat from March to May 2009: (a) March; (b) April; (c) May.
Remotesensing 11 01088 g004
Figure 5. Spatial distribution of winter wheat net primary production (NPP) from March to May, 2009: (a) March; (b) April; (c) May.
Figure 5. Spatial distribution of winter wheat net primary production (NPP) from March to May, 2009: (a) March; (b) April; (c) May.
Remotesensing 11 01088 g005
Figure 6. Mean winter wheat NPP as a function of time.
Figure 6. Mean winter wheat NPP as a function of time.
Remotesensing 11 01088 g006
Figure 7. Spatial distribution of winter wheat NPP averaged for the 2008–2009 growing season.
Figure 7. Spatial distribution of winter wheat NPP averaged for the 2008–2009 growing season.
Remotesensing 11 01088 g007
Figure 8. Spatial distribution of estimated winter wheat yield for the 2008–2009 growing season.
Figure 8. Spatial distribution of estimated winter wheat yield for the 2008–2009 growing season.
Remotesensing 11 01088 g008
Figure 9. Predicted winter wheat yield as a function of measured winter wheat yield for the entire study area for the 2008–2009 growing season.
Figure 9. Predicted winter wheat yield as a function of measured winter wheat yield for the entire study area for the 2008–2009 growing season.
Remotesensing 11 01088 g009
Table 1. NDVI and SR maximum and minimum values for winter wheat for each month of the growing season.
Table 1. NDVI and SR maximum and minimum values for winter wheat for each month of the growing season.
Vegetation IndicesMaximum and MinimumMonths 3
101112123456
NDVI 1MAX0.5570.6840.6870.5480.4300.4930.7570.8540.687
MIN0.2540.2650.2100.2240.1790.1900.2460.4330.239
SR 2MAX3.5195.3375.3863.4262.5122.9437.23612.6655.395
MIN1.6811.7201.5311.5781.4371.4691.6512.5281.629
1 The full name of NDVI is Normalized Difference Vegetation Index; 2 The full name of SR is Ratio Vegetation Index; 3 “Months” refers to the winter wheat growing season from October 2008 to June 2009.

Share and Cite

MDPI and ACS Style

Wang, Y.; Xu, X.; Huang, L.; Yang, G.; Fan, L.; Wei, P.; Chen, G. An Improved CASA Model for Estimating Winter Wheat Yield from Remote Sensing Images. Remote Sens. 2019, 11, 1088. https://doi.org/10.3390/rs11091088

AMA Style

Wang Y, Xu X, Huang L, Yang G, Fan L, Wei P, Chen G. An Improved CASA Model for Estimating Winter Wheat Yield from Remote Sensing Images. Remote Sensing. 2019; 11(9):1088. https://doi.org/10.3390/rs11091088

Chicago/Turabian Style

Wang, Yulong, Xingang Xu, Linsheng Huang, Guijun Yang, Lingling Fan, Pengfei Wei, and Guo Chen. 2019. "An Improved CASA Model for Estimating Winter Wheat Yield from Remote Sensing Images" Remote Sensing 11, no. 9: 1088. https://doi.org/10.3390/rs11091088

APA Style

Wang, Y., Xu, X., Huang, L., Yang, G., Fan, L., Wei, P., & Chen, G. (2019). An Improved CASA Model for Estimating Winter Wheat Yield from Remote Sensing Images. Remote Sensing, 11(9), 1088. https://doi.org/10.3390/rs11091088

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

Article Metrics

Back to TopTop