Changes in the marine environment of Port Kembla Harbour, NSW, Australia, 1975-1995: a review

Data on water quality, sediment quality and aquatic organisms in Port Kembla Harbour from the 1970s to the 1990s are reviewed. In the 1970s, the marine environment of Port Kembla Harbour was in poor condition as a result of pollution from heavy industries. Elevated concentrations of pollutants were found in water, sediment and fish in the harbour; aquatic biodiversity was limited and many fish kills were reported. With the implementation of pollution reduction programs (required by legislation changes) by the industries surrounding the harbour since the 1970s, pollution in the harbour has been reduced dramatically, and the quality of the marine environment of the harbour has noticeably improved. Large reductions in the concentrations of certain toxic wastes and heavy metals in water have occurred. Marine life has returned to the whole harbour (parts were described in 1977 as abiotic). Contaminants in fish have decreased. Despite this achievement, however, there is still considerable room for improvement in the quality of the marine environment of the harbour.     

Stepwise use of GFLOW and MODFLOW to determine relative importance of shallow and deep receptors

A stepwise modeling approach is implemented in which a regional one-layer analytic element model is used to simulate the flow system and to furnish boundary conditions for an extracted local three-dimensional model. In this case study the stepwise approach is used to evaluate the fate of recharge in the Menomonee Valley adjacent to Lake Michigan. Two major receptors exist for recharge that flows through contaminated valley fill: the surface water estuary and a tunnel system constructed approximately 75 to 110 m below land surface to store storm runoff. The primary objective of the modeling is to delineate the contributing areas of recharge to each receptor. Of interest is the ability of the one-layer regional model to furnish flux boundary conditions to the local three-dimensional model despite the presence of vertical flow conditions at the boundaries of the local model. Sensitivity analysis suggests that the local model was insensitive to the vertical distribution of the flux. Each step of the modeling approach demonstrates that both receptors play an important role in capturing valley recharge. The pattern of capture of the one-layer model differed in shape from that delineated by the multi-layer local model in the presence of a flow system with pronounced vertical anisotropy and with sinks drawing water from different elevations.     

Scale effects of hydrostratigraphy and recharge zonation on base flow

Uncertainty regarding spatial variations of model parameters often results in the simplifying assumption that parameters are spatially uniform. However, spatial variability may be important in resource assessment and model calibration. In this paper, a methodology is presented for estimating a critical basin size, above which base flows appear to be relatively less sensitive to the spatial distribution of recharge and hydraulic conductivity, and below which base flows are relatively more sensitive to this spatial variability. Application of the method is illustrated for a watershed that exhibits distinct infiltration patterns and hydrostratigraphic layering. A ground water flow model (MODFLOW) and a parameter estimation code (UCODE) were used to evaluate the influence of recharge zonation and hydrostratigraphic layering on base flow distribution. Optimization after removing spatial recharge variability from the calibrated model altered base flow simulations up to 53% in watersheds smaller than 40 km(2). Merging six hydrostratigraphic units into one unit with average properties increased base flow residuals up to 83% in basins smaller than 50 km(2). Base flow residuals changed <5% in watersheds larger than 40 and 50 km(2) when recharge and hydrostratigraphy were simplified, respectively; thus, the critical basin size for the example area is approximately 40 to 50 km(2). Once identified for an area, a critical basin size could be used to guide the scale of future investigations. By ensuring that parameter discretization needed to capture base flow distribution is commensurate with the scope of the investigation, uncertainty caused by overextending uniform parameterization or by estimating extra parameter values is reduced.     

Investigation of recreational boats as a source of copper at anchorage sites using time-integrated diffusive gradients in thin film and sediment measurements

Antifouling paints on small to medium recreational vessels were first recognized as an important source of pollution in the 1970s. One of the principle biocides in these paints is Cu. Results from a field program in the Gold Coast Broadwater, Queensland, demonstrate a clear correlation between recreational boat numbers at anchorage sites and water column Cu concentrations, for both time-integrated 24 h DGT measurements (n=14, r=0.815, p<0.001) and measurements on composite 0.45 mum-filtered grab samples (n=14, r=0.698, p<0.01) collected every 4 h over 24 h. At boat numbers above 30 the 0.45 microm-filterable Cu concentrations were mostly above the relevant guideline value (1.3 microgl(-1)) but the DGT-reactive Cu concentrations were well below this value at all boat numbers studied. For three-day DGT deployments in Moreton Bay, Queensland, correlations between Cu levels and vessel numbers were not observed, possibly because of uncertain estimates of boat numbers. However, using a multi-factorial ANOVA, DGT-reactive Cu concentrations showed a significant effect for 'sites with vessels vs. sites without vessels' (alpha=0.10, p=0.077) and for 'sampling period' (non-holiday weekdays, weekends, holiday weekdays, holiday weekends) at alpha=0.10, p=0.02. Cu levels in sediments at the same sites were strongly influenced by fraction of clay (<63 microm) material but also by whether the samples were collected at an anchorage or control site. Results from this study further support the view that Cu emissions from antifouling paints may become an important source in waters with high boat numbers and should be taken into account when designing management instruments for coastal waterways.     

Simulating ground water-lake interactions: approaches and insights

Approaches for modeling lake-ground water interactions have evolved significantly from early simulations that used fixed lake stages specified as constant head to sophisticated LAK packages for MODFLOW. Although model input can be complex, the LAK package capabilities and output are superior to methods that rely on a fixed lake stage and compare well to other simple methods where lake stage can be calculated. Regardless of the approach, guidelines presented here for model grid size, location of three-dimensional flow, and extent of vertical capture can facilitate the construction of appropriately detailed models that simulate important lake-ground water interactions without adding unnecessary complexity. In addition to MODFLOW approaches, lake simulation has been formulated in terms of analytic elements. The analytic element lake package had acceptable agreement with a published LAKI problem, even though there were differences in the total lake conductance and number of layers used in the two models. The grid size used in the original LAKI problem, however, violated a grid size guideline presented in this paper. Grid sensitivity analyses demonstrated that an appreciable discrepancy in the distribution of stream and lake flux was related to the large grid size used in the original LAKI problem. This artifact is expected regardless of MODFLOW LAK package used. When the grid size was reduced, a finite-difference formulation approached the analytic element results. These insights and guidelines can help ensure that the proper lake simulation tool is being selected and applied.