A GPP assimilation model for the southeastern Tibetan Plateau based on CO2 eddy covariance flux tower and remote sensing data

The gross primary production (GPP) at individual CO₂ eddy covariance flux tower sites (GPP_Tower) in Dali (DL), Wenjiang (WJ) and Linzhi (LZ) around the southeastern Tibetan Plateau were determined by the net ecosystem exchange of CO₂ (NEE) and ecosystem respiration (R_e). The satellite remote sensing-VPM model estimates of GPP values (GPP_MODIS) used the satellite-derived 8-day surface reflectance product (MOD09A1), including satellite-derived enhanced vegetation index (EVI) and land surface water index (LSWI). In this paper, we assembled a subset of flux tower data at these three sites to calibrate and test satellite-VPM model estimated GPP_MODIS, and introduced the satellite data and site-level environmental factors to develop four new assimilation models. The new assimilation models’ estimates of GPP values were compared with GPP_MODIS and GPP_Tower, and the final optimum model among the four assimilation models was determined and used to calibrate GPP_MODIS. The results showed that GPP_MODIS had similar temporal variations to the GPP_Tower, but GPP_MODlS were commonly higher in absolute magnitude than GPP_Tower with relative error (RE) about 58.85%. While, the assimilation models’ estimates of GPP values (GPP_MODEL) were much more closer to GPP_Tower with RE approximately 6.98%, indicating that the capacity of the simulation in the new assimilation model was greatly improved, the R² and root mean square error (RMSE) of the new assimilation model were 0.57–4.90% higher and 0.74–2.47 g C m⁻² s⁻¹ lower than those of the GPP_MODIS, respectively. The assimilation model was used to predicted GPP dynamics around the Tibetan Plateau and showed a reliable result compared with other researches. This study demonstrated the potential of the new assimilation model for estimating GPP around the Tibetan Plateau and the performances of site-level biophysical parameters in related to satellite-VPM model GPP.     

Yield Estimation Model and Water Productivity of Wheat Crop (Triticum aestivum) in an Irrigation Command Using Remote Sensing and GIS

Crop yield estimation has an important role on economy development and its accuracy and speed influence yield price and helps in deciding the excess or deficit production conditions. The water productivity evaluates the irrigation command through water use efficiency (WUE). Remote sensing (RS) and geographical information system (GIS) techniques were used for crop yield and water productivity estimation of wheat crop (Triticum aestivum) grown in Tarafeni South Main Canal (TSMC) irrigation command of West Bengal State in India. One IRS P6 image and four wide field sensor (WiFS) images for different months of winter season were used to determine the Normalized Difference Vegetation Index (NDVI) and Soil Adjusted Vegetation Index (SAVI) for area under wheat crop. The temporally and spatially distributed spectral growth profile and AREASUM of NDVI (A_NDVI) and SAVI (A_SAVI) with time after sowing of wheat crop were developed and correlated with actual crop yield of wheat (Y_act). The developed relationships between A_SAVI and Y_act resulted high correlation in comparison to that of A_NDVI. Using the developed model the RS based wheat yield (Y_RS) predicted from A_SAVI varied on entire TSMC irrigation command from 22.67 to 33.13 q ha⁻¹ respectively, which gave an average yield of 26.50 q ha⁻¹. The RS generated yield based water use efficiency (WUE_YRS) for water supplied from canal of TSMC irrigation command was found to be 6.69 kg ha⁻¹ mm⁻¹.     

Spatial and Temporal Variability of Net Primary Productivity (NPP) over Terrestrial Biosphere of India Using NOAA-AVHRR Based GloPEM Model

The monitoring of terrestrial carbon dynamics is important in studies related with global climate change. This paper presents results of the inter-annual variability of Net Primary Productivity (NPP) from 1981 to 2000 derived using observations from NOAA-AVHRR data using Global Production Efficiency Model (GloPEM). The GloPEM model is based on physiological principles and uses the production efficiency concept, in which the canopy absorption of photosynthetically active radiation (APAR) is used with a conversion “efficiency” to estimate Gross Primary Production (GPP). NPP derived from GloPEM model over India showed maximum NPP about 3,000 gCm⁻²year⁻¹ in west Bengal and lowest up to 500 gCm⁻²year⁻¹ in Rajasthan. The India averaged NPP varied from 1,084.7 gCm⁻²year⁻¹ to 1,390.8 gCm⁻²year⁻¹ in the corresponding years of 1983 and 1998 respectively. The regression analysis of the 20 year NPP variability showed significant increase in NPP over India (r = 0.7, F = 17.53, p < 0.001). The mean rate of increase was observed as 10.43 gCm⁻²year⁻¹. Carbon fixation ability of terrestrial ecosystem of India is increasing with rate of 34.3 TgC annually (t = 4.18, p < 0.001). The estimated net carbon fixation over Indian landmass ranged from 3.56 PgC (in 1983) to 4.57 PgC (in 1998). Grid level temporal correlation analysis showed that agricultural regions are the source of increase in terrestrial NPP of India. Parts of forest regions (Himalayan in Nepal, north east India) are relatively less influenced over the study period and showed lower or negative correlation (trend). Finding of the study would provide valuable input in understanding the global change associated with vegetation activities as a sink for atmospheric carbon dioxide.     

Estimation of gross primary production in wheat from in situ measurements

Gross primary production (GPP) is a parameter of significant importance for carbon cycle and climate change research. Remote sensing combined with other climate and meteorological data offers a convenient tool for large scale GPP estimation. This paper presents a study of GPP estimation using three methods with in situ measurements of canopy reflectance, LAI, and the photosynthetically active radiation (PAR). First, because LAI is considered as an indicator of the factor of absorbed PAR (f_APAR), it provides reasonable estimates of GPP for all types of wheat with coefficient of determination R² of 0.7353. The second method uses four kinds of vegetation indices (VIs) to estimate GPP because these indices are suggested to be reliable candidates in the estimation of light use efficiency (LUE). Good determination coefficients were acquired in estimating GPP with R² ranging from the lowest of 0.7604 for NDVI to the highest of 0.8505 for EVI. A new method was proposed for the estimation of GPP following the Monteith logic, which considering GPP as a product of VI × VI × PAR. Results indicated that this method can provide the best estimates of GPP as determination coefficient R² increased largely compared to the other two methods. EVI × EVI × PAR was demonstrated to be the most suitable for the estimation of GPP with the highest R² of 0.9207, which was about 10% larger as compared to GPP estimated from the single EVI. These results will be helpful for the development of new models of GPP estimation with all remote sensing inputs.     

Analysis of Seasonal Variability of Vegetation and Methane Concentration over India using SPOT-VEGETATION and ENVISAT-SCIAMACHY Data

This paper highlights the spatial and temporal variability of atmospheric columnar methane (CH₄) concentration over India and its correlation with the terrestrial vegetation dynamics. SCanning IMaging Absorption spectrometer for Atmospheric CHartographY (SCIAMACHY) on board ENVIronmental SATellite (ENVISAT) data product (0.5° × 0.5°) was used to analyze the atmospheric CH₄ concentration. Satellite Pour l'Observation de la Terre (SPOT)-VEGETATION sensor’s Normalized Difference Vegetation Index (NDVI) product, aggregated at 0.5° × 0.5° grid level for the same period (2004 and 2005), was used to correlate the with CH₄ concentration. Analysis showed mean monthly CH₄ concentration during the Kharif season varied from 1,704 parts per billion volume (ppbv) to 1,780 ppbv with the lowest value in May and the highest value in September. Correspondingly, mean NDVI varied from 0.28 (May) to 0.53 (September). Analysis of correlation between CH₄ concentration and NDVI values over India showed positive correlation (r = 0.76; n = 6) in Kharif season. Further analysis using land cover information showed characteristic low correlation in natural vegetation region and high correlation in agricultural area. Grids, particularly falling in the Indo-Gangetic Plains showed positive correlation. This could be attributed to the rice crop which is grown as a predominant crop during this period. The CH₄ concentration pattern matched well with growth pattern of rice with the highest concentration coinciding with the peak growth period of crop in the September. Characteristically low correlation was observed (r = 0.1; n = 6) in deserts of Rajasthan and forested Himalayan ecosystem. Thus, the paper emphasizes the synergistic use of different satellite based data in understanding the variability of atmospheric CH₄ concentration in relation to vegetation.     

Modelling the Gross Primary Productivity of West Africa with the Regional Biomass Model RBM+, using optimized 250 m MODIS FPAR and fractional vegetation cover information

Global warming associated with climate change is one of the greatest challenges of today’s world. Increasing emissions of the greenhouse gas CO₂ are considered as a major contributing factor to global warming. One regulating factor of CO₂ exchange between atmosphere and land surface is vegetation. Measurements of land cover changes in combination with modelling the Gross Primary Productivity (GPP) can contribute to determine important sources and sinks of CO₂.The aim of this study is to accurately model the GPP for a region in West Africa with a spatial resolution of 250 m, and the differentiation of GPP based on woody and herbaceous vegetation. For this purpose, the Regional Biomass Model (RBM) was applied, which is based on a Light Use Efficiency (LUE) approach. The focus was on the spatial enhancement of the RBM from the original 1000–250 m spatial resolution (RBM+). The adaptation to the 250 m scale included the modification of two main input parameters: (1) the fraction of absorbed Photosynthetically Active Radiation (FPAR) based on the 1000 m MODIS MOD15A2 FPAR product which was downscaled to 250 m using MODIS NDVI time series; (2) the fractional cover of woody and herbaceous vegetation, which was improved by using a multi-scale approach. For validation and regional adjustments of GPP and the input parameters, in situ data from a climate station and eddy covariance measurements were integrated.The results of this approach show that the input parameters could be improved significantly: downscaling considerably reduces data gaps of the original FPAR product and the improved dataset differed less than 5.0% from the original data for cloud free regions. The RMSE of the fractional vegetation cover varied between 5.1 and 12.7%. Modelled GPP showed a slight overestimation in comparison to eddy covariance measurements. The in situ data was exceeded by 8.8% for 2005 and by 2.0% for 2006. The model results were converted to NPP and also agreed well with previous NPP measurements reported from different studies. Altogether a high accuracy and suitability of the regionally adjusted and downscaled model RBM+ can be concluded. The differentiation between vegetation growth forms allows a separation of long-term and short-term carbon storage based on woody and herbaceous vegetation, respectively.     

The geopotential value <Emphasis Type="Italic">W</Emphasis><Subscript>0</Subscript> for specifying the relativistic atomic time scale and a global vertical reference system

The TOPEX/Poseidon (T/P) satellite alti- meter mission marked a new era in determining the geopotential constant W ₀. On the basis of T/P data during 1993–2003 (cycles 11–414), long-term variations in W ₀ have been investigated. The rounded value W ₀ = 62636856.0 ± 0.5) m ² s ⁻² has already been adopted by the International Astronomical Union for the definition of the constant L _G = W ₀/c ² = 6.969290134 × 10⁻¹⁰ (where c is the speed of light), which is required for the realization of the relativistic atomic time scale. The constant L _G , based on the above value of W ₀, is also included in the 2003 International Earth Rotation and Reference Frames Service conventions. It has also been suggested that W ₀ is used to specify a global vertical reference system (GVRS). W ₀ ensures the consistency with the International Terrestrial Reference System, i.e. after adopting W ₀, along with the geocentric gravitational constant (GM), the Earth’s rotational velocity (ω) and the second zonal geopotential coefficient (J ₂) as primary constants (parameters), then the ellipsoidal parameters (a,α) can be computed and adopted as derived parameters. The scale of the International Terrestrial Reference Frame 2000 (ITRF2000) has also been specified with the use of W ₀ to be consistent with the geocentric coordinate time. As an example of using W ₀ for a GVRS realization, the geopotential difference between the adopted W ₀ and the geopotential at the Rimouski tide-gauge point, specifying the North American Vertical Datum 1988 (NAVD88), has been estimated.     

Seasonal variation of carbon fluxes in a sparse savanna in semi arid Sudan

Background

Large spatial, seasonal and annual variability of major drivers of the carbon cycle (precipitation, temperature, fire regime and nutrient availability) are common in the Sahel region. This causes large variability in net ecosystem exchange and in vegetation productivity, the subsistence basis for a major part of the rural population in Sahel. This study compares the 2005 dry and wet season fluxes of CO₂ for a grass land/sparse savanna site in semi arid Sudan and relates these fluxes to water availability and incoming photosynthetic photon flux density (PPFD). Data from this site could complement the current sparse observation network in Africa, a continent where climatic change could significantly impact the future and which constitute a weak link in our understanding of the global carbon cycle.

Results

The dry season (represented by Julian day 35–46, February 2005) was characterized by low soil moisture availability, low evapotranspiration and a high vapor pressure deficit. The mean daily NEE (net ecosystem exchange, Eq. 1) was -14.7 mmol d^-1 for the 12 day period (negative numbers denote sinks, i.e. flux from the atmosphere to the biosphere). The water use efficiency (WUE) was 1.6 mmol CO₂ mol H₂O^-1 and the light use efficiency (LUE) was 0.95 mmol CO₂ mol PPFD^-1. Photosynthesis is a weak, but linear function of PPFD. The wet season (represented by Julian day 266–273, September 2005) was, compared to the dry season, characterized by slightly higher soil moisture availability, higher evapotranspiration and a slightly lower vapor pressure deficit. The mean daily NEE was -152 mmol d^-1 for the 8 day period. The WUE was lower, 0.97 mmol CO₂ mol H₂O^-1 and the LUE was higher, 7.2 μmol CO₂ mmol PPFD^-1 during the wet season compared to the dry season. During the wet season photosynthesis increases with PPFD to about 1600 μmol m^-2s^-1 and then levels off.

Conclusion

Based on data collected during two short periods, the studied ecosystem was a sink of carbon both during the dry and wet season 2005. The small sink during the dry season is surprising and similar dry season sinks have not to our knowledge been reported from other similar savanna ecosystems and could have potential management implications for agroforestry. A strong response of NEE versus small changes in plant available soil water content was found. Collection and analysis of flux data for several consecutive years including variations in precipitation, available soil moisture and labile soil carbon are needed for understanding the year to year variation of the carbon budget of this grass land/sparse savanna site in semi arid Sudan.     

Quantification of LST and CO2 levels using Landsat-8 thermal bands on urban environment

Land surface temperature (LST) is an important aspect in global to regional change studies, for control of climate change and balancing of high temperature. Urbanization is one of the influencing factors increasing land surface and atmospheric temperature, by the emission of greenhouse gases (e.g. CO₂, NO and methane). In the present study, LST was derived from Landsat-8 of multitemporal data sets to analyse the spatial structure of the urban thermal environment in relation to the urban surface characteristics and land use–land cover (LULC). LST is influenced by the greenhouse gases i.e. CO₂ plays an important role in increasing the earth’s surface temperature. In order to provide the evidence of influence of CO₂ on LST, the relationship between LST, air temperature and CO₂ was analysed. Landsat-8 satellite has two thermal bands, 10 and 11. These bands were used to accurately to calculate the temperature over the study area. Results showed that the strength of correlation between ground monitoring data and satellite data was high. Based on correlation values of each month April (R² = 0.994), May (R² = 0.297) and June (R² = 0.934), observed results show that band 10 was significantly correlating with air temperature. Relationship between LST and CO₂ levels were obtained from linear regression analysis. band 11 was correlating significantly with CO₂ values in each of the months April (R² = 0.217), May (R² = 0.914) and June, (R² = 0.934), because band 11 is closer to the 15-micron band of CO₂. From the results, it was observed that band 10 can be used for calculating air temperature and band 11 can be used for estimation of greenhouse gases.     

Study of the optimum time span for distinguishing <Emphasis Type="Italic">Avena ludoviciana</Emphasis> from wheat crop based on their spectral reflectance characteristics

The study to establish the optimum time span for distinguishing Avena ludoviciana from wheat crop based on their spectral signatures was carried out at Student’s Research Farm, Department of Agronomy during 2006–07 and 2007–08. The experimental sites during both the seasons were sandy loam in texture, with normal soil reaction and electrical conductivity, low in organic carbon and available nitrogen and medium in available phosphorus and potassium. The experiment was laid out in randomized block design with four replications and consisting of twelve treatments comprising 0, 10, 15, 25, 50, 75, 100, 125, 150, 200, 250 plants m⁻² and a pure Avena ludoviciana plot (T_max). The results revealed that in all the treatments irrespective of wheat and weeds, the red reflectance (%) value decreased from 34 to 95 DAS (days after sowing) in 2006–07 and 45 DAS to 100 DAS during 2007–08, and thereafter a sharp increase was observed in all the treatments. This trend might be due to increased chlorophyll index after 34 DAS as red reflectance was reduced by chlorophyll absorption. Among all the treatments, T_max (Pure Avena ludoviciana plot) had the highest red reflectance and T₀ (Pure wheat plot) had a lowest value of red reflectance during both the years. The highest value of IR reflectance was obtained at 95 DAS (2006–07) and 70 DAS (2007–08) in all the treatments. IR reflectance of wheat crop ranged between 24.61 and 61.21 per cent during 2006–07 and 27.33 and 67.3 per cent during 2007–08. However, IR reflectance values declined after 95 DAS and 70 DAS up to harvesting during 2006–07 and 2007–08. This lower reflectance may have been due to the onset of senescence. The highest RR and NDVI values were recorded under pure wheat treatment and minimum under pure weed plots. This may be due to dark green colour and better vigor of the wheat as compared to Avena ludoviciana. It was observed that by using RR and NDVI, pure wheat can be distinguished from pure populations of Avena ludoviciana after 34 DAS and different levels of weed populations can be discriminated amongst themselves from 68 DAS up to 107 DAS during both the years of investigation.     

Effects of the partitioning of diffuse and direct solar radiation on satellite-based modeling of crop gross primary production

Modeling crop gross primary production (GPP) is critical to understanding the carbon dynamics of agro-ecosystems. Satellite-based studies have widely used production efficiency models (PEM) to estimate cropland GPP, wherein light use efficiency (LUE) is a key model parameter. One factor that has not been well considered in many PEMs is that canopy LUE could vary with illumination conditions. This study investigates how the partitioning of diffuse and direct solar radiation influences cropland GPP using both flux tower and satellite data. The field-measured hourly LUE under cloudy conditions was 1.50 and 1.70 times higher than that under near clear-sky conditions for irrigated corn and soybean, respectively. We applied a two-leaf model to simulate the canopy radiative transfer process, where modeled photosynthetically active radiation (PAR) absorbed by canopy agreed with tower measurements (R² = 0.959 and 0.914 for corn and soybean, respectively). Derived canopy LUE became similar after accounting for the impact of light saturation on leaf photosynthetic capacity under varied illumination conditions. The impacts of solar radiation partitioning on satellite-based modeling of crop GPP was examined using vegetation indices (VI) derived from MODIS data. Consistent with the field modeling results, the relationship between daily GPP and PAR × VI under varied illumination conditions showed different patterns in terms of regression slope and intercept. We proposed a function to correct the influences of direct and diffuse radiation partitioning and the explained variance of flux tower GPP increased in all experiments. Our results suggest that the non-linear response of leaf photosynthesis to light absorption contributes to higher canopy LUE on cloudy days than on clear days. We conclude that accounting for the impacts of solar radiation partitioning is necessary for modeling crop GPP on a daily or shorter basis.     

Camera derived vegetation greenness index as proxy for gross primary production in a low Arctic wetland area

The Arctic is experiencing disproportionate warming relative to the global average, and the Arctic ecosystems are as a result undergoing considerable changes. Continued monitoring of ecosystem productivity and phenology across temporal and spatial scales is a central part of assessing the magnitude of these changes. This study investigates the ability to use automatic digital camera images (DCIs) as proxy data for gross primary production (GPP) in a complex low Arctic wetland site. Vegetation greenness computed from DCIs was found to correlate significantly (R² = 0.62, p < 0.001) with a normalized difference vegetation index (NDVI) product derived from the WorldView-2 satellite. An object-based classification based on a bi-temporal image composite was used to classify the study area into heath, copse, fen, and bedrock. Temporal evolution of vegetation greenness was evaluated and modeled with double sigmoid functions for each plant community. GPP at light saturation modeled from eddy covariance (EC) flux measurements were found to correlate significantly with vegetation greenness for all plant communities in the studied year (i.e., 2010), and the highest correlation was found between modeled fen greenness and GPP (R² = 0.85, p < 0.001). Finally, greenness computed within modeled EC footprints were used to evaluate the influence of individual plant communities on the flux measurements. The study concludes that digital cameras may be used as a cost-effective proxy for potential GPP in remote Arctic regions.     

Validation of a Crop Yield and CO<Subscript>2</Subscript> Fixation Model over Southeast Asia by Carbon Partitioning in Grain Plants

A photosynthetic-sterility model for grain production monitoring has been developed and validated under the background of climate change and Asian economic growth in developing countries. This paper presents an application of the model to evaluate carbon-fixation rates in yields of paddy rice, winter wheat, and maize in Asia. The validation of the model is based on carbon partitioning in grain plants. The carbon hydrate in grains has the same chemical formula as that of cellulose in grain vegetation. The partitioning of carbon in plants can validate fixation amounts of computed carbon using a satellite-based photosynthesis model. The model estimates the photosynthesis fixation of rice reasonably in Japan and China. Results were validated through examination of carbon in grains, but the model tends to underestimate results for winter wheat and maize. This study also provides daily distributions of the PSN, which is the CO₂ fixation in Asian areas combined with a land-cover distribution classified from MODIS data, NDVI from SPOT VEGETATION, and meteorological re-analysis data by European Centre for Medium-Range Forecasts (ECMWF). The mean CO₂ and carbon fixation rates in paddy areas were 25.92 (t CO₂/ha) and 5.28 (t C/ha) in Japan, respectively. Comparisons between the model’s values and MODIS seasonal PSNs show similar trends. The writers are preparing to compare computed photosynthesis rates with observed AsiaFlux data for the validation of this model at field sites of paddy, grassland and forests in Japan and Asian countries.     

Assessment of Vegetation Response to Drought in Nebraska Using Terra-MODIS Land Surface Temperature and Normalized Difference Vegetation Index

Eight-day composite Terra-MODIS cumulative LST and NDVI timeseries data were used to analyze the responses of crop and grassland cover types to drought in Nebraska. Four hundred ninety 1 km pixels that included irrigated and non-irrigated corn and soybeans and three grassland cover types were selected across the state of Nebraska. Statistical analyses revealed that the majority of the land cover pixels experienced significantly higher daytime and nighttime LSTs and lower NDVI during the drought-year growing season (p < 0.01). Among the land cover types analyzed, grassland experienced the highest increase in daytime LST and decrease in NDVI.     

Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops

Leaf area index (LAI) and biomass are important indicators of crop development and the availability of this information during the growing season can support farmer decision making processes. This study demonstrates the applicability of RapidEye multi-spectral data for estimation of LAI and biomass of two crop types (corn and soybean) with different canopy structure, leaf structure and photosynthetic pathways. The advantages of Rapid Eye in terms of increased temporal resolution (∼daily), high spatial resolution (∼5 m) and enhanced spectral information (includes red-edge band) are explored as an individual sensor and as part of a multi-sensor constellation. Seven vegetation indices based on combinations of reflectance in green, red, red-edge and near infrared bands were derived from RapidEye imagery between 2011 and 2013. LAI and biomass data were collected during the same period for calibration and validation of the relationships between vegetation indices and LAI and dry above-ground biomass. Most indices showed sensitivity to LAI from emergence to 8 m²/m². The normalized difference vegetation index (NDVI), the red-edge NDVI and the green NDVI were insensitive to crop type and had coefficients of variations (CV) ranging between 19 and 27%; and coefficients of determination ranging between 86 and 88%. The NDVI performed best for the estimation of dry leaf biomass (CV = 27% and r² = 090) and was also insensitive to crop type. The red-edge indices did not show any significant improvement in LAI and biomass estimation over traditional multispectral indices. Cumulative vegetation indices showed strong performance for estimation of total dry above-ground biomass, especially for corn (CV ≤ 20%). This study demonstrated that continuous crop LAI monitoring over time and space at the field level can be achieved using a combination of RapidEye, Landsat and SPOT data and sensor-dependant best-fit functions. This approach eliminates/reduces the need for reflectance resampling, VIs inter-calibration and spatial resampling.     

Seasonal and inter-annual variability of atmosphere CO<Subscript>2</Subscript> based on NOAA Carbon Tracker analysis and satellite observations

A study on seasonal and inter-annual variability of the atmospheric CO₂ is carried out based on National Oceanic and Atmospheric Administration Carbon Tracker (NOAACT) re-analysis and satellite measurements of mid-troposphere CO₂ by Atmosphere Infrared Sounder on board NASA’s Aqua and lower troposphere CO₂ by Greenhouse-gas Observing Satellite. Seasonal and non-seasonal components of each time series were extracted by means of least square based harmonic analysis procedure. The data of surface CO₂ fluxes used in the NOAACT are also analyzed to examine its relationship with the atmosphere CO₂ variability at different time scales. There exists good consistency between NOAACT analysis and satellite observations in their respective seasonal harmonics and climatology. Surface layer CO₂ exhibits large climatological mean over the regions of major anthropogenic sources together with strong seasonal cycle over the humid and cold climatic terrestrial regions especially over the northern hemisphere. Existence of high coherency with the different components of the surface fluxes shows that surface layer atmosphere CO₂ seasonality is primarily contributed from the terrestrial ecosystem exchanges and secondarily by anthropogenic and oceanic exchanges. The mid-troposphere CO₂ exhibits large values associated with climatology and amplitudes of semi-annual and annual cycles over the northern extra tropics and Polar Regions along with a gradual decreasing trend from northern to southern hemisphere. Inter-annual variability of atmospheric CO₂ in the NOAACT in some extent is consistent with the satellite observations. Large scale circulation patterns, its fluctuations associated with ENSO events and large scale ecosystem disturbances have significant influence on the inter-annual variability.     

Improved maize cultivated area estimation over a large scale combining MODIS–EVI time series data and crop phenological information

The accurate and timely information of crop area is vital for crop production and food security. In this study, the Enhanced Vegetation Index (EVI) data from MODerate resolution Imaging Spectroradiometer (MODIS) integrated crop phenological information was used to estimate the maize cultivated area over a large scale in Northeast China. The fine spatial resolution China’s Environment Satellite (HJ-1 satellite) images and the support vector machine (SVM) algorithm were employed to discriminate distribution of maize in the reference area. The mean MODIS–EVI time series curve of maize was extracted in the reference area by using multiple periods MODIS–EVI data. By analysing the temporal shift of crop calendars from northern to southern parts in Northeast China, the lag value was derived from phenological data of twenty-one agro-meteorological stations; here integrating with the mean MODIS–EVI time series image of maize, a standard MODIS–EVI time series image of maize was obtained in the whole study area. By calculating mean absolute distances (MAD) map between standard MODIS–EVI image and mean MODIS–EVI time series images, and setting appropriate thresholds in three provinces, the maize cultivated area was extracted in Northeast China. The results showed that the overall classification accuracy of maize cultivated area was approximately 79%. At the county level, the MODIS-derived maize cultivated area and statistical data were well correlated (R² = 0.82, RMSE = 283.98) over whole Northeast China. It demonstrated that MODIS–EVI time series data integrated with crop phenological information can be used to improve the extraction accuracy of crop cultivated area over a large scale.     

Evaluating a satellite-based seasonal evapotranspiration product and identifying its relationship with other satellite-derived products and crop yield: A case study for Ethiopia

Satellite-derived evapotranspiration anomalies and normalized difference vegetation index (NDVI) products from Moderate Resolution Imaging Spectroradiometer (MODIS) data are currently used for African agricultural drought monitoring and food security status assessment. In this study, a process to evaluate satellite-derived evapotranspiration (ETa) products with a geospatial statistical exploratory technique that uses NDVI, satellite-derived rainfall estimate (RFE), and crop yield data has been developed. The main goal of this study was to evaluate the ETa using the NDVI and RFE, and identify a relationship between the ETa and Ethiopia’s cereal crop (i.e., teff, sorghum, corn/maize, barley, and wheat) yields during the main rainy season. Since crop production is one of the main factors affecting food security, the evaluation of remote sensing-based seasonal ETa was done to identify the appropriateness of this tool as a proxy for monitoring vegetation condition in drought vulnerable and food insecure areas to support decision makers. The results of this study showed that the comparison between seasonal ETa and RFE produced strong correlation (R² > 0.99) for all 41 crop growing zones in Ethiopia. The results of the spatial regression analyses of seasonal ETa and NDVI using Ordinary Least Squares and Geographically Weighted Regression showed relatively weak yearly spatial relationships (R² < 0.7) for all cropping zones. However, for each individual crop zones, the correlation between NDVI and ETa ranged between 0.3 and 0.84 for about 44% of the cropping zones. Similarly, for each individual crop zones, the correlation (R²) between the seasonal ETa anomaly and de-trended cereal crop yield was between 0.4 and 0.82 for 76% (31 out of 41) of the crop growing zones. The preliminary results indicated that the ETa products have a good predictive potential for these 31 identified zones in Ethiopia. Decision makers may potentially use ETa products for monitoring cereal crop yields and early warning of food insecurity during drought years for these identified zones.     

Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data

A sufficient number of satellite acquisitions in a growing season are essential for deriving agronomic indicators, such as green leaf area index (GLAI), to be assimilated into crop models for crop productivity estimation. However, for most high resolution orbital optical satellites, it is often difficult to obtain images frequently due to their long revisit cycles and unfavorable weather conditions. Data fusion algorithms, such as the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Enhanced STARFM (ESTARFM), have been developed to generate synthetic data with high spatial and temporal resolution to address this issue. In this study, we evaluated the approach of assimilating GLAI into the Simple Algorithm for Yield Estimation model (SAFY) for winter wheat biomass estimation. GLAI was estimated using the two-band Enhanced Vegetation Index (EVI2) derived from data acquired by the Operational Land Imager (OLI) onboard the Landsat-8 and a fusion dataset generated by blending the Moderate-Resolution Imaging Spectroradiometer (MODIS) data and the OLI data using the STARFM and ESTARFM models. The fusion dataset had the temporal resolution of the MODIS data and the spatial resolution of the OLI data. Key parameters of the SAFY model were optimised through assimilation of the estimated GLAI into the crop model using the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm. A good agreement was achieved between the estimated and field measured biomass by assimilating the GLAI derived from the OLI data (GLAI_L) alone (R² = 0.77 and RMSE = 231 g m⁻²). Assimilation of GLAI derived from the fusion dataset (GLAI_F) resulted in a R² of 0.71 and RMSE of 193 g m⁻² while assimilating the combination of GLAI_L and GLAI_F led to further improvements (R² = 0.76 and RMSE = 176 g m⁻²). Our results demonstrated the potential of using the fusion algorithms to improve crop growth monitoring and crop productivity estimation when the number of high resolution remote sensing data acquisitions is limited.