Estuarine modification of nutrient and sediment exports to the Great Barrier Reef Marine Park from the Daintree and Annan River catchments

Nutrient and suspended sediment concentrations were measured in the dry season and during the rising and falling stages of flood events in the Annan and Daintree rivers to estimate catchment exports. These flood events were also sampled along the salinity gradient in the estuary and nearshore shelf to quantify the modification of terrestrial sediment and nutrient loads as they pass through estuaries into the Great Barrier Reef lagoon. In the Daintree River TSS concentrations were found to increase between the catchment and the estuary plume. The source of TSS may have been scour of the estuarine channel or from land use in the catchment of the lower estuary. In the dry season nitrogen enters the Annan and Daintree estuaries predominantly in the form of PON and DON in roughly equal proportions. Nitrogen exports to the GBR are mostly in the form of DON. In the wet season the majority of nitrogen enters the estuaries as DON and leaves as PON. Nitrogen removal in the estuaries and plumes appears to be biologically mediated once suspended sediment concentrations decrease to a point where phytoplankton growth is not light limited. In the dry season phosphorus enters and leaves the estuaries primarily in organic form. PIP is the dominant form of phosphorus in river water, but leaves the estuary more evenly distributed between all forms. These estuarine processes result in less nitrogen and phosphorus being delivered to the GBR lagoon than is exported from the catchment. The differences between these estuaries highlights the need for further work to explore modifications in estuaries that drain into the Great Barrier Reef lagoon.     

Towards ecologically relevant targets for river pollutant loads to the Great Barrier Reef

Degradation of coastal ecosystems in the Great Barrier Reef (GBR), Australia, has been linked with a decline in water quality from land-based runoff. This paper examines the reduction in current end-of-catchment loads required for total suspended solids (TSS) and dissolved inorganic nitrogen (DIN) to achieve GBR water quality guidelines. Based on first-order estimates of sustainable pollutant loads, current TSS and DIN loads would need to be reduced by approximately 7000ktons/y (41%) and 6000tons/y (38%), respectively. Next, these estimated reductions for TSS and DIN are compared with Reef Plan targets for anthropogenic sediment (-20% by 2020) and nitrogen (-50% by 2013) loads. If successful, these targets will accomplish approximately 40% of TSS and 92% of DIN load reductions required to achieve sustainable loads to the GBR lagoon. These first-order estimates elucidate the need to establish ecologically relevant targets for river pollutant loads to the GBR for management and policy.     

Dynamics of herbicide transport and partitioning under event flow conditions in the lower Burdekin region, Australia

This study examined the temporal variability in herbicide delivery to the Great Barrier Reef (GBR) lagoon (Australia) from one of the GBR catchment's major sugarcane growing regions. Annual loads of measured herbicides were consistently in the order of 200+kg. Atrazine, it's degradate desethylatrazine, and diuron contributed approximately 90% of annual herbicide load, with early 'first-flush' events accounting for the majority of herbicide loads leaving the catchment. Assessment of herbicide water-sediment partitioning in flood runoff highlighted the majority of herbicides were transported in predominantly dissolved form, although a considerable fraction of diuron was transported in particulate-bound form (ca. 33%). Diuron was also the herbicide demonstrating the highest concentrations and frequency of detection in sediments collected from catchment waterways and adjacent estuarine-marine environments, an outcome aligning with previous research. Herbicide physico-chemical properties appear to play a crucial role in partitioning between water column and sediment habitat types in GBR receiving ecosystems.     

Exposure of inner-shelf reefs to nutrient enriched runoff entering the Great Barrier Reef Lagoon: post-European changes and the design of water quality targets

We used historical flood plume extent data (modelled) to quantify the typical spatial extent of the summer runoff-seawater mixing zone of the Great Barrier Reef (GBR) lagoon. Spatially explicit analysis of the variability of in situ chlorophyll a concentrations (observed) across the runoff-seawater mixing zone, then allowed us to explore regional differences in the nutrient enrichment impact of runoff events from the various river systems that drain the GBR catchment. We demonstrate the existence of a discernable north-south gradient along the length of the GBR, such that for equivalent runoff:seawater dilutions ratios, lower levels of nutrient enrichment (as indicated by chlorophyll alpha observations) result from the river systems that drain the relatively undisturbed northern areas of the GBR catchment, compared to more human-impacted central and south areas. We identify a strong correlation between this north-south enrichment gradient and the flood concentration of dissolved inorganic nitrogen (DIN) entrained by the various river systems. By substituting the nutrient enrichment characteristics of the human-impacted river discharges with those of the undisturbed northern rivers, we provide a means to compare the short-term enriching 'footprint' for existing runoff intrusions with those that are likely to have occurred under pre-European catchment conditions. We demonstrate that under pre-European conditions, the nutrient enriching impact from river runoff was likely to have been largely constrained within 1-2 km of the coast, whereas existing conditions support the impact of reefs some 20-30 km off the coast. By using the developed spatial relations, we show that for the heavily human-impacted river systems, reductions in the end-of-river concentrations of DIN in the order 50-80% are needed in order to restore parity with pre-European conditions. We discuss these results in regard to developing end-of-catchment water quality targets for the region.     

A multi-criteria approach to Great Barrier Reef catchment (Queensland, Australia) diffuse-source pollution problem

This paper presents a multi-criteria based tool for assessing the relative impact of diffuse-source pollution to the Great Barrier Reef (GBR) from the river basins draining into the GBR lagoon. The assessment integrates biophysical and ecological data of water quality and pollutant concentrations with socio-economic information pertaining to non-point source pollution and (potential) pollutant impact. The tool generates scores for each river basin against four criteria, thus profiling the basins and enabling prioritization of management alternatives between and within basins. The results support policy development for pollution control through community participation, scientific data integration and expert knowledge contributed by people from across the catchment. The results specifically provided support for the Reef Water Quality Protection Plan, released in October 2003. The aim of the plan is to provide a framework for reducing discharge of sediment, nutrient and other diffuse-source loads and (potential) impact of that discharge and for prioritising management actions both between and within river basins.     

The Pollution Load from the River Odra in Comparison to That in Other Polish Rivers in 1988–1997

The discharge of organic matter, nutrients, and heavy metals from the river Odra into the Baltic Sea is compared with the combined discharge of the same substances from other Polish rivers into the Baltic during 1988–1997. The impact of the summer flood of 1997 is also discussed. The discharges of organic matter, nitrogen and phosphorus compounds were closely correlated with the flow rate. The annual loads of these nutrients measured in 1988 and 1996 under similar flow conditions decreased by 15 to 40% in these rivers as a result of the reduction of point sources. A steady decrease in the heavy metal load, independent of weather conditions, has been recorded during the past decade. Zinc decreased by up to 85% in the Odra and 99% in the Wisła. Although the water outflow from the Odra was half that of the Wisła, it transported a proportionally higher load of chemical compounds. Despite the differences in flow rates, the contribution of both rivers to the total riverine phosphorus discharge was approximately the same. The total water outflow during the flood (from 1 July to 28 August 1997) increased by 65% near the mouth of the Odra in comparison to 1996. The phosphate discharge increased by 34%, that of nitrite nitrogen by 88%     

Improvement of surface run‐off in the hydrological model ParFlow by a scale‐consistent river parameterization

We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.     

Hydrodynamic characterization of the Arcachon Bay,using model-derived descriptors

A numerical model (MARS-2D) was developed, with the aim of describing the hydrodynamics that prevail in Arcachon Bay. Direct model results as well as derived mixing and transport time-scales (tidal prism, local and integrated flushing times, age of water masses), were used to understand the behaviour of water masses and exchanges between the Bay and its frontiers. Particular attention was paid to the processes that drive the hydrodynamics (tides, wind and rivers), in order to understand their respective influence.The Arcachon Bay hydrodynamic system appears primarily to be highly influenced by tides; secondarily, by winds. About two third of the lagoon total volume is flushed in and out at each tidal cycle, which represent a mean tidal prism of 384 millions of cubic meters. The percentage of seawater flushed out during the ebb, that returns into the lagoon during the following flood flow is very high (return flow factor=0.95). This pattern leads to calculated integrated flushing times (IFT) ranging from 12.8 to 15.9 days, respectively, for the winter 2001 and summer 2005 simulations (two contrasting climatological situations: in summer, light northwesterly winds and low discharges in the rivers and, in winter, stronger southwesterly winds and higher river flows). Moreover, it has been found that northerly and westerly winds tend to reduce the flushing time, whilst southerly and easterly winds tend to hinder the renewal of the water in the Bay. The behaviour of the waters originating from the two main rivers of the lagoon, was studied also by means of the mean age assessment, under varying conditions of river flow and wind regime.     

River loads of suspended solids, nitrogen, phosphorus and herbicides delivered to the Great Barrier Reef lagoon

Degradation of coastal ecosystems in the Great Barrier Reef (GBR) lagoon, Australia, has been linked with increased land-based runoff of suspended solids, nutrients and pesticides since European settlement. This study estimated the increase in river loads for all 35 GBR basins, using the best available estimates of pre-European and current loads derived from catchment modelling and monitoring. The mean-annual load to the GBR lagoon for (i) total suspended solids has increased by 5.5 times to 17,000ktonnes/year, (ii) total nitrogen by 5.7 times to 80,000tonnes/year, (iii) total phosphorus by 8.9 times to 16,000tonnes/year, and (iv) PSII herbicides is 30,000kg/year. The increases in river loads differ across the 10 pollutants and 35 basins examined, reflecting differences in surface runoff, urbanisation, deforestation, agricultural practices, mining and retention by reservoirs. These estimates will facilitate target setting for water quality and desired ecosystem states, and enable prioritisation of critical sources for management.     

Comparison of impacts of dams on the annual maximum flow characteristics in three regulated hydrologic regimes in Québec (Canada)

Despite the presence of numerous dams in Québec, no study has yet been devoted to their impacts on flood levels. To compensate for this deficiency, we have compared the impacts of dams on the five characteristics (magnitude and its interannual variability, timing and its interannual variability, and asymmetry) of the maximum annual flows between natural rivers and regulated rivers by means of several statistical approaches (analysis of variance, chi‐square test, nonparametric tests, etc.). In the course of this study, we analysed 88 stations on pristine rivers and 60 stations on regulated rivers. The latter group was subdivided into three regulated hydrologic regimes, i.e. inversed flow regimes (25 stations), homogenization flow regimes (15 stations) and natural‐type flow regimes (20 stations). The following observations emerge from this study. (1) In inversed and homogenization flow regimes, generally associated with reservoirs, all the flow characteristics are modified. These modifications notably entrain a decrease in magnitude, a significant reduction in the frequency of the maximum annual spring flows when the snow is melting and an increase in skewness of the distribution and interannual variability of the magnitude and dates of occurrence of the annual maximum flows. We also observed the disappearance of most flows with a recurrence of over 10 years. All these changes particularly affect watersheds larger than 10 000 km². (2) In natural‐type flow regimes, often associated with run‐of‐river dams, very few changes were observed compared with pristine rivers. These changes primarily affected watersheds smaller than 1000 km². Copyright © 2006 John Wiley & Sons, Ltd.     

Estimates of sediment and nutrient loads in 10 major catchments draining to the Great Barrier Reef during 2006-2009

The Great Barrier Reef (GBR) catchment area has been monitored simultaneously for sediment and nutrient exports from 10 priority catchments discharging into the GBR lagoon between 2006 and 2009. This allows GBR catchment-wide exports to be estimated and spatially compared within a discrete time-frame. Elevated levels of sediment and nutrient exports were recorded in all monitored catchments as compared to pre-European estimates, but vary around previous estimates of mean annual loads. During the period of monitoring, the Burdekin and Fitzroy catchments contributed the highest sediment and nutrient exports, however when loads were normalised for area, these catchments produced the lowest unit yields. In contrast, the highest yields were produced in the wetter and proportionately more intensively cultivated Johnstone, O'Connell, and Pioneer catchments particularly for dissolved nitrogens. This assessment offers the necessary scientific foundation for future monitoring, assessment, and management of sediment and nutrient loads entering the GBR.     

Analysis of channel transmission losses in a dryland river reach in north‐eastern Brazil using streamflow series,groundwater level series and multi‐temporal satellite data

Scarcity of hydrological data, especially streamflow discharge and groundwater level series, restricts the understanding of channel transmission losses (TL) in drylands. Furthermore, the lack of information on spatial river dynamics encompasses high uncertainty on TL analysis in large rivers. The objective of this study was to combine the information from streamflow and groundwater level series with multi‐temporal satellite data to derive a hydrological concept of TL for a reach of the Middle Jaguaribe River (MJR) in semi‐arid north‐eastern Brazil. Based on this analysis, we proposed strategies for its modelling and simulation. TL take place in an alluvium, where river and groundwater can be considered to be hydraulically connected. Most losses certainly infiltrated only through streambed and levees and not through the flood plains, as could be shown by satellite image analysis. TL events whose input river flows were smaller than a threshold did not reach the outlet of the MJR. TL events whose input flows were higher than this threshold reached the outlet losing on average 30% of their input. During the dry seasons (DS) and at the beginning of rainy seasons (DS/BRS), no river flow is expected for pre‐events, and events have vertical infiltration into the alluvium. At the middle and the end of the rainy seasons (MRS/ERS), river flow sustained by base flow occurs before/after events, and lateral infiltration into the alluvium plays a major role. Thus, the MJR shifts from being a losing river at DS/BRS to become a losing/gaining (mostly losing) river at MRS/ERS. A model of this system has to include the coupling of river and groundwater flow processes linked by a leakage approach. Copyright © 2012 John Wiley & Sons, Ltd.     

Flow transmission along an arid zone anastomosing river,cooper creek,australia

Cooper Creek is characterized by an exceptionally wide floodplain and a complex system of anastomosing channels which together can accommodate a large range of highly variable flows. In consequence flood frequency curves are unusually steep. With minimal tributary contribution over the two study lengths — Currareva-Nappa Merrie (> 400 km) and Nappa Merrie-Innamincka (32 km) — downstream output is largely controlled by upstream input, so that similar flow regimes and hydrograph forms can be maintained despite transmission losses and the variable activation of channels. Input-output relationships defined in terms of total flow volume and peak discharge suggest that, above a threshold flow of about 25% duration, transmission losses exceed 75% on average over the Currareva-Nappa Merrie length. The large difference in absorption between primary channels on the one hand and subsidiary channels and the floodplain surface on the other means that transmission losses vary non-linearly with stage. They are relatively low when flow is confined to the primary channels, but increase at higher stages as lesser channels and the floodplain are activated. Early levels do not recur until flows of < 0.5% duration are attained when output/input ratios increase rapidly with discharge. The influences of seasonality and background flow on output discharge are small relative to that of input discharge. Despite the long river distance, peak flows at Currareva and Nappa Merrie are well correlated. However, over the shorter Nappa Merrie-Innamincka length with its more confined system of channels, correlations are better and transmission losses less. They vary from 60% just above the threshold discharge to <10% at flows below 20% duration. Over this length outflow hydrographs can be effectively estimated by the three parameter Muskingum procedure, at least on an event by event basis, but further application must await additional records to understand fully how the parameters vary in an environment where streamflow is transmitted through a highly variable system of channels.     

The influence of a season of extreme wet weather events on exposure of the World Heritage Area Great Barrier Reef to pesticides

The 2010-2011 wet season was one of extreme weather for the State of Queensland, Australia. Major rivers adjacent to the Great Barrier Reef (GBR) were discharging at rates 1.5 to >3 times higher than their long term median. Exposure to photosystem II herbicides has been routinely monitored over a period of up to 5 years at 12 inshore GBR sites. The influence of this wet season on exposure to photosystem II herbicides was examined in the context of this long-term monitoring record and during flood plume events in specific regions. Median exposures expressed as diuron equivalent concentration were an average factor of 2.3 times higher but mostly not significantly different (p<0.05) to the median for the long-term monitoring record. The herbicides metolachlor and tebuthiuron were frequently detected in flood plume waters at concentrations that reached or exceeded relevant water quality guidelines (by up to 4.5 times).     

Hydrologic impact of regional climate change for the snowfed and glacierfed river basins in the Republic of Tajikistan: hydrological response of flow to climate change

The warming of the Earth's atmosphere system is likely to change temperature and precipitation, which may affect the climate, hydrology and water resources at the river basins over the world. The importance of temperature change becomes even greater in snow or glacier dominated basins where it controls the snowmelt processes during the late‐winter, spring and summer months. In this study hydrologic responses of streamflow in the Pyanj and Vaksh River basins to climate change are analysed with a watershed hydrology model, based on the downscaled atmospheric data as input, in order to assess the regional climate change impact for the snowfed and glacierfed river basins in the Republic of Tajikistan. As a result of this analysis, it was found that the annual mean river discharge is increasing in the future at snow and glacier dominated areas due to the air temperature increase and the consequent increase in snow/ice melt rates until about 2060. Then the annual mean flow discharge starts to decrease from about 2080 onward because the small glaciers start to disappear in the glacier areas. It was also found that there is a gradual change in the hydrologic flow regime throughout a year, with the high flows occuring earlier in the hydrologic year, due to the warmer climate in the future. Furthermore, significant increases in annual maximum daily flows, including the 100‐year return period flows, at the Pyanj and Vaksh River basins toward the end of the 21st century can be inferred from flood frequency analysis results. Copyright © 2012 John Wiley & Sons, Ltd.     

An assessment of residence times of land-sourced contaminants in the Great Barrier Reef lagoon and the implications for management and reef recovery

We argue that the residence times of key pollutants exported to the Great Barrier Reef (GBR) are greater in the GBR lagoon than those of the water itself, in contradiction to some previous assumptions. Adverse effects of the pollutant discharge will be greater and longer lasting than previously considered, in turn requiring stronger or more urgent action to remediate land practices. Residence times of fine sediments, nitrogen and phosphorus, pesticides and trace metals are suggested to be from years to decades in the GBR lagoon and highly likely to be greater than the residence time of water, estimated at around 15-365days. The recovery of corals and seagrass in the central region of the GBR following current land-use remediation in the catchment depends on the residence time of these contaminants. Ecohydrological modeling suggests that this recovery may take decades even with adequate levels of improved land management practices.     

Circulation in the Great Barrier Reef Lagoon using numerical tracers and in situ data

Numerical hydrodynamic models of the northeastern Queensland shelf, forced by regional winds and modelled boundary currents in the northern Coral Sea, are used to provide improved estimates of general flow trajectories and water residence times within the Great Barrier Reef (GBR) shelf system. Model performance was checked against a limited set of current metre records obtained at Lark Reef (16°S) and the Ribbon Reefs (15.5°S). Estimates of water parcel trajectories are derived from a series of numerical tracer experiments, with daily releases of neutrally buoyant, un-reactive particles at 320 sites along the coast between Cape York (10.7°S) and Hervey Bay (25°S). Flow trajectories and residence times for tracer particles introduced to the GBR lagoon in the southern—ca. 22°S, central—19°S, and northern reef—14°S are emphasised. For purposes of the analysis, the year was divided into two seasons based on mean alongshore current direction. Most coastal sourced tracers entering the central GBR lagoon between 16° and 20°S during the northward-current season (January–August) primarily encounter the outer-shelf reef matrix after exiting the lagoon at its northern “head” (nominally 16°S), after 50–150 days. Up to 70% of tracer particles entering in the southward-current season (August–December) eventually crossed the lagoon to the outer-shelf reef matrix, with median crossing times between 20 and 330 days. During favourable wind conditions, tracers introduced at the coast may move rapidly across the lagoon into the reef matrix. The tracer experiments indicate that most coastal-sourced tracers entering the GBR lagoon remain near the coast for extended periods of time, moving north and south in a coastal band. Residence times for conservative tracer particles (and implied residence times for water-borne materials) within the GBR shelf system ranged from ca. 1 month to 1 year—time frames that are very long relative to development times of planktonic larvae and cycling times for nutrient materials in the water column, implying they are transformed long before reaching the outer reef matrix.     

Analysis of trends in annual streamflow to the Arctic Ocean

There is increasing interest in the magnitude of the flow of freshwater to the Arctic Ocean due to its impacts on the biogeophysical and socio‐economic systems in the north and its influence on global climate. This study examines freshwater flow based on a dataset of 72 rivers that either directly or indirectly contribute flow to the Arctic Ocean or reflect the hydrologic regime of areas contributing flow to the Arctic Ocean. Annual streamflow for the 72 rivers is categorized as to the nature and location of the contribution to the Arctic Ocean, and composite series of annual flows are determined for each category for the period 1975 to 2015. A trend analysis is then conducted for the annual discharge series assembled for each category. The results reveal a general increase in freshwater flow to the Arctic Ocean with this increase being more prominent from the Eurasian rivers than from the North American rivers. A comparison with trends obtained from an earlier study ending in 2000 indicates similar trend response from the Eurasian rivers, but dramatic differences from some of the North American rivers. A total annual discharge increase of 8.7 km³/y/y is found, with an annual discharge increase of 5.8 km³/y/y observed for the rivers directly flowing to the Arctic Ocean. The influence of annual or seasonal climate oscillation indices on annual discharge series is also assessed. Several river categories are found to have significant correlations with the Arctic Oscillation, the North Atlantic Oscillation, or the Pacific Decadal Oscillation. However, no significant association with climate indices is found for the river categories leading to the largest freshwater contribution to the Arctic Ocean.     

Water quality in the inshore Great Barrier Reef lagoon: Implications for long-term monitoring and management

Coastal and inshore areas of the Great Barrier Reef lagoon receive substantial amounts of material from adjacent developed catchments, which can affect the ecological integrity of coral reefs and other inshore ecosystems. A 5-year water quality monitoring dataset provides a 'base range' of water quality conditions for the inshore GBR lagoon and illustrates the considerable temporal and spatial variability in this system. Typical at many sites were high turbidity levels and elevated chlorophyll a and phosphorus concentrations, especially close to river mouths. Water quality variability was mainly driven by seasonal processes such as river floods and sporadic wind-driven resuspension as well as by regional differences such as land use. Extreme events, such as floods, caused large and sustained increases in water quality variables. Given the highly variable climate in the GBR region, long-term monitoring of marine water quality will be essential to detect future changes due to improved catchment management.     

ENSO and the natural variability in the flow of tropical rivers

This paper examines the relationship between the annual discharges of the Amazon, Congo, Paran á, and Nile rivers and the sea surface temperature (SST) anomalies of the eastern and central equatorial Pacific Ocean, an index of El Niño-Southern Oscillation (ENSO). Since river systems are comprehensive integrators of rainfall over large areas, accurate characterization of the flow regimes in major rivers will increase our understanding of large-scale global atmospheric dynamics. Results of this study reveal that the annual discharges of two large equatorial tropical rivers, the Amazon and the Congo, are weakly and negatively correlated with the equatorial Pacific SST anomalies with 10% of the variance in annual discharge explained by ENSO. Two smaller subtropical rivers, the Nile and the Paraná, show a correlation that is stronger by about a factor of 2. The Nile discharge is negatively correlated with the SST anomaly, whereas the Paraná river discharge shows a positive relation. The tendency for reduced rainfall/discharge over large tropical convection zones in the ENSO warm phase is attributed to global scale subsidence associated with major upwelling in the eastern Pacific Ocean.