Taxonomy for rainfall events based on pollutant wash-off potential in urban areas

Conventional rainfall classification for modelling and prediction is quantity based. This approach can lead to inaccuracies in stormwater quality modelling due to the assignment of stochastic pollutant parameters to a rainfall event. A taxonomy for natural rainfall events in the context of stormwater quality is presented based on an in-depth investigation of the influence of rainfall characteristics on stormwater quality. In the research study, the natural rainfall events were classified into three types based on average rainfall intensity and rainfall duration and the classification was found to be independent of the catchment characteristics. The proposed taxonomy provides an innovative concept in stormwater quality modelling and prediction and will contribute to enhancing treatment design for stormwater quality mitigation.

Impact of roof surface runoff on urban water quality

The pollutant impacts of urban stormwater runoff on receiving waters are well documented in research literature. However, it is road surfaces that are commonly identified as the significant pollutant source. This paper presents the outcomes of an extensive program of research into the role of roof surfaces in urban water quality with particular focus on solids, nutrients and organic carbon. The outcomes confirmed that roof surfaces play an important role in influencing the pollutant characteristics of urban stormwater runoff. Pollutant build-up and wash-off characteristics for roads and roof surfaces were found to be appreciably different. The pollutant wash-off characteristics exhibited by roof surfaces show that it influences the first flush phenomenon more significantly than road surfaces. In most urban catchments, as roof surfaces constitutes a higher fraction of impervious area compared to road surfaces, it is important that the pollutant generation role of roof surfaces is specifically taken into consideration in stormwater quality mitigation strategies.

Inherent errors in pollutant build-up estimation in considering urban land use as a lumped parameter

Stormwater quality modelling results is subject to uncertainty. The variability of input parameters is an important source of overall model error. An in-depth understanding of the variability associated with input parameters can provide knowledge on the uncertainty associated with these parameters and consequently assist in uncertainty analysis of stormwater quality models and the decision making based on modelling outcomes. This paper discusses the outcomes of a research study undertaken to analyse the variability related to pollutant build-up parameters in stormwater quality modelling. The study was based on the analysis of pollutant build-up samples collected from 12 road surfaces in residential, commercial and industrial land uses. It was found that build-up characteristics vary appreciably even within the same land use. Therefore, using land use as a lumped parameter would contribute significant uncertainties in stormwater quality modelling. Additionally, it was also found that the variability in pollutant build-up can also be significant depending on the pollutant type. This underlines the importance of taking into account specific land use characteristics and targeted pollutant species when undertaking uncertainty analysis of stormwater quality models or in interpreting the modelling outcomes.

Performance characterisation of a constructed wetland

Performance of a constructed wetland is commonly reported as variable due to the site specific nature of influential factors. This paper discusses outcomes from an in-depth study which characterised treatment performance of a wetland based on the variation in runoff regime. The study included a comprehensive field monitoring of a well established constructed wetland in Gold Coast, Australia. Samples collected at the inlet and outlet was tested for Total Suspended Solids (TSS), Total Nitrogen (TN) and Total Phosphorus (TP). Pollutant concentrations in the outflow were found to be consistent irrespective of the variation in inflow water quality. The analysis revealed two different treatment characteristics for events with different rainfall depths. TSS and TN load reduction is strongly influenced by hydraulic retention time where performance is higher for rainfall events below the design event. For small events, treatment performance is higher at the beginning of the event and gradually decreased during the course of the event. For large events, the treatment performance is comparatively poor at the beginning and improved during the course of the event. The analysis also confirmed the variable treatment trends for different pollutant types.

Influence of Rainfall and Catchment Characteristics on Urban Stormwater Quality

The accuracy and reliability of urban stormwater quality modelling outcomes are important for stormwater management decision making. The commonly adopted approach where only a limited number of factors are used to predict urban stormwater quality may not adequately represent the complexity of the quality response to a rainfall event or site-to-site differences to support efficient treatment design. This paper discusses an investigation into the influence of rainfall and catchment characteristics on urban stormwater quality in order to investigate the potential areas for errors in current stormwater quality modelling practices. It was found that the influence of rainfall characteristics on pollutant wash-off is step-wise based on specific thresholds. This means that a modelling approach where the wash-off process is predicted as a continuous function of rainfall intensity and duration is not appropriate. Additionally, other than conventional catchment characteristics, namely, land use and impervious surface fraction, other catchment characteristics such as impervious area layout, urban form and site specific characteristics have an important influence on both, pollutant build-up and wash-off processes. Finally, the use of solids as a surrogate to estimate other pollutant species was found to be inappropriate. Individually considering build-up and wash-off processes for each pollutant species should be the preferred option.

Climate variability over the last 35,000 years recorded in marine and terrestrial archives in the Australian region : an OZ-INTIMATE compilation

The Australian region spans some 60° of latitude and 50° of longitude and displays considerable regional climate variability both today and during the Late Quaternary. A synthesis of marine and terrestrial climate records, combining findings from the Southern Ocean, temperate, tropical and arid zones, identifies a complex response of climate proxies to a background of changing boundary conditions over the last 35,000 years. Climate drivers include the seasonal timing of insolation, greenhouse gas content of the atmosphere, sea level rise and ocean and atmospheric circulation changes. Our compilation finds few climatic events that could be used to construct a climate event stratigraphy for the entire region, limiting the usefulness of this approach. Instead we have taken a spatial approach, looking to discern the patterns of change across the continent. The data identify the clearest and most synchronous climatic response at the time of the Last Glacial Maximum (LGM) (21 ± 3 ka), with unambiguous cooling recorded in the ocean, and evidence of glaciation in the highlands of tropical New Guinea, southeast Australia and Tasmania. Many terrestrial records suggest drier conditions, but with the timing of inferred snowmelt, and changes to the rainfall/runoff relationships, driving higher river discharge at the LGM. In contrast, the deglaciation is a time of considerable south-east to north-west variation across the region. Warming was underway in all regions by 17 ka. Post-glacial sea level rise and its associated regional impacts have played an important role in determining the magnitude and timing of climate response in the north-west of the continent in contrast to the southern latitudes. No evidence for cooling during the Younger Dryas chronozone is evident in the region, but the Antarctic cold reversal clearly occurs south of Australia. The Holocene period is a time of considerable climate variability associated with an intense monsoon in the tropics early in the Holocene, giving way to a weakened monsoon and an increasingly El Niño-dominated ENSO to the present. The influence of ENSO is evident throughout the southeast of Australia, but not the southwest. This climate history provides a template from which to assess the regionality of climate events across Australia and make comparisons beyond our region. The data identify the clearest and most synchronous climatic response at the time of the Last Glacial Maximum (LGM) (21 ± 3 ka), with unambiguous cooling recorded in the ocean, and evidence of glaciation in the highlands of tropical New Guinea, southeast Australia and Tasmania. Many terrestrial records suggest drier conditions, but with the timing of inferred snowmelt, and changes to the rainfall/runoff relationships, driving higher river discharge at the LGM. In contrast, the deglaciation is a time of considerable south-east to north-west variation across the region. Warming was underway in all regions by 17 ka. Post-glacial sea level rise and its associated regional impacts have played an important role in determining the magnitude and timing of climate response in the north-west of the continent in contrast to the southern latitudes. No evidence for cooling during the Younger Dryas chronozone is evident in the region, but the Antarctic cold reversal clearly occurs south of Australia. The Holocene period is a time of considerable climate variability associated with an intense monsoon in the tropics early in the Holocene, giving way to a weakened monsoon and an increasingly El Niño-dominated ENSO to the present. The influence of ENSO is evident throughout the southeast of Australia, but not the southwest. This climate history provides a template from which to assess the regionality of climate events across Australia and make comparisons beyond our region.

Quantitative genetic parameter estimates for body and carcass traits in a cultured stock of giant freshwater prawn (Macrobrachium rosenbergii) selected for harvest weight in Vietnam

We estimated the heritability and correlations between body and carcass weight traits in a cultured stock of giant freshwater prawn (GFP) (Macrobrachium rosenbergii) selected for harvest body weight in Vietnam. The data set consisted of 18,387 body and 1,730 carcass records, as well as full pedigree information collected over four generations. Variance and covariance components were estimated by restricted maximum likelihood fitting a multi-trait animal model. Across generations, estimates of heritability for body and carcass weight traits were moderate and ranged from 0.14 to 0.19 and 0.17 to 0.21, respectively. Body trait heritabilities estimated for females were significantly higher than for males whereas carcass weight trait heritabilities estimated for females and males were not significantly different (P>. 0.05). Maternal effects for body traits accounted for 4 to 5% of the total variance and were greater in females than in males. Genetic correlations among body traits were generally high in the mixed sexes. Genetic correlations between body and carcass weight traits were also high. Although some issues remain regarding the best statistical model to be fitted to GFP data, our results suggest that selection for high harvest body weight based on breeding values estimated by fitting an animal model to the data can significantly improve mean body and carcass weight in GFP.

Sunshine, rainfall, humidity and child pneumonia in the tropics : time-series analyses

Few studies have formally examined the relationship between meteorological factors and the incidence of child pneumonia in the tropics, despite the fact that most child pneumonia deaths occur there. We examined the association between four meteorological exposures (rainy days, sunshine, relative humidity, temperature) and the incidence of clinical pneumonia in young children in the Philippines using three time-series methods: correlation of seasonal patterns, distributed lag regression, and case-crossover. Lack of sunshine was most strongly associated with pneumonia in both lagged regression [overall relative risk over the following 60 days for a 1-h increase in sunshine per day was 0·67 (95% confidence interval (CI) 0·51–0·87)] and case-crossover analysis [odds ratio for a 1-h increase in mean daily sunshine 8–14 days earlier was 0·95 (95% CI 0·91–1·00)]. This association is well known in temperate settings but has not been noted previously in the tropics. Further research to assess causality is needed.

First flush analysis in urban catchments

The effectiveness of structural elements employed for stormwater mitigation such as bioretention basins and constructed wetlands depend on the compatibility between their design specifications and actual stormwater quality and quantity characteristics. These structural elements are commonly designed to accommodate the initial portion of runoff considering the occurrence of first flush. Therefore, the effectiveness of stormwater quality treatment primarily depends on the in-depth knowledge of the first flush phenomenon and the ability to provide appropriate treatment. The current scientific knowledge relating to first flush is limited primarily due to research investigations being undertaken based on lumped rainfall and runoff parameters. This paper presents the outcomes of an in-depth study undertaken of the first flush phenomenon using a set of indicators which are not only innovative, but is also able to accurately represent the characteristics of the different sectors in a runoff hydrograph. The analysis undertaken confirmed that pollutant wash-off during the initial 10% of runoff volume was critical for the occurrence of first flush. Typically first flush was found to last up to 40% of the runoff volume. The study outcomes provide new knowledge to enhance the effectiveness of structural stormwater treatment measures.

Real - time prediction of rainfall induced instability of residual soil slopes associated with deep cracks

The early warning based on real-time prediction of rain-induced instability of natural residual slopes helps to minimise human casualties due to such slope failures. Slope instability prediction is complicated, as it is influenced by many factors, including soil properties, soil behaviour, slope geometry, and the location and size of deep cracks in the slope. These deep cracks can facilitate rainwater infiltration into the deep soil layers and reduce the unsaturated shear strength of residual soil. Subsequently, it can form a slip surface, triggering a landslide even in partially saturated soil slopes. Although past research has shown the effects of surface-cracks on soil stability, research examining the influence of deep-cracks on soil stability is very limited. This study aimed to develop methodologies for predicting the real-time rain-induced instability of natural residual soil slopes with deep cracks. The results can be used to warn against potential rain-induced slope failures. The literature review conducted on rain induced slope instability of unsaturated residual soil associated with soil crack, reveals that only limited studies have been done in the following areas related to this topic: - Methods for detecting deep cracks in residual soil slopes. - Practical application of unsaturated soil theory in slope stability analysis. - Mechanistic methods for real-time prediction of rain induced residual soil slope instability in critical slopes with deep cracks. Two natural residual soil slopes at Jombok Village, Ngantang City, Indonesia, which are located near a residential area, were investigated to obtain the parameters required for the stability analysis of the slope. A survey first identified all related field geometrical information including slope, roads, rivers, buildings, and boundaries of the slope. Second, the electrical resistivity tomography (ERT) method was used on the slope to identify the location and geometrical characteristics of deep cracks. The two ERT array models employed in this research are: Dipole-dipole and Azimuthal. Next, bore-hole tests were conducted at different locations in the slope to identify soil layers and to collect undisturbed soil samples for laboratory measurement of the soil parameters required for the stability analysis. At the same bore hole locations, Standard Penetration Test (SPT) was undertaken. Undisturbed soil samples taken from the bore-holes were tested in a laboratory to determine the variation of the following soil properties with the depth: - Classification and physical properties such as grain size distribution, atterberg limits, water content, dry density and specific gravity. - Saturated and unsaturated shear strength properties using direct shear apparatus. - Soil water characteristic curves (SWCC) using filter paper method. - Saturated hydraulic conductivity. The following three methods were used to detect and simulate the location and orientation of cracks in the investigated slope: (1) The electrical resistivity distribution of sub-soil obtained from ERT. (2) The profile of classification and physical properties of the soil, based on laboratory testing of soil samples collected from bore-holes and visual observations of the cracks on the slope surface. (3) The results of stress distribution obtained from 2D dynamic analysis of the slope using QUAKE/W software, together with the laboratory measured soil parameters and earthquake records of the area. It was assumed that the deep crack in the slope under investigation was generated by earthquakes. A good agreement was obtained when comparing the location and the orientation of the cracks detected by Method-1 and Method-2. However, the simulated cracks in Method-3 were not in good agreement with the output of Method-1 and Method-2. This may have been due to the material properties used and the assumptions made, for the analysis. From Method-1 and Method-2, it can be concluded that the ERT method can be used to detect the location and orientation of a crack in a soil slope, when the ERT is conducted in very dry or very wet soil conditions. In this study, the cracks detected by the ERT were used for stability analysis of the slope. The stability of the slope was determined using the factor of safety (FOS) of a critical slip surface obtained by SLOPE/W using the limit equilibrium method. Pore-water pressure values for the stability analysis were obtained by coupling the transient seepage analysis of the slope using finite element based software, called SEEP/W. A parametric study conducted on the stability of an investigated slope revealed that the existence of deep cracks and their location in the soil slope are critical for its stability. The following two steps are proposed to predict the rain-induced instability of a residual soil slope with cracks. (a) Step-1: The transient stability analysis of the slope is conducted from the date of the investigation (initial conditions are based on the investigation) to the preferred date (current date), using measured rainfall data. Then, the stability analyses are continued for the next 12 months using the predicted annual rainfall that will be based on the previous five years rainfall data for the area. (b) Step-2: The stability of the slope is calculated in real-time using real-time measured rainfall. In this calculation, rainfall is predicted for the next hour or 24 hours and the stability of the slope is calculated one hour or 24 hours in advance using real time rainfall data. If Step-1 analysis shows critical stability for the forthcoming year, it is recommended that Step-2 be used for more accurate warning against the future failure of the slope. In this research, the results of the application of the Step-1 on an investigated slope (Slope-1) showed that its stability was not approaching a critical value for year 2012 (until 31st December 2012) and therefore, the application of Step-2 was not necessary for the year 2012. A case study (Slope-2) was used to verify the applicability of the complete proposed predictive method. A landslide event at Slope-2 occurred on 31st October 2010. The transient seepage and stability analyses of the slope using data obtained from field tests such as Bore-hole, SPT, ERT and Laboratory tests, were conducted on 12th June 2010 following the Step-1 and found that the slope in critical condition on that current date. It was then showing that the application of the Step-2 could have predicted this failure by giving sufficient warning time.

First flush behaviour in urban residential catchments

The current state of knowledge in relation to first flush does not provide a clear understanding of the role of rainfall and catchment characteristics in influencing this phenomenon. This is attributed to the inconsistent findings from research studies due to the unsatisfactory selection of first flush indicators and how first flush is defined. The research study discussed in this thesis provides the outcomes of a comprehensive analysis on the influence of rainfall and catchment characteristics on first flush behaviour in residential catchments. Two sets of first flush indicators are introduced in this study. These indicators were selected such that they are representative in explaining in a systematic manner the characteristics associated with first flush. Stormwater samples and rainfall-runoff data were collected and recorded from stormwater monitoring stations established at three urban catchments at Coomera Waters, Gold Coast, Australia. In addition, historical data were also used to support the data analysis. Three water quality parameters were analysed, namely, total suspended solids (TSS), total phosphorus (TP) and total nitrogen (TN). The data analyses were primarily undertaken using multi criteria decision making methods, PROMETHEE and GAIA. Based on the data obtained, the pollutant load distribution curve (LV) was determined for the individual rainfall events and pollutant types. Accordingly, two sets of first flush indicators were derived from the curve, namely, cumulative load wash-off for every 10% of runoff volume interval (interval first flush indicators or LV) from the beginning of the event and the actual pollutant load wash-off during a 10% increment in runoff volume (section first flush indicators or P). First flush behaviour showed significant variation with pollutant types. TSS and TP showed consistent first flush behaviour. However, the dissolved fraction of TN showed significant differences to TSS and TP first flush while particulate TN showed similarities. Wash-off of TSS, TP and particulate TN during the first 10% of the runoff volume showed no influence from corresponding rainfall intensity. This was attributed to the wash-off of weakly adhered solids on the catchment surface referred to as "short term pollutants" or "weakly adhered solids" load. However, wash-off after 10% of the runoff volume showed dependency on the rainfall intensity. This is attributed to the wash-off of strongly adhered solids being exposed when the weakly adhered solids diminish. The wash-off process was also found to depend on rainfall depth at the end part of the event as the strongly adhered solids are loosened due to impact of rainfall in the earlier part of the event. Events with high intensity rainfall bursts after 70% of the runoff volume did not demonstrate first flush behaviour. This suggests that rainfall pattern plays a critical role in the occurrence of first flush. Rainfall intensity (with respect to the rest of the event) that produces 10% to 20% runoff volume play an important role in defining the magnitude of the first flush. Events can demonstrate high magnitude first flush when the rainfall intensity occurring between 10% and 20% of the runoff volume is comparatively high while low rainfall intensities during this period produces low magnitude first flush. For events with first flush, the phenomenon is clearly visible up to 40% of the runoff volume. This contradicts the common definition that first flush only exists, if for example, 80% of the pollutant mass is transported in the first 30% of runoff volume. First flush behaviour for TN is different compared to TSS and TP. Apart from rainfall characteristics, the composition and the availability of TN on the catchment also play an important role in first flush. The analysis confirmed that events with low rainfall intensity can produce high magnitude first flush for the dissolved fraction of TN, while high rainfall intensity produce low dissolved TN first flush. This is attributed to the source limiting behaviour of dissolved TN wash-off where there is high wash-off during the initial part of a rainfall event irrespective of the intensity. However, for particulate TN, the influence of rainfall intensity on first flush characteristics is similar to TSS and TP. The data analysis also confirmed that first flush can occur as high magnitude first flush, low magnitude first flush or non existence of first flush. Investigation of the influence of catchment characteristics on first flush found that the key factors that influence the phenomenon are the location of the pollutant source, spatial distribution of the pervious and impervious surfaces in the catchment, drainage network layout and slope of the catchment. This confirms that first flush phenomenon cannot be evaluated based on a single or a limited set of parameters as a number of catchment characteristics should be taken into account. Catchments where the pollutant source is located close to the outlet, a high fraction of road surfaces, short travel time to the outlet, with steep slopes can produce high wash-off load during the first 50% of the runoff volume. Rainfall characteristics have a comparatively dominant impact on the wash-off process compared to the catchment characteristics. In addition, the pollutant characteristics also should be taken into account in designing stormwater treatment systems due to different wash-off behaviour. Analysis outcomes confirmed that there is a high TSS load during the first 20% of the runoff volume followed by TN which can extend up to 30% of the runoff volume. In contrast, high TP load can exist during the initial and at the end part of a rainfall event. This is related to the composition of TP available for the wash-off.

Consistent HMM parameter estimation using Kerridge inaccuracy rates

In this paper, we propose a novel online hidden Markov model (HMM) parameter estimator based on Kerridge inaccuracy rate (KIR) concepts. Under mild identifiability conditions, we prove that our online KIR-based estimator is strongly consistent. In simulation studies, we illustrate the convergence behaviour of our proposed online KIR-based estimator and provide a counter-example illustrating the local convergence properties of the well known recursive maximum likelihood estimator (arguably the best existing solution).

Plasmodium vivax malaria incidence over time and its association with temperature and rainfall in four counties of Yunnan Province, China

Background Transmission of Plasmodium vivax malaria is dependent on vector availability, biting rates and parasite development. In turn, each of these is influenced by climatic conditions. Correlations have previously been detected between seasonal rainfall, temperature and malaria incidence patterns in various settings. An understanding of seasonal patterns of malaria, and their weather drivers, can provide vital information for control and elimination activities. This research aimed to describe temporal patterns in malaria, rainfall and temperature, and to examine the relationships between these variables within four counties of Yunnan Province, China. Methods Plasmodium vivax malaria surveillance data (1991–2006), and average monthly temperature and rainfall were acquired. Seasonal trend decomposition was used to examine secular trends and seasonal patterns in malaria. Distributed lag non-linear models were used to estimate the weather drivers of malaria seasonality, including the lag periods between weather conditions and malaria incidence. Results There was a declining trend in malaria incidence in all four counties. Increasing temperature resulted in increased malaria risk in all four areas and increasing rainfall resulted in increased malaria risk in one area and decreased malaria risk in one area. The lag times for these associations varied between areas. Conclusions The differences detected between the four counties highlight the need for local understanding of seasonal patterns of malaria and its climatic drivers.

A Bayesian regression approach to assess uncertainty in pollutant wash-off modelling

Due to knowledge gaps in relation to urban stormwater quality processes, an in-depth understanding of model uncertainty can enhance decision making. Uncertainty in stormwater quality models can originate from a range of sources such as the complexity of urban rainfall-runoff-stormwater pollutant processes and the paucity of observed data. Unfortunately, studies relating to epistemic uncertainty, which arises from the simplification of reality are limited and often deemed mostly unquantifiable. This paper presents a statistical modelling framework for ascertaining epistemic uncertainty associated with pollutant wash-off under a regression modelling paradigm using Ordinary Least Squares Regression (OLSR) and Weighted Least Squares Regression (WLSR) methods with a Bayesian/Gibbs sampling statistical approach. The study results confirmed that WLSR assuming probability distributed data provides more realistic uncertainty estimates of the observed and predicted wash-off values compared to OLSR modelling. It was also noted that the Bayesian/Gibbs sampling approach is superior compared to the most commonly adopted classical statistical and deterministic approaches commonly used in water quality modelling. The study outcomes confirmed that the predication error associated with wash-off replication is relatively higher due to limited data availability. The uncertainty analysis also highlighted the variability of the wash-off modelling coefficient k as a function of complex physical processes, which is primarily influenced by surface characteristics and rainfall intensity.