Sheet and Rill Erosion and Sediment Delivery to Streams: A Basin Wide Estimation at hillslope to Medium Catchment Scale Report E to Project D10012 of Murray Darling Basin Commission: Basin-wide Mapping of Sediment and Nutrient Exports in Dryland Regions of the MDB
Hua Lu, Chris J. Moran, Ian P. Prosser, Michael R. Raupach, Jon Olley, and Cuan Petheram
CSIRO Land and Water, Canberra Technical Report 15/03, June 2003
C S I R O L A N D a n d W AT E R
Sheet and Rill Erosion and Sediment Delivery to Streams: A Basin Wide Estimation at hillslope to Medium Catchment Scale Report E to Project D10012 of Murray Darling Basin Commission: Basin-wide Mapping of Sediment and Nutrient Exports in Dryland Regions of the MDB
Hua Lu, Chris J. Moran, Ian P. Prosser, Michael R. Raupach, Jon Olley, and Cuan Petheram
CSIRO Land and Water, Canberra Technical Report 15/03, June 2003
Copyright © 2003 CSIRO and the Murray-Darling Basin Commission. This work is copyright. It may be reproduced subject to the inclusion of an acknowledgement of the source.
Authors Hua Lu, Chris J. Moran, Ian P. Prosser, Michael R. Raupach, Jon Olley and Cuan Petheram CSIRO Land and Water, PO Box 1666, Canberra, 2601, Australia. E-mail:
[email protected] Phone: 61-2-6246-5923 For bibliographic purposes, this document may be cited as: Lu, H., Moran, C.J., Prosser, I.P., Raupach, M.R., Olley, J. and Petheram, C. (2003) Hillslope erosion and sediment delivery: A basin wide estimation at medium catchment scale, Technical Report 15/03, CSIRO Land and Water. A PDF version is available at: http://www.clw.csiro.au/publications/technical2003/tr15-03.pdf
ISSN 1446-6163
Important Disclaimer CSIRO Land and Water and the Murray-Darling Basin Commission advise that the information contained in this publication comprises general statements based on scientific research. The reader is advised and needs to be aware that such information may be incomplete or unable to be used in any specific situation. No reliance or actions must therefore be made on that information without seeking prior expert professional, scientific and technical advice. To the extent permitted by law, CSIRO Land and Water and the Murray-Darling Basin Commission (including their employees and consultants) excludes all liability to any person for any consequences, including but not limited to all losses, damages, costs, expenses and any other compensation, arising directly or indirectly from using this publication (in part or in whole) and any information or material contained in it.
2
Table of Contents CSIRO Land and Water........................................................................................................ 1 Executive Summary.............................................................................................................. 7 1 Introduction .................................................................................................................... 9 2 Background .................................................................................................................. 11 2.1 Characteristics of the MDB ..................................................................................... 11 2.1.1 Climate............................................................................................................. 11 2.1.2 Topography ...................................................................................................... 11 2.1.3 Geology and Soil.............................................................................................. 12 2.1.4 Land Use and Management Practices .............................................................. 12 2.2 Spatial Settings and Terminology............................................................................ 13 3 Methodology................................................................................................................. 14 3.1 Hillslope Sheetwash and Rill Erosion...................................................................... 14 3.2 Hillslope Sheet and Rill Erosion under Natural Condition...................................... 15 3.3 Sediment Delivery Ratio (SDR) .............................................................................. 16 3.3.1 Background of SDR ......................................................................................... 16 3.3.2 A New SDR Theory......................................................................................... 18 3.4 Statistical Analysis of Effective Rainfall duration and Intensity............................. 22 3.5 Estimations of Residence Time ............................................................................... 30 3.5.1 Sediment Residence Time as a Function of Particle Size................................ 30 3.5.2 Estimating Travel Time of Water Particles th0 and tn0 ..................................... 31 4 Results ........................................................................................................................... 36 4.1 Hillslope Erosion under Current Land Use.............................................................. 36 4.2 Hillslope Erosion under Natural Conditions............................................................ 38 4.3 Spatial Characteristics of Effective Rainfall Duration and Intensity....................... 40 4.4 Spatial Distribution of Sediment Delivery Ratio ..................................................... 41 5 Discussions and Conclusions................................................................................. 44 Acknowledgments .............................................................................................................. 45 References ............................................................................................................................ 47 Appendix I: Land Use Data .................................................................................................. 50 Appendix II: Pluviograph Rainfall Data ............................................................................. 52
3
List of Figures Figure 1: Major landuse groups in the MDB. .......................................................................... 13 Figure 2: SDR vs catchment area relationships obtained from different areas around the world. ............................................................................................................................... 17 Figure 3: Diagram of a two storage lumped linear model of SDR at catchment scale (after Sivapalan et al. 2001, modified). See text for detail........................................................ 19 Figure 4. Comparison of SDR (%) measurements (Roehl 1962), modeled average SDR and flow response (Robinson and Sivapalan 1997). It shows that flow response represents the upper envelope of the SDR. ............................................................................................. 20 Figure 5: SDR as a function of channel residence time for different values of ter and th (upper panel); SDR as a function of catchment area for different values of ter and th. SDR measurements from USA catchments (Roehl 1962) are also shown as red dots (lower panel)................................................................................................................................ 21 Figure 6: Site locations of pluviograph rainfall data and their relative position to MDB. ...... 22 Figure 7: All rainfall events characterised by their 30 intensity and duration (upper panel); Fit probability density functions to Gamma and exponential distributions for both duration and intensity (second and lower panels). ......................................................................... 24 Figure 8: Rainfall events which have depth equal or greater than 12.7 mm (upper panel); Fit probability density functions to Gamma and exponential distributions for both duration and intensity (second and lower panels). ......................................................................... 25 Figure 9: Effective rainfall events which have depth equal or greater than 12.7 mm (upper panel); Fit probability density function of effective duration to Gamma and exponential distributions (lower panel). .............................................................................................. 26 Figure 10: Relationships between effective 30-min. rainfall intensity and the ratio between mean annual R-factor and mean annual rainfall. ............................................................. 27 Figure 11: Relationships between rainfall duration (tr) to mean annual rainfall (MAR), effective 30-min intensity (MI30), MAR/MI30, and MAR2/R. ................................ 27 Figure 12: Relationships of effective rainfall duration and it relative errors. ................. 28 Figure 13: Error estimations of rainfall duration. Upper Panel: Comparison between rainfall duration estimated using site specific pluviograph data and that estimated using regionalised relationships. Middle Panel: Absolute error [hrs] plotted against number of year with complete data. Lower Panel: Relative error plotted against number of year with complete data. The crosses are the sites have shorter records and relatively larger errors. They are not used in the final relationships that are applied across the MDB................. 29 Figure 14: Diagram of the particle size effect on sediment travel time in relation to the travel time of water particles...................................................................................................... 30
4
Figure 15: Flow chart for the calculation of travel time of water particles. ............................ 35 Figure 16: Estimated annual average sheet and rill erosion rate. ..................................... 37 Figure 17: Monthly distribution of total soil loss rate for the Basin. ................................ 38 Figure 18: Comparison between natural C-factor values modeled using Cubist and those extracted from C-factor map at the locations of minimum cover disturbance. There are 9916 points in total. The line of best fit and 1:1 line are shown. ......................... 39 Figure 19: Estimated annual erosion rate under natural conditions (pre-European settlement conditions). ................................................................................................... 39 Figure 20: Estimated Ratio between erosion rate under current landuse and that under natural conditions......................................................................................................................... 40 Figure 21: Estimated effective 30-min. rainfall intensity for south-eastern Australia. The boundary of MDB is shown........................................................................................... 40 Figure 22: Estimated rainfall duration (left panel) and effective storm duration (right panel for south-eastern Australia. The boundary of MDB is shown. ....................... 41 Figure 23: Estimated travel time of water particles th0 and tn0 for each sub-catchment element in MDB. ........................................................................................................................... 42 Figure 24: Estimated Sediment delivery ratio fro clay, silt and sand particles........................ 42 Figure 25: Estimated overall sediment delivery ratio from each sub-catchment elements. .......................................................................................................................................... 43 Figure 26: Estimated specific sediment yield [t/ha/yr] for each sub-catchment element. ....... 43 Figure AI.1: Data sources of land use used in this project. ..................................................... 51
5
List of Tables Table 1: Landuse groups used to calculate sheet and rill erosion rate..................................... 14 Table 2: Typical values of CN for some land use group. ........................................................ 33 Table 3: Values of Manning’s n used in this study for common land use and vegetation cover groups for overland flow.................................................................................................. 34 Table 4: channel roughness parameter a values used in this study.......................................... 36 Table 5: Three erosion groups (high, medium and low) and their relation to percentage of agricultural lands.............................................................................................................. 37 Table 6: Soil loss rate from land use categories. ..................................................................... 38 Table AI.1: Summary of locally supplied land use data used in this study. ............................ 50 Table AII.1: Details of the pluviograph rainfall sites. ............................................................. 52
6
Executive Summary This report presents a scientific and technical description of the modelling framework and main results for the long-term average hillslope erosion and sediment delivery to streams at hillslope to medium scale catchment over the Murray Darling Basin. The work was a part of project D10012 of the Murray-Darling Basin Commission (MDBC), "Basin-wide mapping of sediment and nutrient exports in dryland regions". Gully and stream bank erosion, sediment transport at larger scale with intervening deposition to flood plain, and associated nutrient exports are dealt with in separated reports (DeRose et al. 2003; Hughes and Prosser 2003). The specific objectives of this part of the work are basin-wide mapping by: (1) Quantifying the hillslope sheetwash and rill erosion under current land use condition; (2) Quantifying the inherent natural hazard of hillslope sheetwash and rill erosion; (3) Determining the amount of sediment generated by sheetwash and rill erosion delivered to the stream network from the sub-catchment elements with contributing area around 50 - 100 km2; (4) Interpreting results in terms of comparison with pre-European land use conditions. The modelling frameworks are described as follows. We undertook new assessments of hillslope sheetwash and rill erosion across the MDB, building upon our previous work for the National Land and Water Resources Audit (NLWRA) (Lu et al. 2001; Lu et al. 2003b). The mean annual hillslope sheetwash and rill erosion was modelled using the Universal Soil Loss Equation (USLE) (Wischmeier and Smith 1978; Renard et al. 1997), which is a model of surface wash and rill erosion based upon factors of rainfall erosivity, terrain, soil erodibility and vegetation cover. USLE factors were calculated from digital elevation models (DEMs), soil attribute maps, land use maps, remote sensing imagery and daily rainfall surfaces. Time series of remote sensing imagery and daily rainfall were used to incorporate the effects of seasonally varying cover and rainfall intensity. Further, we used new digital maps of soil and terrain properties. In this project, improvements to the assessment of sheetwash and rill erosion were made by compiling higher resolution land use data for the MDB from a range of sources and by incorporating a database on crop rotation, tillage and other land management practices. These new data, together with improved analysis of remote sensing data, enabled a more accurate prediction of the effect of vegetation cover and cover management on hillslope erosion. To relate hillslope erosion estimates to riverine water quality, it is necessary to estimate the properties of eroded soil that is delivered to the waterways for further transport. A theory that relates long-term averaged sediment delivery to the statistics of rainfall and catchment parameters was proposed. The derived flood frequency approach was adapted to investigate the problem of regionalization of the sediment delivery ratio (SDR) across the Basin. SDR, a measure of catchment response to the upland erosion rate, was modeled by two lumped linear stores arranged in series: hillslope transport to the nearest streams and flow routing in the channel network. The theory shows that the ratio of catchment sediment residence time (SRT)
7
to average effective rainfall duration is the most important control in the sediment delivery processes. In this study, catchment SRTs were estimated using time of the concentration for overland flow multiplied by an enlargement factor which is a function of particle size. Rainfall intensity and effective duration statistics were regionalized by using long-term measurements from 195 pluviograph sites within and around the Basin. The major findings are as follows. Hillslope Erosion under Current Land Use: Our spatial modelling of sheetwash and rill erosion estimated that a total 2.2 × 108 t yr-1 of sediment moves locally across the Basin, at a mean rate of 2.1 t ha-1 yr-1. Erosion rate increases from south to north and from arid areas to temperate regions. About two-thirds of erosion occurs in the summer period. Agricultural lands have slightly higher erosion rates, at a mean rate of 2.3 t ha-1 yr-1 as most of the agricultural lands are located in the flood plains. Inherent Natural Hillslope Erosion: It was estimated that soil erosion rates are low under preEuropean natural vegetation conditions. The rates are 3 to 10 times on average and up to 100 times smaller than that under current land use. Spatial Characteristics of Erosive Rainfall: Rainfall intensity increases from south to north and from west to east. Coinciding with the effect of topography, it defines the broad pattern of hillslope erosion. Rainfall duration is greater in temperate than arid regions and decreases from uplands to flat inlands. This dissipates sediment transport energy and the whole system is inefficient for sediment transport from erosion sources to basin outlet. Sediment Delivery Ratio and Sediment Yield: The averaged SDR is about 5%, which is lower than the average estimated in other countries for catchments with similar contributing area. Most sub-catchment elements have SDR smaller than 5%, suggesting inefficiency of sediment transport in the broad areas of the Basin. Larger SDRs are obtained at the eastern edge of the Basin, with the Australian Alps having the highest SDR values, followed by the central south of the Murrumbidgee and Bathurst regions. Sediment yield is low for the majority of subcatchment elements with area-specific sediment yield around 0.13 t ha-1 yr-1, which also represents the Basin average. Problem areas are located mainly in the eastern edge of the Basin.
8
1
Introduction
Soil erosion and sediment transport are recognised as major environmental hazards in the Murray Darling Basin (MDB). It is governed by topography, climate, soil, vegetation cover, land use and management factors, through mechanisms including, particle detachment by raindrop impact, hydrology, flow hydraulics and other processes. Ability to estimate erosion rate across the whole basin is significant for three reasons. Firstly, soil erosion has a range of environmental impacts, including loss of organic matter and nutrients, reduction of crop productivity, and downstream water quality degradation (Newcombe and MacDonald 1991). The integrated impacts are often revealed and of importance at catchment or even larger scales. Secondly, effective control of soil erosion is a critical component of natural resource management when the aim is to achieve sustainable agriculture and acceptable ecosystem integrity (Pimentel et al. 1995; Rutherfurd et al. 1998). With limited resources, national scale erosion maps are useful for guiding investment prioritization in effective remediation programs. Thirdly, to aid estimations of soil erosion contributions and their impacts, the effects of changes in climatic conditions, vegetation, and land use on soil erosion rates need to be assessed at regional to continental scales (Pimentel et al. 1995). Due to the prevalence of high-value commodities in the Basin, comprehensive data on full areal extent and severity of the Basin’s soil erosion and sediment delivery is of both economic and environmental importance. Information on spatially distributed sediment delivery is useful in identifying relative importance between sediment sources and the effectiveness of sediment delivery. It helps to establish strategies in effective erosion control, rehabilitation planning, and achieving long-term sustainable productivity in the Basin. In the past, there have been several attempts to estimate soil erosion rates at regional to continental scale, i.e., reviews of erosion data (Edwards 1993); synthesis of hillslope erosion rates and sediment transport (Wasson, et al. 1996), reconnaissance survey using caesium-137 (Loughran and Elliott 1996), and quantitative spatial modelling using USLE (Rosewell 1997). Variations and uncertainties exist in all the previous estimations. The major discrepancies of previous studies are largely due to lack of high quality consistent spatial data and our inability to model the complex systems which involve subsystem interactions both in time and in space. In the late 1990s, Australia launched the National Land and Water Resources Audit (NLWRA 2001) to assess the condition of its land and water resources. The continent-wide assessment of sheetwash and rill erosion was conducted as part of a broader assessment of the conditions of Australian agricultural land (NLWRA 2001). Hillslope sheetwash and rill erosion estimation in this project was building upon our work in the NLWRA project (Lu et al. 2001; Lu et al. 2003b). Improvements were made by compiling higher resolution land use data for the MDB from a range of sources and by incorporating a database on crop rotation, tillage and other land management practices. These new data, together with improved analysis of remote sensing data, enabled a more accurate prediction of the effect of vegetation cover and cover management on hillslope erosion. Only a small fraction of the soil moving on hillslopes is actually delivered to streams (Edwards 1993; Wasson et al. 1996). This implies that most of the sediment travels only a short distance (Parsons and Stromberg 1998) and is deposited before leaving the hillslope. In general, the amount of sediment deposited is intimately related to the topography, climate, soil, vegetation cover, and land use conditions, which are all closely related to the hydrological processes. The travel time for transport of sediment across a field or hillslope is 9
often longer than the duration of runoff-generating events so that runoff infiltrates and is not delivered to the stream, along with the sediment it carries. In some environments there is also patchy generation of runoff on impermeable areas which then infiltrates on other patches of high infiltration, often at sites with better cover. Topography can induce deposition through its influence on the capacity of overland flow to transport sediment. Reductions in gradient and the dispersion of overland flow can both cause deposition. Farm structures, such as contour banks and dams, can have similar effects, altering flow paths or trapping runoff. Deposition also results from abrupt changes to vegetation cover as runoff travels downslope. This causes deposition in backwaters and reduces the sediment transport capacity of flow. In large-scale modelling of sediment transport, this phenomenon of different sediment transport rates between hillslope and catchment scale is usually modelled using a scaling factor called the hillslope sediment delivery ratio (HSDR). This avoids the need to explicitly model patterns of deposition on hillslopes which is not possible across such large areas as the MDB. Prosser et al. (2001) developed a spatially distributed model of mean annual sediment budgets for river basins. The model, SedNet (Sediment River Network Model), used spatial modelling of the erosion, deposition, and transport processes that move sediment and nutrients within landscapes and streams to produce regional budgets for the Murray Darling Basin. The sources of sediment considered are soil erosion by surface (hillslope) processes (Lu et al. 2001), gully erosion and riverbank erosion (Hughes and Prosser 2003). These sediment sources were routed through the river network using a simple conceptual model of the primary controls on sediment export and deposition. The results demonstrate that there is a reasonable correlation between observed and predicted specific sediment yields. The hillslope delivery in SedNet is modelled by USLE-SDR approach (Lu et al. 2001). In the NLWRA project, uniform hillslope sediment delivery ratio was applied to the whole MDB by using HSDR as a calibration factor to obtain the best results (Prosser et al. 2001). This neglects the environmental variation across the MDB that would give varying sediment delivery potential. For example, short steep hillslopes experiencing long storms in the east of the Basin will have a greater delivery ratio than the long, flat hillslopes of the western regions. In this project, a new approach to model HSDR was proposed and implemented to estimate the long-term averaged spatially-distributed HSDR over the entire MDB. HSDR, a measure of catchment response to the upland erosion rate, is modeled by two lumped linear stores arranged in series: hillslope transport to the nearest streams and flow routing in the channel network. A theory developed in hydrologic scaling is adapted here to relate long-term averaged sediment delivery to the effective rainfall duration and catchment sediment residence time (SRT). Average rainfall intensity and effective duration were regionalized by using long-term measurements from 195 pluviograph sites within and around the Basin. SRT is estimated using time of the concentration for overland flow multiplied by an enhancement factor which is a function of particle size. The model was implemented across the MDB by using spatially distributed soil, vegetation, topographical and land use properties under a GIS environment. The report is organised as follows: Section 2 briefly describes some general characteristics of the MDB in terms of climate, topography and soils which have major impact on erosion and sediment transport processes. A set of terminology is given for clarity of later spatial description and interpretation of results. Section 3 presents the methods used for the modelling with emphasis on hillslope sediment delivery ratio (HSDR) and analysis of rainfall
10
intensity data (6-min. interval pluviograph data). The main results are presented in Section 4. Discussions and conclusions are given in Section 5.
2 Background 2.1
Characteristics of the MDB
The Murray-Darling Basin covers an area of 1.1 × 106 km2 or about 14% of Australia. The Basin includes the three longest rivers in Australia. The Darling is 2,740 km long from its source in the north to its confluence with the Murray at Wentworth, the Murray is 2,530 km long from its source in the Australian Alps to its mouth on Encounter Bay in South Australia, and the Murrumbidgee is 1,690 km long. As a semi-arid country with relatively high economic dependence on agricultural revenue, the MDB is of national economic importance with rich irrigation, farming and grazing land. The Basin accounts for 40% of Australia’s agricultural production, utilizing about 70% of all water used for agriculture across the nation. The 1,500,000 hectares under irrigation for crops and pastures represents 70% of the total area under irrigation in Australia. More than 80% of the divertible surface water resource is consumed in the Basin. The Basin holds a population of 2 million people, which is about 10% of the national population. 2.1.1
Climate
There are a range of climatic conditions across the Basin, with cool humid conditions on the eastern uplands, and sub-tropical conditions in the northeast. The climate to the southeast is temperate, while the large western plains are semi-arid and arid areas. Annual precipitation in the Basin ranges from 185 mm to 2,500 mm. The potential evaporation rate is more than twice the precipitation rate. Mean annual evapotranspiration generally increases as rainfall decreases. Less than 10% of stream flow reaches the major rivers (Murray and Darling) and less than 5% of total rainfall is exported to the sea (Crabb 1997). The Basin has large inter-annual variability of the rainfall, mainly due to the impact of the El Nino - Southern Oscillation (ENSO) on the climate of southeastern Australia. This variability in rainfall is amplified in the annual runoff, which is more variable than runoff elsewhere in the world (except for parts of Southern Africa that experience a similar climate). These variations have profound effect on sediment delivery and transport in the basin. 2.1.2
Topography
Combined with the mountains of the Australian Alps and steep hills and colluvial slopes of the Great Dividing Range, much of the Basin consists of the Murray-Murrumbidgee Riverine plain, the Darling floodplain and alluvial floodplains of other tributaries. The low relief over most of the Murray Basin occupies most of the area towards the arid west. Due to this topographic setting, the stream flow and sediment generated from the high rainfall areas are often dispersed or evaporated when the water reaches lowland floodplains.
11
2.1.3
Geology and Soil
The spatial distribution and the properties of soils reflect the effect of climate, topography, flora and fauna acting on parent material over time. Organic soils are found at high altitude in alpine areas. Soils on steep mountain slopes, upper valleys and their terraces reflect the sequence of periodic erosion-deposition driven by tectonic activity and/or climate change (Butler et al. 1983). Soil thickness and horizon development are a function of the age of the deposits; buried soils are widespread (Rowe et al. 1978). Gradational soils of various levels of differentiation reflect the age of their parent material. Red and yellow duplex soils are widespread on rounded spurs, ridges, and hills, and dissected colluvial deposits (Rowe et al. 1978). On the alluvial deposits of the Riverine Plain, sediment texture and drainage control soil profile colour and development. A characteristic soil catena is associated with leveefloodplain transects of prior streams (Butler et al. 1983). Red-brown earths are found on levees, sandy on the crest of the levee and loamy on the backslope. Grey, brown, and red clays are extensive on the floodplains, with gilgai and soluble salt content increasing downstream and with distance transverse from the levee (Butler et al. 1983). Wind-blown parna mantles much of the Riverine Plain but is most prominent on foothills and hilly inliers (Butler et al. 1983). In the Darling alluvial plain, grey self-mulching clay soils derived from basalt are extensive (Butler and Hubble 1978). 2.1.4
Land Use and Management Practices
The major land use types in the Basin are dryland grazing (native pasture), cropping dominated by winter cereals, improved pasture, open forest, and agroforestry. The native vegetation is diverse with grassland, open woodland, woodland and shrubland environments and a very small area of dense vegetation growth in the eastern part of the Basin. The predominant land use is grazing but due to the economic benefits there is a shift towards cropping. As the high resolution land use and land management data sets currently only cover selected dryland areas as part of the Landmark project, 1-km resolution snap-shot land use data derived from 1996 NOAA LAC remote sensed images by Bureau of Rural Sciences (BRS 2001) was used for the other areas (see Figure AI.1 and Table AI.1 in Appendix I for detail). Figure 1 shows spatially distributed current land use classes in the Basin. In terms of management practices, burning was a widespread practice in the early 1970s throughout the Basin. In the 1990s, general observations suggest that stubble retention is much more common north of about Parkes, relative to southern areas. In the northern part of the NSW south-western slopes, a survey carried out by Vanclay (1997) suggested that about 50% of farmers burn stubble, and that burning tends to be associated with cultivation rather than direct drilling. In northern NSW, a survey by Martin et al. (1988) showed that the incidence of stubble burning ranged from 0% in dry years to 28% in years with disease problems. Burning generally occurred soon after harvest – it removes approximately 90% of the stubble cover. The trend in tillage practice is for greater retention of crop residues and fewer tillage operations, especially in the northern part of the Basin. This helped not only to reduce erosion but also increased profits of the farmers. Those features of tillage and crop rotation were considered in our hillslope sheetwash and rill erosion modelling.
12
Typical pasture dry matter production in the Basin is around 8–10 t/ha DM in a good season down to 3 t/ha DM in a poor season. In severe droughts, dry matter production is negligible. Stocking rates tend to be in the range 1.0 – 1.5 DSE/ha. In marginal western areas, suitable pasture systems for rotation with cereals have not yet been developed. In the Cobar-NynganWalgett region, where most of the new development has occurred since 1970, both forest and open grassland have been converted to crop and pasture production (Swift and Skjemstad 2002).
Figure 1: Major landuse groups in the MDB.
2.2
Spatial Settings and Terminology
For modelling purposes, SedNet (Prosser et al. 2001; DeRose et al. 2003) spatially divided the MDB into around 10,000 sub-areas according to its topography using ESRI ArcInfo software (ESRI 2003) and 9’’ digital elevation model (DEM) derived by the Australian National University (Hutchinson et al. 2001). The sub-areas, which are constituted by many grid cells, are the basic constituent elements used to compute hillslope sheet and rill erosion, hillslope sediment delivery ratio, gully erosion, and bank erosion. Those sub-areas are called subcatchment elements and have contributing area around 50 - 100 km2. Grid cells are the basic constituent element for hillslope sheetwash and rill erosion modelling and the results are presented as the same raster GIS formant. For HSDR modelling, grid cells remain the basic constituent element but the results are presented at the sub-catchment element level. For clarity, in this report, we use the following terminology: Contribution point: This refers to the river export point for evaluation of suspended sediment contribution. Sometimes, it is called sediment control location. Sub-catchment element: It is the basic constituent element of SetNet model. Normally, each has a contributing area around 50 - 100 km2. There are nearly 10,000 subcatchment elements covering the MDB. 13
Sub-catchment: A group of sub-catchment elements. These equate to tributary rivers of catchments. E.g., the Cotter sub-catchment is a tributary of the Murrumbidgee catchment. Catchment: Refers to the major upland catchment areas and associated rivers, e.g., the Murrumbidgee Catchment. Basin: Refers to the Murray Darling Basin as a whole. SDR: The mean value of sediment delivery ratio from a sub-catchment element. It is also called HSDR in the other reports produced by this project. Hillslope Erosion: Refers to hillslope sheetwash and rill erosion only.
3 Methodology 3.1
Hillslope Sheetwash and Rill Erosion
Mean annual soil erosion under current land use was predicted using the Universal Soil Loss Equation (USLE), a model of surface wash and rill erosion based upon factors of rainfall erosivity, terrain, soil erodibility, and vegetation cover. We mapped USLE factors from digital elevation models (DEMs), soil property maps, remotely sensed images and climate surfaces. Innovations were made in obtaining high resolution terrain properties from coarse resolution DEMs (Gallant 2001), and seasonal vegetation cover mapping from 14 yr of imagery (Lu et al. 2003c), and seasonal rainfall erosivity estimation using 20 yr of daily rainfall (Lu and Yu 2002). Technique details of the modelling process can be found in Lu et al. (2001). Further improvement on the estimation of erodibility (USLE K-factor), cover and management factor (USLE C-factor) and model validation can be found in Lu et al. (2003b). Table 1: Landuse groups used to calculate sheet and rill erosion rate. Groups 10 11 12
Land use Descriptions Built-up area Perennial watercourse and lake, Mangrove, Reservoir, Saline, Coastal flat, Swamp Non-perennial watercourse and lake
21
Closed forest
22
Open forest
23
Woodland
24 25
Commercial native forest production, Plantation fruit, Agroforestry, Apples, Citrus, Grapes, Stone fruit, Pears, Plantation National Park
31
Cereals excluding rice
32
Legumes
33
Other non-cereal crops
34
Oilseeds
35
Non-cereal forage crops
36
Rice
37
Cotton
38
Potatoes
39
Sugar cane
14
40 41
Other vegetables Nuts
42
Improved Pastures
51
Residual/Native Pasture
This report provides an update on the previous assessment by Lu et al. (2001) by including improved vegetation cover estimation (Lu et al. 2003c) and high resolution land use data supplied by local agencies. Vegetation cover and land use data were used to calculate USLE C-factor. Compared with the snap-shot 1 km resolution land use coverage produced by the Bureau of Rural Sciences (BRS), we found that the locally supplied land use data are more consistent with long-term averaged vegetation cover estimated from remote sensing images. As locally supplied land use data only covers part of the Basin, the BRS 1 km land use map (BRS 2001) was used for the rest of the areas. Details of land use data sources are given in the Appendix I. All the landuse data are re-classified into 22 classes given in the Table 1 for sheetwash and rill erosion rate calculations. 3.2
Hillslope Sheet and Rill Erosion under Natural Condition
To better understand the relative impact of land use and management practices on hillslope erosion, the predicted sheetwash and rill erosion needs to be put in the context of erosion under natural vegetation cover. We predicted natural erosion using a similar procedure as for modelling soil erosion under current land use conditions with a cover factor for native vegetation, keeping the other factors as for the present day. An empirical modelling framework to predict the pre-European settlement (undisturbed) USLE C-factor was implemented. Instead of directly using the results from NLWRA sediment delivery and transport project Theme 5.4b, the modelling work was redone for this project using an updated current C factor and an improved technique for sampling remnant native vegetation. There are two basic assumptions of this modelling framework: 1) climate, soil type, geology and terrain conditions remain unchanged since European settlement; 2) the natural vegetation and soil surface conditions remains similar to pre-settlement condition for those areas with limited disturbance. Based on those two assumptions, we sampled the C-factor from those areas with limited disturbance, built statistical models using climate, soil, geological and terrain variables as predictor variables and used the models to extrapolate to those areas with substantial disturbance by human intervention, especially agricultural activities such as cropping, grazing and tree clearing. Statistical models were constructed using the Cubist data mining tool (Rulequest Research 2001) in a similar way as we used for the predictions of hillslope length and slope (Lu et al. 2003b). In this study, we reserved a proportion of the sample set to test the model, calculating statistics of model performance for both the model-built data and test data sets. 50% of the total sampling points were used for model building and the other 50% of points for model testing. The sampling and modelling were carried out at 0.05 degree resolution. A stepwise model building approach was used. For the first step, each predictive variable was used independently and the best variable was identified on the basis of correlation coefficient and relative error. This variable was then combined with every other variable, to find the second variable that most improved the model. This procedure was repeated until all variables were included. Final selection of the model was based on the statistical diagnostics, and visual comparisons of predicted and measured maps.
15
The predictive variables for modelling the C-factor under pre-European conditions using Cubist were selected to represent the major factors presumed to determine vegetation cover and soil distribution across the continent. They can be broadly grouped into four categories: natural vegetation; soil parent material; climate; and geomorphology. Specifically, the following nineteen predictive variables were selected: (1) Australia - Natural Vegetation (Carnahan, J.A. and AUSLIG (1989) 1:5 M scale); (2) aggregated geology classifications derived from the 1:2.5M scale geology map of Australia; (3) the Australian Soil Classification derived from the Atlas of Australian Soils; (4) mean annual temperature, mean diurnal change, isothermality, temperature seasonality and diurnal temperature range; (5) mean annual rainfall, rainfall seasonality index, annual moisture index and moisture index seasonality; (6) mean annual radiation and radiation seasonality; and (7) 9” DEM, averaged slope and slope length derived from 9” DEM and relief, and their scaled estimations (Gallant 2001). 3.3 3.3.1
Sediment Delivery Ratio (SDR) Background of SDR
Soil erosion models, such as the Universal Soil Loss Equation (USLE) (Wischmeier and Smith 1978) estimate gross soil erosion rate at plot-scale. Erosion rates estimated by USLE are often higher than those measured at catchment outlets. Sediment delivery ratio (SDR) is used to correct for this reduction effect. It is defined as the fraction of gross erosion that is transported from a given area in a given time interval and it is a measure of sediment transport efficiency which accounts for the amount of sediment that is actually transported from the eroding sources to a measurement point or catchment outlet compared to the total amount of soil that is detached over the same area above that point. Mathematically, it is expressed as SDR =
Y E
(3.1)
where Y is average annual sediment yield per unit area and E is average annual erosion over that same area. It compensates for areas of sediment deposition that become increasingly important with increasing catchment area, and therefore, determines the relative significance of sediment sources and their delivery. Factors influence SDR including hydrological inputs (mainly rainfall), landscape properties (e.g., vegetation, topography, and soil properties) and their complex interactions at the land surface. The multitude of such interactions makes it difficult to identify the dominant controls on catchment sediment response and on catchment-to-catchment variability. In reality, erosion is not normally measured directly. It is measured as sediment yield at a small scale, such as a hillslope plot. Thus, SDR is a scaling factor used to accommodate differences in arealaveraged sediment yields between measurement scales. Physically, it stands as a mechanism for compensating for areas of sediment deposition that becoming increasingly important with increasing catchment area. Therefore, transport and storage lie in the heart of SDR. At regional scale, the most widely used method to estimate SDR is through a SDR-area power function: SDR = α Aβ
(3.2)
16
where A is the catchment area (in km2), the constant α and a scaling exponent β are empirical parameters (Maner 1958; Roehl 1962). Field measurements suggest that β is in the range 0.01 to -0.025 (Walling 1983; Richards 1993), which means that SDR decreases with increasing catchment area. The scaling exponent β contains key physical information about catchment sediment transport processes and its close linkage to rainfall-runoff processes. It seems that β decreases with increasing aridity (Richards 1993). Lower value of β (up to –0.7) were found in the Sicilian region and in former USSR catchments (Ferro and Minacapilli 1995). Field data (Figure 2) show that the relationships between SDR and drainage area changes considerably between different catchments over the world. Extrapolation of those empirical relationships can be misleading and results in SDR exceeding 100%.
Figure 2: SDR vs catchment area relationships obtained from different areas around the world.
For catchments with similar area, field data show the values of α and β in equation (3.2) are also different in different regions (Walling 1983; Roehl 1962). It is because the SDR-area relationship does not take into account local descriptors, such as rainfall, topography, vegetation, land use and soil characteristics. There are other empirical relationships which show that SDR varies with various physiographic attributes but the data that went into these relationships are few and of only local extent (Khanbilvardi and Rogowski 1984). This limits the usefulness of such a lumped empirical approach. Williams (1977) developed a procedure for determining SDR based on runoff models for small catchments. Recent development in this direction is towards the spatially distributed modelling using GIS techniques (Ferro 1995). There are other methods to predict sediment delivery and deposition through calculation of sediment transport capacity, avoiding the need for a lumped SDR (Morgan et al. 1998; Van Rompaey et al. 2001). Although those methods were based on improved physical understanding of sediment transport processes, they require high resolution DEMs to route the flow and sediment. They also rely on detailed sediment transport or runoff data to calibrate parameters, such as the sediment transport capacity coefficient. However, such methods often require many parameters which are generally too expensive or even impossible to determine reliably and the input data such as hydraulic resistance, infiltration rate, and soil properties including particle size distributions are not commonly available over large spatial extents. The traditional SDR methods are often data-driven. They depend on the existence of long periods of sediment yield records at the stream gauging stations and a sensible measure or estimation of hillslope erosion rate. However, there are few consistent long periods of
17
sediment yield data available in the MDB to allow such an analysis to be carried out. In addition, approaches based on analyzing sediment yield records cannot identify the separate effects of changing climate, land use and management practices on sediment delivery as catchment response to change is often longer than the record length. It is known that there are some limitations of SDR methods (Walling 1983; Richards 1993). One is that SDR methods cannot explicitly predict the locations and rates of sediment deposition in the lowland phases, and another is the problem of temporal and spatial lumping and lack of physical basis. However, SDR is a very useful concept to model regional scale sediment delivery processes. It avoids the need to explicitly model patterns of deposition on hillslopes which is not possible across such large areas as the MDB. There is little quantitative SDR information available within the Basin for the scale we are interested here. Existing measurements are either at much smaller or larger catchment scales. Studies based on sediment budgets carried out in forest areas of south-eastern New South Wales (NSW) and the East Gippsland show that the values of SDR are in the range of 10% 45% for catchment areas of about 2 km2 (Croke et al. 1999) and range from 2% - 95% for those sub-catchments (with areas of around 100 km2) within the Bega Catchment (Fryirs and Brierley 2001). A SDR of 70% was found for the Upper Wolumlar Creek (area = 18 km2), located in the South Coast of New South Wales (Brierley and Fryirs 1998). For the catchment area in which we are interested, most of the measurements and studies were carried out in humid areas outside of the Basin. Little has been done for the arid and semiarid regions. In summary, SDR is the result of numerous complex interactions among hydrological inputs (mainly rainfall) and landscape properties (e.g., vegetation, topography, and soil properties) through a number of hydrological processes at the land surface. The multitude of such interactions makes it difficult to identify the dominant controls on catchment sediment response and on catchment-to-catchment variability within the MDB. In addition, field measurement of SDR is severely limited. Therefore, it is difficult to model spatially distributed SDR accurately. 3.3.2
A New SDR Theory
One important aim of this study is to develop a SDR model that incorporates the key elements of the catchment storm response and sediment delivery process. Sivapalan et al. (2001) showed that the interactions between time scales, namely between rainfall duration and catchment response lay at the heart of the regional flood frequency estimations. The way that catchment response time varies with catchment area depends on the relative dominance of hillslope response, channel hydraulic response, and network geomorphology. A simple linear model of catchment response (Sivapalan et al. 2001) is used in this study. Instead of using the model for studying catchment response of flood, we use the same concept to model SDR. The model consists of two independent components: sediment transport on hillslopes and sediment routing in the channel network. As shown in Figure 3, these are represented through two linear stores, arranged in series. The hillslope store is supplied with sediment by soil eroison at a rate e [mass/area/time] over an effective storm duration ter (erosion only occurs during this time period). The hillslope stores part of the eroded sediment and delivers the rest to the channel network store, located downstream of it, at a rate yh [mass/area/time]. yh is assumed to be a linear function of the mass of sediment stored in the hillslope per unit area, denoted by Sh [mass/area]. The area specific sediment yield from the 18
network store, y [mass/area/time], which is the same as the area specific sediment yield from the catchment outlet, is assumed to be a linear function of the sediment stored in the channel network, denoted by Sn [mass/area]. The continuity equation of sediment for the two stores can be expressed as: dS h (t ) = e(t ) − yh (t ) dt yh (t ) = S h (t ) / th
(3.3)
dS n (t ) = yh (t ) − y (t ) dt y (t ) = S n (t ) / tn
where th is the mean hillslope residence time and tn is the mean channel residence time. e(t)
Sh(t)
yh (t) =
Sh (t) th
Hillslope Storage Sn(t)
y(t) =
Sn (t) tn
Channel Storage Figure 3: Diagram of a two storage lumped linear model of SDR at catchment scale (after Sivapalan et al. 2001, modified). See text for detail.
For simplicity, we assume that the upland erosion rate e is constant during ter. Equation (3.3) can then be solved analytically. The final expressions for the ratio between the peak of the resulting sedigraph, denoted by yp [mass/area/time] (which is equal to max(y)), and upland erosion rate e can be written as follows: ter 1 − exp − e tn t t − h 1 − exp − er tn − th th
yp
=
tn t n − th
1 ter2 1 ter3 = − + ... e 2 tn2 3 tn3
yp
th ≥ 0
tn ≠ th
(3.4)
t n = th
On an event basis, we assume SDR = y p / e . The peak sediment yield Yp [mass/time] can be estimated by multiplying area specific sediment yield yp [mass/area/time] by the catchment area A. Equations were firstly derived by Sivapalan et al. (2001) for studying the scaling effects on regional flood frequency under different rainfall and catchment conditions.
19
Sivapalan et al. (2001) showed that equation (3.4) is capable of explaining the power law relationship between flow response and catchment area and changing value of the scaling exponent which is caused by a change of hydrological processes. Similarly, equation (3.4) can be used to explain the obtained SDR vs area relationships. As shown in Figure 4, SDR measurements gathered by Roehl (1962) in several American catchments including Blackland Prairies, the Red Hills of Texas and Oklahoma, the Missouri Basin Loess Hills, the Mississippi Sand Clay Hill, and the South-eastern Piedmont (shown in dots) suggested that, in general, SDR decreases with catchment area. The solid line, which is the average flow response (the scaling factor of mean flood discharge defined as the ratio between average rainfall input rate and runoff at the catchment outlet during flood events) calculated using the equation (28a) of Robinson and Sivapalan (1997), represents the upper envelope of SDR. The averaged modeled SDR estimated by equation (3.4) is shown as the dashed line. The reason that SDR is often smaller than flow response is due to the settling velocity of soil particles (compared with water particles) and other effects such as sediment transport capacity. For a given catchment area, the variations in SDR measurements (by up to two orders of magnitude) are due to heterogeneity in catchment properties (e.g. rainfall, catchment slope and curvature, soil texture, etc). The combination of the above physical properties results in differences in the time variables ter, tn and th in eq. (2). Therefore, eq. (2) can be used to model spatially distributed SDR if the time variables ter, tn and th can be spatially differentiated. 1000 SDR (Roehl 1962) SDR (modelled) Flow Response
SDR
100
10
1 0.01
0.1
1
10 Area (km 2)
100
1000
Figure 4. Comparison of SDR (%) measurements (Roehl 1962), modeled average SDR and flow response (Robinson and Sivapalan 1997). It shows that flow response represents the upper envelope of the SDR.
Equations (3.3) and (3.4) can be used to compute the magnitudes of the SDR for different values of the timescales ter, th and tn. The results are presented in the Figure 5 (upper panel) as families of curves relating SDR to tn for different values of ter and th. They show that the SDR remains constant for small values of tn, while for larger values of tn they decrease linearly with increasing values of tn (scaling exponent is -1). The effect of th is to smooth and reduce the magnitude of SDR, without changing the scaling exponents with respect to tn .
20
1
SDR
0.1
0.01
0.001
0.0001 0.1
1
10
100
1000
10000
tn (hrs) ter= 1, th = 50 ter= 5, th = 50 ter= 10, th = 50
ter= 1, th = 5 ter= 5, th = 5 ter= 10, th = 5
ter= 1, th = 0 ter= 5, th = 0 ter= 10, th = 0
1
Measurements
SDR
0.1
0.01
0.001
0.0001 1.E-03
1.E-01
1.E+01
1.E+03
1.E+05
1.E+07
1.E+09
A (km2)
Figure 5: SDR as a function of channel residence time for different values of ter and th (upper panel); SDR as a function of catchment area for different values of ter and th. SDR measurements from USA catchments (Roehl 1962) are also shown as red dots (lower panel).
Equations (3.3) and (3.4) can also be used to explain the SDR and area relationship observed from measurements. Assume that th is independent of catchment area for the scale of catchment shown in Figure 5 and tn can be expressed as a function of catchment size A of the form: tn = α Aβ
(3.5)
where tn is in hours, A is in km2. Assuming that parameters α = 0.76 and β = 0.38 (ARR 1987) Figure 5 (lower panel) shows SDR versus catchment area relationships. It shows that for a given catchment area, SDR values vary by three orders of magnitude due to the differences in ter and th. While the SDR, in the sense of equations (3.3) and its solution (3.4) remains linear for all catchment sizes, the observed change in scaling exponent is cause by a change of hydrological 21
processes. In small catchments, effective rainfall duration is long compared to the catchment residence time, and consequently the sediment transport reaches steady state, with the whole catchment area contributing to the sediment yield. As long as this remains true, the sediment yield increases linearly with catchment area with the exponent at unity. On the other hand, in large catchments, effective storm duration is smaller than the catchment residence time. The fraction of the catchment area contributing to the sediment yield is proportional to the ratio of effective storm duration to catchment residence time. This ratio decreases at the rate of A-β with an increase of catchment area A. Thus the partial area contributing to sediment yield increases only at the rate of A1-β , with exponent 1-β less than unity. The above analysis was based on a single storm. The derived flood frequency method can be used to deal with multiple storms. In this study, for simplicity, we treat effective storm duration ter as a random variable and calculate th and tn as catchment averaged values. According to equation (3.4), SDR becomes a random variable. By knowing the probability distribution of ter, we can derive a probability distribution of SDR. According to standard statistical procedures, we can calculate the mean values of SDR. To apply the model to the whole MDB, we need to estimate two groups of input variables: 1) statistical properties of effective rainfall duration and intensity using pluviograph rainfall data; and 2) hillslope and channel residence time th and tn as a function of particle size and other catchment properties. The following sections describe the procedures for estimating the two groups of input variables. 3.4
Statistical Analysis of Effective Rainfall duration and Intensity
Rainfall varies both in space and time. Both rainfall duration and intensity are important and interrelated, and have significant impacts on sediment generation at the point scale (source) and sediment transport at catchment scale. The main aim of analyzing high temporal resolution rainfall data is to discover the possible controls on the spatial variability of sediment delivery due to temporal variability of rainfall intensity. $ $
$ $
$
$
$
$
$
$ $ $
$ $
$
$
$
$
$
$
$ $ $$ $$$ $ $ $$ $ $$ $ $$ $ $#$ $ $ $ $ $$ $$$$Brisbane $ $ $ $$ $ $ $ $ $ $ $ $ $ $$ $ $ $ $ $$ $ $ $
$ $$ $ $ $ $ $$ $$ $ $ $$ $$$ $$ $ $ $ $$ $ $ $ $ # Sydney $$ $ $ $$ $ $ $#$ $ $ Canberra $ $ $ $ $ $ $ $$ $ $ $$ $ $ $ $ $ $ $ $$$ $$ $ $$ $ $ $$ $ $ $$$$ $ $ $$ $ $ $$$$$# Melbourne $$$ $
$
$
Adelaide $ #$ $
$
$
$
$
$ $ $
$ $$
Figure 6: Site locations of pluviograph rainfall data and their relative position to MDB.
22
Pluviograph rainfall data with 6-min sampling interval was collected for 195 sites from the Bureau of Meteorology (BoM). The site selection is based on two conditions: 1) that sites are within or nearby the MDB; and 2) they have at least 10 years of rainfall record. Among those 195 sites, four sites which have a large number of missing values were eliminated. The locations of the gauges and their relative position to the MDB are shown Figure 6. More detailed site information can be found in Table AII.1 of Appendix II. The analyses of rainfall data is divided in two parts: 1) a statistically analysis of rainfall data at a single site to search for a suitable probability distribution function; and 2) regionalization of the parameters of the suitable probability distribution function. Search for a suitable probability distribution function
The storm events, characterized by their intensity and duration, are estimated by statistical analyzing pluviograph rainfall data, with temporal resolution of 6 minute interval. In this study, the storm events are defined as rain periods separated by dry periods of at least 6 hours or longer. Once a storm is defined, the 30 minute rainfall intensity and storm duration can be estimated. The events which have total rainfall depth equal or larger than 12.7 mm are considered to be potentially harmful in terms of erosion and sediment transport. Therefore, only those events are included in the calculation. The value of 12.7 mm is chosen here to be consistent with that used in the Universal Soil Loss Equation (USLE) (Wischmeier and Smith 1978). In addition, an event which has total rainfall depth equal or larger than 12.7 mm does not necessarily cause erosion during its whole rainfall duration. For those events with low intensity but long duration, the effective duration in terms of causing erosion is shorter than the rainfall duration. To consider this effect, we calculate the effective duration for an event for a given site using the following equation
ter ,i
R = min N i tr ,i , tr ,i ∑ Rj j =1
(3.6)
N
where Ri is the R factor for i-th event,
∑R j =1
j
is the sum of R factor for all the events, and tr,i is
the duration for i-th event. Equation (3.6) simply states that the effective duration is shorter for events with smaller erosivity compared to events with larger erosivity, and the effective duration cannot exceed the actual rainfall duration. The computer program of Yu and Rosewell (1998) was modified for the calculations of rainfall duration and intensity. Figure 7 to Figure 9 show the event 30-intensity and duration for the Wagga Wagga site (upper panels). Figure 7 shows all the rainfall events. It shows that many events have small 30-min intensity and duration. Those events have little effect on sediment generation and transport and can be excluded from the analysis. This exclusion is done by only including the events which have total rainfall depth equal or greater than 12.7 mm as suggested by USLE. Figure 8 shows the events after excluding those small events. The events in which the effective duration is calculated using equation (3.6) are shown in Figure 9.
23
Figure 7: All rainfall events characterised by their 30 intensity and duration (upper panel); Fit probability density functions to Gamma and exponential distributions for both duration and intensity (second and lower panels).
24
Figure 8: Rainfall events which have depth equal or greater than 12.7 mm (upper panel); Fit probability density functions to Gamma and exponential distributions for both duration and intensity (second and lower panels).
25
Figure 9: Effective rainfall events which have depth equal or greater than 12.7 mm (upper panel); Fit probability density function of effective duration to Gamma and exponential distributions (lower panel).
Figure 7 to Figure 9 show that the Gamma distribution fits the data better than the exponential distribution for both intensity and duration (lower two panels). However, due to the awkwardness of operating with the Gamma distribution in analytical form and the reasonable fit of the exponential distribution for effective rainfall for larger events, we decided to use the exponential distribution in this study. The exponential distribution density function is written as f ( x) =
x exp − λ λ 1
(3.7)
where x is the random variable, λ is the mean value of the random variable x (in this study, x can be 30-min rainfall intensity, rainfall duration or effective rainfall duration). For each rainfall site, the mean value λ is obtained by averaging values across all the events considered. To apply the statistical estimation of SDR using the exponential distribution (3.7) to the whole Basin, we need to regionalize the mean value λ for both intensity and duration. Regionalization of the Mean Values of Intensity and Duration
To apply the SDR model cross the Basin, we need to regionalise the mean values of rainfall intensity and duration by linking them to the existing climatic surfaces. The regionalization is done as follows.
26
As shown in Figure 10, a regression relationship with r2 = 0.98 is obtained between MI30 and the ratio between mean annual R factor (R) and mean annual rainfall (MAR). Good relationships are obtained for the mean values of rainfall duration (tr) and the effective duration (ter) in relation to the ratio between mean annual rainfall (MAR) and mean 30min maximum intensity (MI30) (Figure 11 and Figure 12, respectively). As we have already regionalized R and MAR, those relationships shown in Figures 10 to 12 are used to regionalise tr, ter and MI30. Note that only the events with rainfall depth equal or greater than 12.7 mm are considered in the calculation.
Figure 10: Relationships between effective 30-min. rainfall intensity and the ratio between mean annual R-factor and mean annual rainfall.
Figure 11: Relationships between rainfall duration (tr) to mean annual rainfall (MAR), effective 30-min intensity (MI30), MAR/MI30, and MAR2/R.
27
30
25
25
20
20
15
15
10
10
5
5
t er
(hr)
30
0
0
1000
500
1500
0
2000
0
MAR (mm) 30
10
5
20
15
MI30 (mm) 1
y = 0.057 x (r 2 = 0.84)
2
25 0.5
(hr)
Relative Error
20
t er
15
0
10 -0.5 5
100
200
300
400
MAR/MI30
500
600
-1
0
5
10
20
15
t er
25
30
(hr)
Figure 12: Relationships of effective rainfall duration and it relative errors.
Figure 13 shows the errors between rainfall duration and rainfall duration estimated using regionalisation equations derived in this study for all the rainfall sites. It suggests that the prediction accuracy increases as the number of years with complete records increases. This is consistent with the USLE which suggests at least 22 years pluviograph records are appropriate for long term erosion estimation. Relative larger errors occur mainly at the sites which have shorter rainfall record. Larger errors are observed at some sites within the high rainfall regions (Australian Alps), which might suggest more complex and non-linear rainfall patterns in those regions.
28
t r (Simulated)
40 30 20 10 0
0
10
5
15
20
30
25
35
40
t r (Calculated from Pluviograph data)
Absolute Error (hr)
10 5 0 -5 -10 0
10
20
30
40
50
60
50
60
Number of Year with Complete Data
Relative Error
1 0.5 0 -0.5 -1
0
10
20
30
40
Number of Year with Complete Data Figure 13: Error estimations of rainfall duration. Upper Panel: Comparison between rainfall duration estimated using site specific pluviograph data and that estimated using regionalised relationships. Middle Panel: Absolute error [hrs] plotted against number of year with complete data. Lower Panel: Relative error plotted against number of year with complete data. The crosses are the sites have shorter records and relatively larger errors. They are not used in the final relationships that are applied across the MDB.
29
3.5
Estimations of Residence Time
3.5.1 Sediment Residence Time as a Function of Particle Size The residence time of sediment can be estimated as a function of sediment particle size and the travel time of water particles.
Trajectory of Water Particle
Fine Particle Suspension Silt Particle Suspension and Saltation Sand Saltation
Figure 14: Diagram of the particle size effect on sediment travel time in relation to the travel time of water particles.
Suppose we can estimate the travel time of water particles as a function of local slope, roughness, rainfall intensity, etc. The effect of sediment particle size can be reflected as shown in the diagram in Figure 14. Very small clay particles, which are characterized by their slower settling velocity, remain suspended in the water most of the time and their trajectories of travel differ little to that of water particles. Silt particles with faster settling velocity travel with water particles during high velocity flows and settle to the soil bed during low flows. Large sand particles saltate near the soil bed with slow overall velocity. The travel time for different size particles within flowing water was modelled as follows:
th (d ) = th 0 Fh (d ) tn (d ) = tn 0 Fn (d )
(3.8)
where th(d) and tn(d) are the hillslope and channel residence time for particles with diameter d, respectively, and th0 and tn0 are the hillslope and channel travel time of water particles, respectively. Fh(d) and Fn(d) are the enlargement functions describing the influence of particle size d. The mathematical forms of Fh(d) and Fn(d) were modelled as: Fh = exp ( γ h wt (d ) ) Fn = exp ( γ n wt (d ) )
(3.9)
where wt(d) is the settling velocity for particles with diameter equal to d, and γh and γn are the parameters inversely relating to water depth. In general, γh is larger than γn as the typical water depth in overland flow is of order of millimetres and the water depth in small channels is of order of centimetres. The settling velocity was calculated as:
30
1/ 2
4 ρ p gd wt (d ) = 3ρ C (Re ) D p
(3.10)
where ρp is the particle density, ρ is the water density, g is acceleration due to gravity, Rep = wtd/ν is the particle Reynolds number at the settling velocity, and CD is the drag coefficient modelled as a function of the particle Reynolds number Rep:
CD (Re p ) =
24 1 + 0.15 Re0.687 ( ) p Re p
(3.11)
(Durst et al. 1984). Finally, SDR was calculated for each particle size group and then weighted by the particle size distribution to get an overall SDR as follows: N
SDR = ∑ wi SDRi i =1
N
∑w i =1
i
(3.12)
=1
where N is the total number of particle groups, wi and SDRi are the mass percentage and SDR for particle group i, respectively. Three particle size groups are considered in this study. These are: d ≤ 2 µm (clay), 2 ≤ d ≤ 20 µm (silt), and 20 ≤ d ≤ 1000 µm (sand). Particles with diameter larger than 1000 µm are considered too large to be transport far away from their source areas. The mass percentage of each particle size group was estimated using the Australian Soil Resource Information System (ASRIS) product (Carlile et al. 2001). 3.5.2 Estimating Travel Time of Water Particles th0 and tn0 Novotny and Olem (1994) pointed out that land cover and slope are the key factors in affecting sediment delivery rates. Additionally, they stated the importance of factors specific to storm events, such as rainfall intensity, infiltration, ponding, and overland flow energy. However, because this research used average annual erosion rate, consideration of detailed infiltration, ponding and storm specific factors was not feasible. The travel time of water particles is calculated separately for overland flow and stream flow. The travel time is inversely related to flow velocity. During a storm event when overland flow occurs, the flow carries sediment from surface runoff until it reaches a stream. In the stream component, the runoff water is influenced by a different set of factors affecting travel-time compared to that of the overland component. To capture this, the travel time of channel flow is calculated from each cell in the catchment to the outlet by aggregating stream segments. Along each path, travel-time is calculated by aggregating time taken within each cell using procedures described below.
31
Overland flow component:
For the hillslope cells, the overland flow velocity is estimated by combining a kinematic wave approximation with Manning’s equation. The depth of flow at equilibrium (m) is given by (Overton and Meadows 1976): ni L y = 0e.5 s
0.6
(3.13)
where L is the travel distance along the flow path (m), n is Manning’s roughness coefficient, ie is the rainfall excess rate (mm/s), and s is the decimal slope. By substituting the depth of flow at equilibrium in Manning’s equation, the velocity of overland flow (m/s) can be calculated as: (ie L) 0.4 s 0.3 n 0.6
Vo =
(3.14)
The travel time (s) through each cell can be estimated as: to =
D Vo
(3.15)
where D is the distance travelled through that cell (m). For orthogonal flow, the flow distance is the cell width, while for diagonal flow, it is equal to 2 D.
To calculate travel time by implementing the above procedure, four input parameters are needed in the overland component: rainfall excess rate ie, Manning’s coefficient n, flow travel length L, and slope s. Estimations of those input parameters are made as follows. Estimation of Excess Rainfall Rate ie: Excess rainfall generated in a catchment is known to vary spatially. The variation in excess rainfall follows that of land use, land cover, and soil type. Typically, the way to account for this variation is to divide the catchment into smaller areas of “uniform” land use, land cover, and soil type combinations. An average curve number (CN) for the whole catchment determined using the area weighting method is then given by: CN =
CN1 A1 + CN 2 A2 + .... + CN m Am m
∑A
(3.16)
i
i =1
where CNi is the curve number of the sub-area i (with area equal to Ai ). m is the total number of sub-areas. This procedure is the standard procedure used in the USDA SCS rainfall-runoff relationship (SCS 1983). It gives an average excess rainfall depth for the entire catchment, Pe that corresponds to an average rainfall depth, P. The equations used to calculate Pe are: Pe =
( P − 0.2 S ) 2 P + 0.8S
(3.17)
where S is the storage term (in mm) which can be obtained using the formula:
32
S=
254 CN
× 100 − 254
(3.18)
where CN is the curve number that can be obtained from standard tables for different combinations of land use and land cover, soil hydrologic group, treatment, and conditions. The hydrologic soil group reflects soil permeability and surface runoff potential. Following is a description of the four different hydrologic soil groups: Group A are soils with low total surface runoff potential due to their high infiltration rates. They consist mainly of excessively drained sands and gravels. Group B are soils with low to moderate surface runoff potential. They have moderate infiltration rates and moderately fine to moderately coarse texture. Group C are soils with moderate to high surface runoff potential. They have slow infiltration rates and moderately fine to fine textures. Group D are soils with high surface runoff potential. They have very slow infiltration rates and consist chiefly of clay soils.
Typical values of CN for certain land use groups are given in Table 2. Table 2: Typical values of CN for some land use group.
Sources of CN: SCS (1983; 1986), Novotny and Olem (1994)
33
Once the spatially-distributed CN map is developed, the total storage can be obtained by equation (3.18). The excess rainfall equation (3.20) gives the accumulated depth of excess rainfall from the start of the storm to the current time. For an unsteady rainfall/flow event, the incremental value of excess rainfall of a time interval ∆t, ie can be calculated as the difference between the accumulated excess rainfall at the end of that time interval and the accumulated excess rainfall at the beginning of the that same interval as follows: ie (t ) = Pe (t ) − Pe (t − 1)
(3.19)
In this study, we assume steady-state rainfall, we calculate ie as: ie =
Q tr
(3.20)
where tr is the rainfall event duration and Q = Pe is the total excess rainfall for the event. The procedure of estimating tr has been given in Section 3.4. Estimation of Manning’s Roughness Coefficient, n: For simplicity, Manning’s n roughness coefficient is estimated using available land use and landcover data. Table 3 shows estimated typical values of n for overland flow. Table 3: Values of Manning’s n used in this study for common land use and vegetation cover groups for overland flow.
Veg. cover (cv)
cv ≤ 30%
30% < cv ≤ 70%
cv > 70%
Annual (not managed) Pasture
0.15
0.4
0.6
Sow (improved) Pasture
0.15
0.4
0.6
Crop
0.15
0.25
0.4
Forest
0.2
0.6
0.8
Built-up areas
0.1
0.3
0.5
Wetland and ponds
0.125
0.125
0.125
Land use
Estimation of Travel Length and Slope
These parameters can be extracted from a digital elevation model (DEM) by using a geographic information system (GIS). As the 9” DEM is used in this study, the resolution of the DEM is not capable of capturing the overland flow path in detail. In the implementation of equations (3.14) and (3.15), flow length L and the distance of travel D are approximately equal to hillslope length, which is a product of NLWRA sediment transport and delivery project (Gallant 2001). Like hillslope length, slope grid s is also a product of NLWRA sediment transport and delivery project. Both grids were statistically derived using higher resolution DEM, 9” DEM and other climatic, geology, and soil attributes (Gallant 2001).
34
Channel component:
The travel time in the channels can be calculated based on the SCS flow velocity equation (Haan et al. 1994) Vch = as1/ 2
(3.21)
where Vch is flow velocity [m/s], s is the slope [m/m], and a is a coefficient relating to stream roughness condition. Landuse
Soil hydrologic Group
Manning’s n
Representative Rainfall Event
DEM
Slope
Curve Number Flow direction Rainfall Excess Volume Flow Accumulation Rainfall excess intensity Delineate Channel Network
Flow length
Flow Velocity Overland component Channel Component
Calculate travel time for each cell by dividing the travel distance by the flow velocity
Calculate the cumulative travel time Figure 15: Flow chart for the calculation of travel time of water particles.
Similar to overland flow, the travel time (tc) through each channel cell can be estimated as: tc =
D Vch
(3.22)
where D is the distance travelled through that cell (equal to horizontal, vertical or diagonal distance across a cell flow direction). For a given cell i, the cumulative travel time was estimated by summing the travel time along its flow path. More specifically, if a sediment particle in cell i travels through mo cells overland and mc cells in the stream to reach the catchment outlet, equations (3.14) and (3.15) were used in each of the mo upland cells to calculate the concentrated shallow flow travel time and equation (3.22) was used in each of the mc stream cells and aggregated to estimated total
35
stream flow time (Tic). Figure 15 shows the overall procedure for calculating travel-time th0 and tn0. Two input parameters are needed for the channel component: slope s and channel roughness parameter a. The channel roughness parameter a is parameterised as in Table 4. Table 4: channel roughness parameter a values used in this study.
Channel Section
Upstream Area (ha)
A
Concentrated shallow flow
1.8 – 18
4.0
Intermittent stream (grass waterway)
18 – 360
4.5
Permanent Stream (little cover)
360 and up
5.0
4
Results
4.1
Hillslope Erosion under Current Land Use
Figure 16 shows the predicted sheet and rill erosion across the Basin. In general, it is predicted that erosion rate increases from south to north and from west to east. The major source areas are: Brigalow Belt, New South Wales South West Slopes and north part of Darling Riverine Plains. It was predicted that about 2.1 × 108 tones of soil is moved annually on hillslopes over the Basin. The average erosion rate across the basin is 2.1 t ha-1 yr-1. If we denote that a pixel with soil loss rate below 0.5 t ha-1 yr-1 as low erosion, larger than 10 t ha-1 yr-1 as high erosion, and in between as medium, it is estimated that about 40% of the Basin experiences low erosion, 4% faces high erosion and 56% of the Basin experiences medium hillslope erosion. Agricultural lands in steep and higher rainfall intensity areas experience higher erosion rate than other land use groups, showing the potential to target erosion control. Table 5 shows the percentage erosion for those three groups and in relation to percentage agricultural lands in each group. Table 6 divides hillslope erosion into land use classes. In general, the average erosion rate is higher for agricultural lands compared with other land use groups though many of agricultural lands are located in floodplains where the slope is low. The rates are higher compared with surrounding non-cropping areas where other conditions are similar. This confirms that land use and management practices have a major impact on soil erosion.
36
Figure 16: Estimated annual average sheet and rill erosion rate.
Table 5: Three erosion groups (high, medium and low) and their relation to percentage of agricultural lands.
High Erosion Rate (> 10 t/ha/year) Percentage in Basin Area (%) Area (103 km2) 4 42 Area (103 km2) Percentage of Agr.Lands (%) 21 8.7 Medium Erosion Rate (0.5 - 10 t/ha/year) Percentage in Basin Area (%) Area (103 km2) 566 56 Percentage of Agr.Lands (%) Area (103 km2) 19 110 Low Erosion Rate ( < 0.5 t/ha/year) Area (103 km2) Percentage in Basin Area (%) 400 40 Percentage of Agr.Lands (%) Area (103 km2) 14 54 *Total Basin area: 108 million ha. Lakes and reservoirs are not included for erosion statistics. Agricultural lands in the Basin: 17.3 million ha (around 16% of the Basin area).
Figure 17 shows the monthly distribution of total soil loss. It is found that over 75% of the erosion occurs in the summer period, especially in the north part of the Basin. However, the high erosion zone detected within the east part of the Basin, shows weaker summer dominance due to its temperate climate condition.
37
Table 6: Soil loss rate from land use categories. Landuse
Approx. Total Area
Total Erosion
Ave. Erosion Rate
Rate of acceleration
Group
(km2 * 10^3)
(Mt yr-1)
(t ha-1 yr-1)
ratio of current and natural rates
196
15.4
0.79
2.5
Woodland
2178
156.8
0.72
2.8
Plantation
139
12.9
0.92
3.5
Forest
267
14.9
0.56
3.7
Residual/Native Pastures
4239
1017.6
2.40
4.2
Improved Pastures/Legumes
201
22.1
1.10
5.9
Cereals excluding Rice
199
66.7
3.35
9.8
Other agricultural lands
17
6.6
3.88
8.1
45 40 35 30 25 20 15 10 5
De c
No v
O ct
Se p
Au g
Ju ly
Ap ril M ay Ju ne
M ar ch
Fe b
0 Ja n
Total Soil Loss (Mt/month)
National park
Month
Figure 17: Monthly distribution of total soil loss rate for the Basin.
4.2
Hillslope Erosion under Natural Conditions
The best predicting variable is the mean annual soil moisture index, which explains 72% of the variance of the sampled data, followed by the annual mean radiation explaining 59% of the variance of the sampled data. The correlation to polygon based data, such soil, and geology are lower (around 9% to 20%) and they improve little to the overall needed correlation coefficient or spatial patterns of predicted maps. The final predicting variables are clim1, clim2, clim3, clim4, clim7, clim12, clim15, clim20, clim23, clim28, clim31 and austdem. Figure 18 shows the comparison between samples of C values extracted from the current C map for those undisturbed (or minimum disturbed) points and modelled C values using Cubist for the testing data from the final model.
38
5 Discussions and Conclusions Hillslope sheetwash and rill erosion:
There are three dominant forms of water-borne erosion in the Murray-Darling Basin. These are sheetwash and rill erosion (sometimes termed hillslope erosion in the reports produced by this project), the formation and erosion of gullies, and the erosion of riverbanks. In this study, new assessments of hillslope erosion across the MDB were reported, building upon our previous work for the National Land and Water Resources Audit (NLWRA) (Lu et al. 2001; Lu et al. 2003b). Improvements to the assessment of sheetwash and rill erosion were made by compiling higher resolution land use data for the MDB from a range of sources and by incorporating a database on crop rotation, tillage and other land management practices. These new data, together with improved analysis of remote sensing data, enabled a more accurate prediction of the effect of vegetation cover and cover management on hillslope erosion. It is estimated that 2.2 × 108 tonnes of sediment were moved in the MDB annually as hillslope erosion at a mean rate of 2.1 t ha-1 yr-1. Erosion rate increases from south to north and from arid areas to temperate regions with most of the erosion generated from the east and north part of the MDB. Under any given rainfall regime, the reduction of protective ground cover increases the risk of high soil losses. About two-thirds of erosion occurs in the summer period. Agricultural lands have relatively high erosion rates and higher increment of soil erosion rates. Very low soil erosion rates are estimated under pre-European natural vegetation conditions. The rates are 3 – 10 times on average and up to 100 times smaller than that under current land use. A New Theory for Modelling Spatially Distributed Sediment Delivery Ratio:
In this report, we have proposed a theory for sediment delivery ratio and implemented the theory across the MDB in a spatially distributed manner. Spatially, sediments are produced from different sources distributed throughout the Basin. Each source is characterized by its sediment detachment, transport and storage. The SDR model argues that sediment delivery can be closely linked to temporal hydrological control. For each source area, SDR is characterized by two important time variables, namely, its travel time, i.e., the time that particles eroded from the source area and transported through the hillslope conveyance system take to arrive at the channel network and eventually to the catchment outlet, and the typical rainfall duration, which is the primary driving force of sediment transport. For instance, for the same rainfall event, we expect that a source area with a shorter travel time would have a higher SDR. Alternatively, for the same source area, a rainfall event with shorter duration would have a lower SDR as less eroded particles would make their way to the catchment outlet. Those particles will be stored (or deposited) somewhere in the system. These interactions between rainfall attributes (including intensity, duration and intermittency) and catchment characteristics are important factors for understanding spatially distributed sediment delivery. For the arid part of the Basin, rainfall events are often smaller in size spatially with shorter duration but more intense than in humid temperate climates. For a given slope steepness and slope length, local erosion rate in the arid areas is relatively high due to insufficient vegetation cover and relatively more intensive rainfall. However, the sediment delivery follows a different pattern. The shorter rainfall duration and larger variations in interannual rainfall also cause a greater variation in sediment transport. The sediment yield differs from one catchment to another depending upon whether the storm duration is larger or smaller than the sediment residence time (SRT) of the catchments. The SDR model proposed
44
in this study is able to differentiate the catchments for which storms usually last longer than the SRT or for those for which residence time is seldom met. The SDR model allows quantitative estimates of the non-linear effects on sediment delivery due to changes in climate and land use. It expresses the spatial variability of catchmentaveraged SDR in terms of the statistical time variables and particle size distributions. It relies on rainfall intensity (6-min interval) and daily rainfall records (which cover larger areas) instead of stream flow records. It offers a means to understand the dominate processes which control sediment delivery. The model has a simple analytical form which can be implemented in a GIS environment. Applying the model to the MDB, we found: 1) sediment delivery ratio and sediment yield are low for most parts of the Basin except some upland areas in the east and north part of the Basin; 2) the sediment transport can be very effective at sub-catchment level, especially in the areas of the Australia Alps, South West Slopes, Brigalow Belt South, and Darling Downs; 3) only about 5% of sheet and rill erosion are transported from sub-catchment elements in to the streams. The average area specific sediment yield at sub-catchment element level is around 0.1 t ha-1 yr-1. About 14 million tones in total of sediment generated from sheet and rill erosion is delivered from the sub-catchment elements to the major streams. The quantitative, spatially distributed estimations of SDR have important implications not only for the study of off-site environment impact due to exported sediment but also to on-site erosion control. It has been demonstrated that there is economic advantage from identifying the areas that have a higher potential to deliver sediment and prioritizing control implementation in those areas (Dickinson et al. 1990). The spatially distributed SDR map contributes to the development of cost-effective strategies for erosion control (Lu et al. 2003a). SDR and Sediment Yield due to Hillslope Erosion:
In summary, it is found that sediment delivery ratio and sediment yield are low for most part of the Basin except some sloping land in the eastern part of the Basin. Estimated at Basin outlet, spatial patterns of topography, rainfall intensity and rainfall duration suggests the system is not effective in terms of sediment transport. However, the sediment transport can be very effective at sub-catchment scale, such as in the areas of South West Slopes, Brigalow Belt South, and Darling Down regions. Average area specific sediment yield from subcatchment is 0.13 t ha-1 yr-1. On average, about 5% of sheet and rill erosion is transported from sub-catchment elements in to the streams. In total, around 14 million tonnes per year of sediment generated from hillslope sheet and rill erosion is delivered from the sub-catchment elements to the major streams. The Australian Alps have relatively high sediment delivery ratio though the local erosion rate is medium. Caution is recommended for any vegetation clearance in those ranges.
Acknowledgments We acknowledge with appreciation the financial support for this work from the Murray Darling Basin Commission, and the personal support from the MDBC Office particularly from Ms. Lisa Robins. The project steering committee led by Dr. Pat Feehan are thanked for their efforts and time to oversee and to guide the progress of the work. 45
We thank Bofu Yu for supplying the source code of RECS, which was modified to calculate 30-min rainfall intensity and duration for each site with rainfall record. We also thank the local agencies who kindly supplied us the land use data. The interactions and discussions with colleagues within and outside the team are gratefully acknowledged. Individuals include Elisabeth Bui, Greg Cannon, Francis Chiew, Barry Croke, Mick Fleming, John Gallant, Tony Jakeman, Russell Mein, Neil McKenzie, David Simon, David Smiles, and Bill Young.
46
References ARR (1987) Australian Rainfall and Runoff. A Guide to Flood Estimation Volume 1, D.H. Pilgrim (ed.), Institution of Engineers, Australia. Brierley, G.J. and Fryirs, K. (1998) A fluvial sediment budget for upper Wolumla Creek, South Coast, New South Wales, Australia, Australian Geographer, 29, 107-124. BRS (2001) Land use mapping at catchment scale: Principles, procedures and definitions. Bureau of rural Sciences, Canberra. Butler, B.E., Blackburn, G. and Hubble, G.D. (1983) Murray-Darling Plains (VII). In: (Eds.), Soils: An Australian viewpoint. pp. 231-240. Butler, B.E. and Hubble, G.D. (1978) The general distribution and character of the soils in the Murray-Darling River system. Proceedings of the Royal Society of Victoria. 90, 149-156. Carlile P, Bui, E., Moran, C., Minasny, B., and McBratney A.B. (2001) Estimating soil particle size distributions and percent sand, silt and clay for six texture classes using the Australian Soil Resource Information System point database. Technical Report 29/01, CSIRO Land and Water, Canberra. Carnahan, J.A. and Australian Surveying and Land Information Group (1989) Australia - Natural Vegetation. Crabb, P. (1997) Murray-Darling Basin Resources. Murray-Darling Basin Commission, Canberra. Croke, J., Wallbrink, P., Fogarty, P., Hairsine, P., Mockler, S., McCormack, B. and Brophy, J. (1999) Managing Sediment Sources and Movement in Forests: The Forest Industry and Water Quality, Cooperative Research Centre for Catchment Hydrology Industry Report 99/11. DeRose R.C., Prosser I.P., Weisse M., and Hughes A.O. (2003) Summary of sediment and nutrient budgets for the Murray-Darling Basin. Technical Report K to Murray-Darling Basin Commission, Basin-wide mapping of sediment and nutrient exports in dryland regions of the MDB (Project D10012). Dickinson, W.T., Rudra, R.P. and Wall, G.J. (1990) Targeting remedial measures to control nonpoint source pollution. Water Resour. Bull. 26, 499-507. Durst, F., Milojevic, D. and Schönung, B. (1984) Eulerian and Lagrangian predictions of particulate two-phase flows: a numerical study. Appl. Math. Modelling, 8, 101-115. Eagleson, P.S. (1972) Dynamics of flood frequency, Water Resour. Res. 8(4), 878-898. Edwards. K. (1987) Runoff and soil loss studies in New South Wales. Technical Handbook No. 10, Soil Conservation Service of NSW, Sydney, NSW. Edwards, K. (1988) How much soil loss is acceptable? Search 19, 136-140. Edwards, K. (1993) Soil erosion and conservation in Australia. In: World soil erosion and conservation, Ed: Pimentel D, Cambridge, 147-169. ESRI (2003) ESRI GIS software, World Wide Web http://www.esri.com/software/index.html. ESRI, Redlands, California. Ferro, V. and Minacapilli, M. (1995) Sediment delivery processes at basin scale. Hydrological Sciences Journal 40, 703-716. Fryirs, K. and Brierley, G.J. (2001) Variability in sediment delivery and storage along river courses in Bega catchment, NSW, Australia: implications for geomorphic river recovery, Geomorphology 38, 237-265. Freebairn, D.M. and Wockner, G.H. (1986). A study of soil erosion on vertisols of the eastern Darling Downs, Queensland. 1. Effect of surface conditions on soil movement within contour bays. Australian Journal of Soil Research 24, 135-58. Gallant, J.C. (2001) Topographic scaling for the NLWRA sediment project. Technical Report 27/01, CSIRO Land and Water, Canberra. Haan, C.T., Barfield, B.J. and Hays. J.C. (1994) Design Hydrology and Sedimentology for Small Catchments. Academic Press, New York. Hughes, A.O., and Prosser, I.P. (2003) Gully and Riverbank Erosion Mapping for the Murray-Darling Basin. CSIRO Technical Report 3/03, CSIRO Land and Water, Canberra. 20pp. http://www.clw.csiro.au/publications/technical2003/tr3-03.pdf Hutchinson, M.F., Stein, J.A. and Stein, J.L. (2001) Upgrade of the 9 Second Australian Digital Elevation Model. http://cres.anu.edu.au/dem. Khanbilvardi, R.M. and Rogowski, A.S. (1984) Quantitative evaluation of sediment delivery ratios. Water Resour. Bull. 20, 865-874. Loughran, R.J. and Elliott, G.L. (1996) Rates of soil erosion in Australia determined by the caesium-137 technique: a national reconnaissance survey. In: Erosion and Sediment Yield: Global and Regional Perspectives. IAHS Publication. 236, 275-282.
47
Lu, H., Gallant, J., Prosser, I.P., Moran, C., and Priestley, G. (2001) Prediction of Sheet and Rill Erosion Over the Australian Continent, Incorporating Monthly Soil Loss Distribution. Technical Report 13/01, CSIRO Land and Water, Canberra, Australia. Lu, H., Moran, C. J., DeRose, R. and Cannon, G. (2003a) Spatially Distributed Investment Prioritization for Sediment Control in the Murray Darling Basin. Technical Report, (In press)/03. Lu, H., Prosser, I.P., Moran, C.J., Gallant, J. Priestley, G. and Stevenson, J.G (2003b) Predicting sheetwash and rill erosion over the Australian continent, Australian Journal of Soil Research. (revised) Lu, H., Raupach, M.R., McVicar, T.R. Barrett, D.J. (2003c) Decomposition of Vegetation Cover into Woody and Herbaceous Components Using AVHRR NDVI Time Series. Remote Sensing of Environment 5850, 1-18. Lu, H. and Yu, B. (2002) Spatial and seasonal distribution of rainfall erosivity in Australia, Australian Journal of Soil Research 40, 887-901. Maner, S.B. (1958) Factors affecting sediment delivery rates in the Red Hills physiographic area. Trans. Am. Geophys. 39, 669-675. Martin, R.J., MCMillan, M.G. and Cook, J.B. (1988) Survey of farm management practices of the northern wheat belt of New South Wales. Australian Journal of Experimental Agriculture 28, 499-509. Morgan, R.P.C., Quinton, J. N., Smith, R.E., Govers, G., Poesen, J.W.A., Auerswald, K., Chisci, G., Torri., D. and Styczen, M. E. (1998) The European soil erosion model (EUROSEM): A dynamic approach for predicting sediment transport form fields and small catchments. Earth surf. Process. Landforms 23, 527544. Newcombe, C.P. and MacDonald D.D. (1991) Effects of suspended sediments on aquatic ecosystems. North American Journal of Fisheries Management 11, 72-82. NLWRA (2001) Australian Agriculture Assessment. National Land and Water Resources Audit, Canberra. Novotny, V. and Olem, H. (1994) Water Quality-Prevention, Identification, and Management of Diffuse Pollution. Van Nostrand Reinhold, New York. Overton, D.E. and Meadows, M.E. (1976) Stormwater modeling, Academic Press, New York. Parsons, A.J. and Stromberg, S.G.L. (1998) Experimental analysis of size and distance of travel of unconstrained particles in interrill flow Water Resour. Res. 34, 2377-2381. Pimentel, D., Harvey, C., Resosudarmo, P., Sinclair, K., Kurz, D., McNair, M., Crist, S., Shpritz, L., Fitton, L., Saffouri, R. and Blair, R. (1995) Environmental and economic costs of soil erosion and conservation benefits. Science 267, 1117-1123. Prosser, I.P., Rustomji P., Young, W.J., Moran, C.J. and Hughes, A.O. (2001) Constructing river basin sediment budgets for the National Land and Water Resources Audit. CSIRO Land and Water Technical Report 15/01, CSIRO Land and Water, Canberra, Australia. Quinlan, J.R. (1993) C4.5: Programs for Machine Learning, Morgan Kaufmann Publishers, San Mateo, California. . For further information see: http://www.rulequest.com. Renard, K.G., Foster, G.A., Weesies, D.K., McCool, D.K., and Yoder, D.C. (1997) Predicting soil erosion by water: A guide to conservation planning with the Revised Universal Soil Loss Equation. Agriculture Handbook 703, United States Department of Agriculture, Washington DC. Richards, K. (1993) Sediment delivery and the drainage network, In: Channel Network Hydrology, Eds. Beven, K. and Kirkby, M.J., 221-254. Robinson, J.S. and Sivapalan, M. (1997) An investigation into the physical causes of scaling and heterogeneity in regional flood frequency, Water Resour. Res., 33, 1045-1059. Roehl, J.E. (1962) Sediment source areas, and delivery ratios influencing morphological factors. Int. Assoc. Hydro. Sci. 59, 202-213. Rosewell, C.J. (1997) Potential source of sediments and nutrients: Sheet and rill erosion and phosphorus sources. Australia: State of the Environment Technical Paper Series (Inland Waters), Department of the Environment, Sport and Territories, Canberra. Rowe, R.K., Crouch, R.J. and van Dyke, D.C. (1978) Soils in the upper valleys of the Murray River basin. Proceedings of the Royal Society of Victoria 90, 157-166. Rulequest Research (2001) Rulequest Research Data Mining Tools. World Wide Web http://rulequest.com. Rutherfurd, I.D., Ladson, A., Tilleard, J., Stewardson, M., Ewing, S., Brierley, G. and Fryirs, K. (1998) Research and development needs for river restoration in Australia. LWRRDC Occasional Paper No. 15/98, NSW. SCS (1983) TR-20 Computer Program for Project Formulation Hydrology, U. S. Department of Agriculture, Soil Conservation Service, Washington, D. C., May, 1983. SCS (1986) Urban Hydrology for Small Watersheds, Technical Release No. 55, U. S.Department of Agriculture, Soil Conservation Service, Washington, D. C., June,1986. Sivapalan, M., Jothityangkoon, C. and Menabde, M. (2001) Linearity and non-linearity of basin response as a function of scale: Discussion of alternative definitions. Water Resour. Res. 24, 1001-1014.
48
Swift, R. and Skjemstad, J. (2002) Agricultural land use and management information. National Carbon Accounting System, Technical Report No. 13. Vanclay, F. (1997) Land degradation and land management in central NSW. Charles Sturt University, Wagga Wagga). Van Rompaey, A.J.J., Verstraeten, G., Van Oost, K., Govers, G. and Poesen, J. (2001) Modelling mean annual sediment yield using a distributed approach, Earth surf. Process. Landforms 26, 1221-1236. Walling, D.E. (1983) The sediment delivery problem, Journal of Hydrology 65, 209-37. Walling, D.E. (1988) Erosion and sediment yield research - some recent perspectives. Journal of Hydrology 100, 113-141. Wasson, R.J., Olive, L.J. and Rosewell, C. (1996) Rates of Erosion and Sediment Transport in Australia. In: Erosion and Sediment Yield: Global and Regional Perspectives. D.E. Walling and R.Webb (eds) IAHS Publ. Williams. J. R. (1977). Sediment delivery ratios determined with sediment and runof models. In: Erosion and solid matter transport in in land waters. IAHS-AISH publication, 122: 168-179. Wischmeier, W.H. and Smith, D.D. (1978) Predicting rainfall erosion losses: A guide to conservation planning. US Department of Agriculture. Agriculture Handbook No. 537. (US Government Printing Office: Washington, DC). Yu, B. and Rosewell, C.J. (1998) A program to calculate the R-factor for the USLE/RUSLE using BoM/AWS pluviograph data. ENS Working Paper 8/98.
49
Appendix I: Land Use Data
Table AI.1 shows the details of the land use data supplied status and contact details of the agencies. Figure AI.1 shows the land use extent used in this study. Areas shown in blue have land use supplied by local agencies. Table AI.1: Summary of locally supplied land use data used in this study.
Region Murrumbidgee (New South Wales)
Barwon region (New South Wales)
Murray region (New South Wales) South-east South Australia (New South Wales) Golbourn region (New South Wales)
Condamine (Queensland) Upper Billabong* (New South Wales) Bendigo Region (Victoria)
Christian Writte
Data Supply Status Data supplied by a CD with 1:100 000 map sheets of the Murubidgee region. CD includes Arcview shape files, metadata and readme documents. A CD of Landuse mapping was supplied. It contains 1:100 000 landuse mapping of Barwon region done during late 1980s, Showing timber, pasture, and cropping lands. It also contains 1: 50 000 landuse mapping of eastern part of Walgett shire and all of Moree shire (19982000). Data supplied for Upper Murray and Billabong regions. A CD with landuse data was supplied.
Email address Rob Brownbill
[email protected]
Sally Keane Resource Officer (GIS) (02) 69230437 Angela McCormack
[email protected]
Stuart Lucas
[email protected]
David Tonkin
[email protected]
Michael Htun GPO Box 409 Canberra, ACT 2601 Ph: (02) 6279 0100 Fax: (02) 6248 8053
A CD was supplied by MDBC. It contains land use data for Golbourn, Condamine and Upper Billabong regions. Products of MDBC “Project D2006” Land mark Task 6a. See above.
See above.
See above.
See above.
Maree Platt Current landuse mapping for
[email protected] Victoria is as same as BRS 1:250 000 landuse. Just starting a landuse project for NC CMA but won’t be finished until next financial year. No data supplied.
David Burton
50
(Queensland)
GIS/Drafting Officer(Graphics Unit) Department of Natural Resources and Mines
[email protected]
Far West of New South Wales
We were told that detailed landuse data does exist for the region but they were a bit hesitant to give it to us because the information was quite sensitive. Despite assuring that the data would be used in the strictest of confidence, no data supplied to us. Central West Region No data supplied. of New South Wales
Aaron Colbran
[email protected]
Michael Casey
[email protected]
* The Murrumbidgee data includes 90% of the Billabong data. The Murrumbidgee landuse data is used for the overlapped part.
Figure AI.1: Data sources of land use used in this project.
51
Appendix II: Pluviograph Rainfall Data Table AII.1: Details of the pluviograph rainfall sites. SA Pluviograph Stations Site
Name
Lat
Lon
Start
End
Years
% A
21060
JAMESTOWN DPI
-33.2042
138.6011
Oct 1951
Aug 1998
45.3
92 N
23034
ADELAIDE AIRPORT
-34.9581
138.5342
Jan 1967
Apr 2001
30.2
83 Y
23090
ADELAIDE (KENT TOWN)
-34.9211
138.6216
Feb 1977
Jan 2001
24
96 Y
23763
MOUNT CRAWFORD FOREST
-34.7139
138.9453
Jan 1971
Jan 2001
29.3
86 N
24023
LOXTON RESEARCH CENTRE
-34.4333
140.6
Nov1972
May 1984
11.6
94 N
24515
LANGHORNE CREEK
-35.2958
139.033
Jan 1973
Jun 2000
26.5
92 N
Lat
Lon
Start
End
Years
% A
QLD Pluviograph Stations Site
Name
35000
ALPHA POST OFFICE
-23.6497
146.6411
May 1963
Jun 1981
14
68 N
35025
DINGO POST OFFICE
-23.645
149.3303
Apr 1963
Jul 1981
18.2
96 N
35069
TAMBO POST OFFICE
-24.8819
146.2564
Aug 1963
Jan 1999
33.9
90 N
35098
EMERALD DPI TOWN SITE
-23.5
148.15
Aug 1962
Oct 1985
23.3
96 N
35104
KILMACOLM
36031
LONGREACH AERO
-22.4
147.5333
Sep 1963
Jun 1981
17.2
94 N
-23.4372
144.2769
Mar 1966
Oct 2000
34.5
96 Y 85 N
36047
TWIN HILLS POST OFFICE
-21.95
146.9517
Dec 1965
Jun 1981
14.9
39006
BILOELA DPI
-24.3789
150.5164
Nov 1937
May 1996
55.3
88 N
39061
MANNERSLEY
-24.0428
150.8211
Apr 1975
Dec 1993
17.7
86 N
39083
ROCKHAMPTON AERO
-23.3769
150.4761
Nov 1939
Oct 2000
61
94 Y
39090
THEODORE DPI
-24.9503
150.0725
Aug 1951
Jun 1981
29.4
93 N
39104
MONTO POST OFFICE
-24.8667
151.1239
Apr 1963
Jun 1992
29.2
93 N
39123
GLADSTONE RADAR
-23.8558
151.2617
Dec1961
Jun 1993
31.6
92 Y
39128
BUNDABERG AERO
-24.8885
152.3235
Dec 1963
Jan 1999
17.8
47 Y
39171
NARAYEN RES STN
-25.6875
150.8689
Nov 1969
Jun 1981
11.1
90 N
40004
AMBERLEY AMO
-27.6294
152.7114
Oct 1961
Jul 2000
38.3
93 Y
40019
BENARKIN FOREST STATION
-26.9
152.15
Jul 1961
Jun 1981
19.4
87 N
40059
COOROY COMPOSITE
-26.4181
152.9128
Nov 1971
Apr 2001
28.2
86 N
40062
CROHAMHURST
-26.8094
152.87
Oct 1960
Jun 1981
12.8
56 N
40082
UNI. OF QUEENSLAND GATTON
-27.5508
152.3358
Jun 1956
May 2000
41
83 N
40093
GYMPIE
-26.1831
152.6414
May 1963
May 1981
17
90 Y
40102
JIMNA COMPOSITE
-26.6656
152.4594
Feb 1972
Jan 2000
12.3
40 N
40112
KINGAROY PRINCE STREET
-26.5544
151.8456
May 1963
Oct 2000
21
52 N
40126
MARYBOROUGH
-25.5181
152.7111
Aug 1964
Jun 1981
16.5
87 N
40133
MONSIDALE
-26.7
152.4
Aug 1963
Dec 1977
13.8
91 N
40135
MOOGERAH DAM
-28.0317
152.5517
Oct 1964
Jun 1997
29.9
84 N
40160
NERANG GILSTON RD
-28.0092
153.3175
Nov 1970
Feb 2001
15.3
45 N
40189
SOMERSET DAM
-27.1169
152.555
Nov 1936
Nov 1969
23.7
70 N
40192
SPRINGBROOK FORESTRY
-28.2264
153.2786
Oct 1965
Jun 1983
17.7
95 N
40197
MT TAMBORINE FERN ST
-27.9697
153.195
May 1972
Apr 1999
21.5
69 N
40214
BRISBANE REGIONAL OFFICE
-27.4778
153.0306
Jan 1908
Jun 1994
84.4
94 N
40223
BRISBANE AERO
-27.4178
153.1142
May 1949
Feb 2000
50.7
94 Y
40270
RAVENSBOURNE
-27.3628
152.1594
Sep 1965
Jun 1981
15.1
84 N
40282
NAMBOUR DPI
-26.6431
152.9392
Jan 1954
Mar 1999
45
88 Y
40318
KIRKLEAGH
-27.0258
152.5642
Nov 1959
Jun 1990
30.7
94 N
40382
CROWS NEST
-27.2639
152.055
Jul 1965
Jun 1981
15.9
88 N
52
40386
KENILWORTH BRIDGE
-26.5892
152.7322
Jun 1963
Jun 1981
16.4
40406
BEENLEIGH BOWLS CLUB
-27.7092
153.2011
Dec 1968
Jul 2000
12.4
83 N 35 N
40584
HINZE DAM
-28.0483
153.2878
Sep 1974
Apr 1999
19.2
74 N
40606
UPPER MUDGEERABA WATER
-28.1056
153.3289
Nov 1974
Apr 1999
19
71 N
40608
BENOWA WATER TREAT
-28.0039
153.4008
Nov 1974
Jul 1990
15.3
82 N
40609
ELANORA WATER TREAT
-28.1181
153.4456
Nov 1974
Apr 2000
19.7
72 N
40677
MAROON DAM
-28.1753
152.6553
Dec 1977
Apr 2000
20.4
85 N
41044
HERMITAGE
-28.2061
152.1003
Feb 1952
May 2000
47.2
89 N
41060
LEYBURN
-28.0092
151.5861
Mar 1959
Nov 1996
32.2
77 N
41140
WAMBO SHIRE COUNCIL
-27.1864
151.2553
Mar 1959
Dec 1991
32.3
95 N
41175
STANTHORPE (GRANITE BELT HRS)
-28.6214
151.9528
Nov 1965
Apr 1995
28.4
86 Y
41359
OAKEY AERO
-27.4036
151.7414
Dec 1990
Nov 2000
10
96 Y
41467
TOOWOOMBA CITY COUNCIL
-27.5667
151.885
Jan 1957
Dec 1983
24.3
81 N
42016
HANNAFORD POST OFFICE
-27.3408
150.0622
Jul 1969
Mar 1999
20
63 N
43020
MITCHELL POST OFFICE
-26.4908
147.9778
Aug 1963
Jul 1999
32.2
84 N
43044
AMOOLEE FOREST R 238
-26.6383
149.4261
Jan 1967
Jul 2000
18.5
52 N
44021
CHARLEVILLE AERO
-26.4131
146.2611
Jan 1953
Oct 2000
46.7
94 Y
44026
CUNNAMULLA POST OFFICE
-28.0706
145.6808
Jan 1980
Jul 1999
19.5
95 N
45015
QUILPIE AIRPORT
-26.6122
144.2578
Aug 1963
Oct 2000
33.3
84 N
Lat
Lon
Start
End
Years
% A
NSW Pluviograph stations Site
Name
48027
COBAR MO
-31.4853
145.8292
Jun 1962
Feb 2000
37.6
94 Y
50102
CONDOBOLIN SOIL CONSERVATION
-33.0833
147.15
Jul 1957
Dec 1974
17.4
96 N
51049
TRANGIE RESEARCH STATION AWS
-31.9861
147.9489
Aug 1968
Aug 1998
20.1
62 Y
52069
PILLIGA (RIVERVIEW)
-30.2747
148.8222
Dec 1970
Jul 1983
12.6
92 N
53048
MOREE COMPARISON
-29.4819
149.8383
Apr 1964
Jun 1995
31.2
95 N
54036
WALLANGRA (WALLANGRA STATION)
-29.2443
150.8922
Jun 1954
Dec 1969
14.9
90 N
54102
BARRABA (ROSEVALE)
-30.3735
150.6723
Jan 1971
Sep 1999
24.5
81 N
54104
PINDARI DAM
-29.3946
151.2398
Jan 1980
Sep 1999
19.7
88 N
54105
BUNDARRA (GRANITE HEIGHTS)
-30.3367
150.9333
Jan 1975
Sep 1999
24.5
89 N
54138
UPPER HORTON (DUNBEACON)
-30.156
150.3889
Nov 1976
Aug 2000
23.7
92 N
55024
GUNNEDAH SCS
-31.0261
150.2687
Apr 1946
Oct 2000
51.5
89 N
55031
MANILLA POST OFFICE
-30.7478
150.7196
Jan 1953
Dec 1969
11.8
56 N
55054
TAMWORTH AIRPORT
-31.0867
150.8467
Aug 1958
Dec 1992
33.9
96 N
55136
WOOLBROOK (DANGLEMAH ROAD)
-30.9672
151.3451
Jan 1971
Sep 1999
24.3
77 N
55194
GOWRIE NORTH
-31.3365
150.8537
Jan 1971
Mar 2000
25.3
75 N
56016
GUYRA POST OFFICE
-30.2217
151.67
Apr 1959
Nov 1972
10.4
68 N
56018
INVERELL RESEARCH CENTRE
-29.7767
151.0806
Aug 1947
Oct 2000
50.8
89 N
56041
BONSHAW (MONKSTADT)
-29.1333
151.45
May 1954
Dec 1969
15.4
88 N
56059
TENTERFIELD (COOREDULLA)
-29.05
152.1
Jan 1956
Nov 1974
18.3
83 N
56224
GLEN INNES SCS
-29.7
151.7
Oct 1947
Dec 1973
25.8
91 N
57033
WOLLOMOMBI POST OFFICE
-30.5167
152.05
Apr 1959
Dec 1982
21.6
75 N
57058
CARRAI STATE FOREST (DAISY PLAINS)
-30.9367
152.2917
May 1963
Sep 1983
17.8
82 N
57059
STYX R.STATE FOREST
-30.5517
152.2783
Oct 1963
Sep 1983
17
80 N
57104
YARROWITCH (MARETTO)
-31.2739
151.9655
Apr 1959
Aug 1999
19.3
42 N
57105
WALCHA (BULIMBA DOWNS)
-31.0667
151.9167
Apr 1959
Dec 1974
15.2
87 N
61029
KULNURA (WILLIAM ROAD)
-33.2333
151.2
Feb 1969
Aug 1981
12.5
90 N
61078
WILLIAMTOWN RAAF
-32.7939
151.8386
Dec 1952
Oct 2000
45
87 Y
61089
SCONE SCS
-32.0632
150.9272
Jul 1952
May 1999
44.2
85 N
61142
BUCKETTY
-33.1167
151.1333
Jan 1959
Jun 1969
10.3
95 N
61151
CHICHESTER DAM
-32.2426
151.683
Jun 1960
Dec 1980
20.1
91 N
61152
CONGEWAI (GREENOCK)
-32.9995
151.2908
Feb 1959
May 1971
12.2
91 N
53
61158
GLENDON BROOK (LILYVALE)
61171
JERRYS PLAINS (CARRINGTON)
-32.5069 -32.5167
151.3756
Dec 1980
12.4
69 N
Dec 1980
19.7
80 N
61174
MILLFIELD COMPOSITE
-32.9
151.2667
Oct 1958
Jun 1981
20.4
86 N
61178
WOLLOMBI (BIG YENGO LTP)
-32.9333
61181
BROKE (OAKLEY)
-32.75
150.9167
Mar 1958
Nov 1975
15.4
82 N
151.1667
Jan 1959
Jun 1969
10.1
91 N
61193
WOLLOMBI (STOCKYARD CREEK)
61209
PUTTY TEA ROOMS
-32.9
151.0833
Aug 1959
Dec 1969
10.3
97 N
-32.9614
150.6742
Nov 1962
May 1983
17.3
61211
COLO HEIGHTS (THE MILE RIDGE)
73 N
-33.3256
150.7042
Nov 1962
Jun 1999
20.3
49 N
61212 61223
LIDDELL (POWER STATION)
-32.3767
150.96
Jan 1965
Dec 1986
15.4
62 N
MARYVILLE
-32.9131
151.75
Jan 1964
Sep 1991
26.7
80 N
61238
POKOLBIN (SOMERSET)
-32.8126
151.3043
Aug 1965
Jun 1981
15.9
89 N
61240
WOLLOMBI (BLAIR)
-32.9667
151.1333
Sep 1955
Jul 1973
17.8
95 N
61287
MERRIWA (ROSCOMMON)
-32.1897
150.1728
Mar 1969
Dec 1986
12.8
68 N
61288
LOSTOCK DAM
-32.3283
151.4583
Oct 1969
Sep 1999
14.3
43 N
61309
MILBRODALE (HILLSDALE)
-32.6881
150.9728
Jan 1970
Jun 1981
11
86 N
61343
SCONE SCS.2.
-32.0667
150.9333
Oct 1952
Nov 1970
15
78 N
62005
CASSILIS POST OFFICE
-32.0067
149.98
Jan 1975
May 2000
14.7
52 N
62020
BYLONG (MONTORO)
-32.5014
150.0333
Feb 1965
Jun 1989
17.3
67 N
62026
RYLSTONE (ILFORD RD)
-32.8073
149.9768
Sep 1955
Dec 1973
18.3
91 N
63023
COWRA RESEARCH STN
-33.8088
148.7072
Oct 1941
Sep 1999
37.6
61 N
63035
HILL END POST OFFICE
-33.0362
149.4146
Sep 1959
Jan 1975
15.1
85 N
63039
KATOOMBA (NARROW NECK RD)
-33.7135
150.2983
Jun 1965
Jun 1999
13.1
34 N
63108
OBERON DAM
-33.7167
149.8667
Jan 1955
Jun 1989
22.7
59 N
63253
ORANGE (ROSETEAGUE)
-33.3167
149.05
Aug 1955
Jun 1973
17.4
87 N
64009
DUNEDOO POST OFFICE
-32.0163
149.3953
Sep 1959
Jan 1975
14.5
76 N
64046
COONABARABRAN (WESTMOUNT)
-31.2886
149.0687
Jul 1971
Aug 2000
16.3
52 N
65035
WELLINGTON RESEARCH CENTRE
-32.5059
148.9708
Feb 1961
Aug 2000
39.3
94 N
70012
BUNGONIA (INVERARY PARK)
-34.8996
149.9709
May 1965
Aug 2000
22.8
58 N
70014
CANBERRA AIRPORT
-35.3049
149.2014
Dec 1937
Oct 2000
47.7
71 Y
150.9667
Jan 1965 Sep 1958
70015
CANBERRA FORESTRY
-35.3
149.1
Jan 1932
Feb 1971
36.3
82 N
70025
CROOKWELL POST OFFICE
-34.4572
149.469
Feb 1956
Oct 1974
17.8
74 N
70080
TARALGA POST OFFICE
-34.4048
149.8197
Jun 1977
Jul 2000
10.1
39 N
70099
CANBERRA (ACTON)
-35.3
149.1
Jan 1921
Dec 1939
14.8
71 N
70282
CANBERRA CITY
-35.2667
149.1167
Dec 1974
Nov 1988
14
92 N
71010
KIANDRA CHALET
-35.8833
148.5
May 1957
Dec 1967
10.7
92 N
71063
GUTHEGA DAM SMHEA
-36.3833
148.3667
Dec 1957
Jun 1969
11.5
95 N
72023
HUME RESERVOIR
-36.104
147.0329
Mar 1955
Jun 2000
29.2
63 N
72056
BLOWERING DAM
-35.3947
148.247
Mar 1955
Dec 1973
18.1
84 N
72091
CABRAMURRA SMHEA
-35.9383
148.3842
Mar 1956
Oct 1974
18.6
92 N
72112
VALENTINES HUT
-36.2333
148.3667
Jan 1958
Dec 1968
10.3
84 N
72116
JAGUNGAL SMHEA
-36.1333
148.3833
Apr 1959
Sep 1974
15.4
93 N
72150
WAGGA WAGGA AMO
-35.1583
147.4573
Jan 1945
Feb 2000
37.8
63 Y
73007
BURRINJUCK DAM
-35.0008
148.5969
May 1911
Aug 2000
50
51 N
74039
DENILIQUIN FALKINER MEMORIAL
-35.3667
145.05
Jul 1950
Nov 1977
25.9
86 N
74114
WAGGA WAGGA RESEARCH CENTRE
-35.1311
147.3091
Sep 1946
Jul 1990
38.2
81 N
75050
NARADHAN (URALBA)
-33.6104
146.3161
Apr 1970
Aug 2000
14.7
45 N
Lon
Start
End
Years
% A
VIC Pluviograph Stations Site
Name
Lat
76031
MILDURA AIRPORT
-34.2306
142.0839
Apr 1953
Nov 2001
45.3
88 Y
77087
HOPETOUN RWC
-35.7253
142.3656
Mar 1958
Dec 1991
29.7
81 N
79046
WARTOOK RESERVOIR
-37.0944
142.4322
May 1974
Dec 2001
20.9
63 N
79052
ROCKLANDS RESERVOIR
-37.2311
141.9594
Jan 1955
Nov 2001
39.8
75 N
54
79079
ST ARNAUD (TOTTINGTON)
-36.7928
143.1194
Feb 1973
Nov 2001
24.2
75 N
79082
HORSHAM
-36.7022
142.2017
Jan 1958
Nov 1997
33.3
73 N
79086
SUPPLE (AVON NO.3)
-36.8653
143.1214
Jan 1973
Jan 1998
24.8
86 N
80006
CHARLTON POST OFFICE
-36.3
143.4
Sep 1951
Jul 1963
10.1
82 N
80067
CHARLTON
-36.2714
143.3453
Sep 1951
Dec 2001
45.4
82 N
80102
PYRAMID HILL
-36.0597
144.115
Jan 1973
Oct 1986
13.8
95 N
80109
COBRAM (GOULBURN MURRAY)
-35.9133
145.6431
Sep 1957
Nov 2001
43
90 N
80110
GOULBURN-MURRAY WATER (KERANG)
-35.7372
143.9286
Oct 1957
Nov 2001
33.3
68 N
81003
BENDIGO PRISON
-36.7533
144.2825
Feb 1968
Jul 1992
24.3
89 N
81013
DOOKIE AGRICULTURAL COLLEGE
-36.3714
145.7036
Jan 1950
Dec 2001
50.1
89 N
81026
LAANECOORIE WEIR
-36.835
143.89
Aug 1973
Oct 2001
23
74 N
81038
NATTE YALLOCK
-36.9439
143.4697
May 1974
Nov 2001
23
78 N
81049
TATURA INST SUSTAINABLE AG
-36.4381
145.2672
Jul 1960
Nov 2001
24.7
54 Y
81114
TATURA THEISS ENVIRON SERV
-36.3997
145.3325
Jan 1975
Oct 1993
15.6
76 N
81115
WANALTA DAEN STATION
-36.6281
144.8725
Jul 1974
Nov 2001
22.7
78 N
82011
CORRYONG (PARISH LANE)
-36.2014
147.8942
Feb 1972
Nov 2001
29.7
88 N
82016
EUROA
-36.7564
145.5717
Dec 1967
Dec 2001
30.3
80 N
82039
RUTHERGLEN RESEARCH
-36.1064
146.5083
Feb 1975
Nov 2000
19.4
68 Y
82042
STRATHBOGIE
-36.8469
145.7308
Jan 1974
Nov 2001
23.9
80 N
82076
DARTMOUTH RESERVOIR
-36.5364
147.4967
Jul 1975
Oct 2001
21.7
75 N
82107
GOULBURN (LAKE NILLAHCOOTIE)
-36.8564
146.0031
Jul 1968
Oct 2001
33.1
90 N
82121
OVENS RIVER (WANGARATTA)
-36.35
146.3417
Aug 1957
Oct 1993
28.4
70 N
82138
WANGARATTA AERO
-36.4181
146.3056
May 1987
Jul 2001
10.8
70 Y
83017
JAMIESON
-37.3028
146.1392
Jan 1984
Nov 2001
13.9
74 N
83025
OMEO
-37.1011
147.5981
Mar 1985
Nov 2001
12.1
64 N
83031
WHITFIELD
-36.7531
146.4139
Oct 1962
Oct 1991
28.9
91 N
83033
WOODS POINT
-37.5644
146.2467
Aug 1954
Nov 2001
16.8
28 N
83067
BRIGHT
-36.7331
146.9603
Jun 1969
Dec 1994
20.9
69 N
83074
LAKE WILLIAM HOVELL RESERVOIR
-36.9153
146.3858
Jan 1988
Dec 1997
10
89 N
87017
BLACKWOOD (POST OFFICE)
-37.47
144.3069
Feb 1974
Oct 2001
26.4
77 N
87029
LANCEFIELD
-37.2708
144.7217
Jan 1929
Jun 1975
44.8
86 N
87031
LAVERTON RAAF
-37.8636
144.7458
Mar 1965
Aug 1999
31.9
82 Y
87033
LITTLE RIVER
-37.9911
144.4944
Mar 1965
Dec 1999
25.5
61 N
87036
MACEDON FORESTRY
-37.4172
144.5556
Jan 1929
Dec 1974
37.3
70 N
87061
SUNBURY (SALESIAN COLLEGE )
-37.5706
144.7403
Jan 1929
Dec 1967
35.1
88 N
87065
WERRIBEE RESEARCH FARM
82 N
87075
BULLENGAROOK EAST (SESKINORE)
87097
PARWAN
-37.9
144.6833
Jul 1968
Jun 1980
11.6
-37.4964
144.5028
Oct 1966
Dec 2001
34.3
83 N
-37.7
144.3167
Apr 1954
Nov 1973
19.4
91 N
87104
WERRIBEE CATTLE YARDS
-37.9767
144.6317
Apr 1965
Mar 1991
25.6
87 N
87105
WERRIBEE SEWERAGE FARM
-37.95
144.6167
Mar 1965
Dec 1980
15.7
86 N
87107
POINT COOK RAAF ACADEMY
-37.9333
144.75
Apr 1965
May 1976
11.2
94 N
87133
GEELONG NORTH
-38.1164
144.3667
Jan 1972
Oct 1993
19.3
78 N
88023
LAKE EILDON GOULBURN
-37.2325
145.9108
Oct 1957
Oct 2001
39.6
81 N
88029
HEATHCOTE
-36.9578
144.6933
Apr 1968
Oct 2001
24
65 N
88037
LAURISTON RESERVOIR
-37.255
144.3811
Apr 1958
Dec 2001
37.3
75 N
88049
PUCKAPUNYAL
-37
145
Apr 1968
Dec 1988
19.2
85 N
89002
BALLARAT AERODROME
-37.5128
143.7914
Aug 1954
Dec 2001
46.8
90 Y
89016
LAKE BOLAC (POST OFFICE)
-37.7125
142.8381
Apr 1968
Dec 2001
29.2
74 N
89019
MIRRANATWA (BOWACKA)
-37.4061
142.3828
May 1974
Nov 2001
22.9
76 N
89025
SKIPTON
-37.6842
143.3597
Jan 1974
Dec 1991
17.5
87 N
89082
BEAUFORT (SHEEPWASH)
-37.5053
143.2778
May 1974
Dec 2001
27.5
85 N
89085
ARARAT PRISON
-37.2789
142.9797
Dec 1983
Nov 2001
15.1
72 N
-37.7
144.85
Aug 1928
Aug 1973
32
64 N
-37.8333
144.9
Jul 1931
Jun 1976
43.5
93 N
587008
KEILOR
587012
SPOTSWOOD PUMPING STATION
55