Port of Brisbane stormwater quality : stage 2 development of a port specific water quality model

Background: The quality of stormwater runoff from ports is significant as it can be an important source of pollution to the marine environment. This is also a significant issue for the Port of Brisbane as it is located in an area of high environmental values. Therefore, it is imperative to develop an in-depth understanding of stormwater runoff quality to ensure that appropriate strategies are in place for quality improvement, where necessary. To this end, the Port of Brisbane Corporation aimed to develop a port specific stormwater model for the Fisherman Islands facility. The need has to be considered in the context of the proposed future developments of the Port area. ----------------- The Project: The research project is an outcome of the collaborative Partnership between the Port of Brisbane Corporation (POBC) and Queensland University of Technology (QUT). A key feature of this Partnership is that it seeks to undertake research to assist the Port in strengthening the environmental custodianship of the Port area through ‘cutting edge’ research and its translation into practical application. ------------------ The project was separated into two stages. The first stage developed a quantitative understanding of the generation potential of pollutant loads in the existing land uses. This knowledge was then used as input for the stormwater quality model developed in the subsequent stage. The aim is to expand this model across the yet to be developed port expansion area. This is in order to predict pollutant loads associated with stormwater flows from this area with the longer term objective of contributing to the development of ecological risk mitigation strategies for future expansion scenarios. ----------------- Study approach: Stage 1 of the overall study confirmed that Port land uses are unique in terms of the anthropogenic activities occurring on them. This uniqueness in land use results in distinctive stormwater quality characteristics different to other conventional urban land uses. Therefore, it was not scientifically valid to consider the Port as belonging to a single land use category or to consider as being similar to any typical urban land use. The approach adopted in this study was very different to conventional modelling studies where modelling parameters are developed using calibration. The field investigations undertaken in Stage 1 of the overall study helped to create fundamental knowledge on pollutant build-up and wash-off in different Port land uses. This knowledge was then used in computer modelling so that the specific characteristics of pollutant build-up and wash-off can be replicated. This meant that no calibration processes were involved due to the use of measured parameters for build-up and wash-off. ---------------- Conclusions: Stage 2 of the study was primarily undertaken using the SWMM stormwater quality model. It is a physically based model which replicates natural processes as closely as possible. The time step used and catchment variability considered was adequate to accommodate the temporal and spatial variability of input parameters and the parameters used in the modelling reflect the true nature of rainfall-runoff and pollutant processes to the best of currently available knowledge. In this study, the initial loss values adopted for the impervious surfaces are relatively high compared to values noted in research literature. However, given the scientifically valid approach used for the field investigations, it is appropriate to adopt the initial losses derived from this study for future modelling of Port land uses. The relatively high initial losses will reduce the runoff volume generated as well as the frequency of runoff events significantly. Apart from initial losses, most of the other parameters used in SWMM modelling are generic to most modelling studies. Development of parameters for MUSIC model source nodes was one of the primary objectives of this study. MUSIC, uses the mean and standard deviation of pollutant parameters based on a normal distribution. However, based on the values generated in this study, the variation of Event Mean Concentrations (EMCs) for Port land uses within the given investigation period does not fit a normal distribution. This is possibly due to the fact that only one specific location was considered, namely the Port of Brisbane unlike in the case of the MUSIC model where a range of areas with different geographic and climatic conditions were investigated. Consequently, the assumptions used in MUSIC are not totally applicable for the analysis of water quality in Port land uses. Therefore, in using the parameters included in this report for MUSIC modelling, it is important to note that it may result in under or over estimations of annual pollutant loads. It is recommended that the annual pollutant load values given in the report should be used as a guide to assess the accuracy of the modelling outcomes. A step by step guide for using the knowledge generated from this study for MUSIC modelling is given in Table 4.6. ------------------ Recommendations: The following recommendations are provided to further strengthen the cutting edge nature of the work undertaken: * It is important to further validate the approach recommended for stormwater quality modelling at the Port. Validation will require data collection in relation to rainfall, runoff and water quality from the selected Port land uses. Additionally, the recommended modelling approach could be applied to a soon-to-be-developed area to assess ‘before’ and ‘after’ scenarios. * In the modelling study, TSS was adopted as the surrogate parameter for other pollutants. This approach was based on other urban water quality research undertaken at QUT. The validity of this approach should be further assessed for Port land uses. * The adoption of TSS as a surrogate parameter for other pollutants and the confirmation that the <150 m particle size range was predominant in suspended solids for pollutant wash-off gives rise to a number of important considerations. The ability of the existing structural stormwater mitigation measures to remove the <150 m particle size range need to be assessed. The feasibility of introducing source control measures as opposed to end-of-pipe measures for stormwater quality improvement may also need to be considered.

A new geological model of Cenozoic New Zealand suggests a complex phylogeographic history for New Zealand terrestrial taxa

The new model of North Island Cenozoic palaeogeography developed by Kamp et al. has a range of important implications for the evolution of New Zealand terrestrial taxa over the past 30 Ma. Key aspects include the prolonged isolation of the biota on the North Island landmass from the larger and more diverse greater South Island, and the founding of North Island taxa from the potentially unusual ecosystem of a small island around Northland. The prolonged period of isolation is expected to have generated deep phylogenetic splits within taxa present on both islands, and an important current aim should be to identify such signals in surviving endemics to start building a picture of the historical phylogeography, and inferred ecology of both islands through the Cenozoic. Given the potential differences in founding terrestrial species and climatic conditions, it seems likely that the ecology may have been very diferent between the North and South Islands. New genetic data from the 10 or so species of extinct moa suggest that the radiation of moa was much more recent than previously suggested, and reveals a complex pattern that is inferred to result from the interplay of the Cenozoic biogeography, marine barriers, and glacial cycles.

The essential parameters of a resource-based carrying capacity assessment model : an Australian case study

Carrying capacity assessments model a population’s potential self-sufficiency. A crucial first step in the development of such modelling is to examine the basic resource-based parameters defining the population’s production and consumption habits. These parameters include basic human needs such as food, water, shelter and energy together with climatic, environmental and behavioural characteristics. Each of these parameters imparts land-usage requirements in different ways and varied degrees so their incorporation into carrying capacity modelling also differs. Given that the availability and values of production parameters may differ between locations, no two carrying capacity models are likely to be exactly alike. However, the essential parameters themselves can remain consistent so one example, the Carrying Capacity Dashboard, is offered as a case study to highlight one way in which these parameters are utilised. While examples exist of findings made from carrying capacity assessment modelling, to date, guidelines for replication of such studies in other regions and scales have largely been overlooked. This paper addresses such shortcomings by describing a process for the inclusion and calibration of the most important resource-based parameters in a way that could be repeated elsewhere.

Courtyard housing as a subtropical urban design model

Through both theory and practice, this practice-based research develops and tests the idea that courtyard housing can deliver sustainable, compact housing in rapidly growing subtropical cities. It proposes a contemporary urban design model that incorporates urban design courtyard housing prototypes. These prototypes can be further developed by architects and urban designers for similar climatic conditions across the world.

Climate response to physiological forcing of carbon dioxide simulated by the coupled Community Atmosphere Model (CAM3.1) and Community Land Model (CLM3.0)

Increasing concentrations of atmospheric CO2 decrease stomatal conductance of plants and thus suppress canopy transpiration. The climate response to this CO2-physiological forcing is investigated using the Community Atmosphere Model version 3.1 coupled to Community Land Model version 3.0. In response to the physiological effect of doubling CO2, simulations show a decrease in canopy transpiration of 8%, a mean warming of 0.1K over the land surface, and negligible changes in the hydrological cycle. These climate responses are much smaller than what were found in previous modeling studies. This is largely a result of unrealistic partitioning of evapotranspiration in our model control simulation with a greatly underestimated contribution from canopy transpiration and overestimated contributions from canopy and soil evaporation. This study highlights the importance of a realistic simulation of the hydrological cycle, especially the individual components of evapotranspiration, in reducing the uncertainty in our estimation of climatic response to CO2-physiological forcing. Citation: Cao, L., G. Bala, K. Caldeira, R. Nemani, and G.Ban-Weiss (2009), Climate response to physiological forcing of carbon dioxide simulated by the coupled Community Atmosphere Model (CAM3.1) and Community Land Model (CLM3.0).

Modelling climatic risks of aflatoxin contamination in maize

Aflatoxins are highly carcinogenic mycotoxins produced by two fungi, Aspergillus flavus and A. parasiticus, under specific moisture and temperature conditions before harvest and/or during storage of a wide range of crops including maize. Modelling of interactions between host plant and environment during the season can enable quantification of preharvest aflatoxin risk and its potential management. A model was developed to quantify climatic risks of aflatoxin contamination in maize using principles previously used for peanuts. The model outputs an aflatoxin risk index in response to seasonal temperature and soil moisture during the maize grain filling period using the APSIM's maize module. The model performed well in simulating climatic risk of aflatoxin contamination in maize as indicated by a significant R2 (P ≤ 0.01) between aflatoxin risk index and the measured aflatoxin B1 in crop samples, which was 0.69 for a range of rainfed Australian locations and 0.62 when irrigated locations were also included in the analysis. The model was further applied to determine probabilities of exceeding a given aflatoxin risk in four non-irrigated maize growing locations of Queensland using 106 years of historical climatic data. Locations with both dry and hot climates had a much higher probability of higher aflatoxin risk compared with locations having either dry or hot conditions alone. Scenario analysis suggested that under non-irrigated conditions the risk of aflatoxin contamination could be minimised by adjusting sowing time or selecting an appropriate hybrid to better match the grain filling period to coincide with lower temperature and water stress conditions.

Application of a model to assess aflatoxin risk in peanuts.

When exposed to hot (22-35 degrees C) and dry climatic conditions in the field during the final 4-6 weeks of pod filling, peanuts (Arachis hypogaea L.) can accumulate highly carcinogenic and immuno-suppressing aflatoxins. Forecasting of the risk posed by these conditions can assist in minimizing pre-harvest contamination. A model was therefore developed as part of the Agricultural Production Systems Simulator (APSIM) peanut module, which calculated an aflatoxin risk index (ARI) using four temperature response functions when fractional available soil water was <0.20 and the crop was in the last 0.40 of the pod-filling phase. ARI explained 0.95 (P <= 0.05) of the variation in aflatoxin contamination, which varied from 0 to c. 800 mu g/kg in 17 large-scale sowings in tropical and four sowings in sub-tropical environments carried out in Australia between 13 November and 16 December 2007. ARI also explained 0.96 (P <= 0.01) of the variation in the proportion of aflatoxin-contaminated loads (>15 mu g/kg) of peanuts in the Kingaroy region of Australia during the period between the 1998/99 and 2007/08 seasons. Simulation of ARI using historical climatic data from 1890 to 2007 indicated a three-fold increase in its value since 1980 compared to the entire previous period. The increase was associated with increases in ambient temperature and decreases in rainfall. To facilitate routine monitoring of aflatoxin risk by growers in near real time, a web interface of the model was also developed. The ARI predicted using this interface for eight growers correlated significantly with the level of contamination in crops (r=095, P <= 0.01). These results suggest that ARI simulated by the model is a reliable indicator of aflatoxin contamination that can be used in aflatoxin research as well as a decision-support tool to monitor pre-harvest aflatoxin risk in peanuts.

Climatic potential of delta O-18 of Abies spectabilis from the Nepal Himalaya

A 50-year tree-ring delta O-18 chronology of Abies spectabilis growing close to the tree line (3850 m asl) in the Nepal Himalaya is established to explore its dendroclimatic potential. Response function analysis with ambient climatic records revealed that tree-ring delta O-18 is primarily governed by rainfall during the monsoon season (June September), and the regression model accounts for 35% of the variance in rainfall. Extreme dry years identified in instrumental weather data are detected in the delta O-18 chronology. Further, tree-ring delta O-18 is much more sensitive to rainfall fluctuations than other tree-ring parameters such as width and density typically used in dendroclimatology. Correlation analyses with Nino 3.4 SST reveal time-dependent behavior of ENSO-monsoon relationships. (C) 2009 Elsevier GmbH. All rights reserved.

Causality analysis of groundwater dynamics based on a Vector Autoregressive model in the semi-arid basin of Gundal (South India)

Causal relationships existing between observed levels of groundwater in a semi-arid sub-basin of the Kabini River basin (Karnataka state, India) are investigated in this study. A Vector Auto Regressive model is used for this purpose. Its structure is built on an upstream/downstream interaction network based on observed hydro-physical properties. Exogenous climatic forcing is used as an input based on cumulated rainfall departure. Optimal models are obtained thanks to a trial approach and are used as a proxy of the dynamics to derive causal networks. It appears to be an interesting tool for analysing the causal relationships existing inside the basin. The causal network reveals 3 main regions: the Northeastern part of the Gundal basin is closely coupled to the outlet dynamics. The Northwestern part is mainly controlled by the climatic forcing and only marginally linked to the outlet dynamic. Finally, the upper part of the basin plays as a forcing rather than a coupling with the lower part of the basin allowing for a separate analysis of this local behaviour. The analysis also reveals differential time scales at work inside the basin when comparing upstream oriented with downstream oriented causalities. In the upper part of the basin, time delays are close to 2 months in the upward direction and lower than 1 month in the downward direction. These time scales are likely to be good indicators of the hydraulic response time of the basin which is a parameter usually difficult to estimate practically. This suggests that, at the sub-basin scale, intra-annual time scales would be more relevant scales for analysing or modelling tropical basin dynamics in hard rock (granitic and gneissic) aquifers ubiquitous in south India. (c) 2012 Elsevier B.V. All rights reserved.

A simple model for spatial disaggregation of evaporative fraction: Comparative study with thermal sharpened land surface temperature data over India

[1] Evaporative fraction (EF) is a measure of the amount of available energy at the earth surface that is partitioned into latent heat flux. The currently operational thermal sensors like the Moderate Resolution Imaging Spectroradiometer (MODIS) on satellite platforms provide data only at 1000 m, which constraints the spatial resolution of EF estimates. A simple model (disaggregation of evaporative fraction (DEFrac)) based on the observed relationship between EF and the normalized difference vegetation index is proposed to spatially disaggregate EF. The DEFrac model was tested with EF estimated from the triangle method using 113 clear sky data sets from the MODIS sensor aboard Terra and Aqua satellites. Validation was done using the data at four micrometeorological tower sites across varied agro-climatic zones possessing different land cover conditions in India using Bowen ratio energy balance method. The root-mean-square error (RMSE) of EF estimated at 1000 m resolution using the triangle method was 0.09 for all the four sites put together. The RMSE of DEFrac disaggregated EF was 0.09 for 250 m resolution. Two models of input disaggregation were also tried with thermal data sharpened using two thermal sharpening models DisTrad and TsHARP. The RMSE of disaggregated EF was 0.14 for both the input disaggregation models for 250 m resolution. Moreover, spatial analysis of disaggregation was performed using Landsat-7 (Enhanced Thematic Mapper) ETM+ data over four grids in India for contrasted seasons. It was observed that the DEFrac model performed better than the input disaggregation models under cropped conditions while they were marginally similar under non-cropped conditions.

Glacial mass balance changes in the Karakoram and Himalaya based on CMIP5 multi-model climate projections

The impact of future climate change on the glaciers in the Karakoram and Himalaya (KH) is investigated using CMIP5 multi-model temperature and precipitation projections, and a relationship between glacial accumulation-area ratio and mass balance developed for the region based on the last 30 to 40 years of observational data. We estimate that the current glacial mass balance (year 2000) for the entire KH region is -6.6 +/- 1 Gta(-1), which decreases about sixfold to -35 +/- 2 Gta(-1) by the 2080s under the high emission scenario of RCP8.5. However, under the low emission scenario of RCP2.6 the glacial mass loss only doubles to -12 +/- 2 Gta(-1) by the 2080s. We also find that 10.6 and 27 % of the glaciers could face `eventual disappearance' by the end of the century under RCP2.6 and RCP8.5 respectively, underscoring the threat to water resources under high emission scenarios.

The variations and climatic significance of D/H ratios in the carbon-bound hydrogen of cellulose in trees

The σD values of nitrated cellulose from a variety of trees covering a wide geographic range have been measured. These measurements have been used to ascertain which factors are likely to cause σD variations in cellulose C-H hydrogen.

It is found that a primary source of tree σD variation is the σD variation of the environmental precipitation. Superimposed on this are isotopic variations caused by the transpiration of the leaf water incorporated by the tree. The magnitude of this transpiration effect appears to be related to relative humidity.

Within a single tree, it is found that the hydrogen isotope variations which occur for a ring sequence in one radial direction may not be exactly the same as those which occur in a different direction. Such heterogeneities appear most likely to occur in trees with asymmetric ring patterns that contain reaction wood. In the absence of reaction wood such heterogeneities do not seem to occur. Thus, hydrogen isotope analyses of tree ring sequences should be performed on trees which do not contain reaction wood.

Comparisons of tree σD variations with variations in local climate are performed on two levels: spatial and temporal. It is found that the σD values of 20 North American trees from a wide geographic range are reasonably well-correlated with the corresponding average annual temperature. The correlation is similar to that observed for a comparison of the σD values of annual precipitation of 11 North American sites with annual temperature. However, it appears that this correlation is significantly disrupted by trees which grew on poorly drained sites such as those in stagnant marshes. Therefore, site selection may be important in choosing trees for climatic interpretation of σD values, although proper sites do not seem to be uncommon.

The measurement of σD values in 5-year samples from the tree ring sequences of 13 trees from 11 North American sites reveals a variety of relationships with local climate. As it was for the spatial σD vs climate comparison, site selection is also apparently important for temporal tree σD vs climate comparisons. Again, it seems that poorly-drained sites are to be avoided. For nine trees from different "well-behaved" sites, it was found that the local climatic variable best related to the σD variations was not the same for all sites.

Two of these trees showed a strong negative correlation with the amount of local summer precipitation. Consideration of factors likely to influence the isotopic composition of summer rain suggests that rainfall intensity may be important. The higher the intensity, the lower the σD value. Such an effect might explain the negative correlation of σD vs summer precipitation amount for these two trees. A third tree also exhibited a strong correlation with summer climate, but in this instance it was a positive correlation of σD with summer temperature.

The remaining six trees exhibited the best correlation between σD values and local annual climate. However, in none of these six cases was it annual temperature that was the most important variable. In fact annual temperature commonly showed no relationship at all with tree σD values. Instead, it was found that a simple mass balance model incorporating two basic assumptions yielded parameters which produced the best relationships with tree σD values. First, it was assumed that the σD values of these six trees reflected the σD values of annual precipitation incorporated by these trees. Second, it was assumed that the σD value of the annual precipitation was a weighted average of two seasonal isotopic components: summer and winter. Mass balance equations derived from these assumptions yielded combinations of variables that commonly showed a relationship with tree σD values where none had previously been discerned.

It was found for these "well-behaved" trees that not all sample intervals in a σD vs local climate plot fell along a well-defined trend. These departures from the local σD VS climate norm were defined as "anomalous". Some of these anomalous intervals were common to trees from different locales. When such widespread commonalty of an anomalous interval occurred, it was observed that the interval corresponded to an interval in which drought had existed in the North American Great Plains.

Consequently, there appears to be a combination of both local and large scale climatic information in the σD variations of tree cellulose C-H hydrogen.

Growth and survival of sockeye salmon (Oncorhynchus nerka) from Karluk Lake and River, Alaska, in relation to climatic and oceanic regimes and indices, 1922–2000

We examined whether the relationship between climate and salmon production was linked through the effect of climate on the growth of sockeye salmon (Oncorhynchus nerka) at sea. Smolt length and juvenile, immature, and maturing growth rates were estimated from increments on scales of adult sockeye salmon that returned to the Karluk River and Lake system on Kodiak Island, Alaska, over 77 years, 1924–2000. Survival was higher during the warm climate regimes and lower during the cool regime. Growth was not correlated with survival, as estimated from the residuals of the Ricker stock-recruitment model. Juvenile growth was correlated with an atmospheric forcing index and immature growth was correlated with the amount of coastal precipitation, but the magnitude of winter and spring coastal downwelling in the Gulf of Alaska and the Pacific Northwest atmospheric patterns that influence the directional bifurcation of the Pacific Current were not related to the growth of Karluk sockeye salmon. However, indices of sea surface temperature, coastal precipitation, and atmospheric circulation in the eastern North Pacific were correlated with the survival of Karluk sockeye salmon. Winter and spring precipitation and atmospheric circulation are possible processes linking survival to climate variation in Karluk sockeye salmon.

High-resolution natural recording systems: intercomparisons of the response to interannual climatic variability [abstract]

The goal of this work is to examine the properties of recording mechanisms which are common to continuously recording high-resolution natural systems in which climatic signals are imprinted and preserved as proxy records. These systems produce seasonal structures as an indirect response to climatic variability over the annual cycle. We compare the proxy records from four different high-resolution systems: the Quelccaya ice cap of the Peruvian Andes; composite tree ring growth from southern California and the southwestern United States; and the marine varve sedimentation systems in the Santa Barbara basin (off California, United States) and in the Gulf of California, Mexico. An important focus of this work is to indicate how the interannual climatic signal is recorded in a variety of different natural systems with vastly different recording mechanisms and widely separated in space. These high-resolution records are the products of natural processes which should be comparable, to some degree, to human-engineered systems developed to transmit and record physical quantities. We therefore present a simple analogy of a data recording system as a heuristic model to provide some unifying concepts with which we may better understand the formation of the records. This analogy assumes special significance when we consider that natural proxy records are the principal means to extend our knowledge of climatic variability into the past, beyond the limits of instrumentally recorded data.

Atmospheric carbon dioxide and the climatic record

This paper is an attempt to provide a summary review of conclusions from previous studies on this subject. They have been organized under the following subject headings: Conceptualization of the greenhouse effect; The climatic effect of doubled carbon dioxide; Interpretation of the climatic record; Diagnosis of apparent and possible model deficiencies; The paleoclimatic record.