Agrometeorology
Permanent URI for this communityhttps://hdl.handle.net/10413/6537
Browse
Browsing Agrometeorology by Title
Now showing 1 - 20 of 40
- Results Per Page
- Sort Options
Item Analysis and integration of regional scale temperature datasets into a seasonal crop monitoring system.(2018) Magadzire, Tinomutenda Tamuka.; Savage, Michael John.; Funk, Christopher C.Populations in excess of 20 million people in southern Africa annually face food insecurity. This number increases appreciably when detrimental seasonal climate conditions lead to widespread reductions in crop harvests. This situation has led to the development of regionalscale crop monitoring systems that incorporate crop-specific water balance (CSWB) models for early detection and warning of impending weather-related crop production shortfalls. Early warning of anticipated reductions in crop harvests facilitates early action in responding to potential crises. One such system, used by the Famine Early Warning Systems Network (FEWS NET) in southern Africa for crop monitoring, calculates the water requirements satisfaction index (WRSI) using a CSWB model. Operationally, CSWB models for calculating WRSI have used a static length of growing period (LGP) to bracket the period over which rainfall and evapotranspiration variations can affect crop yields. In the long term, concerns have been raised by some studies on the impact of rising air temperatures on crop production. There is therefore a need to incorporate the impacts of air temperature on crops directly into food security monitoring systems, in order to improve the accuracy of these monitoring systems in identifying and locating weather-related crop production shortfalls. This study sought to assess the potential improvements that can be introduced to the crop monitoring system in general and the WRSI in particular, by incorporating air temperature data into the CSWB model. To address this objective, daily maximum and minimum air temperature grids derived from a general circulation model reanalysis were used to generate thermal time estimates, expressed as growing degree day (GDD) grids for a maize crop. The GDDs were used to estimate the LGP of maize for each pixel of each summer season (which typically runs between around October and March) from 1982/1983 to 2016/2017 in southern Africa. The variable, temperature-driven LGP estimates compared favourably with LGP values obtained from literature for a few sample locations. The variable LGP was used to calculate the WRSI for 35 seasons, and the resultant WRSI showed improved correlation with historical yield estimates compared to the static-LGP WRSI, particularly after the farming practice of planting on multiple dates was taken into consideration. Various expressions of WRSI were considered in the analysis, including WRSI calculated assuming planting at the onset of rains, WRSI aggregated from varying number of separate planting dates, including three and six planting dates as test examples, and WRSI calculated using a modified soil water holding capacity to better capture local soil management practices. Historical maize yields for sub-national administrative units from seven southern African countries were correlated with the various WRSI expressions. Gridded GDD data that reflect the accumulated severity of extreme warm temperatures experienced during the crop growth period, referred to as extreme growing degree days, or eGDD, were also noted to have significant correlations with historical maize yields in several southern African countries. A number of variants of the eGDDs were tested, including eGDDs accumulated throughout the crop’s growth period, eGDDs that only occurred simultaneously with periods of crop water deficit, eGDDs that occurred during the crop flowering stage, and eGDDs scaled by the severity of crop water deficit. In several areas, the various eGDD expressions indicated higher correlations with yield than any of the WRSI variants indicated. The eGDD parameter showed strong correlations with WRSI, suggesting that the accumulated high temperatures were a reflection of the influence of low rainfall and low soil moisture during episodes of high temperature. More work is required to calibrate and refine the temperature-based monitoring parameters that were developed in this study, at local, sub-national scales. In particular, assumptions of the linearity of maize yield response for the various parameters should be tested. Potential improvements of the combined eGDD-water deficit parameter through the incorporation of prediction coefficients and constants should also be tested. A secondary aim of the study was to explore how readily available temperature-related datasets can be utilized to derive air-temperature metrics. To this end, satellite-derived thermal infrared (TIR) brightness temperature data were analysed, and a method was developed for identifying cloud cover, while simultaneously estimating cloud-free diurnal brightness temperature curves, using a single TIR satellite channel. The diurnal brightness temperature curves were developed using a sinusoidal and exponential model for daytime and nighttime respectively, utilizing modifications that enabled the curves to be estimated from two known temperatures at any two given times with cloud-free brightness temperature scenes. Comparison of the cloud mask developed in this study with an existing operational cloud mask based on a methodology developed by the EUMETSAT Satellite Applications Facility for Nowcasting gave an accuracy of 85.4%, when the operational method was considered as truth in a confusion matrix analysis. Situations were identified in which the different cloud detection methods showed superior performance, and could therefore complement each other. A statistical method was also developed for calibrating the cloud-free brightness temperatures to station-observed 2-m air temperatures using relationships between the means and diurnal temperature ranges of the two datasets. This enabled the identification of periods of occurrence of extreme warm air temperatures with a coefficient of determination of 0.91, and demonstrated the potential for the usage of TIR data for generating estimates of useful air temperature metrics. The efficiency of the algorithms that were used for simultaneous cloud masking and generation of cloud-free brightness temperature should be improved, in order to enable the methodology to be scaled up to a regional or global gridded level of analysis. Further work for improving operational gridded air temperature datasets by combining station-observed temperature data, modelled data from global circulation models, satellite-derived modelled cloud-free brightness temperature data and cloud masks is recommended.Item Application of crop models in assessing impact of climate change and identifying adaptation options in a tropical environment in Ethiopia.(2021) Feleke, Hirut Getachew.; Savage, Michael John.; Tesfaye, Kindie.Abstract available in PDF.Item Assessment of maize crop health and water stress based on multispectral and thermal infrared unmanned aerial vehicle phenotyping in smallholder farms.(2021) Brewer, Kiara Raynise.; Clulow, Alistair David.; Sibanda, Mbulisi.; Mabhaudhi, Tafadzwanashe.Abstract available in PDF.Item Assessment of vulnerability of cattle farming to climate variability and change in South Africa.(2021) Masemola, Makgethwa Jillie.; Savage, Michael John.; Tongwane, Mphethe Isaac.Livestock are dependent upon weather for their comfort and food supplies. Sometimes, adverse weather conditions can cause production losses, especially if experienced during critical stages of growth. Heat stress is a major cause of production losses in the dairy and beef industries. Heat stress occurs when the temperature of the environment increases above the comfort zone of cattle as a result of solar irradiance. Heat stress decreases grazing and feed intake in cattle, while drought can limit pasture availability for grazing cattle. The temperature-humidity index (THI), a combination of air temperature and relative humidity, was used to determine the influence of heat stress on the productivity of cattle. The aim of the study was to investigate the air temperature and relative humidity conditions over South Africa accountable for high THI values for cattle farming for the period 1985 to 2015. The standard precipitation index (SPI) and normalised difference vegetation index (NDVI) at three months were computed to assess the soil moisture conditions and vegetation greenness for the season with high THI averages. The THI data analysis was performed seasonally, using a 15-year average and daily values from 75 weather stations in South Africa. Monthly rainfall data from 192 weather stations were used to compute SPI at three months. The NDVI used MODIS satellite information to create vegetation images for the three summer months. Results indicated summer as a season when cattle are vulnerable to heat stress. The periods (2005/06, 2007/08, 2012/13, 2013/14) experienced high seasonal averages (THI > 80) compared to the remaining years. Daily THI extremes were prevalent in February in South Africa. The SPI results indicated that the North West, the western Free State and east of the Eastern Cape provinces were vulnerable to dry conditions for the four summer periods. The NDVI results indicated that the eastern and coastal parts of South Africa were areas of high vegetation activity and greenness for all the summers. The Northern Cape, Limpopo, North West and Free State provinces had vegetation that was vulnerable to dry conditions. Heat stress and a healthy vegetation activity were a problem for the Northern Cape and Limpopo provinces, while KwaZulu-Natal and Mpumalanga had heat stress as their only challenge. The northern and eastern parts of South Africa were heat stress areas with air temperature as the main driver in THI. New smaller areas were developing as heat stress zones, and areas that were identified as heat stress zones in the past were increasing in size for South Africa. Natural grazing areas indicated no significant relationship with THI but showed a relationship with SPI. Cattle farmers situated in the heat stressed regions of South Africa and new heat stress developing areas need to take the necessary precautions to make profitable decisions for their cattle production.Item Comparative evaporation measurements above commercial forestry and sugarcane canopies in the KwaZulu-Natal Midlands.(1999) Burger, Caren.; Savage, Michael John.An understanding of the water use of different crops commonly grown in an area is essential for the implementation of integrated catchment management in South Africa. With increasing pressure on water resources, mainly due to the recent changes in the Water Act, it has become important to determine the actual water demands of agricultural and other crops. Policy makers require knowledge of whether forestry canopies use more water than grassland and other agricultural crops. The Bowen ratio and Penman-Monteith methods were used in a comparative study of the evaporation from Saccharum, Acacia and Eucalyptus. All of the research was conducted at marginal sites located in the KwaZulu-Natal Midlands of South Africa over a period of two years. The Bowen ratio energy balance (BREB) technique combines the Bowen ratio (J3) (the ratio between the sensible, H and latent heat flux density, XE), with the net irradiance (Rn) and soil heat flux densities (G) to calculate evaporation. A comparative study of the sitespecific energy balance components (Rn, G, H and AE), general climatic conditions (rainfall, solar irradiance and air temperature) and other site-specific parameters (leaf area index and average canopy height) was conducted on Saccharum and young commercial forests consisting of Acacia and Eucalyptus. The energy balance highlighted important differences in the energy balance components between the different canopies. The differences between the reflection coefficients at the three sites contributed mainly to the differences in the evaporation rates. The low reflection coefficients of the forest canopies (Acacia and Eucalyptus) (0.1 and 0.08 respectively) were smaller than of the sugarcane canopy (0.2). This resulted in more energy available (« 6 %) for partitioning between the sensible and latent heat flux densities and higher evaporation rates for the forestry canopies. Where low leaf area indices existed (Acacia and Eucalyptus sites) (LAI < 2), the soil heat flux density contributed up to 40 % of the net irradiance (G = 0.4 Rn). The evaporation rates for Saccharum, Acacia and Eucalyptus averaged 2 mm day"1 in winter and 5 mm day"1 in summer. The slightly higher summer evaporation rate for Eucalyptus (5.6 mm day"1), compared to Acacia (4.9 mm day"1), resulted from the lower reflection coefficients and canopy resistance (rc) for Eucalyptus (ocs = 0.08, rc = 35 s m"1) compared to Acacia (ocs = 0.1, rc = 45 s m"1). Automatic weather station data (solar irradiance, air temperature, water vapour pressure and windspeed) were applied to site-specific Penman-Monteith equations to predict evaporation for all three sites. Statistically significant relationships (slope, m « 1, r2 > 0.8) were found between the measured (Bowen ratio) and simulated (site-specific Penman-Monteith) evaporation estimates. The current study has demonstrated the effectiveness of applying the Penman-Monteith equation to forest and sugarcane canopies to predict evaporation, provided accurate net irradiance, soil heat flux densities and canopy resistances are used.Item Development of an air quality early warning system for the eThekwini and Msunduzi municipalities.(2023) Sithole, Siphesihle Fortunate.; Mabhaudhi, Tafadzwanashe.; Botai, Joel.The aim of this study was to assess the source and concentration of criteria air pollutants, as conditioned by the weather and climate conditions, in the eThekwini and Msunduzi Local Municipalities, in the KwaZulu-Natal Province. The outcomes of this assessment will contribute towards the establishment of an appropriate early warning system for the two municipalities, in order to reduce the risk of exposure to air pollution and to avoid the associated health risks. For this purpose, the study adopted various methods and tools to determine the status of air pollution within the municipalities, in terms of the exceedance of particles with an aerodynamic of less, or equal to, 2.5 micrometres (PM2.5) and less than 10 micrometres (PM10) of sulphur dioxide (SO2), ozone (O3) and nitrogen dioxide (NO2). The National Ambient Air Quality Standards (NAAQS) and World Health Organisation (WHO) guidelines were used as a threshold for determining these exceedances. In addition, GIS software was used to determine the spatial distribution of the air pollutants, while the IBM SPSS Statistics Version 27 was used to determine the Pearson correlation between the meteorological parameters and the air pollutant concentration levels within the municipalities. The Openair package in R-Studio was used to determine the wind direction and wind speed within the areas, while the US EPA method was adopted to determine the Air Quality Index (AQI). The results of the exceedance study indicated that most of the air pollutants that were emitted exceeded those of the WHO and NAAQS thresholds, especially PM2.5 and PM10, mostly during the winter period, when temperature inversion is at its strongest. A positive relationship was determined between the air pollutants and meteorological conditions. It is worth noting that the meteorological conditions were not the only factors influencing the concentration levels, but the emissions from industries and vehicles, the burning of fossil fuels for domestic purposes and landfill sites were also observed to be contributing factors to the air pollutant concentration within the municipalities. The results of the proposed AQI showed that most of the people living in eThekwini, especially those living south of Durban, were exposed during winter season to high concentrations of PM10, PM2.5 and SO2, while those of the proposed AQI in the Msunduzi Municipality were, in most cases, observed to be moderate throughout the seasons. In two cases, the Pietermaritzburg CBD and at the Pietermaritzburg Airport, the AQI was recorded as hazardous during spring and summer, respectively. According to the proposed AQI, the most dominant air pollutants are NO2 and SO2. The findings of the study suggest that people living in eThekwini and Msunduzi communities may breathe air harmful to their health and well-being, especially those living near industries and national highways. The overall findings of the study showed that the constitutional right of those living within these municipalities, “to live in an environment that is not iv harmful to their health and wellbeing”, is being violated, since they are being exposed to high air pollutant concentrations.Item Estimating solar radiation for water-use and yield simulations under present and projected future climate using Cropsyst.(2003) Abraha, Michael Ghebrekristos.; Savage, Michael John.Agricultural scientists are faced with the challenge of producing enough food for the increasing world population. Hence the need to develop tools for managing soil and plant systems to increase food production in order to meet the world food demand in the future. Crop simulation models have become promising tools in predicting yield and related components fi'om a set of weather, soil, plant and management data inputs. This study describes the estimation of solar radiant density, a crucial input in crop simulation models; calibration and validation of a soil-plant growth simulator, CropSyst, for management purposes; and generation of weather data for assessment of crop production under possible climate changes in the future. Daily solar radiant density, an input required by most crop simulation models, is infiequently observed in many stations. This may prevent application of crop simulation models for specific locations. Long-term data records of daily minimum and maximum air temperatures, precipitation, sunshine hours and/or solar radiant density were obtained for Cedara, Durban, Seven Oaks and Ukulinga in KwaZulu-Natal, South Africa. Solar radiant density was estimated fi'om sunshine hours using the Angstrom equation and ten other models that involved daily minimum and maximum air temperatures and/or precipitation along with extratelTestrial radiant density. Coefficients for the Angstrom equation and one of the other ten models were specifically developed for South African conditions; the remaining models required fitting coefficients using the available data for all locations. The models were evaluated using (i) conventional statistics that involved, root mean square elTor (RMSE) along with its systematic and unsystematic components, slope, intercept, index of agreement (d), and coefficient of determination (R\ and (ii) a fuzzy expert system that involved a single modular indicator (Ira d) aggregated from the modules of accuracy (aggregation of the indices relative RMSE, model efficiency and I-student probability), con'elation (Pearson's correlation coefficient) and pattem (aggregation of pattem index vs day of year and pattem index vs minimum air temperature). For each index, two functions describing membership to the fuzzy subsets Favourable (F) and Unfavourable (V) were defined. The expelt system calculates the modules according to both the degree of membership and a set of decision rules. Solar radiant density estimated from sunshine hours for the Durban station resulted in R2 , RMSE (MJ m,2) and d index of 0.90, 2.32 and 0.97 respectively. In the absence of observed solar radiant density data, estimations from sunshine hours were used for derivation of coefficients as well as evaluation of the models. For Durban, the performance of the models was generally poor. For Cedara, Seven Oakes and Ukulinga two of the models resulted in a high d index and smallest systematic RMSE. The solar radiant density estimated from each model was also used as an input to simulate maize grain yields using the soil-plant growth simulator, CropSyst. The models were ranked according to their ability to simulate grain yields that match those obtained from using the observed solar radiant density. The rankings according to crop simulation, conventional statistics and expert system were compared. The CropSyst model was also evaluated for its ability to simulate crop water-use of fallow and cropped (oats, Italian ryegrass, rye and maize) plots at Cedara, KwaZulu-Natal, South Africa. Soil characteristics, initial soil water conditions, irrigation and weather data were inputted to CropSyst. Crop input parameters for oats, Italian ryegrass and rye were used, with little modifications, as determined from field experiments conducted at Kromdraai open cast mine, Mpumalanga province, South Africa. Crop input parameters for maIze were either determined fi'om field experiments or taken from CropSyst crop input parameters documentation and adjusted within a narrow specification range of values as dictated by CropSyst. The findings indicated that CropSyst was generally able to simulate reasonably well the water-use of fallow and cropped (oats, Italian ryegI°ass, rye and maize) plots; leaf area index and crop evapotranspiration of rye; and grain yield and developmental stages of maize. The validated CropSyst model was also used to simulate timing and amount of irrigation water, and investigate incipient water stress in oats, Italian ryegrass and rye. The CropSyst model was used to investigate potential effects of future climate changes on the productivity of maize grain yields at Cedara, KwaZulu-Natal, South Africa. The effect of planting date (local planting date, a fortnight earlier and a fortnight later) was also included in the study. A 30-year baseline weather data input series were generated by a stochastic weather generator, ClimGen, using 30 years of observed weather data (l971 to 2000). The generated weather data series was compared with the observed for its distributions of daily rainfall and wet and dry series, monthly total rainfall and its variances, daily and monthly mean and variance of precipitation, minimum and maximum air temperature, and solar radiant density. Four months of the year failed to reproduce distributions of wet and dry series, daily precipitation, and monthly variances of precipitation of the observed weather data series. In addition, Penman-Monteith reference evaporation (ETa) was calculated using the observed and generated data series. Cumulative probability function of ETa calculated using the generated weather data series followed the observed distribution well. Moreover, maize grain yields were simulated using the generated and observed weather data series with local, a fortnight earlier and a fortnight later planting dates. The mean simulated grain yields for the respective planting dates were not statistically different from each other; the grain yields simulated using the generated weather data had significantly smaller variance than the grain yields simulated using the observed weather data series. When the generated weather data series was used an input, the early planting date as compared to the locally practiced and late planting dates resulted in significantly greater simulated grain yields. The grain yields simulated using the observed weather data for the early and local planting dates were not statistically different from each other. The baseline period was modified by synthesized climate projections to create future climatic scenarios. The climate changes considered corresponded to doubling of [C02] from 350 to 700 ~t1 ,-I without air temperature and water regime changes, and doubling of [C02] accompanied by increases in mean air temperature and precipitation changes of 2 (lC and 10%, 2 (le and 20%>, 4 °c and 10%, and 4 (lC and 20% respectively. Solar radiant density was also estimated from daily air temperature range for all scenarios that involved change in mean air temperature. In addition, input crop parameters of radiation-use and biomass transpiration efficiencies were modified for maize, in CropSyst, to accommodate changes in elevated levels of [C02]. Equivalent doubling of [C02], without air temperature or water regime changes, resulted in increased simulated grain yields as compared to the baseline period. Adding 2 QC to the mean daily temperature and 10% to the daily precipitation of a [C02] elevated atmosphere reduced the grain yield but still kept it above the level of the baseline period grain yield. Adding 4 QC to the mean daily temperature and 10% to the daily precipitation fLllther decreased the yield. Increasing the daily precipitation by 20% instead of 10% did not change the simulated grain yield as compared to the 10% increments. Early planting date, for all scenarios, also resulted in higher yields, but the relative increment in grain yield was higher for the late planting dates with scenarios that involved increment in mean air temperature. In general, this study confi1l11ed that doubling of [C02] increases yield but the accompanied increase in mean air temperature reduces yield.Item Estimation of reference evaporation and comparison with ET-gage evaporimeter(2003) Abezghi, Tekeste Weldegabrial.; Savage, Michael John.Accurate estimation of reference evaporation is necessary for the estimation of actual evaporation for irrigation and water resource management purposes. Estimation of reference evaporati~n using the Penman-Monteith method using automatic weather station (AWS)measurements requires the available energy to be accurately estimated. The available energy of short grass of 0.12 m was measured using a component net radiometer and soil heat flux plate measurements at the Faculty of Sciences and Agricultural (Agrometeorological station, University of Natal, Pietermaritzburg, latitude ~29.79 oS, longitude ~ 30.95 °E, altitude ~ 650 m). In an attempt to evaluate the accuracy of commonly used procedures of estimating available energy, estimates of net irradiance (from net long wave irradiance and reflection coefficient estimate) and soil heat flux density were compared to the actual measurements. The linear approximation of atmosphere minus crop surface emittance based on air temperature was compared with measured net long wave irradiance and similar empirical formulations. The underestimation of the measured net long wave irradiance was observed using the linear approach. Furthermore, a plot of measured clear sky net long wave irradiance and air temperature showed a logarithmic relation. The estimated reflected solar irradiance was overestimated for the reference crop. The measured soil heat flux density was observed to vary not only with net irradiance but also with cloudiness, wind speed and soil water content. The soil heat flux density measured with plates was noticed to follow the measured net irradiance. The sensitivity of Penman-Monteith latent heat estimate was investigated for the use of estimated reflection coefficient and soil heat flux density as well as ignored soil heat flux density. Results showed the latent heat estimate to be greater when soil heat flux density was ignored. Reduced set assumptions of Penman-Monteith were assessed usmg the microclimatic measurements. The grass reference evaporation estimate using estimated water vapour pressure from the pervious day minimum air temperature and approximated wind speed were found to be seasonal and procedure dependent. The hourly-reduced set estimate of reference evaporation was in good agreement with the grass Penman-Monteith estimate. The estimated daily water vapour pressure underestimated the daily grass Penman-Monteith estimate. The sensitivity of the reduced set reference evaporation estimate was compared for the two values of approximated wind speeds. The assumption of 2 m S-1 wind speed gave a relatively better result. The sensitivity of the surface temperature energy balance (STEB) estimate of reference evaporation was investigated using two different atmospheric stability procedures. The evaporation estimate agreement and performance of the technique were found to vary depending on the stability correction procedure. The Monteith (1973) correction procedure was observed to be more sensitive to a higher surface-air temperature difference. The Monteith (1973) procedure was found to underestimate the reference evaporation and this resulted in a lower correlation coefficient. The uncorrected and Campbell and Norman (1998) stability corrected procedure of STEB estimate overestimated the reference evaporation but resulted in good agreement with actual reference evaporation. The use of estimated available energy using the STEB method resulted in a 7 % overestimate of measured available energy. Different designs of atmometers have been used to measure evaporation. The less expensive and simple ET_gageR (Model A and E) atmometer for daily evaporation measures were compared to grass-based and alfalfa-based Penman-Monteith and STEB estimate of reference evaporation. Two different evaporation surface covers used with the device allowed for the comparison to be made. Measurements using the canvas 30 ET-gage cover for grass reference evaporation were compared to grass based Penman-Monteith and STEB reference evaporation estimates. Correlation between the canvas 30 measures and Penman-Monteith estimates were good compared to the STEB estimate. The ET-gage canvas 54 measures were in a good agreement with alfalfa based Penman-Monteith reference evaporation estimate. There was, however, a slight time lag in ET-gage evaporation with ET-gage evaporation continuing accumulation when the reference evaporation was zero. The sensitivity of the ET-gage for microclimate variation was tested using the measurements made for two levels and three different microclimates. A shade measurement of reference evaporation was overestimated. The response of the ET-gage to one and two meter microclimate measures was similar to the short grass measurement. Furthermore, the ET-gage surface evaporation estimate using the STEB method showed equal response to the ET-gage surface for the microclimate measure and explained the possible cause of the lag of the ET-gage response. Accurate microclimate measurements is a requirement for the performance of the PenmanMonteith approach for the estimation of reference evaporation. The investment cost required for an AWS set up is high. Alternative options for gathering information of the microc1imate measurements required for calculating reference evaporation were assessed in terms of cost saving, accuracy and other advantages. A weather station system using a Hobo H8 logger (internal relative humidity and air temperature sensor and two external channels, one which was used for solar irradiance measurements) was found to be a cost-effective method for calculating the necessary microclimatic information for calculating reference evaporation. With this system reference evaporation was estimated with reasonable accuracy, at 16 % of the cost of normal AWS system. The use of an Event Hobo logger and an ET-gage was found to provide a reasonable estimate of reference evaporation. The use of the reduced set evaporation weather station was found to be unreasonable in terms of cost and accuracy. Air temperature and relative humidity were measured from different design of radiation shields and Stevenson screens. The use of home-made seven-plate plastic radiation shields provided a similar shield to radiation and ventilation compared to manufactured shields. At a low solar angle when wind speed was very low, all the radiation shields including the small Stevenson screens showed a higher air temperature difference relative to the standard Stevenson screen. The highest average difference of air temperature measurement was measured within the small Stevenson screen and metal-radiation shield. The home-made plastic radiation shield showed similar averages of air temperature and water vapour pressure difference compared to the six- and twelve-plate Gill radiation shields. The home-made metal radiation shield showed relatively higher deviation from the mean being cold at night time and hot during the day. More research is needed to explore the efficiency of the ET-gage evaporation from variety of microclimates to establish the cause of the overestimate under shade, to develop better relation of clear day net long wave irradiance and air temperature and the use of a wind speed sensor with Hobo H8 weather station system.Item Frequency domain reflectometry for irrigation scheduling of cover crops.(2003) Gebregiorgis, Mussie Fessehaye.; Savage, Michael John.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.Item Heat and energy exchange above different surfaces using surface renewal.(2008) Mengistu, Michael Ghebrekidan.; Savage, Michael John.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.Item Investigating possible impact of climate change on sugarcane production in KwaZulu-Natal, South Africa.(2016) Sithole, Bonga Benson.; Savage, Michael John.The KwaZulu-Natal Department of Agriculture, Environmental Affairs and Rural Development in 2010 undertook to investigate the vulnerability of KwaZulu-Natal to climate change which identified various sectors of the economy in the region that are impacted by climate change. The aim of the current study was to investigate the possible impact of climate change on sugarcane production in KwaZulu-Natal. The main objective of the study was to identify and synthesize current knowledge, scientific literature and data relating to specific aspects of climate change in KwaZulu-Natal, South Africa. In order to achieve the objectives of the study, a questionnaire was developed to ascertain sugarcane farmers’ awareness about climate change. Based on questionnaire data, the study reveals that sugarcane farmers are aware of the effect of climate change on sugarcane production. Cane growers are also aware of their activities that contribute to climate change. Climatological data for the region were collected and analyzed through the Decision Support System Agro-tech Transfer (DSSAT) daily crop model to assess possible climate change impacts on sugarcane production. Daily rainfall, solar radiation, wind-speed, air temperature (minimum and maximum) and dew point temperature data for various sugarcane mill supply areas were collected from the South African Sugar Research Institute. RClimDex (1.0) software was used to determine if climate change did occur for the period 1966 to 2016. The climate data sets were positive to climate change with respect to daily maximum air temperature, daily minimum air temperature, daily rainfall, daily dew point temperature and daily solar radiation. The impact of climate change on sugarcane production in KwaZulu-Natal has been recognized as the main cause for yield reduction. No major decline in sugarcane production has been noted in KwaZulu-Natal for those farmers practicing irrigation and improved management. The study also demonstrates an increase in the amount of trash on sugarcane in the latter years of the study period. In general, the approach presented in this study encompassed and assessed the effect of climate change on sugarcane production with inclusion farmers perception can be considered as a strategic issue on existing climate change concern in KwaZulu-Natal province. Further research about what measures have been implemented by other countries in addressing climate change is recommended.Item Investigation of microclimate for human comfort in the natural environment and in a car parked in the open.(2016) Luthuli, Sithandiwe Ignatia.; Savage, Michael John.Abstract available in PDF file.Item Irrigation control system with a web-based interface for the management of Eucalyptus planting stock.(2015) Kaptein, Nkosinathi David.; Savage, Michael John.Commercial forestry nurseries use large quantities of water to irrigate the planting stock to meet the annual forestry industry planting demands. However, South Africa is a water scarce country and there is high competition for this limited resource with other sectors. Thus, sustainable water management strategies should be put in place to preserve this precious resource. In commercial forestry nurseries, sustainable water use may be achieved by carefully managing irrigation schedules, such as accurately measuring growing media water content and then replenishing the depleted soil water. Improved nursery irrigation management may not only save water, but also has the potential to reduce the prevalence of pests and diseases, excessive leaching of nutrients, irrigation costs and may produce robust planting stock that is better suited to adapt to field conditions. Growing media water content can be directly measured by the gravimetric method. However, this method is laborious, time consuming, costly and does not allow for near real-time monitoring and control. Several indirect methods for estimating soil water content have been developed and are in use today. Some common methods are: frequent domain reflectometry, time domain reflectometry, time domain transmission and the dual-needle heat pulse method. However, each method has advantages and disadvantages. Considerations for choosing the most appropriate method are the ability to automate, accuracy and precision, skills required to use and the costs to purchase. The main objective of this study was to calibrate the low cost commercially available Decagon EC-5 soil water content sensors using nursery growing media. The calibration equation was then used to measure and control irrigation for Eucalyptus grandis x Eucalyptus urophylla and Eucalyptus dunnii planting stock grown in seedling trays in a greenhouse. A web-based data and information system was utilised to share measurements and display greenhouse environmental conditions in near real-time. The data could be viewed or downloaded using the internet1. Decagon EC-5 soil water content sensors were laboratory-calibrated, using nursery growing media, against the standard gravimetric method. Four nursery growing media were used for calibrations: coir/perlite mix (CP), coir/pine bark/vermiculite mix (CPBV), pine bark (PB) and sandy soil. The calibration relationships between gravimetric water content and sensor output for each growing media were established, and the manufacturer supplied calibration equation was evaluated against laboratory calibration equations. The appropriate laboratory calibration equation was programmed in the datalogger to measure and control irrigation. Greenhouse microclimate measurements of air temperature, relative humidity and solar irradiance were conducted. Hourly grass reference evaporation (ETo) was calculated by the datalogger. The greenhouse microclimate measurements were compared against an open area automatic weather station at the University of KwaZulu-Natal Agrometeorology Instrumentation Mast (UKZN AIM) system measurements. The EC-5 soil water content sensors were used to schedule irrigation for E. grandis x E. urophylla hybrid clones (GxU) (Experiment 1) and E. dunnii seedlings (Experiment 2) grown in seedling trays inside a fully air temperature controlled greenhouse. Irrigation was controlled in three treatments: low, medium and high watering. All treatments were treated the same except for differing irrigation application. Total daily irrigation and drainage were measured per treatment. Irrigation scheduling for GxU was programmed at a set point and E. dunnii seedlings at lower and drained upper limits. Seedling measurements conducted were root collar diameter (RCD), heights, stomatal conductance (gs), root-to-shoot ratio and total leaf area. Total drainage and its electrical conductivity (EC) was also measured. The calibration relationship showed a linear relationship between gravimetric water content and sensor output for all four growing media with an R2 greater than 0.92. The manufacturer supplied calibration equation poorly estimated growing media water content compared to the laboratory calibration. Poor estimation exceeded the 5% error specified by the manufacturer. Air temperature was consistently less than 25°C inside the greenhouse. The external air temperature, as measured by the UKZN AIM system, fluctuated and reached a maximum of 36.6°C during the study period. Solar irradiance inside the greenhouse was 60% lower than that measured by the UKZN AIM system. The relative humidity was consistently higher inside the greenhouse compared to that measured by the UKZN AIM system. Greenhouse grass reference evaporation was consistently lower than the UKZN AIM system due to low air temperature and high RH inside the greenhouse. For Experiment 1, the GxU clones were irrigated too frequently for short periods. This led to over- and under- irrigation in high and low watering treatments, respectively. However, these challenges were addressed in Experiment 2 using E. dunnii seedlings irrigated at lower and drained upper irrigation limits. In Experiment 2, variability in sensor measurements within each treatment were observed at drained upper limit and decreased at lower limit. This was likely caused by a change in the pore space volume from dry to wet growing media. The web-based system was successfully used as an early warning system to monitor soil water content measurements and greenhouse microclimate, averting experimental failure due to lack of irrigation on one occasion. Seedlings in the high watering treatment had the highest RCD, heights and gs followed by the medium and low watering treatments. Although the low watering treatment had the lowest growth rates and gs these seedlings were more robust, hardy and resistant to water stress. The root-to-shoot ratio showed no statistically significant differences between treatments. However, seedlings in the high watering treatment had slightly greater root volume. This was probably due to the increased total seedling leaf area for this treatment which facilitated increased photosynthetic activity and carbohydrates production, enabling increased root growth. The highest EC measurements were recorded in the low watering treatment. This was likely due to low irrigation and therefore nutrients were not washed off the growing media. Medium watering treatment EC was 30% lower than the low watering treatment whilst high watering EC was almost equivalent to the irrigation water. The analysis of economics showed that implementing the fully automated system could be costly. However, there are many potential benefits that may be offered by this system such as reduction in water use, pumping costs and management time. The early warning offered by this system could potentially help avoid the loss of planting stock if there is a problem with the irrigation system. The study showed that irrigation may be automatically scheduled for nursery seedlings trays using low cost Decagon EC-5 soil water content sensors with reasonable accuracy. However, medium-specific calibration is important to improve the soil water content measurement accuracy. The study also showed that reducing irrigation may result in reduced growth rates of seedlings. However, other benefits such as seedling resistance to water stress, robust seedlings, irrigation water savings and a reduction in washing off nutrients may be achieved.Item The long-term measurement of total evaporation over Acacia mearnsii using large aperture scintillometry(2007) Clulow, Alistair David.; Savage, Michael John.A large aperture scintillometer (LAS) was operated continuously over a distance of 575 m from 19 August 2006 to 29 September 2007 in the South African KwaZulu-Natal midlands mistbelt area over Wattle (Acacia mearnsii). The LAS measurements of the structure parameter of the refractive index of air ( 2 n C ), were used to calculate the sensible heat flux. The shortened energy balance equation was used to estimate the latent energy flux as a residual from which the total evaporation (ET) was calculated. The LAS estimates of sensible heat flux during the short transition period (1 hour) between stable and unstable conditions were on occasion erroneous and required verification. Advection was also found to affect sensible heat flux estimates. Long-term operation of the LAS was however found to be possible even at remote sites, producing reliable and continuous results. The LAS estimates of sensible heat are sensitive to zero-plane displacement height and wind speed data impact and these should be derived as accurately as possible. Tree heights were measured at monthly intervals and a zero-plane displacement and effective height were calculated every two weeks. The sensible heat flux was thus processed in two week blocks of data corresponding to progressive effective heights. The tree growth rate was consistent over time and was not affected by seasonality, indicating that reduced air temperatures, rainfall and solar irradiance in winter are not limiting growth. The average growth rate was 0.37 m per month or 4.5 m per year. The LAS ET was compared to the American Society of Civil Engineers - Environmental and Water Resources Institute (ASCE-EWRI) short grass reference evaporation (ETsz) for a seven-month period and was found to compare favourably (R2 = 0.78) with outliers caused by advection and rainfall events. Calculations of grass reference evaporation at hourly and daily intervals provide different results. The daily estimates are lower than the hourly estimates by 17 % on average. Where hourly data is summed to calculate a daily ETsz, night-time values should be included. The LAS ET measurements were validated against the Priestley and Taylor (1972) method of estimating ET and found to be in good agreement (R2=0.94). The Priestley and Taylor daily total latent energy flux, from 22 August 2006 to 29 September 2007, was 9 % higher than the LAS results on average. The Bowen ratio for the entire period is less than 1, indicating that the latent energy flux dominates at the site. The ET over the period of measurement (13 months) is 1250 mm and the rainfall is 750 mm. This confirms previous results at the site using the Bowen ratio energy balance method showing that the ET exceeds the rainfall by 45 % and justifies further research into soil water, ground water and root interactions in the deep soil profile.Item Long-term measurements of spatially-averaged sensible heat flux for a mixed grassland community, using surface layer scintillometry.(2007) Odhiambo, George O.; Savage, Michael John.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.Item Measurement of water potential using thermocouple hygrometers.(1982) Savage, Michael John.; Crass, Albert.; De Jager, James M.Theory predicts that the time dependent voltage curve of a thermocouple psychrometer where there is no change in output voltage with time during the evaporation cycle defines the wet bulb temperature T[w] corresponding to the water potential. In practice, a change in voltage with time does occur and it is convenient to define the voltage corresponding to the water potential as the maximum point of- inflection voltage. A predictive model based on calibration data at a few tempertures is used to obtain the psychrometer calibration slope at any temperature. Use of this model indicates that psychrometers differ from each other and therefore must be individually calibrated if accuracy better than ±5 % in the measurement of water potential is required. Dewpoint hygrometers are shown to be less temperature sensitive than psychrometers and have the added advantage of a voltage sensitivity nearly twice that of psychrometers, typically -7,0 x 10¯³ μV/kPa compared to -3,7 x 10¯³ μV/kPa at 25 °C. The accurate temperature correction of hygrometer calibration curve slopes is a necessity if field measurements are undertaken using either psychrometric or dewpoint techniques. In the case of thermocouple psychrometers, two temperature correction models are proposed, each based on measurement of the thermojunction radius and calculation of the theoretical voltage sensitivity to changes in water potential. The first model relies on calibration at a single temperature and the second at two temperatures. Both these models were more accurate than the temperature correction models currently in use for four leaf psychrometers calibrated over a range of temperatures (15 to 38°C). The model based on calibration at two temperatures is superior to that based on only one calibration. The model proposed for dewpoint hygrometers is similar to that for psychrometers. It is based on the theoretical voltage sensitivity to changes in water potential. Comparison with empirical data from three dewpoint hygrometers calibrated at four different temperatures indicates that these instruments need only be calibrated at, say 25°C, if the calibration slopes are corrected for temperature. A model is presented for the calculation of the error in measured thermocouple hygrometric water potential for individual hygrometers used in the dewpoint or psychrometric mode. The model is based on calculation of the relative standard error in measured thermocouple psychrometric water potential as a function of temperature. Sources of error in the psychrometric mode were in calibration of the instrument as a function of water potential and temperature and in voltage (due to electronic noise and zero offsets) and temperature measurement in the field. Total error increased as temperature decreased, approaching a value usually determined by the shape of the thermocouple junction, electronic noise (at low voltages less than 1 μV) and errors in temperature measurement. At higher temperatures, error was a combination of calibration errors, electronic noise and zero offset voltage. Field calibration data for a number of leaf psychrometers contained total errors that ranged between 6 (at a °C) and 2 %(at 45 °C) for the better psychrometers and between 11 (at 0° C) and 5 % (at 45 C) for the worst assuming that the zero offset was 0,5 μV. Zero offset values were less than 0,7 μV at all times. The dewpoint errors arose from calibration of the dewpoint hygrometer as a function of water potential, extrapolation of the calibration slope to other temperatures, setting the dewpoint coefficient and errors in voltage and temperature measurement. The total error also increased as temperature decreased, because of the differences in temperature sensitivity between dewpoint and psychrometric calibration constants. Consequently, the major source of error in the dewpoint mode arose from the difficulty in determining the dewpoint coefficient. This error, which is temperature dependent, contains three subcomponent errors; the temperature dependence, random variation associated with determining the temperature dependence and error in setting the correct value. Calibration and extrapolation errors were smaller than those of the psychrometric technique. Typically, the error in a dewpoint measurement varied between about 6 and 2 % for the best hygrometer and between 10 and 3 % for the worst for temperatures between 0 and 45 °C respectively. At low temperatures, the dewpoint technique often has no advantage over the psychrometric technique, in terms of measurement errors. In a comparative laboratory study, leaf water potentials were measured using the Scholander pressure chamber, psychrometers and hydraulic press. Newly mature trifoliates cut from field grown soybean (Glycine max (L) Merr. cv. Dribi) were turgidified and, after different degrees of dehydration, leaf water potential measured. One leaflet from the trifoliate was used for the thermocouple psychrometer and another for the press while the central leaflet with its petiolule was retained for use in the pressure chamber. Significant correlations between measurements using these instruments were obtained but the slopes for hydraulic press vs psychrometer measurement curve and hydraulic press vs pressure chamber were 0,742 and 0,775 respectively. Plots of pressure-volume curves indicate that the point of incipient plasmolysis was the same (statistically) for the thermocouple psychrometer and the pressure chamber, but much larger for the hydnaulic press. The above-mentioned differences between the three instruments emphasize the need for calib rating the endpoint defined us i ng the press against one or more of the standard techniques, and, limi ting the use of the press to one person. Cuticular resistance to water vapour diffusion between the substomatal cavity and the sensing psychrometer junction is a problem unique to leaf psychrometry and dewpoint hygrometry; this resistance is not encountered in soil or solution psychrometry. The cuticular resistance may introduce error in the leaf water potential measurement. The effect of abraiding the cuticle of Citrus jambhiri to reduce its resistance, on the measured leaf water potential was investigated. Psychrometric measurements of leaf water potential were compared with simultaneous measurements on nearby leaves using the Scholander pressure chamber, in a field situation. Leaf surface damage, due to abrasion, was investigated using scanning electron microscopy. Thermocouple psychrometers are the only instruments which can measure the in situ water potential of intact leaves, and which may be suitable for continuous, non-destructive monitoring of water potential. Unfortunately, their usefulness is limited by a number of difficulties, among them fluctuating temperatures and temperature gradients within the psychrometer, sealing of the psychrometer chamber to the leaf, shading of the leaf by the psychrometer and resistance to water vapour diffusion by the cuticle when the stomates are closed. Using Citrus jambhiri, several psychrometer designs and operational modifications were tested. In situ psychrometric measurements compared favourably with simultaneous Scholander pressure chamber measurements on neighbouring leaves, corrected for the osmotic potential and the apparent effect of "xylem tension relaxation" following petiole excision. It is generally assumed that enclosure of a leaf by an in situ thermocouple psychrometer substantially modifies the leaf environment, possibly altering leaf water potential, the quantity to be measured. Furthermore, the time response of leaf psychrometers to sudden leaf water potential changes has not been tested under field conditions. In a laboratory investigation, we found good linear correlation between in situ leaf psychrometer (sealed over abraided area) and Scholander pressure chamber measurements (using adjacent leaves) of leaf water potential, 2 to 200 minutes after excision of citrus leaves. A field investigation involved psychrometric measurement prior to petiole excision, and 1 min after excision, simultaneous pressure chamber measurements on adjacent citrus leaves immediately prior to the time of excision and then on the psychrometer leaf about 2 min after excision. Statistical comparisons indicated that within the first two minutes after excision, psychrometer measurements compared favourably with pressure chamber measurements. There was no evidence for a psychrometer leaf water potential time lag. For the high evaporative demand conditions, water potential decreased after excision by as much as 700 kPa in the first minute. Psychrometer field measurements indicated that within the first 5 min of leaf petiole excision, the decrease in leaf water potential with time was linear but that within the first 15 s, there was a temporary increase of the order of a few tens of kilopascal. The thermocouple psychrometer can be used to measure dynamic changes in leaf water potential non-destructively, with an accuracy that compares favourably with that of the pressure chamber. Using in situ thermocouple leaf hygrometers (dewpoint and psychrometric techniques employed) attached to Citrus jambhiri leaves, an increase in measured water potential immediately following petiole excision was observed. The increase ranged between 20 to 80 kPa and occurred 30 s after petiole excision and 100 s after midrib excisions. No relationship between the actual leaf water potential and the increase in water potential due to excision, was found.Item Microclimate modification to improve productivity of ‘Carmen®-Hass’ avocado orchards using shadenet under subtropical conditions of Limpopo Province, South Africa.(2016) Malapana, Clarance Kgethego.; Savage, Michael John.The agricultural environment is a complex and dynamic system. Microclimate, the crop, biosphere, and management practices interact to determine the best yield production. South Africa is a water-scarce country, with high variability in annual rainfall. Thus, water quality and quantity are major limiting factors in agriculture. Hence, shadenetting can be used to modify the orchard microclimate to make the environment more conducive for fruit production. The South African avocado industry is export-oriented, so there is a commercial need to optimise the exportable percentage of avocado fruit. Sunburn, wind and hail damage and small fruit size as a result of water stress are the major cull factors for the industry. It is believed that shadenetting, with changes in management practices, can counter these limiting factors. There is no literature on growing avocado fruit under shadenetting. Therefore, the aim of the research was to determine the effects of a 20% white shadenet on ‘Carmen®-Hass’ avocado orchards and productivity. The long term objective is to improve avocado fruit quality and profitability in the Mooketsi Valley, Limpopo province, South Africa, a subtropical environment by reducing abiotic stress, particularly, solar irradiance, heat and wind. The trial was conducted at Goedgelegen Estate in the Mooketsi Valley on ‘Carmen®-Hass’ trees planted in 2007/8 season. A 1-ha shadenet structure (6 m high) was used, with 20% white shadenet over the roof and 40% green shadenet on the sides. Air and canopy temperature, relative humidity, wind speed, solar irradiance and leaf wetness duration (LWD) and sap flow were monitored at a sub-hourly rate. Evapotranspiration was calculated from the above mentioned parameters. Irrigation was monitored five times per week using tensiometers at 300- and 600-mm soil depths. The comparison between open and shadenet leaf areas showed that leaves in the open treatment were reduced as a result of the abiotic stress. Fruit water content under the shadenet compared to the open was greater, such enabled fruit under the shadenet to reach maturity two weeks earlier when compared with open treatment. Air and canopy temperature and relative humidity were slightly reduced under the shadenet, with the greatest difference occurring during the flowering period in mid-winter. The modification in air temperature and relative humidity was beneficial for bee activity and pollination in 2015 compared to the 2014 season. ‘Carmen®-Hass’ flowers in mid-winter when temperature conditions are not conducive for pollination. Canopy temperature was also reduced under shadenetting compared to the open treatment. The reduction was due to differences in tree density and the role that shadenetting plays. The infrared thermometer measurements were uniform with dense canopies compared to sparse tree canopies. The midday incoming solar irradiance was reduced by 18% under the shadenet compared to the open treatment. Calm conditions were experienced under the shadenet. Hence, windspeed was reduced to negligible levels. Also, the shadenet resist air flow to a certain height compared to the open treatment. LWD was extended by 12% under shadenet. An infestation of the insect pest citrus leaf roller (Archips occidentalis) caused severe damage to the fruit during the 2014/15 season due to the high plant density used. Significant results were that evapotranspiration was reduced by 14 and 29% less water was applied under the shadenet to maintain an adequate soil water content compared to the open treatment. Fruit reached minimum maturity two weeks earlier under shadenet compared to the open treatment. Fruit quality and pack-out were improved under the shadenet due to reduction in sunburn, wind damage and small fruit. But poor yields were experienced during the 2014 season due to poor bee activity, pollination and fruit size distribution were reduced under the shadenet compared to the open treatment. But following the improved bee activity in 2015, the 2016 normal season yield is likely to be improved under the shadenet than in the open treatment. Data collected in the Mooketsi Valley showed that 20% white shadenet has modified the microclimate and improved fruit quality. The water use under the shadenet was improved compared to the open treatment. But a thorough investigation on bee management under shadenet is required to optimise pollination in order to obtain greater yields under the shadenet.Item Modeling variety differences in canopy growth and development of sugarcane (Saccharum officinarum L.) using Canegro.(2003) Zhou, Marvellous Mabeza.; Savage, Michael John.Crop models have great potential as research tools, for crop system management and policy analysis. One of the most promising future uses of crop models is in crop improvement. The limitation in the use of models for crop improvement has been the inability of crop models to predict variety differences. Currently, the CANEGRO model, a sugarcane crop model developed the South African Sugar Association Experiment Station (Inman-Bamber, 1991a) can only model the performance of the NC0376 variety. Experiments were undertaken in the South East Lowveld of Zimbabwe, which is a hot and dry environment where sugarcane is grown under irrigation, to examine the canopy growth and development of four commercial varieties, ZN6, ZN7, N14 and NC0376. The study aimed at determining variety differences in canopy (tillers and leaves) development, develop parameters that can be used to model variety differences and test the improved CANEGRO canopy model for its ability to predict variety differences in canopy growth and development. For the late season, the numbers of leaves and tillers produced by each variety were counted every fortnight throughout the crop cycle. The total leaf area of the varieties and the individual leaf area on a stalk were determined using a Delta-T leaf area meter every fortnight. The date of emergence of successive leaves on a stalk was recorded daily. The leaf angles of each variety were measured every fortnight. The amount of photosynthetically active radiation (PAR) intercepted by the varieties was measured using a SunScan Ceptometer. Tillering and tiller senescence rates, phyllochron intervals, extinction coefficients and base temperatures were determined for the growth and development processes of varieties ZN6, ZN7, N14 and NC0376. Tiller and leaf population development was varietal. Tillering and leaf emergence were highly correlated to thermal time while tiller and leaf senescence were less correlated to thermal time. The poor correlation of the senescence phases to thermal time could mean that tiller and leaf senescence was driven by other factors other than thermal time. PAR interception could be one of these factors. The data showed that PAR interception could be a trigger of tiller senescence. The study showed that the tiller and leaf population development could be approximated by two linear equations. Tillering will be the first linear phase and tiller senescence the second linear phase. The first linear phase is driven by thermal time. While the second linear phase is triggered by PAR interception, the major driving factors need to be determined. This study proposed the use of two linear equations to model tiller and leaf population development as opposed to the polynomial equations used in the current CANEGRO model. Polynomial equations assume the factors driving tillering and tiller senescence are the same. The green leaf numbers per plant showed that all varieties experienced a decline in green leaf numbers with crop age. Varieties NC0376 and ZN7 had the greatest decline in green leaf numbers per plant while varieties N14 and ZN6 had the least decline. Variety ZN7 had the highest number of green leaves per plant while NC0376 had the least. The tiller growth and development was divided into three phases: the exponential phase during the initiation of stalks, the first linear phase during a period of rapid stalk elongation and the second linear phase during sucrose accumulation and maturation. The first two phases of development were driven by thermal time while the sucrose accumulation was not. There were variety differences in tiller growth and development. There were variety differences in base temperature for the development of various components of the canopy. Internode formation occurred at lower air temperatures than stalk elongation and tillering while canopy heights were correlated with higher air temperatures. This implies that internode formation could occur under conditions unsuitable for stalk elongation and may explain the short internodes frequently observed in stalks exposed to winter during rapid stalk elongation. The basic requirements for physiological parameters are that they should be stable across different environments, have significant differences between varieties and have physiological meaning. The parameters studied were thermal time requirement for shoot emergence, leaf appearance, to reach peak tiller population and to start of stalk elongation; surface area of the youngest biggest leaf, leaf number of the youngest biggest leaf, PAR transmission at the start of tiller senescence, extinction coefficients, and peak and mature tiller population. The difference between varieties in thermal time to shoot emergence was least using a base temperature of 16 QC compared to using 10 QC and therefore 16 QC could be a more appropriate base temperature for shoot emergence. The accumulated soil temperatures were less variable than accumulated air temperature and could therefore be a more reliable driver of shoot emergence. However, the limitations in the use of soil temperature are that it is not a readily available measurement and that it is not easy to measure. The gradual increase in phyllochron intervals appeared to be a better method of predicting leaf appearance compared to using a broken stick model. The phyllochron gradient was proposed, as it is likely to be a more robust way of modelling leaf appearance. The varieties had different phyllochron gradients. Variety ZN7 had highest rate of leaf appearance and produced the highest number of leaves per stalk while NC0376 had the lowest rate and produced the least number of leaves. There were statistically significant differences between varieties (P = 0,05) in PAR transmission at the start of tiller senescence and a base temperature of 16 QC was best at determining accumulated thermal time to the start of tiller senescence. Varieties with higher peak tiller population had higher final tiller population, lower thermal time per tiller and a higher ratio of final to peak tiller population. There were differences between varieties in the youngest leaf number attaining maximum leaf area and the leaf area of the youngest biggest leaf. Variety Nl4 had the biggest leaves and NC0376 had the smallest. Variety Nl4 had the highest leaf area index (LAI) while ZN7 had the lowest. There were significant differences (P = 0,01) in PAR intercepted by the varieties but there were no significant differences in extinction coefficients. Extinction coefficients increased with crop age. The varieties had significantly different (P = 0,01) leaf angles and ellipsoidal leaf angle distribution parameters. The measurement of LAI using SunScan ceptometer provided a better estimate of extinction coefficients than LAI measured using Delta-T leaf area meter. Model evaluation showed that CANEGRO canopy model version 2 was improved compared to than version 1. The model (version 2) was accurate in predicting tiller heights and dead leaf numbers per stalk. It was fairly accurate in predicting green leaf numbers per plant, stalk population and intercepted PAR but was poor in predicting LA!. Version 2 has proved to be a substantial improvement over version 1 in predicting stalk population. Generally, the version 2 model overestimated tiller heights early and underestimated later, overestimated the tiller population and LAI after peak, underestimated green leaf numbers per stalk for varieties ZN6, ZN7 and N14 and overestimated dead leaf numbers per stalk and intercepted PAR. The version 2 model predicted a constant green leaf numbers per plant and LAI from peak to harvest while observed data showed that green leaf numbers per stalk and LAI decreased towards harvest. Version 2 model predicted the tiller population of NC0376 closely but underestimated tiller senescence in N14 and also underestimated final tiller population in varieties ZN6 and ZN7. Future model refinements may need to focus on the prediction of the sigmoid pattern of tiller heights. The model may need to be calibrated to predict the green leaf numbers per stalk accurately, which should possibly improve the prediction of LAI that in turn could improve the prediction of intercepted PAR. The improvement in the timing and rate of tiller senescence should improve the prediction of tiller population particularly in varieties ZN6, ZN7 and N14. The study showed that the broken stick method IS superior in explaining leaf and tiller population development compared to using polynomial equations. The development of variety parameters helped improve the prediction of variety differences in canopy growth and development. A major weakness of most crop models is modelling variety differences in canopy growth and development. The inability of crop models to predict variety differences has limited their use in plant breeding. This study has resulted in an improved version of CANEGRO version 1 that is an initial attempt at modelling variety differences of sugarcane.Item Modelling the impacts of increased air temperature on maize yields in selected areas of the South African highveld using the cropsyst model.(2014) Pasi, Jonathan M.; Savage, Michael John.Abstract available in PDF file.Item Modelling the soil water balance and applications using a decision support system (DSSAT v3.5).(2003) Ghebreab, Tesfalidet Alem.; Savage, Michael John.Water is a scarce resource used by various stakeholders. Agriculture is one of the users of this resource especially for growing plants. Plants need to take up carbon dioxide to prepare their own food. For this purpose plants have stomatal openings. These same openings are used for transpiration. Quantifying transpiration is important for efficient water resource management and crop production because it is closely related to dry matter production. Transpiration could be measured using a number of methods or calculated indirectly through quantification of the soil water balance components using environmental instruments. The use of models such as the Decision Support System for Agrotechnology Transfer (DSSAT v3.5) is, however, much easier than environmental instruments. Nowadays, with increased capabilities of computers, the use of crop simulation modelling has become a common practice for various applications. But it is important that models, such as DSSAT v3.5, be calibrated and verified before being used for various applications such as long-term risk assessment, evaluation of cultural practices and other applications. In this study the model inputs have been collected first Then the model was calibrated and verified. Next sensitivitY analysis was carried to observe the model behavior to changes in inputs. Finally the model has been applied for long-term risk assessment and evaluation of cultural practices. In this study, the data collected formed the basis forthe minimum dataset needed for running the DSSAT v3.5 model. In addition, the factory given transmission of shading material over a tomato crop was compared to actual measurements. Missing weather data (solar irradiance, minimum and maximum air temperature and rainfall) were completed after checking that it was homogeneous to measurements from nearby automatic weather station. It was found that factory-given transmission value of 0.7 of the shade cloth was different from the actual one of 0.765. So this value was used for conversion of solar irradiance measured outside the shade cloth to solar irradiance inside the shade cloth. Conventional laboratory procedures were used for the analysis of soil physical and chemical properties. Soil water content limits were determined using texture and bulk density regression based equations. Other model inputs were calculated using the DSSAT model. Crop management inputs were also documented for creation of the experimental details file. The DSSATv3.5 soil water balance model was calibrated for soil, plant and weather conditions at Ukulinga by modifying some of its inputs and then simulations of the soil water balance components were evaluated against actual measurements. For this purpose half of the data available was used for calibration and the other half for verification. Model simulations of soil water content (150 to 300 mm and 450 to 600 mm) improved significantly after calibration. In addition, simulations of leaf area index (LA!) were satisfactory. Simulated evapotranspiration (E1) had certain deviations from the measured ET because the latter calculated ET by multiplying the potential ET with constant crop multiplier so-called the crop coefficient. Sensitivity analysis and long-term risk assessments for yield, runoff and drainage and other model outputs were carried out for soil, plant and weather conditions at Ukulinga. For this purpose, some of the input parameters were varied individually to determine the effect on seven model output parameters. In addition, long-term weather data was used to simulate yield, biomass at harvest, runoff and drainage for various initial soil water content values. The sensitivity analysis gave results that conform to the current understanding of the soil-plant atmosphere system. The long-term assessment showed that it is risky to grow tomatoes during the winter season at Ukulinga irrespective of the initial soil water content unless certain measures are taken such as the use of mulching to protect the plants from frost. The CROPGRO-Soya bean model was used to evaluate the soil water balance and gro'W1:h routines for soil, plant and weather conditions at Cedara. In addition, cultural practices such as row spacing, seeding rate and cultivars were also evaluated using longterm weather data. Simulations of soil water content were unsatisfactory even after calibration of some of the model parameters. Other model parameters such as LAI, yield and flowering date had satisfactory agreement with observed values. Results from this study suggest that the model is sensitive to weather and cultural practices such as seeding rates, row spacing and cultivar maturity groups. The general use of decision support systems is limited by various factors. Some of the factors are: unclear definition of clients/end users; no end user input prior to or during the development of the DSS; DSS does not solve the problems that the client is experiencing; DSS do not match their decision-making style; producers see no reason to change the current management practices; DSS does not provide benefit over current decision-making system; limited computer ownership amongst producers; lack of field testing; producers do not trust the output due to the lack of understanding of the underlying theories of the models utilized; cannot access the necessary data inputs; lack of technical support; lack of training in the development ofDSS software; marketing and support constraints; institutional resistances; short shelf-life of DSS software; technical constraints, user constraints and other constraints. For successful use of DSS, the abovementioned constraints have to be solved before their useful impacts on farming systems could be realized. This study has shown that the DSSAT v3.5 model simulations of the soil water balance components such as evapotranspiration and soil water content were unsatisfactory while simulations of plant parameters such as leaf area index, yield and phonological stages were simulate to a satisfactory standard. Sensitivity analysis gave results that conform to the current understanding of the soil-plant -atmosphere system. Model outputs such as yield and phonological stages were found to sensitive to weather and cultural practices such as seeding rates, row spacing and cultivar maturity groups. It ha been further investigated that the model could be used for risk assessment in various crop management practices and evaluation of cultural practices. However, before farmers can use DSSAT v3.5, several constraints have to be solved.