<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ARS</journal-id><journal-title-group><journal-title>Advances in Remote Sensing</journal-title></journal-title-group><issn pub-type="epub">2169-267X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ars.2018.71001</article-id><article-id pub-id-type="publisher-id">ARS-82736</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Temperature Based Spatiotemporal Growth Monitoring of Rice Plant from Germination-Ripening Stage Using Remote Sensing and GIS Techniques
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>S.</surname><given-names>M. H. Raza</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>S.</surname><given-names>A. Mahmood</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>H.</surname><given-names>Batool</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>T.</surname><given-names>J. Shad</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>S.</surname><given-names>Alvi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>F.</surname><given-names>Waseems</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>M.</surname><given-names>Atif Butt</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>S.</surname><given-names>Shehzad Hassan</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ali</surname><given-names>Imam Mirza</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Geography, University of the Punjab, Quaid e Azam Campus, Lahore, Pakistan</addr-line></aff><aff id="aff1"><addr-line>Remote Sensing and GIS Group, Department of Space Science, University of the Punjab, Quaid-e-Azam Campus, Lahore, 
Pakistan</addr-line></aff><pub-date pub-type="epub"><day>28</day><month>02</month><year>2018</year></pub-date><volume>07</volume><issue>01</issue><fpage>1</fpage><lpage>14</lpage><history><date date-type="received"><day>26,</day>	<month>October</month>	<year>2017</year></date><date date-type="rev-recd"><day>24,</day>	<month>February</month>	<year>2018</year>	</date><date date-type="accepted"><day>28,</day>	<month>February</month>	<year>2018</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Rice (
  Oryza sativa) is the second staple food largely grown and widely consumed in Pakistan. About 10% of the total crop area of Pakistan is cultivated by rice that takes a part in value addition of almost 1.3% - 1.6% in the total Gross Domestic Product (GDP). Due to global warming, temperature has a profound impact on rice crop phenology. Low temperature is the main factor of delay in rice plant growth and very high temperature results in stressed and short heighted plant so the crop sown in a region at the same time is not ready to harvest at same hours but a delay is observed. The study area under investigation was district Sheikhupura, Nankana, Lahore, Gujranawala and Hafizabad, which are famous for rice productivity. Landsat 7, 8 freely available thermal dataset are used to calculated pixel based temperature values to evaluate growth using agricultural growth indicators. The total covered area was 13,480 km
  <sup>2</sup> in which 484 km
  <sup>2</sup> area was marked as less growth rate area with low temperature values due to water body and excess of vegetation over there. About 7960 km
  <sup>2</sup> area is marked as good for growth experiencing optimum temperature for rice plant. Approximately 4944 km2 area is marked as stressed rice plant area experiencing high temperature values adjacent to urban population. An attempt is made here to map this effect of temperature-based growth variability of the rice plant across the study area.
 
</p></abstract><kwd-group><kwd>Climate</kwd><kwd> Gross Domestic Product (GDP)</kwd><kwd> Global Warming</kwd><kwd> Phenology</kwd><kwd>  Temperature</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>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% [<xref ref-type="bibr" rid="scirp.82736-ref1">1</xref>] .</p><p>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.</p><p>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 [<xref ref-type="bibr" rid="scirp.82736-ref2">2</xref>] .</p><p>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.</p></sec><sec id="s2"><title>2. Area under Investigation</title><p>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 <xref ref-type="fig" rid="fig1">Figure 1</xref>. 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.</p><sec id="s2_1"><title>2.1. Rice Sequential Stages</title><p>Rice plant undergoes three sequential growth stages that include, vegetative</p><p>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.</p><sec id="s2_1_1"><title>2.1.1. Germination Stage</title><p>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 [<xref ref-type="bibr" rid="scirp.82736-ref3">3</xref>] .</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref> 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.</p><p>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 [<xref ref-type="bibr" rid="scirp.82736-ref4">4</xref>] .</p></sec><sec id="s2_1_2"><title>2.1.2. Seed Establishment Stage</title><p>In the post germination stage the growth is very sensitive to the temperature. <xref ref-type="fig" rid="fig3">Figure 3</xref> 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.</p><p>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 [<xref ref-type="bibr" rid="scirp.82736-ref5">5</xref>] .</p></sec></sec></sec><sec id="s3"><title>3. Material and Methods</title><sec id="s3_1"><title>3.1. Length of Growing Season</title><p>Landsat 8 satellite image of 19-April 2014 is selected for germination of rice plant by verifying the date using [<xref ref-type="bibr" rid="scirp.82736-ref6">6</xref>] 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</p><p>season. This algorithm is formulated as following.</p><p>NDVI Ratio = NDVI − NDVI min NDVI max − NDVI min (1)</p><p>In Equation (1) NDVI is interpolated daily NDVI value throughout the year 2014. NDVI<sub>max</sub> and NDVI<sub>min</sub> are the maximum and minimum values recorded in a year. Start of season is defined as the day of year when NDVI<sub>Ratio</sub> 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.</p></sec><sec id="s3_2"><title>3.2. Procedure of Temperature Calculation Using Landsat 8 Thermal Bands</title><p>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.</p></sec><sec id="s3_3"><title>3.3. Calculation of Irradiance</title><p>Initially the DN values are converted into radiance. Radiance gives the flux of radiant energy per unit area.</p><p>Radiance is calculated here using the following formula.</p><p>Irradiance = ( 3.342 &#215; 10 − 4 &#215; Thermal   Band ) + 0.1 [<xref ref-type="bibr" rid="scirp.82736-ref7">7</xref>] (2)</p><p>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/(m<sup>2</sup>&#215;srad&#215;μm). <xref ref-type="fig" rid="fig4">Figure 4</xref> determines the irradiance based variations in the study site, where Lahore shows the highest irradiance values due to highest soil fulx.</p></sec></sec><sec id="s4"><title>4. Results and Discussion</title><sec id="s4_1"><title>4.1. Temperature Calculation</title><p>For temperature calculation, the dataset of irradiance calculated earlier is used</p><p>and the following algorithm applied</p><p>T = { K 2 / ln ( ε K 1 / Irradiance + 1 ) } − 272.15 (3)</p><p>In Equation (3), K<sub>1</sub> and K<sub>2</sub> are the constants for thermal bands and their values are mentioned in metadata Landsat, in our case we used following values of these constants.</p><p>K1_CONSTANT_BAND_10 = 774.89</p><p>K1_CONSTANT_BAND_11 = 480.89</p><p>K2_CONSTANT_BAND_10 = 1321.08</p><p>K2_CONSTANT_BAND_11 = 1201.14</p><p>ε 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 <xref ref-type="fig" rid="fig5">Figure 5</xref>.</p></sec><sec id="s4_2"><title>4.2. Histogram Evaluation</title><p>The histogram in <xref ref-type="fig" rid="fig6">Figure 6</xref> is generated using temperature based dataset calculated</p><p>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.</p><p>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 <xref ref-type="fig" rid="fig5">Figure 5</xref>. 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.</p><p><xref ref-type="fig" rid="fig7">Figure 7</xref> describes that absorption of moisture content of rice seed increases with increase in temperature.</p><p><xref ref-type="fig" rid="fig8">Figure 8</xref> describes that germination of rice plant is achieved in less time at higher temperatures.</p><p><xref ref-type="fig" rid="fig9">Figure 9</xref> describes that maximum cell division occur at optimum temperature while high and low temperature values cause stressed or slow cell division.</p><p>Lush green areas in Figures 10-13 are near to the populated regions of the</p><p>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.</p><p>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.</p></sec><sec id="s4_3"><title>4.3. Leaf Emergence</title><p>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.</p></sec><sec id="s4_4"><title>4.4. Tillering Stage</title><p>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</p></sec><sec id="s4_5"><title>4.5. Heading Stage</title><p>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</p></sec><sec id="s4_6"><title>4.6. Ripening Period</title><p>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.</p><p>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.</p></sec></sec><sec id="s5"><title>5. Conclusions</title><p>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.</p><p>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 SO<sub>x</sub>, NO<sub>x</sub> 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.</p></sec><sec id="s6"><title>Acknowledgements</title><p>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.</p></sec><sec id="s7"><title>Cite this paper</title><p>Raza, S.M.H., Mahmood, S.A., Batool, H., Shad, T.J., Alvi, S., Waseems, F., Butt, M.A., Hassan, S.S. and Mirza, A.I. (2018) Temperature Based Spatiotemporal Growth Monitoring of Rice Plant from Germination-Ripening Stage Using Remote Sensing and GIS Techniques. Advances in Remote Sensing, 7, 1-14. https://doi.org/10.4236/ars.2018.71001</p></sec></body><back><ref-list><title>References</title><ref id="scirp.82736-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Qureshi, K.A., et al. (1993) Contribution of Area and Yield of Total Rice Production in Pakistan—An Analysis. Pakistan Journal of Agriculture Science, 30.</mixed-citation></ref><ref id="scirp.82736-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Ziska, E., Fraser, D., et al. (1997) Assessing Risks of Climate Variability and Climate Change for Rice. Science, 240, 996-1002.</mixed-citation></ref><ref id="scirp.82736-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Livingston, B.E. and Haasis, F.W. (1932) Relations of Time and Maintained Temperature to Germination Percentage for a Lot of Rice Seed. American Journal of Botany, 596-615.</mixed-citation></ref><ref id="scirp.82736-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Yoshida, S. (2012) Effects of Temperature on Growth of the Rice Plant (Oryza sativa L.) in a Controlled Environment. Soil Science and Plant Nutrition, 302-304.</mixed-citation></ref><ref id="scirp.82736-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Nakagava (2003) Response of Rice to Increasing Temperature and Atmospheric CO2. Environmental Sciences and Engineering, 63-80.</mixed-citation></ref><ref id="scirp.82736-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">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.</mixed-citation></ref><ref id="scirp.82736-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">http://landsat.usgs.gov/Landsat8_Using_Product.php</mixed-citation></ref></ref-list></back></article>