Male Dominance in the New Zealand Longfin Eel Population of a New Zealand River: Probable Causes and Implications for Management

The endemic New Zealand longfin eel Anguilla dieffenbachi (hereafter, longfin eel), is overfished, and in southern South Island, New Zealand, rivers have recently become predominated by males. This study examined length and age at sexual differentiation in male eels in the Aparima River catchment (area, 1,375 km(2); mean flow, 20 m(3.)s(-1)) and the sex ratio and distribution of eels throughout the catchment. Longfin eels differentiated into males mostly at lengths from 300 to 460 mm and ages from 10 to 25+ years. Females were rare: Of 738 eels examined for sexual differentiation, 466 were males and 5 were females, and a few others, not examined, were large enough to be female. These counts suggest a male : female ratio among differentiated longfin eels of 68:1. Of 31 differentiated shortfin eels A. australis, less common in the Aparima River, 26 were females. Male longfin eels were distributed throughout the main stern and tributaries; undifferentiated eels were more prevalent in lower and middle reaches and in the main stem than in upper reaches and tributaries. In other studies, male longfin eels predominated commercial catches in the Aparima and four other southernmost rivers, by 2.4:1 to 13.6:1 males to females. The Aparima River had the most skewed sex ratio. Longfin eel catches from the Aparima River will become more male predominated because few sublegal-size females were present. The length-frequency distributions of eels in the present samples and in the commercial catches were truncated just above minimum legal size (about 460 mm), showing that few females escape the fishery. Historically, females predominated these rivers. The recent change in sex ratio is attributable partly to selective harvest of females, and partly to changes in the structure of the population from fishing, such that differentiation into males has been favored. Longevity, delayed sexual maturity, semel-parity, and endemism with restricted range make the longfin eel particularly vulnerable to overfishing.

Hydrogeological investigations in the upper Ouémé catchment in Benin, West Africa

The scope of this PhD thesis was the hydrogeological conceptualisation of the Upper Ouémé river catchment in Benin. The study area exceeds 14,500 km**2 and is underlain by a crystalline basement. At this setting the typical sequence of aquifers - a regolith aquifer at the top and a fractured bedrock aquifer at the bottom - is encountered, which is found in basement areas all over Africa and elsewhere in the world. The chosen regional approach revealed important information about the hydrochemistry and hydrogeology of this catchment. Based on the regional conceptual model a numerical groundwater flow model was designed. The numerical model was used to estimate the impact of climate change on the regional groundwater resources. This study was realised within the framework of the German interdisciplinary research project IMPETUS (English translation: "Integrated approach to the efficient management of scarce water resources in West Africa"), which is jointly managed by the German universities of Bonn and Cologne. Since the year 2000 the Upper Ouémé catchment was the principal target for investigations into the relevant processes of the regional water cycle. A first study from 2000 to 2003 (Fass, 2004, http://nbn-resolving.de/urn:nbn:de:hbz:5n-03849) focused on the hydrogeology of a small local catchment (~30 km**2). In the course of this thesis five field campaigns were underdone from the year 2004 to 2006. In the beginning of 2004 a groundwater monitoring net was installed based on 12 automatic data loggers. Manual piezometric measurements and the sampling of groundwater and surface water were realised for each campaign throughout the whole study area. Water samples were analysed for major ions, for a choice of heavy metals and for their composition by deuterium, oxygen-18 and tritium. The numerical model was performed with FEFLOW. The hydraulic and hydrochemical characteristics were described for the regolith aquifer and the bedrock aquifer. The regolith aquifer plays the role of the groundwater stock with low conductivity while the fractures of the bedrock may conduct water relatively fast towards extraction points. Flow in fractures of the bedrock depends on the connectivity of the fracture network which might be of local to subregional importance. Stable isotopes in combination with hydrochemistry proved that recharge occurs on catchment scale and exclusively by precipitation. Influx of groundwater from distant areas along dominant structures like the Kandi fault or from the Atacora mountain chain is excluded. The analysis of tritium in groundwater from different depths revealed the interesting fact of the strongly rising groundwater ages. Bedrock groundwater may possibly be much older than 50 years. Equilibrium phases of the silicate weathering products kaolinite and montmorillonite showed that the deeper part of the regolith aquifer and the bedrock aquifer feature either stagnant or less mobile groundwater while the shallow aquifer level is influenced by seasonal groundwater table fluctuations. The hydrochemical data characterised this zone by the progressive change of the hydrochemical facies of recently infiltrated rainwater on its flow path into deeper parts of the aquifers. Surprisingly it was found out that seasonal influences on groundwater hydrochemistry are minor, mainly because they affect only the groundwater levels close to the surface. The transfer of the hydrogeological features of the Upper Ouémé catchment into a regional numerical model demanded a strong simplification. Groundwater tables are a reprint of the general surface morphology. Pumping or other types of groundwater extraction would have only very local impact on the available groundwater resources. It was possible to integrate IMPETUS scenario data into the groundwater model. As a result it was shown that the impact of climate change on the groundwater resources until the year 2025 under the given conditions will be negligible due to the little share of precipitation needed for recharge and the low water needs for domestic use. Reason for concern is the groundwater quality on water points in the vicinity of settlements because of contamination by human activities as shown for the village of Dogué. Nitrate concentrations achieved in many places already alerting levels. Health risks from fluoride or heavy metals were excluded for the Upper Ouémé area.

Organic compounds in snow, ice, and sub-ice waters from the mouth of the Severnaya Dvina River and the Kandalaksha Gulf (White Sea)

The paper reports data on concentrations of organic compounds (organic carbon, lipids; aliphatic hydrocarbons, and polycyclic aromatic hydrocarbons) in snow, ice, and sub-ice waters from the mouth of the Severnaya Dvina River in March 2005-2007 and the Kandalaksha Gulf (Chupa Bay) in March 2004. It was established that organic compounds are accumulated in snow and the upper ice layer near Archangelsk city. Distribution of molecular markers indicates that pollutions were mainly caused by local fallouts. In the Chupa Bay organic compounds are concentrated in the lower ice layer; it is typical for Arctic snow-ice cover. High contents of organic compounds in the snow-ice cover of the White Sea are caused by pollution of air and water during the winter season.

Particulate organic matter (POM) of samples from the Lena River delta

Particulate organic matter (POM) derived from permafrost soils and transported by the Lena River represents a quantitatively important terrestrial carbon pool exported to Laptev Sea sediments (next to POM derived from coastal erosion). Its fate in a future warming Arctic, i.e., its remobilization and remineralization after permafrost thawing as well as its transport pathways to and sequestration in marine sediments, is currently under debate. We present one of the first radiocarbon (14C) data sets for surface water POM within the Lena Delta sampled in the summers of 2009 - 2010 and spring 2011 (n = 30 samples). The bulk D14C values varied from -55 to -391 per mil translating into 14C ages of 395 to 3920 years BP. We further estimated the fraction of soil-derived POM to our samples based on (1) particulate organic carbon to particulate nitrogen ratios (POC : PN) and (2) on the stable carbon isotope (d13C) composition of our samples. Assuming that this phytoplankton POM has a modern 14C concentration, we inferred the 14C concentrations of the soil-derived POM fractions. The results ranged from -322 to -884 per mil (i.e., 3060 to 17 250 14C years BP) for the POC : PN-based scenario and from -261 to -944 per mil (i.e., 2370 to 23 100 14C years BP) for the d13C-based scenario. Despite the limitations of our approach, the estimated D14C values of the soil-derived POM fractions seem to reflect the heterogeneous 14C concentrations of the Lena River catchment soils covering a range from Holocene to Pleistocene ages better than the bulk POM D14C values. We further used a dual-carbon-isotope three-end-member mixing model to distinguish between POM contributions from Holocene soils and Pleistocene Ice Complex (IC) deposits to our soil-derived POM fraction. IC contributions are comparatively low (mean of 0.14) compared to Holocene soils (mean of 0.32) and riverine phytoplankton (mean of 0.55), which could be explained with the restricted spatial distribution of IC deposits within the Lena catchment. Based on our newly calculated soil-derived POM D14C values, we propose an isotopic range for the riverine soil-derived POM end member with D14C of -495 ± 153 per mil deduced from our d13C-based binary mixing model and d13C of -26.6 ± 1 per mil deduced from our data of Lena Delta soils and literature values. These estimates can help to improve the dual-carbon-isotope simulations used to quantify contributions from riverine soil POM, Pleistocene IC POM from coastal erosion, and marine POM in Siberian shelf sediments.

Compositions of paleosols from the Salarevo Formation, Sukhona River valley (Severnaya Dvina basin)

Paleosols crop out in the Sukhona River valley as several members up to 10 m thick embedded into the Salarevo Formation sediments. Principal characteristics of the paleosols include a dense network of root channels, indications of eluvial gley alteration, redistribution and formation of secondary carbonates represented by several generations, and formation of block-prismatic soil structure with specific clayey films at structural jointing faces. The paleosols are divided into a number of genetically interrelated horizons (from top to bottom): presumably organogenic accumulation (AElg), eluvial gley horizon (Elg), illuvial horizons (B1 and B2), illuvial gley horizon (Bg), and transitional horizons (ElBg and BElg). The paleosols formed under conditions of a semiarid climate with sharp seasonal or secular and multisecular oscillations of atmospheric precipitation. Such soils point to specific ecological environments existed in the northern semiarid belt of the Earth before the greatest (in Phanerozoic) biospheric crisis at the Permian-Triassic boundary.

Characteristics of soil profiles from the Ouémé Upper Catchment (HVO), Térou, Benin

Soil degradation threatens agricultural production and food security in Sub-Saharan Africa. In the coming decades, soil degradation, in particular soil erosion, will become worse through the expansion of agriculture into savannah and forest and changes in climate. This study aims to improve the understanding of how land use and climate change affect the hydrological cycle and soil erosion rates at the catchment scale. We used the semi-distributed, time-continuous erosion model SWAT (Soil Water Assessment Tool) to quantify runoff processes and sheet and rill erosion in the Upper Ouémé River catchment (14500 km**2, Central Benin) for the period 1998-2005. We could then evaluate a range of land use and climate change scenarios with the SWAT model for the period 2001-2050 using spatial data from the land use model CLUE-S and the regional climate model REMO. Field investigations were performed to parameterise a soil map, to measure suspended sediment concentrations for model calibration and validation and to characterise erosion forms, degraded agricultural fields and soil conservation practices. Modelling results reveal current "hotspots" of soil erosion in the north-western, eastern and north-eastern parts of the Upper Ouémé catchment. As a consequence of rapid expansion of agricultural areas triggered by high population growth (partially caused by migration) and resulting increases in surface runoff and topsoil erosion, the mean sediment yield in the Upper Ouémé River outlet is expected to increase by 42 to 95% by 2025, depending on the land use scenario. In contrast, changes in climate variables led to decreases in sediment yield of 5 to 14% in 2001-2025 and 17 to 24% in 2026-2050. Combined scenarios showed the dominance of land use change leading to changes in mean sediment yield of -2 to +31% in 2001-2025. Scenario results vary considerably within the catchment. Current "hotspots" of soil erosion will aggravate, and a new "hotspot" will appear in the southern part of the catchment. Although only small parts of the Upper Ouémé catchment belong to the most degraded zones in the country, sustainable soil and plant management practices should be promoted in the entire catchment. The results of this study can support planning of soil conservation activities in Benin.

Dietary ecology of the Australian freshwater turtle (Elseya sp.: Chelonia : Chelidae) in the Burnett River, Queensland

The Burnett River snapping turtle (Elseya sp.) from the Burnett, Mary and Fitzroy river systems is an undescribed Australian freshwater turtle, of which very little ecological information is known. This paper describes the dietary ecology of the species in the Burnett River catchment. Stomach and faecal samples were collected from turtles and an index of relative importance was used to rank food items found in stomach samples. This index indicated that algae and aquatic ribbon weed (Vallisneria) were the dominant food items consumed. No difference in diet was found between males and females. Although the sample size was small, diet appeared to vary slightly seasonally, with Elseya sp. selectively feeding on the flower buds of the Chinese elm tree (Celtis chinensis) and the seeds of the blackbean tree (Castanospermum australe) when these food items were seasonally available. Faecal samples suggest that the most ingested foods ( algae and aquatic ribbon weed) were also the most digestible. Although predominantly herbivorous, Elseya sp. was seen to eat carrion once in the wild.

Sediment Budget for Cane Land on the Lower Herbert River Floodplain, North Queensland, Australia

Land use in the river catchments of tropical North Queensland appears to have increased the export of sediment and nutrients to the coast. Although evidence of harmful effect of sediment on coastal and riverine ecosystems is limited, there is a growing concern about its possible negative impacts. Sugarcane cultivation on the floodplains of the tropical North Queensland river catchments is thought to be an important source of excess sediment in the river drainage systems. Minimum-tillage, trash blanket harvesting has been shown to reduce erosion from sloping sugarcane fields, but in the strongly modified floodplain landscape other elements (e.g. drains, water furrows and headlands) could still be important sediment sources. The main objectives of this thesis are to quantify the amount of sediment coming from low-lying cane land and identify the important sediment sources in the landscape. The results of this thesis enable sugarcane farmers to take targeted measures for further reduction of the export of sediment and nutrients. Sediment budgets provide a useful approach to identify and quantify potential sediment sources. For this study a sediment budget is calculated for a part of the Ripple Creek catchment, which is a sub-catchment of the Lower Herbert River. The input of sediment from all potential sources in cane land and the storage of sediment within the catchment have been quantified and compared with the output of sediment from the catchment. Input from, and storage on headlands, main drains, minor drains and water furrows, was estimated from erosion pin and surface profile measurements. Input from forested upland, input from fields and the output at the outlet of the catchment was estimated with discharge data from gauged streams and flumes. Data for the sediment budget were collected during two ‘wet’-seasons: 1999-2000 and 2000-2001. The results of the sediment budget indicate that this tropical floodplain area is a net source of sediment. Plant cane fields, which do not have a protective trash cover, were the largest net source of sediment during the 1999-2000 season. Sediment input from water furrows was higher, but there was also considerable storage of sediment in this landscape element. Headlands tend to act as sinks. The source or sink function of drains is less clear, but seems to depend on their shape and vegetation cover. An important problem in this study is the high uncertainty in the estimates of the sediment budget components and is, for example, likely to be the cause of the imbalance in the sediment budget. High uncertainties have particularly affected the results from the 20002001 season. The main source of uncertainty is spatial variation in the erosion and deposition processes. Uncertainty has to be taken into consideration when interpreting the budget results. The observation of a floodplain as sediment source contradicts the general understanding that floodplains are areas of sediment storage within river catchments. A second objective of this thesis was therefore to provide an answer to the question: how can floodplains in the tropical North Queensland catchments can be a source of sediment? In geomorphic literature various factors have been pointed out, that could control floodplain erosion processes. However, their importance is not 'uniquely identified'. Among the most apparent factors are the stream power of the floodwater and the resistance of the floodplain surface both through its sedimentary composition and the vegetation cover. If the cultivated floodplains of the North Queensland catchments are considered in the light of these factors, there is a justified reason to expect them to be a sediment source. Cultivation has lowered the resistance of their surface; increased drainage has increased the drainage velocity and flood control structures have altered flooding patterns. For the Ripple Creek floodplain four qualitative scenarios have been developed that describe erosion and deposition under different flow conditions. Two of these scenarios were experienced during the budget study, involving runoff from local hillslopes and heavy rainfall, which caused floodplain erosion. In the longer term larger flood events, involving floodwater from the Herbert River, may lead to different erosion and deposition processes. The present study has shown that the tropical floodplain of the Herbert River catchment can be a source of sediment under particular flow conditions. It has also shown which elements in the sugarcane landscape are the most important sediment sources under these conditions. This understanding will enable sugarcane farmers to further reduce sediment export from cane land and prevent the negative impact this may have on the North Queensland coastal ecosystems.

Dewfish Demonstration Reach : Restoring native fish populations in the Condamine Catchment

The purpose of the project is to demonstrate how the restoration of riverine habitat and connectivity benefits native biodiversity and promote the importance of a healthy river system for native fish and the greater river catchment. The goal is to restore native fish populations to 60% of pre-European settlement levels and improve aquatic health within the Reach.

Influence of physical and chemical properties of solids on heavy metal adsorption

Partition of heavy metals between particulate and dissolve fraction of stormwater primarily depends on the adsorption characteristics of solids particles. Moreover, the bioavailability of heavy metals is also influenced by the adsorption behaviour of solids. However, due to the lack of fundamental knowledge in relation to the heavy metals adsorption processes of road deposited solids, the effectiveness of stormwater management strategies can be limited. The research study focused on the investigation of the physical and chemical parameters of solids on urban road surfaces and, more specifically, on heavy metal adsorption to solids. Due to the complex nature of heavy metal interaction with solids, a substantial database was generated through a series of field investigations and laboratory experiments. The study sites for the build-up pollutant sample collection were selected from four urbanised suburbs located in a major river catchment. Sixteen road sites were selected from these suburbs and represented typical industrial, commercial and residential land uses. Build-up pollutants were collected using a wet and dry vacuum collection technique which was specially designed to improve fine particle collection. Roadside soil samples were also collected from each suburb for comparison with the road surface solids. The collected build-up solids samples were separated into four particle size ranges and tested for a range of physical and chemical parameters. The solids build-up on road surfaces contained a high fraction (70%) of particles smaller than 150ìm, which are favourable for heavy metal adsorption. These solids particles predominantly consist of soil derived minerals which included quartz, albite, microcline, muscovite and chlorite. Additionally, a high percentage of amorphous content was also identified in road deposited solids. In comparing the mineralogical data of surrounding soil and road deposited solids, it was found that about 30% of the solids consisted of particles generated from traffic related activities on road surfaces. Significant difference in mineralogical composition was noted in different particle sizes of build-up solids. Fine solids particles (<150ìm) consisted of a clayey matrix and high amorphous content (in the region of 40%) while coarse particles (>150ìm) consisted of a sandy matrix at all study sites, with about 60% quartz content. Due to these differences in mineralogical components, particles larger than and smaller than 150ìm had significant differences in their specific surface area (SSA) and effective cation exchange capacity (ECEC). These parameters, in turn, exert a significant influence on heavy metal adsorption. Consequently, heavy metal content in >150ìm particles was lower than in the case of fine particles. The particle size range <75ìm had the highest heavy metal content, corresponding with its high clay forming minerals, high organic matter and low quartz content which increased the SSA, ECEC and the presence of Fe, Al and Mn oxides. The clay forming minerals, high organic matter and Fe, Al and Mn oxides create distinct groups of charge sites on solids surfaces and exhibit different adsorption mechanisms and bond strength, between heavy metal elements and charge sites. Therefore, the predominance of these factors in different particle sizes leads to different heavy metal adsorption characteristics. Heavy metals show preference for association with clay forming minerals in fine solids particles, whilst in coarse particles heavy metals preferentially associate with organic matter. Although heavy metal adsorption to amorphous material is very low, the heavy metals embedded in traffic related materials have a potential impact on stormwater quality.Adsorption of heavy metals is not confined to an individual type of charge site in solids, whereas specific heavy metal elements show preference for adsorption to several different types of charge sites in solids. This is attributed to the dearth of preferred binding sites and the inability to reach the preferred binding sites due to competition between different heavy metal species. This confirms that heavy metal adsorption is significantly influenced by the physical and chemical parameters of solids that lead to a heterogeneity of surface charge sites. The research study highlighted the importance of removal of solids particles from stormwater runoff before they enter into receiving waters to reduce the potential risk posed by the bioavailability of heavy metals. The bioavailability of heavy metals not only results from the easily mobile fraction bound to the solids particles, but can also occur as a result of the dissolution of other forms of bonds by chemical changes in stormwater or microbial activity. Due to the diversity in the composition of the different particle sizes of solids and the characteristics and amount of charge sites on the particle surfaces, investigations using bulk solids are not adequate to gain an understanding of the heavy metal adsorption processes of solids particles. Therefore, the investigation of different particle size ranges is recommended for enhancing stormwater quality management practices.

Hydrogeological Features of the Lockyer Valley

The Lockyer Valley is situated 80 km west of Brisbane and is bounded on the sou th and west by the Great Dividing Range. The valley is a major western sub - catchment of the larger Brisbane River drainage system and is drained by the Lockyer Creek. The Lockyer catchment forms approximately 20% of the total Brisbane River catchment and has an area of around 2900 km2. The Lockyer Creek is an ephemeral drainage system, and the stream and associated alluvium are the main source for irrigation water supply in the Lockyer Valley. The catchment is comprised of a number of well -defined, elongate tributaries in the south, and others in the north, which are more meandering in nature.

Assessing the suitability of existing spatial data for disaster planning and mitigation in Queensland

The Council of Australian Governments (COAG) in 2003 gave in-principle approval to a best-practice report recommending a holistic approach to managing natural disasters in Australia incorporating a move from a traditional response-centric approach to a greater focus on mitigation, recovery and resilience with community well-being at the core. Since that time, there have been a range of complementary developments that have supported the COAG recommended approach. Developments have been administrative, legislative and technological, both, in reaction to the COAG initiative and resulting from regular natural disasters. This paper reviews the characteristics of the spatial data that is becoming increasingly available at Federal, state and regional jurisdictions with respect to their being fit for the purpose for disaster planning and mitigation and strengthening community resilience. In particular, Queensland foundation spatial data, which is increasingly accessible by the public under the provisions of the Right to Information Act 2009, Information Privacy Act 2009, and recent open data reform initiatives are evaluated. The Fitzroy River catchment and floodplain is used as a case study for the review undertaken. The catchment covers an area of 142,545 km2, the largest river catchment flowing to the eastern coast of Australia. The Fitzroy River basin experienced extensive flooding during the 2010–2011 Queensland floods. The basin is an area of important economic, environmental and heritage values and contains significant infrastructure critical for the mining and agricultural sectors, the two most important economic sectors for Queensland State. Consequently, the spatial datasets for this area play a critical role in disaster management and for protecting critical infrastructure essential for economic and community well-being. The foundation spatial datasets are assessed for disaster planning and mitigation purposes using data quality indicators such as resolution, accuracy, integrity, validity and audit trail.

Mid-Holocene Aboriginal occupation of offshore islands in northern Australia? A reassessment of Wurdukanhan, Mornington Island, southern Gulf of Carpentaria, Australia

Claims for mid-Holocene Aboriginal occupation at the shell matrix site of Wurdukanhan, Mornington Island, Gulf of Carpentaria, Australia, are reassessed through an analysis of the excavated assemblage coupled with new surveys and an extensive dating program. Memmott et al. (2006, pp. 38, 39) reported basal ages of c.5000–5500 years from Wurdukanhan as 'the oldest date yet obtained for any archaeological site on the coast of the southern Gulf of Carpentaria' and used these dates to argue for 'a relatively lengthy occupation since at least the mid-Holocene'. If substantiated, with the exception of western Torres Strait, these claims make Mornington Island the only offshore island used across northern Australia in the mid-Holocene where it is conventionally thought that Aboriginal people only (re)colonised islands after sea-level maximum was achieved after the mid-Holocene. Our analysis of Wurdukanhan demonstrates high shellfish taxa diversity, high rates of natural shell predation and high densities of foraminifera throughout the deposit demonstrating a natural origin for the assemblage. Results are considered in the context of other dated shell matrix sites in the area and a geomorphological model for landscape development of the Sandalwood River catchment.

An extended regression approach to estimating loads and their uncertainties in Great Barrier Reef catchments

There are numerous load estimation methods available, some of which are captured in various online tools. However, most estimators are subject to large biases statistically, and their associated uncertainties are often not reported. This makes interpretation difficult and the estimation of trends or determination of optimal sampling regimes impossible to assess. In this paper, we first propose two indices for measuring the extent of sampling bias, and then provide steps for obtaining reliable load estimates by minimizing the biases and making use of possible predictive variables. The load estimation procedure can be summarized by the following four steps: - (i) output the flow rates at regular time intervals (e.g. 10 minutes) using a time series model that captures all the peak flows; - (ii) output the predicted flow rates as in (i) at the concentration sampling times, if the corresponding flow rates are not collected; - (iii) establish a predictive model for the concentration data, which incorporates all possible predictor variables and output the predicted concentrations at the regular time intervals as in (i), and; - (iv) obtain the sum of all the products of the predicted flow and the predicted concentration over the regular time intervals to represent an estimate of the load. The key step to this approach is in the development of an appropriate predictive model for concentration. This is achieved using a generalized regression (rating-curve) approach with additional predictors that capture unique features in the flow data, namely the concept of the first flush, the location of the event on the hydrograph (e.g. rise or fall) and cumulative discounted flow. The latter may be thought of as a measure of constituent exhaustion occurring during flood events. The model also has the capacity to accommodate autocorrelation in model errors which are the result of intensive sampling during floods. Incorporating this additional information can significantly improve the predictability of concentration, and ultimately the precision with which the pollutant load is estimated. We also provide a measure of the standard error of the load estimate which incorporates model, spatial and/or temporal errors. This method also has the capacity to incorporate measurement error incurred through the sampling of flow. We illustrate this approach using the concentrations of total suspended sediment (TSS) and nitrogen oxide (NOx) and gauged flow data from the Burdekin River, a catchment delivering to the Great Barrier Reef. The sampling biases for NOx concentrations range from 2 to 10 times indicating severe biases. As we expect, the traditional average and extrapolation methods produce much higher estimates than those when bias in sampling is taken into account.