Developing indices of temporal dispersion and continuity to map natural vegetation

An accurate and updated natural vegetation map is imperative for sustainable environmental management. This paper proposed a novel natural vegetation mapping algorithm based on time series images. Several indices of temporal dispersion and continuity were established for this purpose: low density (LD), medium density (MD), high density (HD) and medium continuity (MC). These indices were developed based on the particular percentiles-determined section of the EVI2 temporal profiles obtained through continuous wavelet transform. The natural vegetation was generally characterized as with lower temporal dispersion and greater temporal continuity compared with agricultural crops. The proposed methodology incorporated the indices of temporal dispersion and continuity and was applied to 13 provinces in central East China based on 500 m 8-day composite Moderate Resolution Imaging Spectroradiometer (MODIS) Enhanced Vegetation Index with two bands (EVI2) in 2013. An overall accuracy of 92.97% was obtained when compared with 2715 ground truth sites. There was also a good agreement (kappa index = 0.8049) on the distribution and areas of different vegetation types between the MODIS-estimated image and the Landsat 8 OLI interpreted data on two test regions. This study demonstrated the efficiency of the transform and metric integrated time series classification approaches in the fields of land and vegetation cover mapping.     

Effects of reclamation of natural wetlands to a rice paddy on dissolved carbon dynamics in the Sanjiang Plain,Northeastern China

Natural wetlands play an important role in the global carbon cycle, and loss of dissolved carbon through water has been indicated as one of the most important carbon sources for riverine ecosystems. During the last century, a large natural wetland area was reported to be converted to other land use types such as rice paddy land around the world. In this study, we explored the dynamics of dissolved organic carbon (DOC) and dissolved inorganic carbon (DIC) in two natural freshwater wetlands and a rice paddy field, which was reclaimed from the natural wetlands in the Sanjiang Plain, Northeastern China, during the growing season (May–October) of 2009. The DOC and DIC concentrations in the two ecosystems were significantly different (P < 0.05). The mean DOC concentrations during the growing season in the surface water of the Deyeuxia angustifolia and Carex lasiocarpa wetlands were 49.88 ± 5.44 and 27.97 ± 1.69 mg/L, respectively, while it was only 8.63 ± 2.54 mg/L in the rice paddy field. Specific ultra-violet light absorption at 254 nm (SUVA₂₅₄) of DOC increased by an average of 19.54% in the surface water from the natural wetlands to rice paddy, suggesting that DOC mobilized in the natural wetlands was more aromatic than that in the rice paddy field. The mean DIC concentration in surface water of the rice paddy was 5.25 and 5.04 times higher than that in the natural D. angustifolia and C. lasiocarpa wetlands, respectively. The average ratio of DIC to dissolved total carbon (DTC) for the water sampled from the artificial drainage ditch in the rice paddy field was 61.82%, while it was 14.75% from the nearby channel of natural wetlands. The significant differences in dissolved carbon concentration in surface water and channels originating from different land use types suggested that reclamation of natural wetlands to rice paddy field would reduce DOC runoff and increase the DIC concentration to adjacent watersheds. Our study results for the changed pattern in dissolved carbon after the natural wetland was transformed to paddy field could have important implications for studying the impacts of the large-scale land use change to carbon cycle and management.     

Effect of different fertilization modes on soil organic carbon sequestration in paddy fields in South China: A meta-analysis

Evidences have shown that fertilizer application could affect soil organic carbon (SOC) content in cropland. South China is the main production region of rice, in which many different fertilization practices have been widely used during the last several decades, but the effect of different fertilization modes on SOC sequestration in paddy fields in this region was scarcely studied. Based on 274 samples taken from 44 paddy field experimental sites in South China, a meta-analysis was performed to quantify the relative annual change of SOC content (RAC) and SOC sequestration duration in paddy fields under five fertilization modes (inorganic nitrogen fertilization, N; inorganic nitrogen and phosphorus fertilization, NP; inorganic nitrogen, phosphorus, and potassium fertilization, NPK; organic fertilization, O; and inorganic combined with organic fertilization, OF). The results showed that the RAC under the five fertilization modes was 0–0.4 g kg⁻¹ yr⁻¹, with increments of 0.19 and 0.23 g kg⁻¹ yr⁻¹ in double cropping systems and triple cropping systems, respectively. The RAC under the O and OF treatments was higher than that under the N, NP and NPK treatments and was the highest (0.32 g kg⁻¹ yr⁻¹) under OF among the five treatments. The SOC accumulation rate decreased with time, and the SOC sequestration duration under N, NP, NPK, O and OF was approximately 23, 28, 46, 64 and 55 years, respectively. The accumulation enhancement rate of SOC over the whole SOC sequestration period under the N, NP, NPK, O and OF treatments was approximately 12.9%, 23.4%, 29.3%, 47.2% and 55.1%, respectively. The OF treatment performed the largest potential for SOC sequestration in paddy fields in South China under long-term conditions.     

Natural parasitism of Romanomermis iyengari (Welch) (Nematoda: Mermithidae) on various species of mosquitoes breeding in rice fields in Pondicherry,India

Natural habitat of Romanomermis iyengari, an entomophilic mermithid nematode, is rice field. Parasitism of this nematode on mosquito larvae breeding in rice fields in Pondicherry, India, has been studied in 1979. The present study was to find out its persistence in the same habitat, after many years of agricultural practices. Evaluation of parasitism on mosquito larvae was carried out during two rice cultivation seasons at a schedule of 1 week prior, and 1 and 2 weeks after transplantation of paddy. Larval stages of mosquito species such as Culex tritaeniorhynchus, Culex vishnui, Anopheles vagus, and Anopheles subpictus showed parasitism of R. iyengari. Parasitism was the highest on C. tritaeniorhynchus with a maximum of 15.38% during season I and 17.85% during season II. Followed by this, it was on A. vagus (10.82% and 9.43%), C. vishnui (7.37% and 8.69%), and A. subpictus (3.70% and 6.36%). The overall density of mosquito larvae and level of parasitism of R. iyengari on them during season II was significantly higher than that of season I. The parasitism was maximum at 1 week before transplantation and showed a decreasing trend towards 2 weeks after transplantation, during both the seasons. The level of natural parasitism of R. iyengari during the present investigation was higher than that observed in the same ecosystem almost 35 years ago. This is clearly an indication towards the capability of the nematode to withstand several years of agricultural practices, including pesticide and agrichemical application. Hence, if applied at sufficient level and frequency, this nematode could act as a successful biocontrol agent against mosquitoes breeding in rice fields.     

Insecticide resistance monitoring and correlation analysis to select appropriate insecticides against Nilaparvata lugens (Stål), a migratory pest in Korea

Nilaparvata lugens Stål is one of the important migratory pests of rice paddy fields in Korea. Resistance levels to nine insecticides were monitored in 12 local strains and correlation analysis was conducted to determine cross-resistance relationships among the tested insecticides. The local strains revealed 1.3- to 28.0-, 1.6- to 6.0-, 2.8- to 237.0-, 0.6- to 0.9-, and 0.7- to 1.3-fold resistance to carbamates, organophosphates, neonicotinoids, fipronil and etofenprox, respectively. Organophosphate insecticides revealed moderate correlations with benzofuranyl methylcarbamate (r = 0.566–0.614, p > 0.01). Three neonicotinoids were not correlated with each other, but imidacloprid and clothianidin were moderately correlated with several benzofuranyl methylcarbamate and organophosphate insecticides (r = 0.590–0.705, p > 0.05), indicating that unknown common factors (such as detoxification enzymes) might contribute to resistance to both insecticides. Fipronil and etofenprox exhibited low levels of resistance and cross-resistance with other insecticides, suggesting their potential as an effective insecticide for field application. Resistance level monitoring and correlation analysis would be valuable for the selection of appropriate insecticides to control insecticide-resistant N. lugens, a typical migratory pest in Korea.     

Estimating root zone soil moisture at distant sites using MODIS NDVI and EVI in a semi-arid region of southwestern USA

This study investigates the potential of using Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) to estimate root zone soil moisture at native in-situ measured sites, and at distant sites under the same climatic setting. We obtained in-situ data from Soil Climate Analysis Network (SCAN) sites near the Texas-New Mexico border area, and NDVI and EVI products from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on board the Terra satellite. Results show that soil moisture values of the same depth are highly correlated (r = 0.53 to 0.85) among sites as far as 150 km apart, and that NDVI and EVI are highly correlated at the same site (r = 0.87 to 0.91). Correlation based on raw time series of NDVI and soil moisture is in general higher than that based on deseasonalized time series at every depth. The correlation reaches maximum value when vegetation index (VI) lags soil moisture by 5 to 10 days. NDVI shows a slightly higher correlation with soil moisture than EVI does by using the deseasonalized time series of NDVI and soil moisture. It is found that deseasonalized time series of NDVI and soil moisture are correlated at native sites (r = 0.33 to 0.77), but not at sites where soil moisture is very low. Regression analysis was conducted using deseasonalized time series soil moisture and deseasonalized time series NDVI with a 5-day time lag. Regression models developed at one site and applied to a similar distant site can estimate soil moistures, accounting for 50–88% of the variation in observed soil moistures.     

Carbon footprints of rice production in five typical rice districts in China

Five provinces located in the five main rice-growing regions in China were selected as study areas, which were Jiangsu, Heilongjiang, Sichuan, Guangdong and Hunan province respectively in the middle and lower reaches of the Yangtze River, northern, southwest, southern and central rice districts. Carbon footprints of rice production in these five provinces were calculated through the life cycle assessment method using governmental statistical data, industrial standards and relevant technical data separately. Material and energy consumptions were estimated, key stages of energy consumptions and carbon emissions were identified as well. Moreover, improving measurements had been suggested correspondingly. The results indicated that: the energy consumptions of rice production in these five provinces ranked as following (high to low): Guangdong, Heilongjiang, Hunan, Sichuan and Jiangsu. The carbon footprints of rice production were 2504.20 kg carbon dioxide equation per ton rice (kgCO₂-eq./t) (Guangdong province), 2326.47 kgCO₂-eq./t (Hunan province), 1889.97 kgCO₂-eq./t (Heilongjiang province), 1538.90 kgCO₂-eq./t (Sichuan province) and 1344.92 kgCO₂-eq./t (Jiangsu province) respectively. Reducing the quantities of urea and using the intermittent irrigation method could decrease energy consumption as well as carbon footprint.     

Morphological and physiological responses of nine southern U.S. rice cultivars differing in their tolerance to enhanced ultraviolet-B radiation

The impact of climatic change on crop production is a major global concern. One of the climatic factors, ultraviolet-B radiation (UV-B; 280–320 nm), which is increasing as a result of depletion of the global stratospheric ozone layer, can alter crop productivity. As the initial step in development of UV-B tolerant rice cultivars for the southern U.S., in this study we screened popular southern U.S. rice cultivars for variation in tolerance to elevated UV-B radiation with respect to morphological, phenological and physiological parameters. Plants grown in the greenhouse at the Texas AgriLife Research and Extension Center in Beaumont, Texas, U.S. were exposed to 0, 8 or 16 kJ m⁻² day⁻¹ UV-B radiation for 90 days. Our results showed differences among southern US rice cultivars in response to UV-B treatments with respect to leaf photosynthetic rate (P_n), leaf phenolic concentration, pollen germination (PG), spikelet fertility (SF), leaf number, leaf area, and yield. For most of the cultivars, plants exposed to enhanced UV-B radiation showed decreased P_n, PG, SF and yield and increased spikelet abortion and leaf phenolic concentration compared to the plants grown in a UV-B-free environment. In this study, cultivar ‘Clearfield XL729’ performed better than the other cultivars under enhanced UV-B radiation.     

Mapping riparian condition indicators in a sub-tropical savanna environment from discrete return LiDAR data using object-based image analysis

Mapping, monitoring and managing the environmental condition of riparian zones are major focus areas for local and state governments in Australia. New remotely sensed data techniques that can provide the required mapping accuracies, complete spatial coverage and processing and mapping transferability are currently being developed for use over large spatial extents. The research objective was to develop and apply an approach for mapping riparian condition indicators using object-based image analysis of airborne Light Detection and Ranging (LiDAR) data. The indicators assessed were: streambed width; riparian zone width; plant projective cover (PPC); longitudinal continuity; coverage of large trees; vegetation overhang; and stream bank stability. LiDAR data were captured on 15 July 2007 for two 5 km stretches along Mimosa Creek in Central Queensland, Australia. Field measurements of riparian vegetation structural and landform parameters were obtained between 28 May and 5 June 2007. Object-based approaches were developed for mapping each riparian condition indicator from the LiDAR data. The validation and empirical modelling results showed that the object-based approach could be used to accurately map the riparian condition indicators (R² = 0.99 for streambed width, R² = 0.82 for riparian zone width, R² = 0.89 for PPC, R² = 0.40 for bank stability). These research findings will be used in a 26,000 km mapping project assessing riparian vegetation and physical form indicators from LiDAR data in Victoria, Australia.     

An economical approach for d-lactic acid production utilizing unpolished rice from aging paddy as major nutrient source

In order to reduce the raw material cost of d-lactic acid fermentation, the unpolished rice from aging paddy was used as major nutrient source in this study. The unpolished rice saccharificate, wheat bran powder and yeast extract were employed as carbon source, nitrogen source and growth factors, respectively. Response surface methodology (RSM) was applied to optimize the dosages of medium compositions. As a result, when the fermentation was carried out under the optimal conditions for wheat bran powder (29.10 g/l) and yeast extract (2.50 g/l), the d-lactic acid yield reached 731.50 g/kg unpolished rice with a volumetric production rate of 1.50 g/(l h). In comparison with fresh corn and polished rice, the d-lactic acid yield increased by 5.79% and 8.71%, and the raw material cost decreased by 65% and 52%, respectively, when the unpolished rice was used as a major nutrient source. These results might provide a reference for the industrial production of d-lactic acid.     

Carbon and nitrogen footprint of double rice production in Southern China

Agriculture plays an important role in greenhouse gases (GHGs) emissions and reactive nitrogen (Nr) loss. Therefore, carbon (C) and nitrogen (N) footprint reductions in agro-ecosystem have become an increasingly hot topic in global climate change and agricultural adaptation. The objective of this study was to assess the C footprint (CF) and N footprint (NF) of double rice (Oryza sativa L.) production using life cycle assessment method in Southern China. The results showed that fertilizer application and farm machinery operation contributed the most to both GHGs and Nr emissions from agricultural inputs in the double rice production process. The CF for the early, late, and double rice was 0.86, 0.83, and 0.85 kg CO₂-eq kg⁻¹ year⁻¹ at yield-scale, respectively. In addition, the NF was 10.47, 10.89, and 10.68 g N-eq kg⁻¹ year⁻¹ at yield-scale for the early, late and double rice, respectively. The largest fraction of CF and NF of double rice was the share of CH₄ emission and NH₃ volatilization from the paddy field, respectively. Higher CF and NF at yield-scale for Guangdong, Guangxi, and Hainan provinces were presented, compared to the average level in double rice cropping for the region, while smaller than those of Jiangxi, Hubei, and Hunan provinces. Some effective solutions would be favorable toward mitigating climate change and eutrophication of the double rice cropping region in Southern China, including reduction of fertilizer application rates, improvements in farm machinery operation efficiencies, and changes in regional allocation of double rice cropping areas.     

Protective effects of pre-germinated brown rice diet on low levels of Pb-induced learning and memory deficits in developing rat

Lead (Pb) is a known neurotoxicant in humans and experimental animals. Numerous studies have provided evidence that humans, especially young children, and animals chronically intoxicated with low levels of Pb show learning and memory impairments. Unfortunately, Pb-poisoning cases continue to occur in many countries. Because the current treatment options are very limited, there is a need for alternative methods to attenuate Pb toxicity. In this study, the weaning (postnatal day 21, PND21) rats were randomly divided into five groups: the control group (AIN-93G diet, de-ionized water), the lead acetate (PbAC) group (AIN-93G diet, 2 g/L PbAC in de-ionized water), the lead acetate + WR group (white rice diet, 2 g/L PbAC in de-ionized water; PbAC + WR), the lead acetate + BR group (brown rice diet, 2 g/L PbAC in de-ionized water; PbAC + BR) and the lead acetate + PR group (pre-germinated brown rice diet, 2 g/L PbAC in de-ionized water; PbAC + PR). The animals received the different diets until PND60, and then the experiments were terminated. The protective effects of pre-germinated brown rice (PR) on Pb-induced learning and memory impairment in weaning rats were assessed by the Morris water maze and one-trial-learning passive avoidance test. The anti-oxidative effects of feeding a PR diet to Pb-exposed rats were evaluated. The levels of reactive oxygen species (ROS) were determined by flow cytometry. The levels of 8-hydroxy-2-deoxyguanosine (8-OHdG), γ-aminobutyric acid (GABA) and glutamate were determined by HPLC. Our data showed that feeding a PR diet decreased the accumulation of lead and decreased Pb-induced learning and memory deficits in developing rats. The mechanisms might be related to the anti-oxidative effects and large amount of GABA in PR. Our study provides a regimen to reduce Pb-induced toxicity, especially future learning and memory deficits in the developing brain.     

Electromyographic mapping of the erector spinae muscle with varying load and during sustained contraction

The purpose of the study was to explore changes in the spatial distribution of erector spinae electromyography amplitude during static, sustained contractions and during contractions of increasing load. Surface electromyographic (EMG) signals were detected from nine healthy subjects using a grid of 13 × 5 electrodes placed unilaterally over the lumbar erector spinae musculature. Subjects stood in a 20° forward flexed position and performed: (1) six 20-s long contractions with loads ranging from 2.5 kg to 12.5 kg (2.5 kg increments) and (2) a 6 min sustained contraction with 7.5 kg load. Root mean square (RMS) and mean power spectral frequency (MNF) were computed from the recorded EMG signals. EMG RMS increased (P < 0.0001) and MNF remained constant during contractions of increased load. During the sustained contraction, MNF decreased (P < 0.0001) and RMS did not change over time. The centroid (center of activity) of the RMS map shifted caudally (P < 0.0001) with time during the sustained contraction but did not change with varying load. These results suggest a change in the distribution of erector spinae muscle activity with fatigue and a uniform distribution of muscle activation across loads.     

Straighthead disease of rice (Oryza sativa L.) induced by arsenic toxicity

Straighthead disease is a physiological disorder of rice (Oryza sativa L.) characterized by sterility of the florates/spikelets leading to reduced grain yield. Though the exact cause of straighthead is unknown, a glass house experiment was conducted to investigate the effect of inorganic arsenic on straighthead disease in rice (Oryza sativa L.). BRRI dhan 29, a popular Bangladeshi rice strain, was grown in soils spiked with arsenic (prepared from sodium arsenate, Na₂HAsO₄·7H₂O) at the rate of 10, 20, 30, 40, 50, 60, 70, 80 and 90 mg of As kg^?1 and one control treatment was also run to compare the results. Although there may be some other soil physico-chemical factors involved, arsenic concentration was found to be closely associated with straighthead of rice. With the increase of soil arsenic concentration, the severity of straighthead increased significantly. Up to the 50 mg of As kg^?1 soil treatments, the severity of straighthead incidences were not prevalent. Straighthead resulted in sterile florets with distorted lemma and palea, reduced plant height, tillering, panicle length and grain yield. Straighthead caused approximately 17–100% sterile florates/spikelets formation and about 16–100% loss of grain yield. Straighthead also causes the reduction of panicle formation and panicle length significantly (p < 0.01). In the present study, panicle formation was found to be reduced by 21–95% by straighthead.     

Mapping soil organic matter in low-relief areas based on land surface diurnal temperature difference and a vegetation index

Accurate estimates of the spatial variability of soil organic matter (SOM) are necessary to properly evaluate soil fertility and soil carbon sequestration potential. In plains and gently undulating terrains, soil spatial variability is not closely related to relief, and thus digital soil mapping (DSM) methods based on soil–landscape relationships often fail in these areas. Therefore, different predictors are needed for DSM in the plains. Time-series remotely sensed data, including thermal imagery and vegetation indices provide possibilities for mapping SOM in such areas. Two low-relief agricultural areas (Peixian County, 28 km × 28 km and Jiangyan County, 38 km × 50 km) in northwest and middle Jiangsu Province, east China, were chosen as case study areas. Land surface diurnal temperature difference (DTD) extracted from moderate resolution imaging spectroradiometer (MODIS) land surface temperature (LST), and soil-adjusted vegetation index (SAVI) at the peak of growing season calculated from Landsat ETM+ image were used as predictors. Regression kriging (RK) with a mixed linear model fitted by residual maximum likelihood (REML) and residuals interpolated by simple kriging (SK) were used to model and map SOM spatial distribution; ordinary kriging (OK) was used as a baseline comparison. The root mean squared error, mean error and mean absolute error calculated from leave-one-out cross-validation were used to assess prediction accuracy. Results showed that the proposed covariates provided added value to the observations. SAVI aggregated to MODIS resolution was able to identify local highs and lows not apparent from the DTD imagery alone. Despite the apparent similarity of the two areas, the spatial structure of residuals from the linear mixed models were quite different; ranges on the order of 3 km in Jiangyan but 16 km in Peixian, and accuracy of best models differed by a factor of two (3.3 g/kg and 6.3 g/kg SOM, respectively). This suggests that time-series remotely sensed data can provide useful auxiliary variable for mapping SOM in low-relief agricultural areas, with three important cautions: (1) image dates must be carefully chosen; (2) vegetation indices should supplement diurnal temperature differences, (3) model structure must be calibrated for each area.     

Estimating biomass and carbon mitigation of temperate coniferous forests using spectral modeling and field inventory data

Realizing the importance of forest carbon monitoring and reporting in climate change, the present study was conducted to derive spectrally modeled aboveground biomass and mitigation using Landsat data in combination with sampled field inventory data in the coniferous forests of Western Himalaya. After conducting preliminary survey in 2009, 90 quadrats (45 each for calibration and validation) of 0.1 ha were laid in six forest types for recording field inventory data viz. diameter at breast height, height, slope and aspect. Biomass carbon (Mg ha^− 1) was worked out for different forest types and crown density classes (open with 10–40% crown density and closed with > 40% crown density) using recommended volume equations, ratios and factors. Biomass carbon map (aboveground + belowground) was generated for the entire region using geospatial techniques. Normalized difference vegetation index (NDVI) was generated and spectral values were extracted to establish relation (R² = 0.72, p < 0.01) with the field inventory data. The model developed was validated (R² = 0.73, p < 0.01) with 45 sample observations not used earlier for predicting and generating biomass carbon map (2009) for the entire region. The data from field based inventory indicates highest total biomass carbon (171.40, σ ± 23.19) Mg ha^− 1 for Fir–Spruce (closed) which has relatively more mature girth classes and low tree density. This value was found to be significantly higher than other forest types. Lowest biomass carbon was observed for Blue Pine (open) (37.15, σ ± 11.82) Mg ha^− 1. The NDVI values for the entire region ranged from 0 to 0.62 and consequently the spectrally derived aboveground biomass carbon varied from 0 to 600 Mg ha^− 1. The study demonstrates the application of mapping, spectral responses and sampled field inventory for type wise assessment of carbon mitigation in temperate coniferous forests of Himalayas.     

Cenozoic long-term terrestrial climatic evolution in Germany tracked by δ18O of rodent tooth phosphate

Even though the Cenozoic has been recognized as a period of important climate change, long-term climatic changes that took place in the continental domain are still questioned. For an area, southern Germany, for which other long-term palaeoclimatic records exist, analysis of oxygen isotope composition of small mammal teeth has been carried out for localities ranging in age from the late Eocene (c.36 Ma) to the middle Miocene (c.10 Ma). Comparison of this long-term continental δ¹⁸O record with the marine record reveals comparable trends. The major Cenozoic climatic events are thus recognized in the continental oxygen isotope record.Through comparison with other quantitative palaeoclimatic reconstructions available for the studied area, temperatures in southern Germany are broadly estimated to have fluctuated between 12 °C and 25 °C from late Eocene to early late Miocene. According to the different available chronological frameworks, either a ～ 2 °C (considering the classical biostratigraphy) or a ～ 6 °C (considering a revised biostratigraphy) decrease in mean air temperature is estimated for the Eocene/Oligocene boundary.     

Characterization of spatial variability of soil physicochemical properties and its impact on Rhodes grass productivity

Characterization of soil properties is a key step in understanding the source of spatial variability in the productivity across agricultural fields. A study on a 16 ha field located in the eastern region of Saudi Arabia was undertaken to investigate the spatial variability of selected soil properties, such as soil compaction ‘SC’, electrical conductivity ‘EC’, pH (acidity or alkalinity of soil) and soil texture and its impact on the productivity of Rhodes grass (Chloris gayana L.). The productivity of Rhodes grass was investigated using the Cumulative Normalized Difference Vegetation Index (CNDVI), which was determined from Landsat-8 (OLI) images. The statistical analysis showed high spatial variability across the experimental field based on SC, clay and silt; indicated by values of the coefficient of variation (CV) of 22.08%, 21.89% and 21.02%, respectively. However, low to very low variability was observed for soil EC, sand and pH; with CV values of 13.94%, 7.20% and 0.53%, respectively. Results of the CNDVI of two successive harvests showed a relatively similar trend of Rhodes grass productivity across the experimental area (r = 0.74, p = 0.0001). Soil physicochemical layers of a considerable spatial variability (SC, clay, silt and EC) were utilized to delineate the experimental field into three management zones (MZ-1, MZ-2 and MZ-3); which covered 30.23%, 33.85% and 35.92% of the total area, respectively. The results of CNDVI indicated that the MZ-1 was the most productive zone, as its major areas of 50.28% and 45.09% were occupied by the highest CNDVI classes of 0.97–1.08 and 4.26–4.72, for the first and second harvests, respectively.