InstructionsThe following steps are necessary for the selection of the “hot” pixel:1. The selection of the “hot” pixel follows the same proc<strong>edu</strong>re as for the “cold” pixel. Itis more difficult, however, because there is a <strong>br</strong>oader range of temperatures for “hot”pixel candidates to select from.2. The “hot” pixel should be located in a dry and bare agricultural field where one canassume there is no evapotranspiration taking place. It is recommended that one notuse a hot desert area, an asphalt parking lot, a roof, a south facing mountain slope,or other such extremely hot areas, since the relationship between dT and T s may notfollow the same linear relation as for agricultural and bare agricultural soils. Inaddition, the prediction of soil heat flux, G, is less certain for man-made structuresand desert. Desert is difficult because the thin canopy creates a two-source systemfor heat transfer (canopy and underlying soil).3. The “hot” pixel should have a surface albedo similar to other dry and bare fields inthe area of interest. It should have a LAI in the range of 0 to 0.4 (corresponding tolittle or no vegetation). In the Mountain Model, the “hot” pixel should have a slope ofless than one percent. As with the “cold” pixel, it should be located near the weatherstation.4. If there has been precipitation within the past three or four days prior to the image,then it is possible that the hot pixel may exhibit some ET. In this case, it cannot beassumed that the H = R n – G at the hot pixel. Instead, H must be predicted as R n –G – ET bare soil , where ET bare soil is predicted by operating a daily soil water balancemodel of the surface soil using ground-based weather measurements. An examplemodel is included in Annex 8 of FAO-56 available at:http://www.fao.org/docrep/X0490E/X0490E00.htmA downloadable spreadsheet from Annex 8 is available from:http://www.kimberly.uidaho.<strong>edu</strong>/water/fao56/index.htmlGenerally, the ET from a bare soil can be roughly estimated to be equal to 0.8 ET r ,0.5 ET r , 0.3 ET r , 0.2 ET r , and 0.1 ET r , 1, 2, 3, 4, and 5 days following a substantialrain event (15 mm or more).Illustrative Example (Landsat 7 ETM+, 4/8/2000)In the following image, a simple daily water balance model predicted that bare soil in thisarea is totally dry so that we can assume zero ET for the “hot” pixel.72
We open the image or subset image and select true and false color. We open the T s imageusing Raster Options-Pseudo Color and link the two viewers. We go to Raster – Attributesin the T s viewer to view the attribute table shown below:The hottest pixel in the image is around 316 K. However, the purpose is not to select anextreme “hot” pixel but to select a “dry” pixel where ET is zero.The problem is that the surface temperatures of “totally dry surfaces” are not necessarilyuniform. For example, T s of a dry black asphalt road might be extremely high while T s of adry desert might be lower but still quite high. T s of a very dry agricultural field might bemoderately high. ET from these three different surfaces could all be zero or very close tozero.The selection of an extremely high temperature pixel for the “hot” pixel could cause <strong>SEBAL</strong>to predict positive ET for other “dry “ pixels. For example, if the “cold” pixel is in a fieldhaving T s = 18 o C and our “hot” pixel is selected in the middle of a desert with T s = 43 o C,<strong>SEBAL</strong> would regard any totally dry and fallow agricultural field having T s = 30 o C as a fieldwith some moisture. This is because the dT for this dry agricultural field would be predictedto be less than that for the hot pixel based on our dT vs. T s function, and therefore H wouldbe under predicted. Therefore, it is best to limit the land surface condition for candidate hotpixels to agricultural fields.Our first step in selecting a “hot” pixel is to find a representative, totally dry agricultural fieldin the area of interest. It is important to check precipitation records in order to get an idea ofhow wet an area is. One can apply a water balance model to predict the residual moisturein non-irrigated bare soil. For our example, many fallow fields are evident in the image and73
- Page 1 and 2:
SEBALSurface Energy Balance Algorit
- Page 3 and 4:
6. Tables for surface albedo comput
- Page 5 and 6:
δ declination of the earth radians
- Page 7 and 8:
INTRODUCTION1. BACKGROUNDLand manag
- Page 9 and 10:
THE THEORETICAL BASIS OF SEBAL1. OV
- Page 11 and 12:
The surface emissivity is the ratio
- Page 13 and 14:
• Reference evapotranspiration, E
- Page 15 and 16:
Figure 3. Flow Chart of the Net Sur
- Page 17 and 18:
2. The reflectivity for each band (
- Page 19 and 20:
and ET are specified for the two ex
- Page 21 and 22: To complete step 2, model F06_surfa
- Page 23 and 24: D. Choosing the “Hot” and “Co
- Page 25 and 26: irrigated crops near Kimberly, Idah
- Page 27 and 28: where; u * is the friction velocity
- Page 29 and 30: Weather stationu, z x, z om, u *H f
- Page 31 and 32: Note: An alternative method used by
- Page 33 and 34: where;⎛ 2 ⎞⎜1+ x(0.1m)Ψ = 2
- Page 35 and 36: where; ET inst is the instantaneous
- Page 37 and 38: Figure 6 shows an example of a dail
- Page 39 and 40: 1. If there is some cloud cover in
- Page 41 and 42: • For all other formats (includin
- Page 43 and 44: On the main page, select Type in Pa
- Page 45 and 46: These values for z om can vary by r
- Page 47 and 48: Appendix 3Reference Evapotranspirat
- Page 49 and 50: and click on “Open”3. Provide t
- Page 51 and 52: 5. Describe weather station and wea
- Page 53 and 54: 8. Save the output file using a nam
- Page 55 and 56: 4. Go to the Viewer tool and view t
- Page 57 and 58: Appendix 5Time Issues Around Weathe
- Page 59 and 60: 2. What time interval does the data
- Page 61 and 62: The wind speed is 3.4 m/s for the 1
- Page 63 and 64: Table 6.3. ESUN λ for Landsat 5 TM
- Page 65 and 66: Table 6.6. Typical albedo values.Fr
- Page 67 and 68: 4. Make a first guess for T cold by
- Page 69 and 70: Fig.3. P40/R30, Landsat 7 ETM+ 4/8/
- Page 71: Fig.5. P40/R30, Landsat 7 ETM+ 4/8/
- Page 75 and 76: Now, we look at the northwest area
- Page 77 and 78: (2) The regional weather or soil co
- Page 79 and 80: We conclude that this “heat islan
- Page 81 and 82: Part 1: using the manual iteration
- Page 83 and 84: Part 2: Using the automatic iterati
- Page 85 and 86: Appendix 9SAVI and LAI for Southern
- Page 87 and 88: equation for G. The following are e
- Page 89 and 90: and for 24-hours:Gwater = 0.9 Rn -
- Page 91 and 92: The Stability of the AtmosphereAppe
- Page 93 and 94: Appendix 12The SEBAL Mountain Model
- Page 95 and 96: where; t is the standard clock time
- Page 97: N. Momentum Roughness LengthThe mom