Pakistan is fortunate to enjoy the rich export quality of Basmati rice that returns three to four times of all the coarse varieties available in international market. The cultivated land area for rice production doubled from 1948 to 1978 due to high demands at national level. The ratio of increment in the area used for rice production increases by roughly 2% annually. Rice cultivation area increased from 790 thousands ha in 1947 to 2017 thousands ha in 1990 with 163.54% increment over the period of 44 years with same trend on the way. Production also increased along the same patterns with 693 thousand tons in 1947 and 3220 thousand tons in 1990 giving an increase of 364%. In context of yield, it went on increasing from 877 kg in 1948 to 1528 kg in 1990 indicating an increase of 73%  .
Punjab Pakistan, is famous for best quality of rice throughout the world. Punjab lies in the temperate and arid zone with extreme variations in weather. It is generally characterized by long hot summer season and cold winters with little rainfall, not more than 255 mm annually. Summer season is getting longer due to anthropogenic activities which may have a negative impact on rice productivity. This abnormal change in temperature is responsible for dramatic shift in rice productivity regardless of rice crop calendar of Punjab province.
Agricultural productivity is largely influenced by the weather and climate in any region. Productivity is regulated by prevailing climatic factors including light intensity, sunshine duration, rainfall, temperature, air pressure, and humidity level across the region. An increase in atmospheric temperature has considerable effect on growth and quality of rice by affecting its phonology, physiology, and yield components. It is also observed that an increase in temperature by 4˚C during the growth period of rice plant results in early maturation by five to six days in both wet and dry seasons  .
Rice plants require a specific amount of heat to develop from seedling to germination and germination to maturation. Crop calendars are used to predict the development stage of rice plant to take management decisions which can often mislead due to unusual climatic variations.
2. Area under Investigation
Five districts, including Sheikhupura, Hafizabad, Gujranwala, Nankana Sahib and Lahore are considered, as largest producers of rice comprise of an area known as Eastern Punjab, mapped in Figure 1. Eastern Punjab is the plain area of Pakistan which is connected to India on one side with an average elevation of about 712ft sea level and adjoining boundary with Sargodha and Faisalabad on the western side. Geospatially it is located between 73E to 75E longitude and 31N to 32.5N latitude with high profile fertile land. Our study area is a subset of Landsat patch with 149 path No and 38 row No.
2.1. Rice Sequential Stages
Rice plant undergoes three sequential growth stages that include, vegetative
Figure 1. Eastern Punjab Pakistan.
stage, reproductive stage and maturity stage. Vegetative stage takes about 50 to 60 days, reproductive stage takes about 20 to 30 days and the final maturity stage is completed in 30 to 45 days. These stages are further divided in Germination, seed establishment, leaf emergence, tillering, heading, and ripening stage. Duration for each stage is dependent on water availability and environmental conditions in which it is grown.
2.1.1. Germination Stage
Germination initializes as white tip appears from surface of land or water, when seed dormancy is broken. A single grain of rice weighs about 12 - 44 mg depending upon the variety at 0% moisture content. Hull weighs about 21% of total grain weight. Germination starts when seed dormancy is broken after absorption of adequate water at temperature ranging between 10˚C - 40˚C  .
Figure 2 describes that water absorption of seed is achieved in less time at higher temperature and vice versa. 60% moisture content is absorbed by the seed in 54 hours at 30˚C and the same rate of absorption is obtained in 75 - 80 hours at 25˚C and 110 hours at 20˚C and so on.
Figure 2. Rate of moisture content absorption at different temperature values.
Temperature has a profound impact on germination. Germination precentage varies at different temperature with a completion period of 2, 3, and 14 days after incubation. 90% - 97% germination is achieved in 2 days at a temperature range of 27˚C - 37˚C, in 3 days at a temperature between 25˚C - 28˚C, in 6 days between temperature range of 18˚C - 24˚C and in 14 days at a temperature range between 14˚C - 17˚C  .
2.1.2. Seed Establishment Stage
In the post germination stage the growth is very sensitive to the temperature. Figure 3 describes that maximum cell division occurs at 25˚C - 30˚C so this temperature range is considered good for high growth rates. Above and below this range, the seed growth rate becomes slow and the seed may die at 40˚C.
In the second week of germination, rice plant enters into the phase of tissue elongation which results in cell division and cell enlargement. Critical temperature for this growth stage is 30˚C which slows down while moving toward 15˚C and it stops below 15˚C  .
3. Material and Methods
3.1. Length of Growing Season
Landsat 8 satellite image of 19-April 2014 is selected for germination of rice plant by verifying the date using  midpoint pixel algorithm. Initially time series NDVI data is produced for the area using MODIS 13A2 algorithm which is meaningful for phonological monitoring. White et al., pixel formula is then applied to identify the date of start of season for rice plantation. This algorithm is sophisticated and widely used to identify the time span of length of growing
Figure 3. Rate cell division at different temperature values.
season. This algorithm is formulated as following.
In Equation (1) NDVI is interpolated daily NDVI value throughout the year 2014. NDVImax and NDVImin are the maximum and minimum values recorded in a year. Start of season is defined as the day of year when NDVIRatio exceeds 0.5 as this is considered to correspond to the timing of most rapid increase in NDVI to the initial leafing of our stormy canopy.
3.2. Procedure of Temperature Calculation Using Landsat 8 Thermal Bands
Landsat 8 datasets used in this research are freely available on earth explorer website. These datasets consist of two thermal bands, B10 and B11 with 100 m spatial resolution in atmospheric window between 10 μm and 12 μm. These datasets are used here to calculate temperature against each pixel on the ground.
3.3. Calculation of Irradiance
Initially the DN values are converted into radiance. Radiance gives the flux of radiant energy per unit area.
Radiance is calculated here using the following formula.
The value of 3.3420E−04 in Equation (2) is extracted from metadata of Landsat 8 Image. This formula helps to calculate Irradiance against each pixel across the image in Watts/(m2×srad×μm). Figure 4 determines the irradiance based variations in the study site, where Lahore shows the highest irradiance values due to highest soil fulx.
4. Results and Discussion
4.1. Temperature Calculation
For temperature calculation, the dataset of irradiance calculated earlier is used
Figure 4. Irradiance map of study area.
and the following algorithm applied
In Equation (3), K1 and K2 are the constants for thermal bands and their values are mentioned in metadata Landsat, in our case we used following values of these constants.
K1_CONSTANT_BAND_10 = 774.89
K1_CONSTANT_BAND_11 = 480.89
K2_CONSTANT_BAND_10 = 1321.08
K2_CONSTANT_BAND_11 = 1201.14
ε is emissivity which is taken as 0.95 generally. However many authors have calculated this factor by using indices such as Normalized Difference Vegetation Index etc. to examine soil emissivity level. The spatial distribution of the temperature values in the study site are mapped in Figure 5.
4.2. Histogram Evaluation
The histogram in Figure 6 is generated using temperature based dataset calculated
Figure 5. Temperature map of study area.
Figure 6. Histogram of temperature and pixel density.
above. In this histogram, temperature values are taken on x-axis and pixel density on y-axis. It clearly describes the actual temperature variation is between 27˚C - 31˚C and maximum number of temperature values are falling in this range. The values other than this range are due to the flux reflected by some manmade features like (Plastic sheds, solar panel, steel pipes etc..) which should be ignored for appropriate results.
Rice seed germination initializes between 15 April to 20 May in our study area as these dates are favorable for rice seed germination. We calculated the date of start of season in Figure 5. and obtained Landsat 8 satellite images for the same date of start of season when rice plant is in dormancy broken and germination phase to investigate the spatial growth rate variations.
Figure 7 describes that absorption of moisture content of rice seed increases with increase in temperature.
Figure 8 describes that germination of rice plant is achieved in less time at higher temperatures.
Figure 9 describes that maximum cell division occur at optimum temperature while high and low temperature values cause stressed or slow cell division.
Lush green areas in Figures 10-13 are near to the populated regions of the
Figure 7. Water absorption capacity map of Study area using Landsat 8 image of date 19-04-2014.
Figure 8. Germination map of study area using Landast8 satellite image of 05-May-2014.
Figure 9. Rice seed cell division map of study area using Landsat 8 image of date 05-May-2014.
Figure 10. Leaf emergence map of rice plant in the study area using Landsat Image of date 21-May-2014.
Figure 11. Tillering map of rice plant in the study area using Landsat Image of the date 22-June-2014.
Figure 12. Heading/Anthesis map of rice plant in the study area using Landsat Image of the date 25-August-2014.
Figure 13. Repining map of rice crop in the study area using Landsat image of the date 26-Sep-2014.
city with higher temperature values that produce short heighted, stressed plants with less yield. The stressed plants are the result of low metabolism with considerable damage to the plant, sometimes leading to cell death. High concentration of environmental pollution in metropolitan areas, cause an artificial green house impact on rice crop which results in an increase of temperature and finally the stressed vegetation.
Rice crop cultivated near to water body is exposed to comparatively low temperatures that are less than optimum range which results in slow growth for each stage of rice plant development and a sharp decrement in rice plant development is observed. Finally the crop is delayed and results a lengthy growing season in the regions adjacent to water body.
4.3. Leaf Emergence
A rice leaf is composed of a sheet, ligules, blade and auricle. The sheet ribbon shaped structure is rolled into a cylinder that includes all the parts of a shoot. The blade is a longer, flat but narrow part of a leaf. Temperature has profound impact on leaf emergence stage. High temperature speeds up leaf emergence process. Leaf emerges every fifth day for the rice plant facing 20˚C. 100 degree days are recommended for first leaf emergence so 5 days are required (100/20 = 5 days) at 20˚C and every fourth day at 25˚C which is (100/25 = 4 days). Succeeding leaf emergence happen after the former leaf has completed its emergence. Upper leaves have longer life span than the lower ones.
4.4. Tillering Stage
Tiller are the branches that emerge from the axil during vegetative stage. The first leaf of tiller emerges from second leaf of main Culm when the fifth leaf on the main Culm emerges. Similarly the second leaf of tiller develops. This process of tiller emergence continues till the emergence of last leaf. This rule is not only applicable on the main Culm but to all tillers
4.5. Heading Stage
The time of initiation of Panicle from the boot is known as heading. This process takes about 10 - 14 days for 100 percent development depending upon environmental conditions in which it is exposed. Temperature conditions play a vital role in this stage, for example, a sharp increase in number of days is observed as temperature drops from 24˚C to 21˚C. A temperature drop of 1˚C causes a delay of 13 days in heading. When temperature is increased, number of days are decreased from 91 at 27˚C to 86 days at 30˚C. Flowering starts as the process of heading ends that continues for up to 7 days. The floral organ of rice including Lemma and Palea is called spikelet. Anthesis is a series of event of opening and closing spikelet which lasts within 1 - 2.5 hours
4.6. Ripening Period
Ripening period is defined by increase in size and weight and change in color of grain from green to yellow periodically. In the early stage, the grains of rice are green that turn to yellow and the texture of grain turns from milky to dough and yellow ripe to maturity. Considerable amount of sugar and starch in the plant leaf sheet and Culm is translocated to the grain during ripening period. Duration of grain filling up to maturity differs with different rice varieties. For example Indica rice variety IR20 takes 13 days at mean temperate of 28˚C and 33 days at 16˚C to reach ripening period. A japonica rice variety Fujiska 5 takes 18 days at a mean temperature of 28˚C and 43 days at 16˚C because water content decreases faster at high temperatures, however above a certain limit temperatures can cause stress in growth which results in low yield. 15˚C - 21˚C is considered optimum temperature in repining stage while 22˚C - 29˚C is critical high. High temperature during grain filling results in increased respiration and more carbohydrates consumption. This result in less filled spikelet and reduced yield, therefore night time is considered good during the grain filling stage.
The objective of this paper is to incorporate the Landsat 8, Thermal band data to calculate temperature values against each pixel in the study area. Spatial variation of these heat units will help to mark the regions with slow or fast growth rates of rice plant in various growth stages e.g. (Germination, Tillering, Anthesis and ripening etc.) described briefly above.
Remote sensing and GIS techniques applied in this research, not only proved to be very useful in spatio-temporal monitoring of rice plant at various growth stages, but also remained time efficient for a large study area. Remote Sensing and GIS techniques are capable of extraction of error free data to develop temperature based thematic maps instead of manual techniques of data collection of wide areas which are always time consuming and more expansive.
Rice plantation near to cities and closer to the road network experienced high temperature which resulted in stressed and short heighted plants due to high concentration of environmental pollution including SOx, NOx and particulate matter. Rice plantation near to water body e.g. Rivers, Canals and Nalas, experienced comparatively low temperature values which resulted in slow growth rate. All the histograms developed in this study show that highest density of temperature values are falling in optimum temperature for particular growth stage so our study area is good for rice plantation. Rice plantation in our study are sown at same time is not ready to harvest at the same time but there is a delay due to heat units experienced by a single rice plant. We could not get thermal bands with 1day temporal resolution. However, if such high temporal resolution data are obtained, best results regarding growth monitoring against each growth stage may be achieved and finally a best yield.
The authors are thankful to the United States Geological Survey (USGS) for provision of a value able data that we used in this research. We are also thankful to worthy Vice Chancellor of the University of the Punjab Pakistan for provding us a good environment to complete this research.
 White, M.A., de Beurs, K., Didan, K., Inouye, D., Richardson, A., Jensen, O., et al. (2009) Inter Comparison, Interpretation, Assessment of Spring Phonology in North America Estimated from Remote Sensing for 1982-2006. Global change Biology, 15, 2335-2359.