Environmental impact of coal mining and coal seam gas production on surface water quality in the Sydney basin,Australia

The extraction of coal and coal seam gas (CSG) will generate produced water that, if not adequately treated, will pollute surface and groundwater systems. In Australia, the discharge of produced water from coal mining and related activities is regulated by the state environment agency through a pollution licence. This licence sets the discharge limits for a range of analytes to protect the environment into which the produced water is discharged. This study reports on the impact of produced water from coal mine activities located within or discharging into high conservation environments, such as National Parks, in the outer region of Sydney, Australia. The water samples upstream and downstream from the discharge points from six mines were taken, and 110 parameters were tested. The results were assessed against a water quality index (WQI) which accounts for pH, turbidity, dissolved oxygen, biochemical oxygen demand, total dissolved solids, total phosphorus, nitrate nitrogen and E .coli. The water quality assessment based on the trace metal contents against various national maximum admissible concentration (MAC) and their corresponding environmental impacts was also included in the study which also established a base value of water quality for further study. The study revealed that impacted water downstream of the mine discharge points contained higher metal content than the upstream reference locations. In many cases, the downstream water was above the Australia and New Zealand Environment Conservation Council and international water quality guidelines for freshwater stream. The major outliers to the guidelines were aluminium (Al), iron (Fe), manganese (Mn), nickel (Ni) and zinc (Zn). The WQI of surface water at and downstream of the discharge point was lower when compared to upstream or reference conditions in the majority of cases. Toxicology indices of metals present in industrial discharges were used as an additional tool to assess water quality, and the newly proposed environmental water quality index (EWQI) lead to better trend in the impact of coal and coal seam gas mining activities on surface water quality when compared to the upstream reference water samples. Metal content limits were based on the impact points assigned by the Agency for Toxic Substances and Disease Registry, USA. For environmental and health impact assessment, the approach used in this study can be applied as a model to provide a basis to assess the anthropogenic contribution from the industrial and mining activities on the environment.     

An Improved Rating System for Assessing Surface Water Contamination Potential from MSW Landfills

Leachate emission from uncontrolled municipal solid waste landfills (referred to as waste sites in the present study) is a major threat to the environment and living beings in its vicinity. Surface water contamination potential resulting from leachate may be used as one of the criteria for prioritization of sites for remediation purposes. The existing hazard rating systems that prioritize waste sites considering surface water contamination potential as one of the criteria are mainly suited for the developed countries where these were developed initially. In developing countries like India, the set of conditions differ from those in developed countries, and therefore the existing systems may not be suitable for developing countries. Thus in the present study, an improved system is proposed to assess surface water contamination potential from MSW sites. The system is based on the concept of Source, Pathway and Receptor. The proposed system employs parameters derived from the review of existing rating systems and selects their best and worst values based on literature review, design standards and field conditions. The importance weights of the system parameters have been decided based on expert judgment using Delphi technique. Sensitivity analysis of the system shows that the improved system is more sensitive than the existing systems for the site conditions encountered in developing countries. Monte Carlo analysis of the proposed system confirms the spread of the scores obtained from the system over the full scale of 0–1000. The improved system is compared with existing systems by applying it to waste sites from metropolitan cities of India and performing clustering analysis on the rating scores. The clustering analysis shows that the rating scores from the improved system are less clustered as compared to the scores from the existing systems. This demonstrates that the improved system makes a better tool to distinctly prioritize the waste sites for remediation purpose.     

Simulating Nitrate Leaching Profiles in a Highly Permeable Vadose Zone

An approach is developed to simulate leaching of a dissolved chemical constituent in the vadose zone of an aquifer. Specifically, nitrate loading at the water table for different water table depths, for a range of aquifer permeability values, and for different cases of heterogeneity of the aquifer, are considered. Models from the literature are first used to derive soil–water characteristic curves (water retention and hydraulic conductivity) from a grain size distribution curve for unsaturated conditions. Given infiltration from the surface, the initial conditions for the chemical concentration, and the water content profile, leaching of the chemical in the vadose zone is simulated as a function of both time and depth. The methodology is illustrated for a permeable aquifer. Simulations are undertaken using a finite element code for saturated and unsaturated flow. Different scenarios are simulated depending on the heterogeneity of the aquifer and the depth of the water table. Modeling results show that in the example case studied, nitrate concentration loading at the water table does not depend strongly on the position of the water table, but rather on the material properties of the aquifer. The contribution of this endeavor resides in the methodology which allows a prediction of nitrate leaching using only the grain size property of the aquifer. It allows practitioners to obtain a first assessment of leaching with limited data.     

Testing the MAGIC acid rain model in highly organic,low-conductivity waters using multiple calibrations

Accurate predictions of acid precipitation effects on water resources are important in order to allow a better understanding of various pollution control strategy outcomes. Dynamic geochemical models have been developed to address this need, but have to be tested under a variety of environmental conditions to provide confidence in their predictions. The most commonly used aquatic acidification model in North America and Europe is the model of acidification of groundwater in catchments (MAGIC). Though extensively used, MAGIC has never been tested in catchments with extremely low ionic strength water and high in natural organic acids (NOAs) from wetlands, two conditions which are common in large parts of Canada. We calibrated the model for two catchments located in Nova Scotia, Canada, which had some of the most dilute freshwaters reported in the literature and very high NOA. We also evaluated the variability inherent in calibration data sets by using five separate data sets collected over a 15-year period at the same sites. We show good model simulations for the main cations and anions in catchment waters. However, modeling pH is more difficult in the highly organic waters and requires modification to the acid dissociation constants. Calculated acid neutralization capacity can also be more difficult to model due to the low ion content making small errors more important. In theory, multiple calibrations of a model at a same site should produce identical hindcasts and predictions. In reality, the multiple calibrations produced a series of similar, but not identical outcomes which give a probable range of past values and future outcomes. We feel that this practical approach to validation is a useful addition to the arsenal of model testing tools.     

A Compartmental Model for Stable Time-Dependent Modeling of Surface Water and Groundwater

Water resource system planners make decisions that guide water management policy. The fundamental tools for assessing management and infrastructure strategies are combined hydro-economic models of river basins (RBHE models). These models have improved the economic efficiency of water use in situations of competition for scarce water resources. In RBHE models, a groundwater model is coupled with surface water models of the various water resources. Today, the groundwater models used in an RBHE model can be of two types: cell models or numerical models. Cell models are easy to use, but they are too simple to realistically describe the geology and hydrology of the area under investigation. Numerical models, in contrast, are closer to the physical behavior of the aquifer. However, the vast quantity of data to be analyzed makes them impractical for many management scenarios. Moreover, the calibrations of these high-resolution models are generally difficult and sensitive to the variation of parameters, especially when boundary conditions are dynamic. This is the case when dynamic river data or dynamic surface lake data are present. In this work, a compartmental cell model is built on the hydrogeology of the aquifer. In this approach, the hydrogeology of the aquifer and the dynamic boundary conditions are treated with separate models. A general mathematical formulation is presented where the calibration stage, the validation stage, and the prediction stage are formulated as a series of sub-model calibrations and solved using a general least squares routine. With this approach, it becomes possible to treat both the water level and the pumping rate in each cell as variables to be predicted. In most of the models, the pumping rates are known and the goals of the computation are to estimate the groundwater level. However, when for political or technical reasons access to some of the wells is difficult, the pumping rates are only partially known. Then, both groundwater levels and pumping rates are variables to be predicted by the groundwater model. A computer program was developed using MATLAB, with a Visual Basic graphical user interface using COM technology to access the advanced mathematical libraries. The approach was implemented with a real case study of the Yarkon–Taninim aquifer in Israel. The results indicate that the method is more stable than the classical approach.     

Optimizing water quality monitoring networks using continuous longitudinal monitoring data: a case study of Wen-Rui Tang River, Wenzhou, China

Identification of representative sampling sites is a critical issue in establishing an effective water quality monitoring program. This is especially important at the urban-agriculture interface where water quality conditions can change rapidly over short distances. The objective of this research was to optimize the spatial allocation of discrete monitoring sites for synoptic water quality monitoring through analysis of continuous longitudinal monitoring data collected by attaching a water quality sonde and GPS to a boat. Sampling was conducted six times from March to October 2009 along a 6.5 km segment of the Wen-Rui Tang River in eastern China that represented an urban-agricultural interface. When travelling at a velocity of ～2.4 km h(-1), this resulted in water quality measurements at ～20 m interval. Ammonia nitrogen (NH(4)(+)-N), electrical conductivity (EC), dissolved oxygen (DO), and turbidity data were collected and analyzed using Cluster Analysis (CA) to identify optimal locations for establishment of long-term monitoring sites. The analysis identified two distinct water quality segments for NH(4)(+)-N and EC and three distinct segments for DO and turbidity. According to our research results, the current fixed-location sampling sites should be adjusted to more effectively capture the distinct differences in the spatial distribution of water quality conditions. In addition, this methodology identified river reaches that require more comprehensive study of the factors leading to the changes in water quality within the identified river segment. The study demonstrates that continuous longitudinal monitoring can be a highly effective method for optimizing monitoring site locations for water quality studies.     

An Analytical Solution to Model Aquaculture Wetlands Subject to Intermittent Loading and Variable Initial Concentrations

Constructed wetlands are increasingly being used to treat intermittent releases of wastewaters from aquaculture ponds. Most constructed wetland design protocols are based on steady-state releases and as such are incapable of characterizing transient pollutant distribution arising from intermittent releases. Given the seasonal nature of aquaculture operations, the constructed wetland has to be developed using water from other nearby sources (termed as fill water in this study). In some instances, the pollutant concentrations in the fill water are noted to be greater than that measured in the aquaculture ponds. Also, the pollutants in the wetland may exhibit spatial variability due to geochemical interactions with wetland soils. These effects cannot be captured with existing steady-state design schemes. A mathematical model based on the tanks-in-series conceptualization is developed to overcome the aforementioned limitations. An analytical solution for multiple intermittent releases and spatially varying initial concentrations is obtained using Laplace transform and superposition techniques. The utility of the developed model is illustrated using an example based after a constructed wetland at the Loma Alta Shrimp Aquaculture Facility (LASAF) in South Texas. The developed model captures the extremely nonlinear relationship between the maximum discharge concentration and the hydraulic residence time. For the conditions assumed in this study, the extent of mixing is not of major concern when the spacing between the loadings is greater than 0.75 times the hydraulic residence time (HRT). Model results also suggest that the wetland should be oriented in a manner such that the internal geochemical production of the pollutant, if any, is near the inlet.     

A simple linear model for estimating ozone AOT40 at forest sites from raw passive sampling data

A rapid, empirical method is described for estimating weekly AOT40 from ozone concentrations measured with passive samplers at forest sites. The method is based on linear regression and was developed after three years of measurements in Trentino (northern Italy). It was tested against an independent set of data from passive sampler sites across Italy. It provides good weekly estimates compared with those measured by conventional monitors (0.85 ≤R(2)≤ 0.970; 97 ≤ RMSE ≤ 302). Estimates obtained using passive sampling at forest sites are comparable to those obtained by another estimation method based on modelling hourly concentrations (R(2) = 0.94; 131 ≤ RMSE ≤ 351). Regression coefficients of passive sampling are similar to those obtained with conventional monitors at forest sites. Testing against an independent dataset generated by passive sampling provided similar results (0.86 ≤R(2)≤ 0.99; 65 ≤ RMSE ≤ 478). Errors tend to accumulate when weekly AOT40 estimates are summed to obtain the total AOT40 over the May-July period, and the median deviation between the two estimation methods based on passive sampling is 11%. The method proposed does not require any assumptions, complex calculation or modelling technique, and can be useful when other estimation methods are not feasible, either in principle or in practice. However, the method is not useful when estimates of hourly concentrations are of interest.     

Monitoring cyclical air/water elemental mercury exchange

Previous experimental work has demonstrated that elemental mercury evasion from natural water displays a diel cycle; evasion rates during the day can be two to three times evasion rates observed at night. A study with polychlorinated biphenyls (PCBs) found that diurnal PCB air/water exchange rates exceeded nocturnal exchange rates by 32%. Given that the exchange rates of both PCBs and elemental mercury are dominated by the resistance in the aqueous thin film at the air/water interface and that water column elemental mercury concentrations in natural water bodies also display a diel cycle (and water column PCB concentrations do not) the findings here suggest that PCBs can serve as a tracer to assess the relative contribution of diel atmospheric temperature variations on elemental mercury air/water exchange rates. Using previously published data describing water column elemental mercury concentrations and the previously published diel mercury evasion model, four evasion scenarios are examined within the context of monitoring air/water toxicant exchange: constant atmospheric temperatures and constant water column elemental mercury concentrations; variable atmospheric temperatures and constant water column elemental mercury concentrations; constant atmospheric temperatures and variable water column elemental mercury concentrations; and variable atmospheric temperatures and variable water column elemental mercury concentrations. A scenario of monthly elemental mercury air/water exchange also is examined (at constant atmospheric and water column elemental mercury concentrations). Some of the findings include: (1) atmospheric temperature variations do have a significant effect on air/water toxicant exchange; (2) diel atmospheric temperature variations become more significant to overall diel toxicant exchange rates the closer the air/water system is to equilibrium conditions; (3) for refractory toxicants, average diel exchange rates are best estimated by averaging datasets obtained over a 24 h period or, at minimum, by measuring exchange rates at average atmospheric temperature values; (4) for elemental mercury, variable diel water column concentrations are likely to be the dominant contributor to variations in diel evasion rates; (5) diel atmospheric temperature variations amplify the magnitudes of both diel mercury evasion and absorption events and can shift maximum evasion rates to later in the day; (6) variations in monthly elemental mercury air/water exchange rates may exceed diel variations: and (7) 24 h and monthly monitoring efforts will likely be required to accurately describe diel and annual elemental mercury air/water exchange in a given system.     

Study on gas diffusion emitted from different height of point source

The flow and dispersion of stack-gas emitted from different elevated point source around flow obstacles in an urban environment have been investigated, using computational fluid dynamics models (CFD). The results were compared with the experimental results obtained from the diffusion wind tunnel under different conditions of thermal stability (stable, neutral or unstable). The flow and dispersion fields in the boundary layer in an urban environment were examined with different flow obstacles. Gaseous pollutant was discharged in the simulated boundary layer over the flat area. The CFD models used for the simulation were based on the steady-state Reynolds-Average Navier-Stoke equations (RANS) with kappa-epsilon turbulence models; standard kappa-epsilon and RNG kappa-epsilon models. The flow and dispersion data measured in the wind tunnel experiments were compared with the results of the CFD models in order to evaluate the prediction accuracy of the pollutant dispersion. The results of the CFD models showed good agreement with the results of the wind tunnel experiments. The results indicate that the turbulent velocity is reduced by the obstacles models. The maximum dispersion appears around the wake region of the obstacles.     

Bias Correction Techniques to Improve Air Quality Ensemble Predictions: Focus on O3 and PM Over Portugal

Five air quality models were applied over Portugal for July 2006 and used as ensemble members. Each model was used, with its original set up in terms of meteorology, parameterizations, boundary conditions and chemical mechanisms, but with the same emission data. The validation of the individual models and the ensemble of ozone (O₃) and particulate matter (PM) is performed using monitoring data from 22 background sites. The ensemble approach, based on the mean and median of the five models, did not improve significantly the skill scores due to large deviations in each ensemble member. Different bias correction techniques, including a subtraction of the mean bias and a multiplicative ratio adjustment, were implemented and analysed. The obtained datasets were compared against the individual modelled outputs using the bias, the root mean square error (RMSE) and the correlation coefficient. The applied bias correction techniques also improved the skill of the individual models and work equally well over the entire range of observed O₃ and PM values. The obtained results revealed that the best bias correction technique was the ratio adjustment with a 4-day training period, demonstrating significant improvements for both analysed pollutants. The increase in the ensemble skill found comprehends a bias reduction of 88 % for O₃, and 92 % for PM10, and also a decrease in 23 % for O₃ and 43 % for PM10 in what concerns the RMSE. In addition, a spatial bias correction approach was also examined with successful skills comparing to the uncorrected ensemble for both pollutants.     

STUMOD—a Tool for Predicting Fate and Transport of Nitrogen in Soil Treatment Units

A typical onsite wastewater treatment system consists of a septic tank and a soil treatment unit to treat wastewater before it is discharged through the vadose zone to an aquifer. A tool was developed for the purpose of predicting the fate and transport of nitrogen in soil treatment units (STUMOD or Soil Treatment Unit Model). STUMOD calculates nitrogen species concentrations and the fraction of total nitrogen reaching the aquifer or a specified soil depth. Input data include parameters for hydraulics and nutrient transport and transformation. An analytical solution is used to calculate the profile of pressure based on Darcy’s equation and the relationships between suction head, unsaturated hydraulic conductivity, and soil moisture. Chemical transport is based on simplification of the advection–dispersion equation. STUMOD is relatively simple to use but accounts for important processes such as ammonium sorption, nitrification, and denitrification. STUMOD accounts for the effect of soil moisture content (a surrogate for redox conditions) on nitrification and denitrification reactions. The model has provisions to handle the influence of temperature and organic carbon content on nitrogen transformation. Model outputs, generated based on input parameters obtained from extensive literature review, were compared to a numerical model and data from laboratory tests and field sites. Both measured data and STUMOD outputs show a relatively higher removal in clayey soils compared to sandy soils. Consistent with literature data for most soils, STUMOD predicted ammonium conversion to nitrate within the first foot below the trench infiltrative surface.     

Use of solvent extraction for the removal of heavy metals from liquid wastes

Solvent extraction is recommended as a suitable method for the removal of heavy metals from the waste waters of the chemical and electronic industries. Common extractants are organic compounds with molecular mass 200–450, almost insoluble in water (5–50 ppm), that selectively extract metals from aqueous solutions. On the basis of data from the literature, the extraction conditions are reviewed for the metals that cause problems in waste waters. The extraction conditions are understood to mean the type of extractant, anion present in the aqueous phase, and pH. As far as the data permit, the metals are ordered according to their extractability by given extractants. Special attention is paid to the main source of contamination of waste waters by heavy metals, namely electroplating plants. Two examples of the design of mixer-settler batteries are presented, which demonstrate continuous removal of heavy metals from the waste waters of electroplating plants to levels below the permissible limits. The first example illustrates removal of cadmium using triisooctylamine in a mixer-settler battery with the pertraction flow-arrangement; the second, removal of zinc and lead from rinse waters using the organophosphoric acid DEHPA (bis (2-ethylhexyl) phosphoric acid) in a battery with counter-current flow and reflux. The extraction equilibrium data were obtained in laboratory experiments using chloride solutions in concentrations close to those in the rinse waters of electroplating plants (erná and Volaufová, 1991).     

Saturated and Unsaturated Water Flow in Inclined Porous Media

This paper considers the two-dimensional saturated and unsaturated flow of water through inclined porous media, namely a waste dump or hill slope. Since the partial differential equation governing this water flow transforms from being parabolic to elliptic as the water flow varies from unsaturated to saturated, an iterative, finite differencing scheme is used to develop a numerical solution. The model can be used to investigate the effects that hill slope angle, depth of soil cover and hilltop width have on water accumulation in the dump and the time required for saturation to occur at different areas in the dump domain. The accuracy and reliability of the computer based solution is tested for two different boundary conditions – (1) no flow on all boundaries (i.e., the internal redistribution of soil moisture to steady state) and (2) a constant rainfall flux on the dump surface. Numerical studies then show the effects of changing the hill slope angle, depth of layer, and dump geometry on the flow characteristics in the dump.     

Inter-comparison of remote sensing sensing-based shoreline mapping techniques at different coastal stretches of India

Many techniques are available for detection of shorelines from multispectral satellite imagery, but the choice of a certain technique for a particular study area can be tough. Hence, for the first time in literature, an inter-comparison of the most widely used shoreline mapping techniques such as Normalized Difference Water Index (NDWI), Modified NDWI (MNDWI), Improved Band Ratio (IBR) Method, and Automatic Water Extraction Index (AWEI) has been done along four different coastal stretches of India using multitemporal Landsat data. The obtained results have been validated with the high-resolution images of Cartosat-2 (panchromatic) and multispectral images from Google Earth. Performance of the above indices has been analyzed based on the statistics, such as overall accuracy, kappa coefficient, user’s accuracy, producer’s accuracy, and the average deviation from the reference line. It is observed that the performance of NDWI and IBR techniques are dependent on the physical characteristics of the sites, and therefore, it varies from one site to another. Results indicate that unlike these two indices, the AWEI algorithm performs consistently well followed by MNDWI irrespective of the land cover types.     

Implementing In-Stream Nutrient Processes in Large-Scale Landscape Modeling for the Impact Assessment on Water Quality

For a long time, watershed models focused on the transport of chemicals from the terrestrial part of the watershed to the surface water bodies by leaching and erosion. After the substances had reached the surface water, they were routed through the channel network often without any further transformation. Today, there is a need to extend watershed models with in-stream processes to bring them closer to natural conditions and to enhance their usability as support tools for water management and water quality policies. This paper presents experience with implementing in-stream processes in a ecohydrological dynamic watershed model and its application on the large scale in the Saale River basin in Germany. Results demonstrate that new implemented water quality parameters like chlorophyll a concentrations or oxygen amount in the reach can be reproduced quite well, although the model results, compared with results achieved without taking into account algal and transformation processes in the river, show obvious improvement only for some of the examined nutrients. Finally, some climate and water management scenarios expected to impact in-stream processes in the Saale basin were run. Their results illustrate the relative importance of physical boundary conditions on the amount and concentration of the phytoplankton, which leads to the conclusion that measures to improve water quality should not only take nutrient inputs into account but also climate influences and river morphology.     

A phycocyanin probe as a tool for monitoring cyanobacteria in freshwater bodies

In many countries, the presence of cyanobacteria in freshwater bodies used for both drinking water and recreational purposes is under increasing public health attention. Water managers are considering how to implement monitoring that leads to a minimization of the risks incurred by the users of potentially contaminated sites. To address this question, this study involved assessing the performance of a submersible probe for measuring phycocyanin-specific fluorescence as a function of cyanobacterial biomass, with the aim of applying it as a tool for surveillance management. Its advantages and limits compared to more traditional analyses are discussed. The monitoring of cyanobacteria in the water bodies of western France was carried out using a minifluorimeter specific to the fluorescence of phycocyanin, a pigment specific to cyanobacteria. The results are compared with the analyses recommended by the World Health Organisation (chlorophyll a and cell counting). This study based on nearly 800 samples shows a significant correlation between the phycocyanin content and the cyanobacterial biomass, expressed as the number of cells per mL (R2 = 0.73). This submersible probe is simple and rapid to use, making it possible to take into account horizontal and vertical heterogeneities in the proliferation growth. In this way, we are able to detect at an early stage the conditions that could potentially lead to a risk, in order to start sampling. Due to its sensitivity, this tool proves suitable for monitoring aimed at reducing the risks incurred by the users of contaminated sites and launching preventative actions. The use of the phycocyanin probe provides an effective tool to complement traditional analyses of cyanobacterial presence. It is suggested that a surveillance protocol based on phycocyanin concentration can significantly improved the accuracy of the extent of cyanobacterial bloom development in the light of spatial and temporal variabilities associated with these occurrences.     

Occurrence of PCBs in raw and finished drinking water at seven public water systems along the Hudson River

Polychlorinated biphenyls (PCBs) were measured in raw and finished drinking water at seven Public Water Systems (PWSs) along the Hudson River as part of a baseline monitoring program prior to the extensive sediment dredging of the Upper Hudson River. Water samples were either analyzed using an Aroclor method (based on USEPA Method 508) or a congener method (Modified Green Bay Mass Balance Method). Using the congener-based method, raw water concentrations ranged from <9.3 to 164.3 ng/L and finished water concentrations ranged from <9.3 to 186.6 ng/L. Using the Aroclor method, finished water concentrations ranged from <5.0 to 200.9 ng/L. Most finished water samples above 73.0 ng/L were from a PWS with wells drilled near the river. Excluding the well data, total PCB concentrations in raw water at systems in the Upper River were similar to concentrations at systems in the Lower River, though the congener patterns differed. Paired comparison of total PCB concentrations using the two analytical methods showed good agreement, although raw water showed a different relationship than finished water.     

滇池水环境质量综合评价方法

研究了2000—2017年滇池水环境质量时空变化和生态环境时空演替规律,结合滇池目前水环境质量监测工作要求,利用以往研究成果筛选出的36个评价指标,采用层次分析法,构建滇池水生态健康评价方法模型,得出滇池水生态健康评价标准,输入滇池湖体2011—2017年4个监测点位相关指标监测数据进行方法验证。验证结果表明,建立的评价方法能较为全面、准确、宏观地反映滇池水环境质量状况,与现有综合评价方法相比,能综合反映现有评价方法结果,且具有宏观、全面等优点,能很好地响应生物多样性变化,服务滇池的治理工作。