Monitoring the dynamic thermal activity in Yellowstone National Park is required by the United States Congress. The continuous monitoring is important to maintain the safety of the visitors and park service personnel, plan and relocate infrastructure, and study potential impact from nearby geothermal development including oil and gas industry. This dissertation is part of a study initiated in the early 2000s to monitor the thermal activity of dynamic areas within the Park, using airborne remote sensing imagery. This study was focused in Norris Geyser Basin, the hottest geyser basin in the park, located near the northwestern rim of the Yellowstone’s caldera. The study is considered the first long-term comprehensive airborne remote sensing study in the basin which took place between August 2008 and October 2013. In this study, at least one 1-meter resolution thermal infrared image and three-band images (multispectral) were acquired and used to estimate year-to-year changes in radiant temperature, radiant flux, and radiant power from the thermal source in Norris. Presence of residual radiant flux in the ground from absorbed solar radiation and atmospheric longwave radiation was the main challenge to compere year-to-year changes in the thermal activity. This residual flux is included in the total radiant flux calculated through the remote sensing images which gives false estimates of the flux generated from the underling thermal source. Two methods were suggested in Chapters 2 and 4 of this dissertation to estimate the residual radiant flux. A method was developed in Chapter 2 to estimate the residual radiant flux in a bare ground area covered with hydrothermal siliceous sinter deposit. The method compared ground-based measurements with high spatial resolution airborne remote sensing measurements to estimate the residual radiant flux. In Chapter 4, a method was developed to estimate the residual radiant flux in the six surface classes in Norris, including bare ground, bare ground with siliceous sinter deposit, lakes and pools, river, forest, and grass. The assumptions and implications of each method were discussed to suggest a reliable method to estimate the geothermal radiant flux after subtracting the absorbed residual radiant flux. Chapter 3 provides an analysis of the four components of heat flux in the ground surface, including conduction of sensible heat, convection of sensible heat by liquid water and water vapor, and convection of latent heat by water vapor. The main purpose from the analysis was to assess the hypothesis that the convection and latent heat flux are negligible which therefore supported the results obtained from the analysis in Chapters 2 and 4.

Capital investment reduction, exergetic efficiency improvement and material compatibility issues have been identified as the primary techno-economic challenges associated, with the near-term development and deployment of thermal energy storage (TES) in commercial-scale concentrating solar power plants. Three TES techniques have gained attention in the solar energy research community as possible candidates to reduce the cost of solar-generated electricity, namely (1) sensible heat storage, (2) latent heat (tank filled with phase change materials (PCMs) or encapsulated PCMs packed in a vessel) and (3) thermochemical storage. Among these the PCM macro-encapsulation approach seems to be one of the most-promising methods because of its potential to develop more effective energy exchange, reduce the cost associated with the tank and increase the exergetic efficiency. However, the technological barriers to this approach arise from the encapsulation techniques used to create a durable capsule, as well as an assessment of the fundamental thermal energy transport mechanisms during the phase change. A comprehensive study of the energy exchange interactions and induced fluid flow during melting and solidification of a confined storage medium is reported in this investigation from a theoretical perspective. Emphasis has been placed on the thermal characterization of a single constituent storage module rather than an entire storage system, in order to, precisely capture the energy exchange contributions of all the fundamental heat transfer mechanisms during the phase change processes. Two-dimensional, axisymmetric, transient equations for mass, momentum and energy conservation have been solved numerically by the finite volume scheme. Initially, the interaction between conduction and natural convection energy transport modes, in the absence of thermal radiation, is investigated for solar power applications at temperatures (300 - 400°). Later, participating thermal radiation within the storage medium has been included in order to extend the conventional natural convection-dominated model and to analyze its influence on the melting and freezing dynamics at elevated temperatures (800 - 850°). A parametric analysis has been performed in order to ascertain the effects of the controlling parameters on the melting/freezing rates and the total and radiative heat transfer rates at the inner surface of the shell. The results show that the presence of thermal radiation enhances the melting and solidification processes. Finally, a simplified model of the packed bed heat exchanger with multiple spherical capsules filled with the storage medium and positioned in a vertical array inside a cylindrical container is analyzed and numerically solved. The influence of the inlet mass flow rate, inner shell surface emissivity and PCM attenuation coefficient on the melting dynamics of the PCM has been analyzed and quantified.

The use of a latent heat storage system using phase change materials (PCMs) is an effective way of storing thermal energy. Latent heat storage enables high-energy storage density which reduces the footprint of the system and the cost. However, PCMs have very low thermal conductivities making them unsuitable for large-scale use without enhancing the effective thermal conductivity. In order to address, the low thermal conductivity of the PCMs, macroencapsulation of PCMs has been adopted as an effective technique. The macroencapsulation not only provides a self-supporting structure of PCM and separates the PCM from thermal fluids but also enhances the heat transfer rate. The current work involves study of various concepts of encapsulation of low cost inorganic PCMs. Sodium nitrate (NaNO3), a low cost PCM, was selected for thermal storage in a temperature range of 300 - 500˚C. Various techniques like electroless coatings, coatings using silicates, coatings with metal oxide (SiO2) and sand encapsulation are discussed. A novel technique of metal oxide coating was developed where firstly a high temperature polymer, such as, polymer (stable > 500˚C) was coated over PCM pellets, and cured, so that the pellet becomes insoluble in water as well as several organic solvents and later the metal oxide is coated over the pellet using self-assembly, hydrolysis, and simultaneous chemical oxidation at various temperatures. The coated PCM pellets were characterized.

Thermal energy storage (TES) is one of the most attractive and cost effective solutions to the intermittent generation systems like solar, wind and other renewable sources, compared to alternatives. It creates a bridge between the power supply and demand during peak hours or at times of emergency to ensure the continuous supply of energy. Among all the TES systems, latent heat thermal energy storage (LHTES) draws lots of interests as it has high energy density and can store or retrieve energy isothermally. Two major technical challenges associated with the LHTES are low thermal conductivity of the phase change materials (PCMs), and corrosion tendency of the containment vessel with the PCMs. Macro-encapsulation of the PCM is one of the techniques to encounter the low thermal conductivity issue as it will maximize the heat transfer area for the given volume of the PCM and restrict the PCMs to get in contact with the containment vessel. However, finding a suitable encapsulation technique that can address the volumetric expansion problem and compatible shell material are significant barriers of this approach. In the present work, an innovative technique to encapsulate PCMs that melt in the 100-350 oC temperature range was developed for industrial and private applications. This technique did not require a sacrificial layer to accommodate the volumetric expansion of the PCMs on melting. The encapsulation consisted of coating a non-reactive polymer over the PCM pellet followed by deposition of a metal layer by a novel non-vacuum metal deposition technique. The fabricated spherical capsules were tested in different heat transfer fluid (HTF) environments like air, oil and molten salt (solar salt). Thermophysical properties of the PCMs were investigated by DSC/TGA, IR and weight change analysis before and after the thermal cycling. Also, the constrained melting and solidification of sodium nitrate PCM inside the spherical capsules of different sizes were compared to explore the charging and discharging time. To accomplish this, three thermocouples were installed vertically inside the capsule at three equidistant positions. Low-density graphene was dispersed in the PCM to increase its conductivity and compared with pure PCM capsules. A laboratory scale packed-bed LHTES system was designed and built to investigate the performance of the capsules. Sodium nitrate (m.p. 306oC) was used as the PCM and air was used as the heat transfer fluid (HTF). The storage system was operated between 286oC and 326oC and the volumetric flow rate of the HTF was varied from 110 m3/h to 151 m3/h. The temperature distribution along the bed (radially and axially) and inside the capsules was monitored continuously during charging and discharging of the system. The effect of the HTF mass flow rate on the charging and discharging time and on the pressure drop across the bed was evaluated. Also, the energy and exergy efficiencies were calculated for three different flow rates. Finally, a step-by-step trial manufacturing process was proposed to produce large number of spherical capsules.

The present work investigates, numerically, the process of melting and solidification in hollow vertical cylinders, filled with air and phase change material (PCM). The PCM used is sodium nitrate, which expands upon melting. Therefore, a void must be present within the cylinder, which is filled with air. The influence of cylinder shape on melting time is determined. The numerical model takes both conductive and convective heat transfer into account during the melting process. The Volume-of-Fluid (VOF) model is used to track the interface between the PCM and air as the PCM melts. Three dimensionless numbers represent the characteristics of the problem, which are the Grashof, Stefan, and Prandtl numbers. The Stefan and Prandtl numbers are held constant, while the Grashof number varies. Inner Aspect Ratio (AR) is used to characterize the shape of the cylinder, which is defined as the ratio of the height to the diameter of the vertical cylinder. In this study, a range of AR values from 0.23 to 10 is investigated. Cylinders with small AR, corresponding to high Grashof numbers, lead to lower melting times compared with cylinders with high AR. The molten PCM velocity was also influenced greatly by this difference between solid PCM shape between high and low AR cases. Cylinders with small AR, corresponding to high Grashof numbers, lead to higher solidification times compared with cylinders with high AR. It was found that the velocity decreased during the solidification process, but the shape of the cylinder had an effect on the decrease. Natural convection velocity was found to decrease during the solidification process and, therefore, its effects diminish as solidification proceeds.

Sensible energy storage for solar domestic hot water (SDHW) systems is space consuming and heavy. Latent heat energy storage systems (LHESSs) offer a solution to this problem. However, the functionality of a LHESS during simultaneous charging/discharging, an operating mode encountered when used with a SDHW, had not been studied experimentally. A small scale vertical cylindrical LHESS, with dodecanoic acid as the phase change material (PCM), was studied during separate and simultaneous charging/discharging. Natural convection was found to have a strong influence during melting, but not during solidification. During simultaneous operation heat transfer was limited by the high thermal resistance of the solid PCM. However, when the PCM was melted, direct heat transfer occurred between the hot and cold heat transfer fluids, indicating the significance of the PCM phase on heat transfer in the system. The results of this research will lead to more optimally designed LHESS for use with SDHW. ?

This thesis work presents a numerical study of the effects of fins and thermal fluid velocities on the storage characteristics of a cylindrical latent heat energy storage system (LHESS). The work consists of two main components: 1. The development of a numerical method to study and solve the phase change heat transfer problems encountered in a LHESS during charging of the system, which results in melting of the phase change material (PCM). The numerical model is based on the finite element method. The commercial software COMSOL Multiphysics was used to implement it. The effective heat capacity method was applied in order to account for the large amount of latent energy stored during melting of a PCM, and the moving interface between the solid and liquid phases. The fluid flow, heat transfer and phase change processes were all validated using known analytical solutions or correlations. 2. Due to the low thermal conductivity of PCMs, the heat transfer characteristics of an enhanced LHESS was studied numerically. The effects of fins and the thermal fluid velocity on the melting rate of the PCM in the LHESS were analyzed. Results obtained for configurations having between 0 and 27 fins show that the heat transfer rate increases with addition of fins and thermal fluid velocity. The effect of the HTF velocity was observed to be small with few fin configurations since the thermal resistance offered by the LHESS system, mostly PCM, is vastly more important under these conditions; while its effect becomes more pronounced with addition of fins, since the overall thermal resistance decreases greatly with the addition of fins. The total energy stored after 12 hours for 0 and 27 fins configurations range between 3.6 MJ and 39.7 MJ for a thermal fluid velocity of 0.05 m/s and between 3.7 MJ and 57 MJ for a thermal fluid velocity of 0.5 m/s. The highest system efficiencies for the 0.05 m/s and 0.5 m/s, obtained with 27 fins configuration are 68.9% and 97.9% respectively.

This work presents an experimental and numerical study of phase change behaviour in a horizontal cylindrical latent heat energy storage system (LHESS). Fins with two orientations, straight fins and angled fins, are added into the PCM to enhance heat transfer. The PCM used in this study is lauric acid which has desirable thermal properties for LHESS. The experimental work concentrates on studying the heat transfer mechanism during phase change, impacts of the HTF inlet temperature and HTF flow rates. Moreover, heat transfer enhancement effectiveness of straight fins and angles fins is compared. Numerical model is simulated using COMSOL Multiphysics software package. It is observed that conduction is the dominant heat transfer mechanism during the initial stage of charging, and natural convection plays a more important role afterwards. Conduction plays a major role during solidification. Complete melting time is affected by the HTF inlet temperature and HTF flow rates.

Evapotransipration by vegetation cover is an important component of the water budget and energy balance in any ecosystem. A key to more improved water management therefore lies in improving our understanding of evapotranspiration, the process that drives water use by plants. Estimations of the turbulent fluxes are required for various applications in micrometeorology, hydrology, environmental studies and agriculture. Numerous methods for estimation of turbulent fluxes have been developed and tested. Direct measurements of fluxes are usually achieved by the eddy covariance (EC) method, which is considered as the most reliable. However, the application of the EC method is often problematic. The necessary sensors for wind, temperature and humidity must respond very fast (resolution of 10 Hz or better) and at the same time must not show noticeable drift. This makes them delicate, expensive and difficult to calibrate among other problems associated with the method. Due to their ability to integrate atmospheric processes along a path length that may range between a few hundred metres to a few kilometres, optical methods based on the analysis of scintillation appear to be an alternative and possible supplement to classical micrometeorological methods such as the EC method, which may provide local fluxes typically at the scale of 100 m. The use of the scintillometry technique in surface flux measurements is therefore gaining in popularity. The accuracy of the measurements obtained by one method is judged by comparison of the measurements obtained by those of another method considered as the standard. For turbulent flux measurements, the EC method is taken as the standard method for the determination of sensible heat fluxes. This research presents the measurement of sensible heat fluxes using the surface layer scintillometer (SLS). The SLS system used has a dual-beam and a recommended path length of between 50 and 250 m. The method was tested against the EC method for different Bowen ratio (f3) values, as required by the theory, under different atmospheric stability conditions, as well as for different wind directions relative to the SLS beam path and slanting beam path orientation. Also presented is an analysis of the different forms of the Monin-Obukhov Similarity (MOST) functions used in micrometeorology and suggested by various authors, done by comparing the resulting sensible heat flux measured by the SLS method with the ones calculated through an iterative determination of the Monin-Obukhov parameters. A comparison of the structure function parameter of temperature (Ci ) corrected for fJ and those measured (using SLS) was carried out, with the results showing very good correspondence between the corrected and uncorrected ci values, indicating that not correcting for fJ for SLS measured ci does not result in significant error in the resulting ci values, and hence sensible heat flux estimates. A comparison of the sensible heat flux Fh obtained using EC and SLS methods for fJ < 0.6 and fJ > 0.6 followed and the results also show good correspondence between the values obtained using the EC and SLS methods, although the agreement is slightly improved for cases when fJ > 0.6. A sensitivity analysis indicates that both the ECand SLS-measurements of Fh are influenced by fJ values. A sensitivity analysis on the influence of fJ on Fh measurements by both the EC and SLS methods further indicates that the influence of fJ on Fh measurements is not large enough to warrant correcting Fh measurements for fJ . The F" measurements by the EC method appears to be influenced more by fJ especially for fJ values less than 0.74. A comparison of the various methods for computing the empirical similarity functions used by MOST was also carried out and the results show a significant difference in the Fh computed following the various methods suggested by different researchers. As for the agreement between the EC and SLS methods determination of Fh for the different atmospheric stability conditions, there seems to be a better agreement in the Fh measurements as noted by correlation coefficients closer to 1 and greater tvalues obtained during unstable atmospheric conditions in the colder months of June and August while reduced agreement in the values is recorded in the warmer summer period from November to December. Also noted is a slight difference in the EC measurements compared to the SLS measurement of F". The difference in the measurements is noticed for unstable atmospheric conditions. Also noted is that EC and SLS measurements of Fh differ slightly when the atmospheric condition is nearneutral. However the agreement between the Fh values measured by the two measurement methods is still good. was set up in an inclined position, with the receiver set at 0.68 m above the ground level and transmitter at 1.68 m, resulting in an effective height difference of 1.00 m. There was generally good agreement in the 2-min measurements of F" by the two methods for the SLS set up in inclined position, with the 30-min data resulting in even better agreements. The findings confirm that the SLS set up does not impair its performance in measuring sensible heat fluxes. This also shows that the SLS would also work well in non-ideal (heterogeneous) conditions which the inclined optical beam path mimics. For those days when wind direction was mainly approximately perpendicular to the beam, the F" values obtained by SLS and EC methods are more in agreement than when the wind direction was either irregular or parallel to the SLS beam path. Wind speed also seems to influence the F" estimates by the two methods since the agreement in the Fh values obtained by the two methods is greater when wind speed is higher compared to times of the day when the wind speed is reduced. The atmospheric stability influences the peak position of footprint with the peak footprint position being further from the measurement point when the atmospheric stability condition is closer to stable as denoted by the Obukhov length of -5 and closer to the measurement point for convectively unstable atmospheric conditions as shown by the Obukhov length of -30. Also shown is that a larger fetch is required when the atmosphere is convectively unstable as indicated by the contours plotted on top of the footprint plots. In general, there seems to be very good agreement in the sensible heat flux values obtained by the two methods, especially since SLS offers areal-averaged sensible heat flux measurements compared to the EC method which basically provides a point measurement. The SLS method therefore offers a better alternative for obtaining sensible heat flux from larger and heterogeneous area - although to a limit of250 m since beyond 250 m, the method suffers from a saturation problem. / Thesis (Ph.D.)-University of KwaZulu-Natal, Pietermaritzburg, 2007.

The high population growth rate and rapid urbanization that the world is experiencing today has aggravated the competition for the already scarce resource ¡V water ¡V between the agricultural sector and the other economic sectors. Moreover, within the agricultural sector, water is increasingly being used for commercial plantations as opposed to growing food crops, threatening food security. Therefore, it is very important that this scarce resource is managed in an efficient and sustainable manner, for now and future use. This requires understanding the process of evaporation for accurate determination of water-use from agricultural lands. In the past, direct measurements of evaporation have proven difficult because of the cost and complexity of the available equipments, and level of expertise involved. This justifies a quest for relatively simple, accurate and inexpensive methods of determining evaporation for routine field applications. Estimation of sensible heat flux (H) from high frequency air temperature measurements and then calculating latent energy flux (ƒÜE) and hence evaporation as a residual of the shortened surface energy balance equation, assuming that closure is met, is appealing in this sense. Concurrent net irradiance (Rn) and soil heat flux (G) measurements can be conducted with relative ease for use in the energy balance equation. Alternately, evaporation can also be mathematically modelled, using single- or multi-layer models depending on vegetation cover, from less expensive routine meteorological observations. Therefore, the ultimate objective of this study is to estimate and model H and ƒÜE, and thereby evaporation, accurately over sparsely vegetated agricultural lands at low cost and effort. Temperature-variance (TV) and surface renewal (SR) methods, which use high-frequency (typically 2 to 10 Hz) air temperature measurements, are employed for estimation of H. The TV method is based on the Monin and Obukhov Similarity Theory (MOST) and uses statistical measures of the high frequency air temperature to estimate H, including adjustments for stability. The SR method is based on the principle that an air parcel near the surface is renewed by an air parcel from above and, to determine H, it uses higher order air temperature differences between two consecutive sample measurements lagged by a certain time interval. Single- and double-layer models that are based on energy and resistance combination theory were also used to estimate evaporation and H from sparse vegetation. Single- and double-layer models that were extended to include inputs of radiometric temperature in order to estimate H were also used. The transmission of solar irradiance to the soil beneath in sparse canopies is variable and depends on the vegetation density, cover and apparent position of the sun. A three-dimensional radiation interception model was developed to estimate this transmission of solar irradiance and was used as a sub-module in the double-layer models. Estimations of H from the TV (HTV), SR (HSR) and double-layer models were compared against H obtained from eddy covariance (HEC), and the modelled ƒÜE (single- and double-layer) were compared with that obtained from the shortened energy balance involving HEC. Besides, long-term ƒÜE calculated from the shortened energy balance using HTV and HSR were compared with those calculated using HEC. Unshielded and naturally-ventilated fine-wire chromel-constantan thermocouples (TCs), 75 ƒÝm in diameter, at different heights above the ground over sparse Jatropha curcas trees, mixed grassland community and bare fallow land were used to measure air temperature. A three-dimensional sonic anemometer mounted at a certain height above the ground surface was also used to measure virtual temperature and wind speed at all three sites. All measurements were done differentially at 10-Hz frequency. Additional measurements of Rn, G and soil water content (upper 60 mm) were also made. The Jatropha trees were planted in a 3-m plant and inter-row spacing in a 50 m ¡Ñ 60 m plot with the surrounding plots planted to a mixture of Jatropha trees and Kikuyu grass. Average tree height and leaf area index measurements were taken on monthly and bimonthly basis respectively. An automatic weather station about 10 m away from the edge of the Jatropha plot was also used to obtain solar irradiance, air temperature and relative humidity, wind speed and direction and precipitation data. Soil water content was measured to a depth of 1000 mm from the surface at 200 mm intervals. Soil and foliage surface temperatures were measured using two nadir-looking infrared thermometers with one mounted directly above bare soil and the other above the trees. The three-dimensional solar irradiance interception model was validated using measurements conducted on different trees and planting patterns. Solar irradiance above and below tree canopies was measured using LI-200 pyranometer and tube solarimeters respectively. Leaf area density (LAD) was estimated from LAI, canopy shape and volume measurements. It was also determined by scanning leaves using either destructive sampling or tracing method. The performance of the TV method over sparse vegetation of J. curcas, mixed grassland community and fallow land was evaluated against HEC. Atmospheric stability conditions were identified using (i) sensor height (z) and Obukhov length (L) obtained from EC and (ii) air temperature difference between two thermocouple measurement heights. The HTV estimations, adjusted and not adjusted for skewness (actual and estimated) of air temperature (sk), for unstable conditions only and for all stability conditions were used. An improved agreement in terms of slope, coefficient of determination (r2) and root mean square error (RMSE), almost over all surfaces, was obtained when the temperature difference rather than the z/L means of identifying stability conditions was used. The agreement between the HTV and HEC was improved for estimations adjusted for actual sk than not adjusted for sk. Improved agreement was also noted when HTV was adjusted using estimated sk compared to not adjusting for sk over J. curcas. The TV method could be used to estimate H for surfaces with varying homogeneity with reasonable accuracy. Long-term water-use of a fetch-limited sparse vegetation of J. curcas was determined as a residual of the shortened surface energy balance involving HTV and HSR and compared with those estimated using HEC. Concurrent measurements of Rn and G were also performed. The long-term water-use of J. curcas trees calculated from the shortened surface energy balance involving HTV and HSR agreed very well when compared with those obtained from HEC. The seasonal HTV and HSR also agreed very well when compared with HEC. Changes in structure of the canopy and environmental conditions appeared to influence partitioning of the available energy into H and ƒÜE. The seasonal total evaporation for the EC, TV and SR methods amounted to 626, 640 and 674 mm respectively with a total rainfall of 690 mm. Footprint analysis also revealed that greater than 80% of the measured flux during the day originates from within the surface of interest. The TV and SR methods, therefore, offer a relatively low-cost means for long-term estimation of H, and ƒÜE, hence the total evaporation, using the shortened surface energy balance along with measurements of Rn and G. Evaporation and biomass production estimations from tree crops requires accurate representation of solar irradiance transmission through the canopy. A relatively simple three-dimensional, hourly time-step tree-canopy radiation interception model was developed and validated using measurements conducted on isolated trees, hedgerows and tree canopies arranged in tramline mode. Measurements were obtained using tube solarimeters placed 0.5 m from each other starting from the base of a tree trunk in four directions, along and perpendicular to the row up to mid-way between trees and rows. Model-simulations of hourly radiant transmittance were in good agreement with measurements with an overall r2 of 0.91; Willmott.s index of agreement of 0.96; and general absolute standard deviation of 17.66%. Agreement between model-estimations and measurements, however, was influenced by distance and direction of the node from the tree trunk, sky conditions, symmetry of the canopy, and uniformity of the stand and leaf distribution of the canopy. The model could be useful in planning and management applications for a wide range of tree crops. Penman-Monteith (PM) equation and the Shuttleworth and Wallace (SW) model, representing single- and dual-source models respectively, were used to determine the total evaporation over a sparse vegetation of J. curcas from routine automatic weather station observations, resistance parameters and vegetation indices. The three-dimensional solar irradiance interception model was used as a sub-module in the SW model. The total evaporation from the sparse vegetation was also determined as a residual of the shortened surface energy balance using measurements of Rn, G and HEC. The PM equation failed to reproduce the .measured. daily total evaporation during periods of low LAI, with improved agreement with increased LAI. The SW model, however, produced total evaporation estimates that agreed very well with the .measured. with a slope of 0.96, r2 of 0.91 and RMSE of 0.45 mm for a LAI ranging from 0 (no leaves) to 1.83 m2 m-2. The SW model also estimated soil evaporation and plant transpiration separately, and about 66 % of the cumulative evaporation was attributed to soil evaporation. These findings suggest that the PM equation should be replaced by the SW model for surfaces that assume a range of LAI values during the growing season. The H was estimated using (i) SW model that was further developed to include surface radiometric temperature measurements; (ii) one-layer model, but linked with a two-layer model for estimation of excess resistance, that uses surface radiometric temperature; and (iii) the SW model (unmodified). The agreement between modelled and measured H, using 10-min data, was in general reasonably good with RMSE (W m-2) of 45.11, 43.77 and 39.86 for the three models respectively. The comparative results that were achieved from (iii) were not translated into the daily data as all models appeared to have a tendency to underestimate H. The resulting RMSEs for the daily H data for the three models were (MJ m-2) 1.16, 1.17 and 1.18 respectively. It appears that similar or better agreement between measured and estimated H can be forged without the need for surface radiometric temperature measurements. The study showed, in general, that high frequency air temperature measurements can be used to estimate H with reasonable accuracy using the simple and relatively low-cost TV and SR methods. Moreover, these methods can be used to calculate ƒÜE, hence ET, as a residual of the shortened surface energy balance equation along with measurements of Rn and G assuming that energy balance closure is met. The simple and low-cost nature of these methods makes replication of measurements easier and their robust nature allows long-term measurements of energy fluxes. The study also showed that H and ƒÜE can be modeled using energy and resistance combination equations with reasonable accuracy. It also reiterated that the SW-type models, which treat the plant canopy and soil components separately, are more appropriate for estimation of H and ƒÜE over sparse vegetation as opposed to the PM-type models. / Thesis (Ph.D.)-University of KwaZulu-Natal, Pietermaritzburg, 2010.