Este trabalho objetivou estabelecer relações entre evaporação de tanques Classe A e de 20 m2, em Piracicaba, SP e testar o desempenho de métodos meteorológicos de estimativa da evaporação de superfícies de água livre nessa localidade e em Jaboticabal, SP. Utilizou-se dados meteorológicos coletados em ambiente de um pequeno lago artificial em Piracicaba para cálculos dos métodos de balanço de energia-razão de Bowen, Penman, Priestley-Taylor e Penman simplificado. Dados meteorológicos observados sobre gramado foram usados para estimativa pelos métodos de Penman, Penman modificado, Dalton adaptado, Linacre e tanque Classe ASnyder. Os métodos foram avaliados utilizando índices de concordância e desempenho a partir de análise de regressão e de desvios por comparação com medidas do evaporímetro de 20 m2. Foram observadas relações lineares significativas entre os valores da evaporação dos tanques Classe A e de 20 m2 em todos os meses do ano em Piracicaba. A razão mensal entre as medidas dos evaporímetros apresentou variação ao longo do ano, com superestimativa pequena (até 3%) de fevereiro a junho e em novembro, e grande (10% a 15%) nos outros meses, em relação aos dados de Oliveira (1971) para a mesma localidade. A maior consistência estatística e representatividade da série de medidas, com 3 repetições de tanque Classe A, indicam que os fatores mensais de relações entre os dois evaporímetros podem substituir os obtidos por Oliveira (1971). Todos os métodos utilizados nos dois ambientes apresentaram boa concordância em termos médios com os observados no tanque de 20 m2, sempre com grande dispersão dos dados, diminuindo com o aumento do número de dias de contabilização da evaporação. O método de balanço de energia apresentou bom desempenho com cálculos realizados em períodos de 24 horas, mostrando desvios acentuados com contabilização no período diurno. O método de Penman mostrou bom desempenho quando utilizado para períodos de 24 horas, com subestimativas médias entre -6,0% e -10,5%, mas com superestimativas de +7,2 e +12,4% no período diurno. Os métodos de Priestley-Taylor (1972) e de Penman simplificado mostraram desvios médios pequenos (+2,7 e +5,7%), sendo considerados de bom desempenho. Para os métodos que usaram observações no posto meteorológico, o do tanque Classe ASnyder mostrou os menores desvios médios em Piracicaba (+3,9% a + 4,6%), mas apresentou maiores superestimativas (+14,0% a +14,4%) em Jaboticabal. Outra discrepância entre localidades foi observada para o método de Penmam modificado, com subestimativas entre 5,5 e 14,1% em Piracicaba, mas desvio de +1,0%, em Jaboticabal. O método original de Penman mostrou desvios similares em ambas localidades (-8,9% a 12,3%), enquanto o de Dalton adaptado apresentou desvios de +2,1% a 8,1%, mas com alta dispersão em todos dias de contabilização. Os métodos de Penman simplificado e o de Priestley-Taylor, no caso do lago, e o tanque Classe ASnyder, no caso do posto meteorológico, foram os que no geral apresentaram resultados mais comparáveis com os do tanque de 20 m2, sendo que par ao último deve-se atentar para o desvio sistemático em Jaboticabal. / This study aimed to establish relations between evaporation Class A pan and 20 m2 tank, in Piracicaba region, State of Sao Paulo, Brazil and to test the performance of meteorological methods to estimate open-water surfaces evaporation in this same region and also in Jaboticabal, SP. Computation of daily evaporation was performed by using weather data, collected in the environment of a small artificial lake in Piracicaba, and the following methods: energy balance-Bowen ratio, Penman, Priestley-Taylor and Penman simplified. The methods of Penman, Penman adapted, Dalton adapted, Linacre and of Class A pan-Snyder were evaluated by using weather station data in both regions. The results of the methods were compared with the measurements of the 20 m2 tank evaporation, in a daily basis, and also for two, three or for the mean values for all days of measurements in every month, using indexes of agreement and performance from the regression analysis and deviations from the estimated and measured values. The linear relations between evaporation measured by the two evaporimeters were significative for all months in Piracicaba. The ratio of the two measurements showed variation throughout the year, with a small overestimation (about 3%) from February to June and November, and greater (10% to 15%) in the other months, compared to results of Oliveira (1971) for the same region. Due to the statistical consistency and representativeness of the series of 9 years and 3 repetitions of Class A pan, it is indicated that the factors of monthly relations between the two evaporímeters observed in this study can advantageously replace those obtained by Oliveira (1971). All methods used in the two environments showed good agreement, in average basis, with measurements by 20 m2 tank, with large dispersion of data, which decreased with the increase of the time step. The energy balance method showed good performance when the calculations were performed in 24-h periods, but large deviations compared to the tank with data integrated over the daylight period. The Penman method showed good performance when used for periods of 24 h, with an average underestimation between - 6.0 and -10.5% when integrated over the 24-h period, but with overestimation of +7.2 and +12.4% for the daylight period. The methods of Priestley-Taylor (1972) and Penman simplified showed small deviations between +2.7 and + 5.7%, being considered as having good performance compared to the evaporimeter. Among the methods that used weather station data, the class A- Snyder was the one that showed the lowest average deviation in Piracicaba (+3.9 to + 4.6%), but the greatest overestimation (+14.0% to 14.4%) in Jaboticabal. Another discrepancy between locations was observed for the Penman method modified, with underestimation between 5.5% and 14.1% in Piracicaba, but deviation of +1.0% in Jaboticabal. The original Penman method showed similar deviations in the two regions (-8.9% to 12.3%), while Dalton adapted method showed deviations of +2.1% to 8.1%. The methods of Penman simplified and Priestley-Taylor, used in the lake environment, and the Class A pan-Snyder method, with the use of wheater station data, were those that showed most comparable results with the 20 m2 tank evaporation, but in Jaboticabal, it is necessary to take to account the observed systematic mean deviation.

Submitted by Erika Demachki (erikademachki@gmail.com) on 2014-10-15T19:25:11Z No. of bitstreams: 2 Dissertação - Diogo Silva Pena - 2013.pdf: 3166468 bytes, checksum: 0b15748ebfd5186f66f8e329d0d57146 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Approved for entry into archive by Jaqueline Silva (jtas29@gmail.com) on 2014-10-16T18:11:42Z (GMT) No. of bitstreams: 2 Dissertação - Diogo Silva Pena - 2013.pdf: 3166468 bytes, checksum: 0b15748ebfd5186f66f8e329d0d57146 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Made available in DSpace on 2014-10-16T18:11:42Z (GMT). No. of bitstreams: 2 Dissertação - Diogo Silva Pena - 2013.pdf: 3166468 bytes, checksum: 0b15748ebfd5186f66f8e329d0d57146 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Previous issue date: 2013-12-20 / Among the plants recommended for the production of biodiesel there is the Jatropha (Jatropha curcas L.), however there is a lack of studies for this crop in particular regarding its management. Climate variables can sometimes restrain the growth and development of certain kinds of plant specious, especially by extremes ranges of temperature and rainfall. This study`s purpose was to define areas suitable in Goiás state for the cultivation of Jatropha by conducting a study based on an agroclimatic zoning using a Geographic Information System. Some other products were also created, as the climate characterization of the state, by generating temperature, precipitation and evapotranspiration maps, as well as defining the water balance of the state. The zoning does not show evidences of any area that is considered unsuitable for cultivation of Jatropha in the state, which much of the state appears to be fully suitable, with some marginal regions, either by water deficiency, or by thermal deficiency. A total of 64,6% of Goiás is suitable for cultivation of Jatropha, while 35,4% have a condition of marginality for the development of the plant. Among all the areas defined as marginal for the cultivation of Jatropha, 28,8% is considered to be marginal by water deficiency, and 6,6% is considered to be marginal by thermal deficiency. / A cultura do pinhão manso destaca-se no cenário do agronegócio brasileiro como uma das mais promissoras para a produção de biodiesel, muito embora sejam poucos os estudos voltados para a implantação, manejo e técnicas de cultivo dessa cultura. As variáveis climáticas podem em algumas regiões restringir o plantio e a condução da cultura, em especial a temperatura e a pluviosidade. Por meio da caracterização do clima do estado através da geração de mapas de temperatura, precipitação e evapotranspiração, bem como a definição do balanço hídrico do estado, objetivou-se neste trabalho, definir áreas aptas no estado de Goiás para o cultivo do pinhão manso por meio de realização de um zoneamento agroclimático. Para isso foram utilizados dados climáticos de precipitação pluviométrica de 114 estações com média de 20 anos de registro, e dados climáticos de temperatura do ar referentes à 34 estações, também com média de 20 anos de registro, que foram usadas para estimar os valores médios de temperatura do ar para as demais localidades do estado, via modelo de regressão linear múltipla. O zoneamento não evidenciou nenhuma área que seja considerada inapta ao cultivo do pinhão no estado, estando grande parte do território sujeito à aptidão plena da cultura, com algumas regiões marginais, ora por deficiência hídrica, ora por deficiência térmica. Do total do estado de Goiás, 64,6% do território goiano é apto ao cultivo do pinhão manso, enquanto que 35,4% apresenta alguma condição de marginalidade para o desenvolvimento da cultura. Deste total de áreas definidas como marginais ao cultivo do pinhão, 28,8% é tida como marginal por deficiência hídrica, e 6,6% é tida como marginal por deficiência térmica.

O amendoim é uma cultura de grande importância econômica que apresenta facilidade na sua condução e é cultivada em diversas regiões do Brasil como também no exterior. Portanto, objetivou-se com a presente pesquisa verificar a influência de ambientes protegido de cultivo e aplicações, via fertirrigação, de íons cálcio e potássio, sobre a cultura do amendoim, em dois ciclos distintos. O experimento foi instalado no posto agrometeorológico da ESALQ/USP, do Departamento de Engenharia de Biossistemas, e dividido em dois ciclos vegetativos, onde foram utilizados doze lisímetros, preenchidos com solo Argissolo Vermelho Amarelo. Foram considerados sobre os lisímetros três ambientes: a céu aberto e utilizando filmes plásticos (75 e 150 micrometros de espessura). Cada ambiente foi constituído por 4 lisímetros. As análises estatísticas foram realizadas separadamente, a saber: ambientes protegidos diferenciados pela espessura do filme plástico; ambiente a céu aberto. Para o manejo da irrigação a fertirrigação foi determinada em função do período das aplicações e do turno de rega. As variáveis climáticas foram coletadas por data logger (Hobbo) (temperatura do ar, umidade relativa do ar, temperatura do ponto de orvalho, densidade de fluxo radiante) e precipitação pluvial. Foram avaliados os parâmetros altura da planta e diâmetro do caule. Os parâmetros de pós-colheita avaliados foram produção, produtividade, número de vagens e umidade do grão. Foram determinadas as concentrações dos íons Ca++ e K+ e o potencial hidrogeniônico (pH). Pôde-se concluir que o íon potássio apresentou, no ciclo mais seco, uma baixa mobilidade no solo para os três ambientes estudados e o íon cálcio também apresentou baixa mobilidade no solo para os três ambientes, porém, para os dois ciclos da cultura. Dos três ambientes avaliados, o mais indicado para o cultivo do amendoim foi o instalado a céu aberto, no ciclo mais chuvoso da cultura. Em termos de ambiente protegido, o mais indicado para o cultivo e condução da cultura foi o de 75 µm. Os parâmetros vegetativos altura da planta e diâmetro de caule obtiveram resultados expressivos somente no primeiro ciclo da cultura nos três ambientes avaliados. Entre os ambientes protegidos, no qual foi instalado o filme plástico de 75 µm, para obtenção de uma boa produtividade no primeiro ciclo da cultura, deve ser aplicada uma faixa média de temperatura de ponto de orvalho de 20,1 0C; no ambiente a céu aberto, a faixa média de temperatura de ponto de orvalho deve estar em torno de 20,79 0C para obtenção de uma boa produtividade no ciclo mais chuvoso, ao passo que uma temperatura média de 6,22 0C no ciclo mais seco torna-se um fator limitante a não realização do plantio. A melhor época para o cultivo e condução do amendoim, seja num ambiente a céu aberto como num ambiente protegido, é a denominada amendoim das águas, compreendida geralmente entre setembro ou outubro até fevereiro do próximo ano. A concentração para uma aplicação ideal dos nutrientes estudados foi a de 100% da demanda requerida pela cultura, no período normal de cultivo do amendoim (ciclo chuvoso), considerando o fator pós-colheita. / The peanut is a crop of great economic importance that is easy to grow and is cultivated in several areas of Brazil and the world. The aim was to verify the influence of different environments and applications of calcium and potassium ions by fertigation on the peanut crop. The research was carried out at ESALQ/USP, of Biosystems Engineering Department. This was divided in two vegetative cycles, in which twelve lysimeters were used and filled with an Arenic Abruptic Paleudult soil. The treatments considered for the lysimeters was three solar environments of open sky and using two thicknesses of plastic films (75 and 150 micrometers of thickness). Each environment was replicated by 4 lysimeters. Irrigation management for the fertigation was determined as a function of the period of applications and the interval of irrigation. Climatic variables were collected by data logger (Hobo) (air temperature, air humidity, dew point temperature, density of radiant flow and precipitation. Crop height and stem diameter were evaluated as growth parameters. The parameters after harvest were total yield, number of nuts, and kernel moisture. Calcium and Potassium concentrations were also obtained as well as hydrogenionic potential. It was determine that the potassium ion presented in driest cycle with low mobility in soil for the three studied environments and the calcium ion also had low mobility in soil for the three environments for both cycles of the crop. Of the three evaluated environments, the best for tillage of peanut was that installed to open sky in the rainiest cycle of the crop. When considering the two protected environments, the best for tillage and growth of the crop was that of 75 µm. The vegetative parameters, plant height and stem diameter, only obtained expressive results in first cycle of the crop in the three evaluated environments. For protected environments, in which plastic film of 75 µm was installed, a good yield in first cycle of the crop required a medium value of dew point temperature of 20.1 oC, in the environment to open sky, the medium value of dew point temperature should be around 20.79 oC for obtaining of a good yield in the rainiest cycle. A medium temperature of 6.22 oC in the driest cycle becomes a limiting factor to the non development of the crop. The best time for tillage and growth of peanut in an environment to open sky and in a protected environment, is denominated \"waters peanut\", generally occurred between September or October through February of the following year. An ideal application of nutrients was 100% of the demand required by the crop, in normal period of peanuts tillage (raining cycle) considering factor after harvest.

Studies of trace gas fluxes have advanced the understanding of bulk interactions between the atmosphere and ecosystems. Micrometeorological instrumentation is currently unable to resolve vertical scalar sources and sinks within plant canopies. Inverted analytical Lagrangian equations provide a non-intrusive method to calculate source distributions. These equations are based on Taylor's (1921) description of scalar dispersion, which requires a measure of the degree of correlation between turbulent motions, defined by the Lagrangian length scale (L). Inverse Lagrangian (IL) analyses can be unstable, and the uncertainty in L leads to uncertainty in source predictions. A review of the literature on studies using IL analysis with various scalars in a multitude of canopy types found that parameterizations where L reduces to zero at the ground produce better results in the IL analysis than those that increase closer to the ground, but no individual L parameterization gives better results than any other does. The review also found that the relationship between L and the measurable Eulerian length scale (Le) may be more complex in plant canopies than the linear scaling investigated in boundary layer flows. The magnitude and profile shape of L was investigated within a corn and a forest canopy using field measurements to constrain an analytical Lagrangian equation. Measurements of net CO2 flux, soil-to-atmosphere CO2 flux, and in-canopy profiles of CO2 concentrations provided the information required to solve for L in a global optimization algorithm for half hour intervals. For dates when the corn was a strong CO2 sink, and for the majority of dates for the forest, the optimization frequently located L profiles that follow a convex shape. A constrained optimization then smoothed the profile shape to a sigmoidal equation. Inputting the optimized L profiles in the forward and inverse Lagrangian equations leads to strong correlations between measured and calculated concentrations (corn canopy: C_{calc} = 1.00C_{meas} +52.41 mumol m^{-3}, r^2 = 0.996; forest canopy: C_{calc} = 0.98C_{meas} +276.5 mumol m^{-3}, r^2 = 0.99) and fluxes (corn canopy: F_{soil} = 0.67F_{calc} - 0.12 mumol m^{-2}s^{-1}, r^2 = 0.71, F_{net} = 1.17F_{calc} + 1.97mumol m^{-2}s^{-1}, r^2 = 0.85; forest canopy: F_{soil} = 0.72F_{calc} - 1.92 mumol m^{-2}s^{-1}, r^2 = 0.18, F_{net} = 1.24F_{calc} + 0.65 mumol m^{-2}s^{-1}, r^2 = 0.88). In the corn canopy, coefficients of the sigmoidal equation were specific to each half hour and did not scale with any measured variable. Coefficients of the optimized L equation in the forest canopy scaled weakly with variables related to the stability above the canopy. Plausible L profiles for both canopies were associated with negative bulk Richardson number values. / Funding from NSERC.

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.

Current mining methods used to extract coal from underground mine workings disturb the natural environment and the existing stable geological structures. As a result, the ingress of water into the mines increases and the quality of the water passing through the mine workings deteriorates, irrespective of the operational status of the mines. Water ingress is generated by regional aquifers, local aquifers, recharge from the surface through rainfall, natural drainage paths on the surface, and surface water bodies. The quality of water in the mines deteriorates as a result of contact with the remaining coal in the mine workings. Mining can therefore cause an increased influx of water into a mine and the degradation of this water. The solution to reducing the impact of mines on the environment is to prevent, or at least reduce, the amount of water entering the mines, and to manage this water to prevent further degradation in water quality. This study focused on afforestation with Eucalyptus viminalis trees to manage or inhibit ingress of water into underground mine workings. The hypothesis of this study was that a change in vegetation, from grassland to fast-growing and potentially high water-using trees like Eucalyptus. could possibly reduce the drainage of water below the root-zone and into the mine workings. The hypothesis was tested by estimating the components of the soil water balance for a grassland site and a Eucalyptus tree site. The research site was situated in Mpumalanga, (260 36' Sand 290 08' E, 1650 m a.m.s.l.), one of South Africa's major coal bearing areas. Although the Secunda area is a treeless environment and conditions are not optimal for forestry, some Eucalyptus species are suited for conditions (frost and periodic droughts) encountered in this area. The soil water balance of grassland and E. viminalis trees were studied through a field experiment and a long-term (30 years) modelling exercise. Total evaporation of the grassland site was estimated using the Bowen ratio energy balance technique. The transpiration of six representative E. viminalis trees were estimated using the heat pulse velocity technique. The soil water storage changes at both sites were determined from the soil water content, estimated using water content reflectometers. Measurements were performed in a smectic clay soil which resulted in measurements difficulties. Vertical cracks were formed under soil drying. To establish the importance of climate and plant growth on the drainage beyond the root-zone, the soil water balance of a grassland and an E. viminalis site were simulated over a 30-year period with the Soil Water Atmosphere Plant (SWAP) model. It was concluded from the comparative field experiment and modelling, that a change in vegetation from grassland to E. viminalis will reduce the drainage of water below the root-zone, especially under above-average rainfall conditions. The reduction in drainage beyond the root-zone at the E. viminalis sites, compared to the grassland site, was demonstrated in the modelling exercise and can be deduced from the total evaporation and soil water storage estimated at both sites. The results from the field experiment confirmed the modelling results and showed that usually there were higher transpiration rates for the E. viminalis tree site, compared to the grassland site. The higher transpiration rates for E. viminalis trees resulted in lower relative saturation of soil layers and lower profile soil water contents at the E. viminalis site, and higher daily soil water storage changes at the E. viminalis site compared to the grassland site. These differences were more pronounced during winter when the grassland was dormant. The results from the modelling exercise showed that an E. viminalis tree stand, with a closed canopy, reduced drainage below the root-zone compared to a grassland. The drainage at the grassland site contributed to up to 54 % of the rainfall, compared to the 43 % at the E. viminalis site. However, under below-average rainfall conditions the annual drainage at both sites, were similar. Further, the absolute magnitude of the drainage was similar to the total evaporation at the grassland site under certain conditions. The results not only suggest that a change in vegetation, from grassland to E. viminalis trees, would reduce the drainage beyond the root-zone, but that it may delay the onset of drainage. Under above-average rainfall conditions, the modelled drainage at the E. viminalis site only exceeded 20 mm, a month later than at the grassland site. The simulation results also showed that under conditions of aboveaverage rainfall, drainage occurs whenever the rainfall exceeds the long-term average rainfall, irrespective of the existing vegetation. However, when the rainfall is belowaverage drainage at both sites are limited to large rainfall events. This simulation showed that over a period of eight years, E. viminalis trees could potentially reduce the drainage by 1235 mm more than grassland, which is equivalent to 1540 m3 ha- I a-I, or 1.54 Me ha- I a-I. The annual average reduction in drainage below the root-zone caused by E. viminalis trees (1.79 Mf ha-1 a-\ is a small reduction when compared to the influx of water into mineworkings. E.g. the influx of water into a bord-and-pillar mine range between 0.5 and 4 Mt d-I per area mined and up to 17000 Mt d-I per area mined under high extraction mining (Hodgson and Krantz, 1998; Hodgson et aI., 2001). This work gave a comprehensive account of the differences in the soil water relations of grassland and E. viminalis trees overlying coal mine working. Few other studies in South Africa compared the total evaporation and soil water relations of grassland and E. viminalis trees in so much detail. State of the art monitoring techniques were used and produced valuable comparison of their use in expansive clay profiles. The work should contribute to management decisions focussed on limiting ingress of water into mine workings. / Thesis (Ph.D.)-University of Natal, Pietermaritzburg, 2003.

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.

Increased pressure on the available limited water resources for agricultural production has a significant impact on sugarcane production. Routine monitoring of evaporation with reliable accuracy is essential for irrigation scheduling, for more efficient use of the available water resources and for management purposes. An indirect method for estimating evaporation involves measuring the sensible heat flux (H) from which latent energy flux and hence total evaporation can be calculated, as a residual using the shortened energy balance from measurements of net irradiance and soil heat flux. Various methods for measuring H may include Bowen ratio energy balance, eddy covariance (EC), flux variance (FV), optical scintillation, surface renewal (SR) and temperature variance (TV). Each method has its own advantages and disadvantages, in terms of method theoretical assumptions, accuracy, complexity, cost, fetch requirements and power consumption. The TV and SR methods are inexpensive and reasonably simple with a reduced power requirement compared to other methods since they require high frequency air temperature data which is obtained by using an unshielded naturally-ventilated type-E fine-wire thermocouple at a single point above the canopy surface. The TV method is based on the Monin-Obukhov similarity theory (MOST) and uses the mean and standard deviation of the air temperature for each averaging period. Currently, there are two TV methods used for estimating sensible heat flux (HTV) at sub-hourly time intervals, one includes adjustment for stability, and a second that includes adjustment for air temperature skewness. Another method used to estimate sensible heat flux from the mean and standard deviation of air temperature is based on MOST and uses spatial second-order air temperature structure function. For the TV method adjusted for stability and the method based on MOST that uses a spatial second-order air temperature structure function, the Monin-Obukhov atmospheric stability parameter () is needed. The parameter  can be estimated from EC measurements or alternatively estimated independently using an iteration process using horizontal wind speed measurements. The TV method including adjustment for air temperature skewness requires the mean and standard deviation of the air temperature and air temperature skewness for each averaging time period as the only input. The SR method is based on the coherent structure concept. Currently, there are various SR models method for estimating sensible heat flux. These include an ideal SR analysis model method based on an air temperature structure function analysis, the SR analysis model with a finite micro-front period, combined SR with K-theory and combined SR model method based on MOST. The ideal SR analysis model based on an air temperature structure function analysis should be calibrated to determine the SR weighting factor (). The other SR approaches require additional measurements such as crop height and horizontal wind speed measurements. In all of the SR approaches, air temperature time lags are used when calculating the air temperature structure functions. In this study, the performance of TV and SR methods were evaluated for estimation of sensible heat and latent energy fluxes at different heights for air temperature time lags of 0.4 and 0.8 s for daytime unstable conditions against EC above a sugarcane canopy at the Baynesfield Estate in KwaZulu-Natal, South Africa. For all methods, latent energy flux (LE) and hence evaporation was estimated as a residual from the shortened energy balance equation using H estimates and net irradiance and soil heat flux density measurements. The ideal SR analysis model method based on an air temperature structure function analysis approach was calibrated and validated against the EC method above the sugarcane canopy using non-overlapping data sets for daytime unstable conditions during 2008. During the calibration period, the SR weighting factor was determined for each height and air temperature time lag. The magnitude of ranged from 0.66 to 0.55 for all measurement heights and an air temperature time lag of 0.8 s. The value increased with a decrease in measurement height and an increase in air temperature time lag. For the validation data set, the SR sensible heat flux (HSR) estimates corresponded well with EC sensible heat flux (HEC) for all heights and both air temperature time lags. The agreement between HSR and HEC improved with a decrease in measurement height for the air temperature time lag of 0.8 s. The best HSR vs HEC comparisons were obtained at a height of 0.20 m above the crop canopy using = 0.66 for an air temperature time lag of 0.8 s. The residual estimates of latent energy flux by SR and EC methods were in good agreement. The LESR at a height of 0.20 m above the canopy yielded the best comparisons with LEEC estimated as a residual. The performance of the TV method, including adjustment for stability, and / Thesis (Ph.D.)-University of KwaZulu-Natal, Pietermaritzburg, 2010.

The demand for the world’s increasingly scarce water supply is rising rapidly, challenging its availability for agriculture and other environmental uses, especially in water scarce countries, such as South Africa, with mean annual rainfall is well below the world’s average. The implementation of effective and sustainable water resources management strategies is then imperative, to meet these increasingly growing demands for water. Accurate assessment of evaporation is therefore crucial in agriculture and water resources management. Evaporation may be estimated using different micrometeorological methods, such as eddy covariance (EC), Bowen ratio energy balance (BR), surface renewal (SR), flux variance (FV), and surface layer scintillometry (SLS) methods. Despite the availability of different methods for estimating evaporation, each method has advantages and disadvantages, in terms of accuracy, simplicity, spatial representation, robustness, fetch, and cost. Invoking the shortened surface energy balance equation for which advection and stored canopy heat fluxes are neglected, the measurement of net irradiance, soil heat flux, and sensible heat flux allows the latent energy flux and hence the total evaporation amount to be estimated. The SR method for estimating sensible heat, latent energy, and other scalars has the advantage over other micrometeorological methods since it requires only measurement of the scalar of interest at one point. The SR analysis for estimating sensible heat flux from canopies involves high frequency air temperature measurements (typically 2 to 10 Hz) using 25 to 75 ìm diameter fine-wire thermocouples. The SR method is based on the idea that parcel of air near a surface is renewed by an air parcel from above. The SR method uses the square, cube, and fifth order of two consecutive air temperature differences from different time lags to determine sensible heat flux. Currently, there are three SR analysis approaches: an ideal SR analysis model based on structure function analysis; an SR analysis model with finite micro-front period; and an empirical SR analysis model based on similarity theory. The SR method based on structure function analysis must be calibrated against another standard method, such as the eddy covariance method to determine a weighting factor á which accounts for unequal heating of air parcels below the air temperature sensor height. The SR analysis model based on the finite micro-front time and the empirical SR analysis model based on similarity theory need the additional measurement of wind speed to estimate friction velocity. The weighting factor á depends on measurement height, canopy structure, thermocouple size, and the structure function air temperature lag. For this study, á for various canopy surfaces is determined by plotting the SR sensible heat flux SR H against eddy covariance EC H estimates with a linear fit forced through the origin. This study presents the use of the SR method, previously untested in South Africa, to estimate sensible heat flux density over a variety of surfaces: grassland; Triffid weed (Chromolaena odorata); Outeniqua Yellow wood (Podocarpus Falcatus) forest; heterogeneous surface (Jatropha curcas); and open water surface. The sensible heat flux estimates from the SR method are compared with measurements of sensible heat flux obtained using eddy covariance, Bowen ratio, flux variance, and surface layer scintillometer methods, to investigate the accuracy of the estimates. For all methods used except the Bowen ratio method, evaporation is estimated as a residual using the shortened energy balance from the measured sensible heat and from the additional measurements of net irradiance and soil heat flux density. Sensible heat flux SR H estimated using the SR analysis method based on air temperature structure functions at a height of 0.5 m above a grass canopy with a time lag r = 0.5 s, and á =1 showed very good agreement with the eddy covariance EC H , surface layer scintillometer SLS H , and Bowen ratio BR H estimates. The half-hourly latent energy flux estimates obtained using the SR method SR ë E at 0.5 m above the grass canopy for a time lag r = 0.5 s also showed very good agreement with EC ë E and SLS ë E . The 20-minute averages of SR ë E compared well with Bowen ratio BR ë E estimates. Sensible heat and latent energy fluxes over an alien invasive plant, Triffid weed (C. odorata) were estimated using SR , EC , FV and SLS methods. The performance of the three SR analysis approaches were evaluated for unstable conditions using four time lags r = 0.1, 0.4, 0.5, and 1.0 s. The best results were obtained using the empirical SR method with regression slopes of 0.89 and root mean square error (RMSE) values less than 30 W m-2 at measurement height z = 2.85 and 3.60 m above the soil surface for time lag r = 1.0 s. Half-hourly SR H estimates using r = 1.0 s showed very good agreement with the FV and SLS estimates. The SR latent energy flux, estimated as a residual of the energy balance ë ESR , using time lag r = 1.0 s provided good estimates of EC ë E , FV ë E , and SLS ë E for z = 2.85 and 3.60 m. The performance of the three SR analysis approaches for estimating sensible heat flux above an Outeniqua Yellow wood stand, were evaluated for stable and unstable conditions. Under stable conditions, the SR analysis approach using the micro-front time produced more accurate estimates of SR H than the other two SR analysis approaches. For unstable conditions, the SR analysis approach based on structure functions, corrected for á using EC comparisons produced superior estimates of SR H . An average value of 0.60 is found for á for this study for measurements made in the roughness sublayer. The SR latent energy flux density estimates SR ë E using SR H based on structure function analysis gave very good estimates compared with eddy covariance ( EC ë E ) estimates, with slopes near 1.0 and RMSE values in the range of 30 W m-2. The SR ë E estimates computed using the SR analysis approach using the micro-front time also gave good estimates comparable to EC ë E . The SR and EC methods were used to estimate long-term sensible heat and latent energy flux over a fetch-limited heterogeneous surface (J. curcas). The results show that it is possible to estimate long-term sensible heat and latent energy fluxes using the SR and EC methods over J. curcas. Continuous measurements of canopy height and leaf area index measurements are needed to determine á . The weighting factor á was approximately 1 for placement heights between 0.2 and 0.6 m above the Jatropha tree canopy. The daily sensible heat and latent energy flux estimates using the SR analysis gave excellent estimates of daily EC sensible heat and latent energy fluxes. Measurements of sensible heat and estimates of the latent energy fluxes were made for a small reservoir, using the SR and EC methods. The SR sensible heat flux SR H estimates were evaluated using two air temperature time lags r = 0.4 and 0.8 s at 1.0, 1.3, 1.9, 2.5 m above the water surface. An average á value of 0.175 for time lag r = 0.4 s and 0.188 for r = 0.8 s was obtained. The SR H and EC H estimates were small (-40 to 40 W m-2). The heat stored in water was larger in magnitude (-200 to 200 W m-2) compared to the sensible heat flux. The SR and EC latent energy fluxes were almost the same in magnitude as the available energy, due to the small values of the sensible heat fluxes. The daily evaporation rate ranged between 2.0 and 3.5 mm during the measurement period. The SR method can be used for routine estimation of sensible heat and latent energy fluxes with a reliable accuracy, over a variety of surfaces: short canopies, tall canopies, heterogeneous surface, and open water surface, if the weighting factor á is determined. Alternatively, the SR method can be used to estimate sensible heat flux which is exempt from calibration using the other two SR analysis approaches, with additional measurement of wind speed for estimating friction velocity iteratively. The advantages of the SR method over other micrometeorological methods are the relatively low cost, easy installation and maintenance, relatively low cost for replicate measurements. These investigations may pave the way for the creation of evaporation stations from which real-time and sub-hourly estimates of total evaporation may be obtained relatively inexpensively. / Thesis (Ph.D.)-University of KwaZulu-Natal, Pietermaritzburg, 2008.

A well-managed irrigation scheduling system needs a rapid, preCIse, simple, costeffective and non-destructive soil water content sensor. The PRl profile probe and Diviner 2000 were used to determine the timing and amount of irrigation of three cover crops (Avena sativa L., Secale cereale L. and Lolium multiflonlm Lam.), which were planted at Cedara, KwaZulu-Natal. The PRl profile probe was first calibrated in the field and also compared with the Diviner 2000. For the calibration of the PRl profile probe the factory-supplied parameters (aJ = 8.4 and ao = 1.6) showed good correlation· compared to the soil-estimated parameters (aJ = 11.04 and ao = 1.02). The factorysupplied parameters gave a linear regression coefficient (r2 ) of 0.822 and root mean square error (RMSE) of 0.062. The soil-estimated parameter showed a linear regression coefficient of 0.820 with RMSE of 0.085. The comparison between the soil water content measured using the PR1 profile probe and Diviner 2000 showed a linear regression coefficient of 0.947 to 0.964 with a range of RMSE of 0.070 to 0.109 respectively for the first 100 to 300 mm soil depths. The deeper depths (400, 600 and 1000 mm) showed a linear regression coefficient ofO.716to 0.810 with a range of 0.058 to 0.150 RMSE. These differences between the shallow and deeper depths could be due to soil variability or lack of good contact between the access tube and the surrounding soil. To undertake irrigation scheduling using the PRl profile probe and Diviner 2000, the soil water content limits were determined using field, laboratory and regression equations. The field method was done by measuring simultaneously the soil water content using the PR1 profile probe and soil water potential using a Watermark sensor and tensiometers at three depths (100, 300 and 600 mm) from a 1 m2 bare plot, while the soil dries after being completely saturated. The retentivity function was developed from these measurements and the drained upper limit was estimated to be 0.355 m3 m-3 when the drainage from the pre-wetted surface was negligible. The lower limit was calculated at -1500 kPa and it was estimated to be 0.316 m3m,3. The available soil water content, which is the difference between the upper and lower limit, was equal to 0.039 m3 m,3. In the laboratory the soil water content and matric potential were measured from the undisturbed soil samples taken from the edge of the 1 m2 bare plot before the sensors were installed. Undisturbed soil samples were taken using a core sampler from 100 to 1000 mm soil depth in three replications in 100 mm increments. These undisturbed soil samples were saturated and subjected to different matric potentials between -1 to -1500 kPa. In the laboratory, the pressure was increased after the cores attained equilibrium and weighed before being subjecting to the next matric potential. The retentivity function was then developed from these measurements. The laboratory method moved the drained upper limit to be 0.390 m3 m,3 at -33 kPa and the lower limit be 0.312 m3m-3 at -1500 kPa. The regression equation, which uses the bulk density, clay and silt percentage to calculate the soil water content at a given soil water potential, estimated the drained upper limit to be 0.295 m3m-3at -33 kPa and the lower limit 0.210 m3 m,3 at -1500 kPa. Comparison was made between the three methods using the soil water content measured at the same soil water potential. The fieldmeasured soil water content was not statistically the same with the laboratory and estimated soil water content. This was shown from the paired-t test, where the probability level (P) for the laboratory and estimated methods were 0.011 and 0.0005 respectively at 95 % level of significance. However, it showed a linear regression coefficient of 0.975 with RMSE of 0.064 when the field method was compared with the laboratory method. The field method showed a linear regression coefficient of 0.995 with RMSE of 0.035 when compared with the estimated method. The timing and amount of irrigation was determined using the PR1 profile probe and Diviner 2000. The laboratory measured retentivity function was used to define the fill (0.39 m3 m-3 ) and high refill point (0.34 m3 m-3 ). The soil water content was measured using both sensors two to three times per week starting from May 29 (149 day of year, 2002) 50 days after planting until September 20 (263 day of year) 11 days before harvesting. There were five irrigations and twenty rainfall events. The next date of irrigation was predicted graphically using, the PRl profile probe measurements, to be on 3 September (246 day of year) after the last rainfall event on 29 August (241 day of year) with 8 mm. When the Diviner 2000 was used, it predicted two days after the PRl profile probe predicted date. This difference appeared since the Diviner 2000-measured soil water content at the rooting depth was slightly higher than the PRl profile probe measurements. The amount of irrigation was estimated using two comparable methods (graphic and mathematical method). The amount of irrigation that should have been applied on 20, September (263 day of year) to bring the soil water content to field capacity was estimated to be 4.5 hand 23 mm graphically and 5.23 hand 20 mm mathematically. The difference between these two methods was caused due to the error encountered while plotting the correct line to represent the average variation in soil water content and cumulative irrigation as a function of time. More research is needed to find the cause for the very low soil water content measurements of the PRI profile probe at some depths. The research should be focused on the factors, which could affect the measurement of the PRl profile probe and Diviner 2000 like salinity, temperature, bulk density and electrical conductivity. Further research is also needed to extend the non-linear relationship between the electrical resistance of the sensor and soil water potential up to -200 kPa. This non-linear equation of the Watermark is only applicable within the range of soil water potential between -10 and -100 kPa. / Thesis (M.Sc.)-University of Natal, Pietermaritzburg, 2003.