Impacts of improved grazing land management on sediment yields, Part 1: Hillslope processes.

Summary Poor land condition resulting from unsustainable grazing practices can reduce enterprise profitability and increase water, sediment and associated nutrient runoff from properties and catchments. This paper presents the results of a 6 year field study that used a series of hillslope flume experiments to evaluate the impact of improved grazing land management (GLM) on hillslope runoff and sediment yields. The study was carried out on a commercial grazing property in a catchment draining to the Burdekin River in northern Australia. During this study average ground cover on hillslopes increased from ~35% to ~75%, although average biomass and litter levels are still relatively low for this landscape type (~60 increasing to 1100 kg of dry matter per hectare). Pasture recovery was greatest on the upper and middle parts of hillslopes. Areas that did not respond to the improved grazing management had <10% cover and were on the lower slopes associated with the location of sodic soil and the initiation of gullies. Comparison of ground cover changes and soil conditions with adjacent properties suggest that grazing management, and not just improved rainfall conditions, were responsible for the improvements in ground cover in this study. The ground cover improvements resulted in progressively lower runoff coefficients for the first event in each wet season, however, runoff coefficients were not reduced at the annual time scale. The hillslope annual sediment yields declined by ~70% on two out of three hillslopes, although where bare patches (with <10% cover) were connected to gullies and streams, annual sediment yields increased in response to higher rainfall in latter years of the study. It appears that bare patches are the primary source areas for both runoff and erosion on these hillslopes. Achieving further reductions in runoff and erosion in these landscapes may require management practices that improve ground cover and biomass in bare areas, particularly when they are located adjacent to concentrated drainage lines.

Use of CLIMEX modelling to identify prospective areas for exploration to find new biological control agents for prickly acacia

Selection of biocontrol agents that are adapted to the climates in areas of intended release demands a thorough analysis of the climates of the source and release sites. We present a case study that demonstrates how use of the CLIMEX software can improve decision making in relation to the identification of prospective areas for exploration for agents to control the woody weed, prickly acacia Acacia nilotica ssp. indica in the arid areas of north Queensland.

Risk management strategies using seasonal climate forecasting in irrigated cotton production: a tale of stochastic dominance.

Decision-making in agriculture is carried out in an uncertain environment with farmers often seeking information to reduce risk. As a result of the extreme variability of rainfall and stream-flows in north-eastern Australia, water supplies for irrigated agriculture are a limiting factor and a source of risk. The present study examined the use of seasonal climate forecasting (SCF) when calculating planting areas for irrigated cotton in the northern Murray Darling Basin. Results show that minimising risk by adjusting plant areas in response to SCF can lead to significant gains in gross margin returns. However, how farmers respond to SCF is dependent on several other factors including irrigators’ attitude towards risk.

An overview of land application of pig effluent-P using soil P chemistry and mass balance calculations.

Attention is directed at land application of piggery effluent (containing urine, faeces, water, and wasted feed) as a potential source of water resource contamination with phosphorus (P). This paper summarises P-related properties of soil from 0-0.05 m depth at 11 piggery effluent application sites, in order to explore the impact that effluent application has had on the potential for run-off transport of P. The sites investigated were situated on Alfisol, Mollisol, Vertisol, and Spodosol soils in areas that received effluent for 1.5-30 years (estimated effluent-P applications of 100-310000 kg P/ha in total). Total (PT), bicarbonate extractable (PB), and soluble P forms were determined for the soil (0-0.05 m) at paired effluent and no-effluent sites, as well as texture, oxalate-extractable Fe and Al, organic carbon, and pH. All forms of soil P at 0-0.05 m depth increased with effluent application (PB at effluent sites was 1.7-15 times that at no-effluent sites) at 10 of the 11 sites. Increases in PB were strongly related to net P applications (regression analysis of log values for 7 sites with complete data sets: 82.6 % of variance accounted for, p <0.01). Effluent irrigation tended to increase the proportion of soil PT in dilute CaCl2-extractable forms (PTC: effluent average 2.0 %; no-effluent average 0.6%). The proportion of PTC in non-molybdate reactive forms (centrifuged supernatant) decreased (no-effluent average, 46.4 %; effluent average, 13.7 %). Anaerobic lagoon effluent did not reliably acidify soil, since no consistent relationship was observed for pH with effluent application. Soil organic carbon was increased in most of the effluent areas relative to the no-effluent areas. The four effluent areas where organic carbon was reduced had undergone intensive cultivation and cropping. Current effluent management at many of the piggeries failed to maximise the potential for waste P recapture. Ten of the case-study effluent application areas have received effluent-P in excess of crop uptake. While this may not represent a significant risk of leaching where sorption retains P, it has increased the risk of transport of P by run-off. Where such sites are close to surface water, run-off P loads should be managed.

Piggery pond sludge as a nitrogen source for crops. 1. Mineral N supply estimated from laboratory incubations and field application of stockpiled and wet sludge

Piggery pond sludge (PPS) was applied, as-collected (Wet PPS) and following stockpiling for 12 months (Stockpiled PPS), to a sandy Sodosol and clay Vertosol at sites on the Darling Downs of Queensland. Laboratory measures of N availability were carried out on unamended and PPS-amended soils to investigate their value in estimating supplementary N needs of crops in Australia's northern grains region. Cumulative net N mineralised from the long-term (30 weeks) leached aerobic incubation was described by a first-order single exponential model. The mineralisation rate constant (0.057/week) was not significantly different between Control and PPS treatments or across soil types, when the amounts of initial mineral N applied in PPS treatments were excluded. Potentially mineralisable N (No) was significantly increased by the application of Wet PPS, and increased with increasing rate of application. Application of Wet PPS significantly increased the total amount of inorganic N leached compared with the Control treatments. Mineral N applied in Wet PPS contributed as much to the total mineral N status of the soil as did that which mineralised over time from organic N. Rates of C02 evolution during 30 weeks of aerobic leached incubation indicated that the Stockpiled PPS was more stabilised (19-28% of applied organic C mineralised) than the WetPPS (35-58% of applied organic C mineralised), due to higher lignin content in the former. Net nitrate-N produced following 12 weeks of aerobic non-leached incubation was highly correlated with net nitrate-N leached during 12 weeks of aerobic incubation (R^2 = 0.96), although it was <60% of the latter in both sandy and clayey soils. Anaerobically mineralisable N determined by waterlogged incubation of laboratory PPS-amended soil samples increased with increasing application rate of Wet PPS. Anaerobically minemlisable N from field-moist soil was well correlated with net N mineralised during 30 weeks of aerobic leached incubation (R^2 =0.90 sandy soil; R^2=0.93 clay soil). In the clay soil, the amount of mineral N produced from all the laboratory incubations was significantly correlated with field-measured nitrate-N in the soil profile (0-1.5 m depth) after 9 months of weed-free fallow following PPS application. In contrast, only anaerobic mineralisable N was significantly correlated with field nitrate-N in the sandy soil. Anaerobic incubation would, therefore, be suitable as a rapid practical test to estimate potentially mineralisable N following applications of different PPS materials in the field.

Robustness of NIR calibrations for assessing fruit quality

Predictive models based on near infra-red spectroscopy for the assessment of fruit internal quality attributes must exhibit a degree of robustness across the parameters of variety, district and time to be of practical use in fruit grading. At the time this thesis was initiated, while there were a number of published reports on the development of near infra-red based calibration models for the assessment of internal quality attributes of intact fruit, there were no reports of the reliability ("robustness") of such models across time, cultivars or growing regions. As existing published reports varied in instrumentation employed, a re-analysis of existing data was not possible. An instrument platform, based on partial transmittance optics, a halogen light source and (Zeiss MMS 1) detector operating in the short wavelength near infra-red region was developed for use in the assessment of intact fruit. This platform was used to assess populations of macadamia kernels, melons and mandarin fruit for total soluble solids, dry matter and oil concentration. Calibration procedures were optimised and robustness assessed across growing areas, time of harvest, season and variety. In general, global modified partial least squares regression (MPLS) calibration models based on derivatised absorbance data were better than either multiple linear regression or `local' MPLS models in the prediction of independent validation populations . Robustness was most affected by growing season, relative to the growing district or variety . Various calibration updating procedures were evaluated in terms of calibration robustness. Random selection of samples from the validation population for addition to the calibration population was equivalent to or better than other methods of sample addition (methods based on the Mahalanobis distance of samples from either the centroid of the population or neighbourhood samples). In these exercises the global Mahalanobis distance (GH) was calculated using the scores and loadings from the calibration population on the independent validation population. In practice, it is recommended that model predictive performance be monitored in terms of predicted sample GH, with model updating using as few as 10 samples from the new population undertaken when the average GH value exceeds 1 .0 .

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Subsoil constraints to grain production in the cropping soils of the north-eastern region of Australia: an overview

In dryland agricultural systems of the subtropical, semi-arid region of north-eastern Australia, water is the most limiting resource. Crop productivity depends on the efficient use of rainfall and available water stored in the soil during fallow. Agronomic management practices including a period of fallow, stubble retention, and reduced tillage enhance reserves of soil water. However, access to stored water in these soils may be restricted by the presence of growth-limiting conditions in the rooting zone of the crop. These have been termed as subsoil constraints. Subsoil constraints may include compacted or gravel layers (physical), sodicity, salinity, acidity, nutrient deficiencies, presence of toxic elements (chemical) and low microbial activity (biological). Several of these constraints may occur together in some soils. Farmers have often not been able to obtain the potential yield determined by their prevailing climatic conditions in the marginal rainfall areas of the northern grains region. In the past, the adoption of soil management practices had been largely restricted to the top 100 mm soil layer. Exploitation of the subsoil as a source of water and nutrients has largely been overlooked. The key towards realising potential yields would be to gain better understanding of subsoils and their limitations, then develop options to manage them practically and economically. Due to the complex nature of the causal factors of these constraints, efforts are required for a combination of management approaches rather than individual options, with the aim to combat these constraints for sustainable crop production, managing natural resources and avoiding environmental damage.

Fungal colonization of sago starch in Papua New Guinea

Sago starch is an important source of dietary carbohydrates in lowland Papua New Guinea. Over the past 30 years there have been sporadic reports of severe illness following consumption of sago starch. A common assumption is that fungal metabolites might be associated with the illness, leading to the need for a more thorough investigation of the mycoflora of sago starch. Sago starch was collected from areas of high sago consumption in Papua New Guinea for fungal analysis (69 samples). Storage methods and duration were recorded at the time of collection and pH on arrival at the laboratory. Yeasts were isolated from all samples except two, ranging from 1.2 × 103 to 8.3 × 107 cfu/g. Moulds were isolated from 65 of the 69 samples, ranging from 1.0 × 102 to 3.0 × 106 cfu/g. Of 44 samples tested for ergosterol content, 42 samples showed the presence of fungal biomass. Statistical analyses indicated that sago starch stored for greater than five weeks yielded significantly higher ergosterol content and higher numbers of moulds than sago stored for less than five weeks. The method of storage was also shown to influence mould numbers with storage in natural woven fibre containers returning significantly greater numbers than present in other storage methods tested. Potentially mycotoxigenic genera of moulds including Aspergillus and Penicillium were commonly isolated from sago starch, and as such storage factors that influence the growth of these and other filamentous fungi might contribute to the safety of traditional sago starch in PNG.

Achievements in forest tree improvement in Australia and New Zealand 6: Genetic improvement and conservation of Araucaria cunninghamii in Queensland

Araucaria cunninghamii (hoop pine) typically occurs as an emergent tree over subtropical and tropical rainforests, in a discontinuous distribution that extends from West Irian Jaya at about 0°30'S, through the highlands of Indonesian New Guinea and Papua New Guinea, along the east coast of Australia from 11°39'S in Queensland to 30°35'S in northern New South Wales. Plantations established in Queensland since the 1920s now total about 44000 ha, and constitute the primary source for the continuing supply of hoop pine quality timber and pulpwood, with a sustainable harvest exceeding 440 000 m3 y-1. Establishment of these managed plantations allowed logging of all native forests of Araucaria species (hoop pine and bunya pine, A. bidwillii) on state-owned lands to cease in the late 1980s, and the preservation of large areas of araucarian forest types within a system of state-owned and managed reserves. The successful plantation program with this species has been strongly supported by genetic improvement activities since the late 1940s - through knowledge of provenance variation and reproductive biology, the provision of reliable sources of improved seed, and the capture of substantial genetic gains in traits of economic importance (for example growth, stem straightness, internode length and spiral grain). As such, hoop pine is one of the few tropical tree species that, for more than half a century, has been the subject of continuous genetic improvement. The history of commercialisation and genetic improvement of hoop pine provides an excellent example of the dual economic and conservation benefits that may be obtained in tropical tree species through the integration of gene conservation and genetic improvement with commercial plantation development. This paper outlines the natural distribution and reproductive biology of hoop pine, describes the major achievements of the genetic improvement program in Queensland over the past 50+ y, summarises current understanding of the genetic variation and control of key selection traits, and outlines the means by which genetic diversity in the species is being conserved.

Climatic-requirements models of cat's claw creeper Macfadyena unguis-cati (Bignoniaceae) to prioritise areas for exploration and release of biological control agents

Climate matching software (CLIMEX) was used to prioritise areas to explore for biological control agents in the native range of cat's claw creeper Macfadyena unguis-cati (Bignoniaceae), and to prioritise areas to release the agents in the introduced ranges of the plant. The native distribution of cat's claw creeper was used to predict the potential range of climatically suitable habitats for cat's claw creeper in its introduced ranges. A Composite Match Index (CMI) of cat's claw creeper was determined with the 'Match Climates' function in order to match the ranges in Australia and South Africa where the plant is introduced with its native range in South and Central America. This information was used to determine which areas might yield climatically-adapted agents. Locations in northern Argentina had CMI values which best matched sites with cat's claw creeper infestations in Australia and South Africa. None of the sites from where three currently prioritised biological control agents for cat's claw creeper were collected had CMI values higher than 0.8. The analysis showed that central and eastern Argentina, south Brazil, Uruguay and parts of Bolivia and Paraguay should be prioritised for exploration for new biological control agents for cat's claw creeper to be used in Australia and South Africa.

Leptospira weilii serovar Topaz, a new member of the Tarassovi serogroup isolated from a bovine source in Queensland, Australia

This paper reports on a Leptospira isolate of bovine origin and its identification as belonging to a previously unknown serovar, for which the name Topaz is proposed. The isolate (94-79970/3) was cultured from bovine urine from a north Queensland dairy farm in Australia. Strain 94-79970/3 grew at 30 °C in Ellinghausen McCullough Johnson Harris (EMJH) medium but failed to grow at 13 °C in EMJH medium or in the presence of 8-azaguanine. Serologically, strain 94-79970/3 produced titres against the Leptospira borgpetersenii serovar Tarassovi, the reference strain for the Tarassovi serogroup; however, no significant titres to any other serovars within the serogroup were obtained. Using 16S rRNA and DNA gyrase subunit B gene analysis, strain 94-79970/3 was identified as a member of the species Leptospira weilii. We propose that the serovar be named Topaz, after the location where the original isolate was obtained. The reference strain for this serovar is 94-79970/ 3 (=KIT 94-79970/35LT722).

Distribution and eradication of an exotic tephritid fruit fly in Australia: Relevance of invasion theory

Data from the eradication of the incursion of Bactrocera papayae Drew and Hancock (Dipt.: Tephritidae) in Australia (1995-1998) are used to assess the significance of various aspects of invasion theory, including the influence of towns on establishment, influence of propagule pressure on the pattern of establishment, and the existence of source-sink dynamics. Because there were no sentinel traps in place, considerable spread had occurred before the eradication campaign started. The distribution of fly density around the epicentre in the town of Cairns and a transect along the main traffic routes to the north and south fitted a Cauchy model with a tail having the same slope as a power model with an exponent of -2.4 extending to 160 km. The Cauchy model indicated that 50% of the flies on the transect would have occurred within 3.2 km of the epicentre, 90% within 13.2 km, and 99% within 60 km. The two major satellites at Mareeba (35 km from the epicentre in Cairns) and Mossman (65 km) were not used for the transect data and had respectively 15 and 30 times the density predicted by the model. The proportion of traps that caught flies (a measure of site occupancy) fell with distance from the epicentre. B. papayae was trapped consistently on only three of the 16 rainforest transects that were surveyed and these were relatively close to urban areas where eradication efforts were intense. Despite there being no eradication effort in the rainforest, the trends to extinction were similar to those in adjacent areas. The strategy of initially concentrating eradication efforts on the core and major satellites while maintaining a quarantine barrier at the airport and the boundaries of the infested area appears to be the key to the containment and rapid eradication of the incursion.

Encroachment of Echinococcus granulosus into urban areas in eastern Queensland, Australia

Objective: To investigate the prevalence of Echinococcus granulosus in wild dogs (dingos and dingo-domestic dog hybrids) living in and around human habitation on Fraser Island and in townships of the Maroochy Shire, on Queensland's Sunshine Coast, Australia. Design: Wild dogs were humanely killed on Fraser Island and in the Maroochy Shire because they were deemed a potential danger to the public. Their intestines were collected and the contents examined for intestinal parasites. Procedure: Intestines were removed as soon after death as possible, packed in plastic bags and kept frozen until examination. The intestinal contents were washed, sieved and examined microscopically for the presence of helminths, which were identified and counted. Results: Intestines from 108 wild dogs, 7 foxes and 18 Fraser Island dingoes were examined. Echinococcus granulosus was only present in the wild dogs from Maroochy Shire (46.3%) with worm burdens of between 30 and 104,000. Other helminths included Spirometra erinacei, Dipylidium caninum, Taenia spp., Ancylostoma caninum and Toxocara canis. Two specimens of a trematode (Haplorchinae sp.) usually found infecting fish and seabirds were recovered from a Fraser Island dingo. Conclusion: Dingoes on Fraser Island are not infected with E. granulosus and do not pose a hydatid disease public health risk to residents or visitors. However, wild dogs examined from the Maroochy Shire do present a potential hydatid disease public health risk.

Rice responses to soil management in a rice-based cropping system in the semi-arid tropics of southern Lombok, Eastern Indonesia

This paper is the first of a series that investigates whether new cropping systems with permanent raised beds (PRBs) or Flat land could be successfully used to increase farmers' incomes from rainfed crops in Lombok in Eastern Indonesia. This paper discusses the rice phase of the cropping system. Low grain yields of dry-seeded rice (Oryza sativa) grown on Flat land on Vertisols in the rainfed region of southern Lombok, Eastern Indonesia, are probably mainly due to (a) erratic rainfall (870-1220 mm/yr), with water often limiting at sensitive growth stages, (b) consistently high temperatures (average maximum - 31 C), and (c) low solar radiation. Farmers are therefore poor, and labour is hard and costly, as all operations are manual. Two replicated field experiments were run at Wakan (annual rainfall = 868 mm) and Kawo (1215 mm) for 3 years (2001/2002 to 2003/2004) on Vertisols in southern Lombok. Dry-seeded rice was grown in 4 treatments with or without manual tillage on (a) PRBs, 1.2 m wide, 200 mm high, separated by furrows 300 mm wide, 200 mill deep, with no rice sown in the well-graded furrows, and (b) well-graded Flat land. Excess surface water was harvested from each treatment and used for irrigation after the vegetative stage of the rice. All operations were manual. There were no differences between treatments in grain yield of rice (mean grain yield = 681 g/m(2)) which could be partly explained by total number of tillers/hill and mean panicle length, but not number of productive tillers/hill, plant height or weight of 1000 grains. When the data from both treatments on PRBs and from both treatments on Flat land, each year at each site were analysed, there were also no differences in grain yield of rice (g/m(2)). When rainfall in the wet season up to harvest was over 1000 mm (Year 2; Wakan, Kawo), or plants were water-stressed during crop establishment (Year 1; Wakan) or during grain-fill (Year 3: Kawo), there were significant differences in grain yield (g/1.5 m(2)) between treatments; generally the grain yield (g/1.5 m(2)) on PRBs with or without tillage was less than that on Flat land with or without tillage. However, when the data from both treatments on PRBs and from both treatments on Flat land, each year at each site, were analysed, the greater grain yield of dry-seeded rice on Flat land (mean yield 1 092 g/1.5 m(2)) than that on PRBs (mean 815 g/1.5 m(2)) was mainly because there were 25% more plants on Flat land. Overall when the data in the 2 outer rows and the 2 inner rows on PRBs were each combined, there was a higher number of productive tillers in the combined outer rows (mean 20.7 tillers/hill) compared with that in the combined inner rows on each PRB (mean 18.2 tillers/hill). However, there were no differences in grain yield between combined rows (mean 142 g/m row). Hence with a gap of 500 mm (the distance between the outer rows of plants on adjacent raised beds), plants did not compensate in grain yield for missing plants in furrows. This suggests that rice (a) also sown in furrows, or (b) sown in 7 rows with narrower row-spacing, or (c) sown in 6 rows with slightly wider row-spacing, and narrower gap between outer rows on adjacent beds, may further increase grain yield (g/1.5 m(2)) in this system of PRBs. The growth and the grain yield (y in g/m(2)) of rainfed rice (with rainfall on-site the only source of water for irrigation) depended mainly on the rainfall (x in mm) in the wet season up to harvest (due either to site or year) with y = 1. 1x -308; r(2) = 0.54; p < 0.005. However, 280 mm (i.e. 32%) of the rainfall was not directly used to produce grain (i.e. when y = 0 g/m(2)). Manual tillage did not affect growth and grain yield of rice (g/m(2); g/1.5 m(2)), either on PRB or on Flat land.